1
2 /*
3 We define the memory operations here. The reason we just do not use
4 the standard memory routines in the PETSc code is that on some machines
5 they are broken.
6
7 */
8 #include <petscsys.h> /*I "petscsys.h" I*/
9 #include <petscbt.h>
10 #include <../src/sys/utils/ftn-kernels/fcopy.h>
11
12 /*@
13 PetscMemcmp - Compares two byte streams in memory.
14
15 Not Collective
16
17 Input Parameters:
18 + str1 - Pointer to the first byte stream
19 . str2 - Pointer to the second byte stream
20 - len - The length of the byte stream
21 (both str1 and str2 are assumed to be of length len)
22
23 Output Parameters:
24 . e - PETSC_TRUE if equal else PETSC_FALSE.
25
26 Level: intermediate
27
28 Note:
29 PetscArraycmp() is preferred
30 This routine is anologous to memcmp()
31
32 .seealso: PetscMemcpy(), PetscMemcmp(), PetscArrayzero(), PetscMemzero(), PetscArraycmp(), PetscArraycpy(), PetscStrallocpy(),
33 PetscArraymove()
34 @*/
PetscMemcmp(const void * str1,const void * str2,size_t len,PetscBool * e)35 PetscErrorCode PetscMemcmp(const void *str1,const void *str2,size_t len,PetscBool *e)
36 {
37 int r;
38
39 PetscFunctionBegin;
40 if (len > 0 && !str1) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_NULL,"Trying to compare at a null pointer");
41 if (len > 0 && !str2) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_NULL,"Trying to compare at a null pointer");
42 r = memcmp((char*)str1,(char*)str2,len);
43 if (!r) *e = PETSC_TRUE;
44 else *e = PETSC_FALSE;
45 PetscFunctionReturn(0);
46 }
47
48 #if defined(PETSC_HAVE_HWLOC)
49 #include <petsc/private/petscimpl.h>
50 #include <hwloc.h>
51
52 /*@C
53 PetscProcessPlacementView - display the MPI process placement by core
54
55 Input Parameter:
56 . viewer - ASCII viewer to display the results on
57
58 Level: intermediate
59
60 Notes:
61 Requires that PETSc be installed with hwloc, for example using --download-hwloc
62 @*/
PetscProcessPlacementView(PetscViewer viewer)63 PetscErrorCode PetscProcessPlacementView(PetscViewer viewer)
64 {
65 PetscErrorCode ierr;
66 PetscBool isascii;
67 PetscMPIInt rank;
68 hwloc_bitmap_t set;
69 hwloc_topology_t topology;
70 int err;
71
72 PetscFunctionBegin;
73 PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,1);
74 ierr = PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii);CHKERRQ(ierr);
75 if (!isascii) SETERRQ(PetscObjectComm((PetscObject)viewer),PETSC_ERR_SUP,"Only ASCII viewer is supported");
76
77 ierr = MPI_Comm_rank(MPI_COMM_WORLD,&rank);CHKERRQ(ierr);
78 hwloc_topology_init ( &topology);
79 hwloc_topology_load ( topology);
80 set = hwloc_bitmap_alloc();
81
82 err = hwloc_get_proc_cpubind(topology, getpid(), set, HWLOC_CPUBIND_PROCESS);
83 if (err) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error %d from hwloc_get_proc_cpubind()",err);
84 ierr = PetscViewerASCIIPushSynchronized(viewer);CHKERRQ(ierr);
85 ierr = PetscViewerASCIISynchronizedPrintf(viewer,"MPI rank %d Process id: %d coreid %d\n",rank,getpid(),hwloc_bitmap_first(set));CHKERRQ(ierr);
86 ierr = PetscViewerFlush(viewer);CHKERRQ(ierr);
87 hwloc_bitmap_free(set);
88 hwloc_topology_destroy(topology);
89 PetscFunctionReturn(0);
90 }
91 #endif
92
93