1e5c89e4eSSatish Balay #define PETSC_DLL 2e5c89e4eSSatish Balay /* 3e5c89e4eSSatish Balay 4e5c89e4eSSatish Balay This file defines part of the initialization of PETSc 5e5c89e4eSSatish Balay 6e5c89e4eSSatish Balay This file uses regular malloc and free because it cannot know 7e5c89e4eSSatish Balay what malloc is being used until it has already processed the input. 8e5c89e4eSSatish Balay */ 9e5c89e4eSSatish Balay 10e5c89e4eSSatish Balay #include "petsc.h" /*I "petsc.h" I*/ 11e5c89e4eSSatish Balay #include "petscsys.h" 12e5c89e4eSSatish Balay #if defined(PETSC_HAVE_STDLIB_H) 13e5c89e4eSSatish Balay #include <stdlib.h> 14e5c89e4eSSatish Balay #endif 15e5c89e4eSSatish Balay #if defined(PETSC_HAVE_MALLOC_H) 16e5c89e4eSSatish Balay #include <malloc.h> 17e5c89e4eSSatish Balay #endif 18e5c89e4eSSatish Balay #include "petscfix.h" 196dc8fec2Sbcordonn #include "zope.h" 20e5c89e4eSSatish Balay 21e5c89e4eSSatish Balay /* ------------------------Nasty global variables -------------------------------*/ 22e5c89e4eSSatish Balay /* 23e5c89e4eSSatish Balay Indicates if PETSc started up MPI, or it was 24e5c89e4eSSatish Balay already started before PETSc was initialized. 25e5c89e4eSSatish Balay */ 26e5c89e4eSSatish Balay PetscTruth PETSC_DLLEXPORT PetscBeganMPI = PETSC_FALSE; 27e5c89e4eSSatish Balay PetscTruth PETSC_DLLEXPORT PetscInitializeCalled = PETSC_FALSE; 28e5c89e4eSSatish Balay PetscTruth PETSC_DLLEXPORT PetscFinalizeCalled = PETSC_FALSE; 29e5c89e4eSSatish Balay PetscMPIInt PETSC_DLLEXPORT PetscGlobalRank = -1; 30e5c89e4eSSatish Balay PetscMPIInt PETSC_DLLEXPORT PetscGlobalSize = -1; 31e5c89e4eSSatish Balay 32e5c89e4eSSatish Balay #if defined(PETSC_USE_COMPLEX) 33e5c89e4eSSatish Balay #if defined(PETSC_COMPLEX_INSTANTIATE) 34e5c89e4eSSatish Balay template <> class std::complex<double>; /* instantiate complex template class */ 35e5c89e4eSSatish Balay #endif 36e5c89e4eSSatish Balay MPI_Datatype PETSC_DLLEXPORT MPIU_COMPLEX; 37e5c89e4eSSatish Balay PetscScalar PETSC_DLLEXPORT PETSC_i; 38e5c89e4eSSatish Balay #else 39e5c89e4eSSatish Balay PetscScalar PETSC_DLLEXPORT PETSC_i = 0.0; 40e5c89e4eSSatish Balay #endif 41e5c89e4eSSatish Balay MPI_Datatype PETSC_DLLEXPORT MPIU_2SCALAR = 0; 42e5c89e4eSSatish Balay MPI_Datatype PETSC_DLLEXPORT MPIU_2INT = 0; 43e5c89e4eSSatish Balay /* 44e5c89e4eSSatish Balay These are needed by petscbt.h 45e5c89e4eSSatish Balay */ 46e5c89e4eSSatish Balay char PETSC_DLLEXPORT _BT_mask = ' '; 47e5c89e4eSSatish Balay char PETSC_DLLEXPORT _BT_c = ' '; 48e5c89e4eSSatish Balay PetscInt PETSC_DLLEXPORT _BT_idx = 0; 49e5c89e4eSSatish Balay 50e5c89e4eSSatish Balay /* 51e5c89e4eSSatish Balay Function that is called to display all error messages 52e5c89e4eSSatish Balay */ 53e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT (*PetscErrorPrintf)(const char [],...) = PetscErrorPrintfDefault; 54e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT (*PetscHelpPrintf)(MPI_Comm,const char [],...) = PetscHelpPrintfDefault; 55e5c89e4eSSatish Balay 56e5c89e4eSSatish Balay /* ------------------------------------------------------------------------------*/ 57e5c89e4eSSatish Balay /* 58e5c89e4eSSatish Balay Optional file where all PETSc output from various prints is saved 59e5c89e4eSSatish Balay */ 60e5c89e4eSSatish Balay FILE *petsc_history = PETSC_NULL; 61e5c89e4eSSatish Balay 62e5c89e4eSSatish Balay #undef __FUNCT__ 63e5c89e4eSSatish Balay #define __FUNCT__ "PetscLogOpenHistoryFile" 64e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscLogOpenHistoryFile(const char filename[],FILE **fd) 65e5c89e4eSSatish Balay { 66e5c89e4eSSatish Balay PetscErrorCode ierr; 67e5c89e4eSSatish Balay PetscMPIInt rank,size; 68e5c89e4eSSatish Balay char pfile[PETSC_MAX_PATH_LEN],pname[PETSC_MAX_PATH_LEN],fname[PETSC_MAX_PATH_LEN],date[64]; 69e5c89e4eSSatish Balay char version[256]; 70e5c89e4eSSatish Balay 71e5c89e4eSSatish Balay PetscFunctionBegin; 72e5c89e4eSSatish Balay ierr = MPI_Comm_rank(PETSC_COMM_WORLD,&rank);CHKERRQ(ierr); 73e5c89e4eSSatish Balay if (!rank) { 74e5c89e4eSSatish Balay char arch[10]; 75e5c89e4eSSatish Balay ierr = PetscGetArchType(arch,10);CHKERRQ(ierr); 76e5c89e4eSSatish Balay ierr = PetscGetDate(date,64);CHKERRQ(ierr); 773f6e4ae9SSatish Balay ierr = PetscGetVersion(&version,256);CHKERRQ(ierr); 78e5c89e4eSSatish Balay ierr = MPI_Comm_size(PETSC_COMM_WORLD,&size);CHKERRQ(ierr); 79e5c89e4eSSatish Balay if (filename) { 80e5c89e4eSSatish Balay ierr = PetscFixFilename(filename,fname);CHKERRQ(ierr); 81e5c89e4eSSatish Balay } else { 82e5c89e4eSSatish Balay ierr = PetscGetHomeDirectory(pfile,240);CHKERRQ(ierr); 83e5c89e4eSSatish Balay ierr = PetscStrcat(pfile,"/.petschistory");CHKERRQ(ierr); 84e5c89e4eSSatish Balay ierr = PetscFixFilename(pfile,fname);CHKERRQ(ierr); 85e5c89e4eSSatish Balay } 86e5c89e4eSSatish Balay 87e5c89e4eSSatish Balay *fd = fopen(fname,"a"); if (!fd) SETERRQ1(PETSC_ERR_FILE_OPEN,"Cannot open file: %s",fname); 88e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"---------------------------------------------------------\n");CHKERRQ(ierr); 89e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"%s %s\n",version,date);CHKERRQ(ierr); 90e5c89e4eSSatish Balay ierr = PetscGetProgramName(pname,PETSC_MAX_PATH_LEN);CHKERRQ(ierr); 91e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"%s on a %s, %d proc. with options:\n",pname,arch,size);CHKERRQ(ierr); 92e5c89e4eSSatish Balay ierr = PetscOptionsPrint(*fd);CHKERRQ(ierr); 93e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"---------------------------------------------------------\n");CHKERRQ(ierr); 94e5c89e4eSSatish Balay fflush(*fd); 95e5c89e4eSSatish Balay } 96e5c89e4eSSatish Balay PetscFunctionReturn(0); 97e5c89e4eSSatish Balay } 98e5c89e4eSSatish Balay 99e5c89e4eSSatish Balay #undef __FUNCT__ 100e5c89e4eSSatish Balay #define __FUNCT__ "PetscLogCloseHistoryFile" 101e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscLogCloseHistoryFile(FILE **fd) 102e5c89e4eSSatish Balay { 103e5c89e4eSSatish Balay PetscErrorCode ierr; 104e5c89e4eSSatish Balay PetscMPIInt rank; 105e5c89e4eSSatish Balay char date[64]; 106e5c89e4eSSatish Balay 107e5c89e4eSSatish Balay PetscFunctionBegin; 108e5c89e4eSSatish Balay ierr = MPI_Comm_rank(PETSC_COMM_WORLD,&rank);CHKERRQ(ierr); 109e5c89e4eSSatish Balay if (!rank) { 110e5c89e4eSSatish Balay ierr = PetscGetDate(date,64);CHKERRQ(ierr); 111e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"---------------------------------------------------------\n");CHKERRQ(ierr); 112e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"Finished at %s\n",date);CHKERRQ(ierr); 113e5c89e4eSSatish Balay ierr = PetscFPrintf(PETSC_COMM_SELF,*fd,"---------------------------------------------------------\n");CHKERRQ(ierr); 114e5c89e4eSSatish Balay fflush(*fd); 115e5c89e4eSSatish Balay fclose(*fd); 116e5c89e4eSSatish Balay } 117e5c89e4eSSatish Balay PetscFunctionReturn(0); 118e5c89e4eSSatish Balay } 119e5c89e4eSSatish Balay 120e5c89e4eSSatish Balay /* ------------------------------------------------------------------------------*/ 121e5c89e4eSSatish Balay 122e5c89e4eSSatish Balay /* 123e5c89e4eSSatish Balay This is ugly and probably belongs somewhere else, but I want to 124e5c89e4eSSatish Balay be able to put a true MPI abort error handler with command line args. 125e5c89e4eSSatish Balay 126e5c89e4eSSatish Balay This is so MPI errors in the debugger will leave all the stack 127e5c89e4eSSatish Balay frames. The default abort cleans up and exits. 128e5c89e4eSSatish Balay */ 129e5c89e4eSSatish Balay 130e5c89e4eSSatish Balay #undef __FUNCT__ 131e5c89e4eSSatish Balay #define __FUNCT__ "Petsc_MPI_AbortOnError" 132e5c89e4eSSatish Balay void Petsc_MPI_AbortOnError(MPI_Comm *comm,PetscMPIInt *flag) 133e5c89e4eSSatish Balay { 134e5c89e4eSSatish Balay PetscFunctionBegin; 135e5c89e4eSSatish Balay (*PetscErrorPrintf)("MPI error %d\n",(int)*flag); 136e5c89e4eSSatish Balay abort(); 137e5c89e4eSSatish Balay } 138e5c89e4eSSatish Balay 139e5c89e4eSSatish Balay #undef __FUNCT__ 140e5c89e4eSSatish Balay #define __FUNCT__ "Petsc_MPI_DebuggerOnError" 141e5c89e4eSSatish Balay void Petsc_MPI_DebuggerOnError(MPI_Comm *comm,PetscMPIInt *flag) 142e5c89e4eSSatish Balay { 143e5c89e4eSSatish Balay PetscErrorCode ierr; 144e5c89e4eSSatish Balay 145e5c89e4eSSatish Balay PetscFunctionBegin; 146e5c89e4eSSatish Balay (*PetscErrorPrintf)("MPI error %d\n",(int)*flag); 147e5c89e4eSSatish Balay ierr = PetscAttachDebugger(); 148e5c89e4eSSatish Balay if (ierr) { /* hopeless so get out */ 149e5c89e4eSSatish Balay MPI_Finalize(); 150e5c89e4eSSatish Balay exit(*flag); 151e5c89e4eSSatish Balay } 152e5c89e4eSSatish Balay } 153e5c89e4eSSatish Balay 154e5c89e4eSSatish Balay #undef __FUNCT__ 155e5c89e4eSSatish Balay #define __FUNCT__ "PetscEnd" 156e5c89e4eSSatish Balay /*@C 157e5c89e4eSSatish Balay PetscEnd - Calls PetscFinalize() and then ends the program. This is useful if one 158e5c89e4eSSatish Balay wishes a clean exit somewhere deep in the program. 159e5c89e4eSSatish Balay 160e5c89e4eSSatish Balay Collective on PETSC_COMM_WORLD 161e5c89e4eSSatish Balay 162e5c89e4eSSatish Balay Options Database Keys are the same as for PetscFinalize() 163e5c89e4eSSatish Balay 164e5c89e4eSSatish Balay Level: advanced 165e5c89e4eSSatish Balay 166e5c89e4eSSatish Balay Note: 167e5c89e4eSSatish Balay See PetscInitialize() for more general runtime options. 168e5c89e4eSSatish Balay 169e5c89e4eSSatish Balay .seealso: PetscInitialize(), PetscOptionsPrint(), PetscMallocDump(), PetscMPIDump(), PetscFinalize() 170e5c89e4eSSatish Balay @*/ 171e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscEnd(void) 172e5c89e4eSSatish Balay { 173e5c89e4eSSatish Balay PetscFunctionBegin; 174e5c89e4eSSatish Balay PetscFinalize(); 175e5c89e4eSSatish Balay exit(0); 176e5c89e4eSSatish Balay return 0; 177e5c89e4eSSatish Balay } 178e5c89e4eSSatish Balay 179e5c89e4eSSatish Balay PetscTruth PetscOptionsPublish = PETSC_FALSE; 180e5c89e4eSSatish Balay EXTERN PetscErrorCode PetscSetUseTrMalloc_Private(void); 181e5c89e4eSSatish Balay extern PetscTruth petscsetmallocvisited; 182e5c89e4eSSatish Balay static char emacsmachinename[256]; 183e5c89e4eSSatish Balay 184e5c89e4eSSatish Balay PetscErrorCode (*PetscExternalVersionFunction)(MPI_Comm) = 0; 185e5c89e4eSSatish Balay PetscErrorCode (*PetscExternalHelpFunction)(MPI_Comm) = 0; 186e5c89e4eSSatish Balay 187e5c89e4eSSatish Balay #undef __FUNCT__ 188e5c89e4eSSatish Balay #define __FUNCT__ "PetscSetHelpVersionFunctions" 189e5c89e4eSSatish Balay /*@C 190e5c89e4eSSatish Balay PetscSetHelpVersionFunctions - Sets functions that print help and version information 191e5c89e4eSSatish Balay before the PETSc help and version information is printed. Must call BEFORE PetscInitialize(). 192e5c89e4eSSatish Balay This routine enables a "higher-level" package that uses PETSc to print its messages first. 193e5c89e4eSSatish Balay 194e5c89e4eSSatish Balay Input Parameter: 195e5c89e4eSSatish Balay + help - the help function (may be PETSC_NULL) 196e5c89e4eSSatish Balay - version - the version function (may be PETSc null) 197e5c89e4eSSatish Balay 198e5c89e4eSSatish Balay Level: developer 199e5c89e4eSSatish Balay 200e5c89e4eSSatish Balay Concepts: package help message 201e5c89e4eSSatish Balay 202e5c89e4eSSatish Balay @*/ 203e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscSetHelpVersionFunctions(PetscErrorCode (*help)(MPI_Comm),PetscErrorCode (*version)(MPI_Comm)) 204e5c89e4eSSatish Balay { 205e5c89e4eSSatish Balay PetscFunctionBegin; 206e5c89e4eSSatish Balay PetscExternalHelpFunction = help; 207e5c89e4eSSatish Balay PetscExternalVersionFunction = version; 208e5c89e4eSSatish Balay PetscFunctionReturn(0); 209e5c89e4eSSatish Balay } 210e5c89e4eSSatish Balay 211e5c89e4eSSatish Balay #undef __FUNCT__ 212e5c89e4eSSatish Balay #define __FUNCT__ "PetscOptionsCheckInitial_Private" 213e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscOptionsCheckInitial_Private(void) 214e5c89e4eSSatish Balay { 215e5c89e4eSSatish Balay char string[64],mname[PETSC_MAX_PATH_LEN],*f; 216e5c89e4eSSatish Balay MPI_Comm comm = PETSC_COMM_WORLD; 2176dc8fec2Sbcordonn PetscTruth flg1,flg2,flg3,flag, flgz; 218e5c89e4eSSatish Balay PetscErrorCode ierr; 219e5c89e4eSSatish Balay PetscInt si; 220e5c89e4eSSatish Balay int i; 221e5c89e4eSSatish Balay PetscMPIInt rank; 222e5c89e4eSSatish Balay char version[256]; 223e5c89e4eSSatish Balay 224e5c89e4eSSatish Balay PetscFunctionBegin; 225e5c89e4eSSatish Balay ierr = MPI_Comm_rank(PETSC_COMM_WORLD,&rank);CHKERRQ(ierr); 226e5c89e4eSSatish Balay 227e5c89e4eSSatish Balay /* 228e5c89e4eSSatish Balay Setup the memory management; support for tracing malloc() usage 229e5c89e4eSSatish Balay */ 230e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-malloc_log",&flg3);CHKERRQ(ierr); 231e5c89e4eSSatish Balay #if defined(PETSC_USE_DEBUG) 232e5c89e4eSSatish Balay ierr = PetscOptionsGetTruth(PETSC_NULL,"-malloc",&flg1,&flg2);CHKERRQ(ierr); 233e5c89e4eSSatish Balay if ((!flg2 || flg1) && !petscsetmallocvisited) { 234e5c89e4eSSatish Balay ierr = PetscSetUseTrMalloc_Private();CHKERRQ(ierr); 235e5c89e4eSSatish Balay } 236e5c89e4eSSatish Balay #else 237e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-malloc_dump",&flg1);CHKERRQ(ierr); 238e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-malloc",&flg2);CHKERRQ(ierr); 239e5c89e4eSSatish Balay if (flg1 || flg2 || flg3) {ierr = PetscSetUseTrMalloc_Private();CHKERRQ(ierr);} 240e5c89e4eSSatish Balay #endif 241e5c89e4eSSatish Balay if (flg3) { 242e5c89e4eSSatish Balay ierr = PetscMallocSetDumpLog();CHKERRQ(ierr); 243e5c89e4eSSatish Balay } 244e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-malloc_debug",&flg1);CHKERRQ(ierr); 245e5c89e4eSSatish Balay if (flg1) { 246e5c89e4eSSatish Balay ierr = PetscSetUseTrMalloc_Private();CHKERRQ(ierr); 247e5c89e4eSSatish Balay ierr = PetscMallocDebug(PETSC_TRUE);CHKERRQ(ierr); 248e5c89e4eSSatish Balay } 249e5c89e4eSSatish Balay 2507783f70dSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-malloc_info",&flg1);CHKERRQ(ierr); 2517783f70dSSatish Balay if (!flg1) { 252e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-memory_info",&flg1);CHKERRQ(ierr); 2537783f70dSSatish Balay } 254e5c89e4eSSatish Balay if (flg1) { 255e5c89e4eSSatish Balay ierr = PetscMemorySetGetMaximumUsage();CHKERRQ(ierr); 256e5c89e4eSSatish Balay } 257e5c89e4eSSatish Balay 258e5c89e4eSSatish Balay /* 259e5c89e4eSSatish Balay Set the display variable for graphics 260e5c89e4eSSatish Balay */ 261e5c89e4eSSatish Balay ierr = PetscSetDisplay();CHKERRQ(ierr); 262e5c89e4eSSatish Balay 263e5c89e4eSSatish Balay /* 264e5c89e4eSSatish Balay Print the PETSc version information 265e5c89e4eSSatish Balay */ 266e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-v",&flg1);CHKERRQ(ierr); 267e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-version",&flg2);CHKERRQ(ierr); 268e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-help",&flg3);CHKERRQ(ierr); 269e5c89e4eSSatish Balay if (flg1 || flg2 || flg3){ 270e5c89e4eSSatish Balay 271e5c89e4eSSatish Balay /* 272e5c89e4eSSatish Balay Print "higher-level" package version message 273e5c89e4eSSatish Balay */ 274e5c89e4eSSatish Balay if (PetscExternalVersionFunction) { 275e5c89e4eSSatish Balay ierr = (*PetscExternalVersionFunction)(comm);CHKERRQ(ierr); 276e5c89e4eSSatish Balay } 277e5c89e4eSSatish Balay 2783f6e4ae9SSatish Balay ierr = PetscGetVersion(&version,256);CHKERRQ(ierr); 279e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"--------------------------------------------\ 280e5c89e4eSSatish Balay ------------------------------\n");CHKERRQ(ierr); 281e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"%s\n",version);CHKERRQ(ierr); 282e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"%s",PETSC_AUTHOR_INFO);CHKERRQ(ierr); 283e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"See docs/copyright.html for copyright information\n");CHKERRQ(ierr); 284e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"See docs/changes/index.html for recent updates.\n");CHKERRQ(ierr); 285e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"See docs/troubleshooting.html for problems.\n");CHKERRQ(ierr); 286e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"See docs/manualpages/index.html for help. \n");CHKERRQ(ierr); 287e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"Libraries linked from %s\n",PETSC_LIB_DIR);CHKERRQ(ierr); 288e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"--------------------------------------------\ 289e5c89e4eSSatish Balay ------------------------------\n");CHKERRQ(ierr); 290e5c89e4eSSatish Balay } 291e5c89e4eSSatish Balay 292e5c89e4eSSatish Balay /* 293e5c89e4eSSatish Balay Print "higher-level" package help message 294e5c89e4eSSatish Balay */ 295e5c89e4eSSatish Balay if (flg3){ 296e5c89e4eSSatish Balay if (PetscExternalHelpFunction) { 297e5c89e4eSSatish Balay ierr = (*PetscExternalHelpFunction)(comm);CHKERRQ(ierr); 298e5c89e4eSSatish Balay } 299e5c89e4eSSatish Balay } 300e5c89e4eSSatish Balay 301e5c89e4eSSatish Balay /* 302e5c89e4eSSatish Balay Setup the error handling 303e5c89e4eSSatish Balay */ 304e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-fp_trap",&flg1);CHKERRQ(ierr); 305e5c89e4eSSatish Balay if (flg1) { ierr = PetscSetFPTrap(PETSC_FP_TRAP_ON);CHKERRQ(ierr); } 306e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-on_error_abort",&flg1);CHKERRQ(ierr); 307e5c89e4eSSatish Balay if (flg1) { ierr = PetscPushErrorHandler(PetscAbortErrorHandler,0);CHKERRQ(ierr)} 308e8fb0fc0SBarry Smith ierr = PetscOptionsHasName(PETSC_NULL,"-on_error_mpiabort",&flg1);CHKERRQ(ierr); 309e8fb0fc0SBarry Smith if (flg1) { ierr = PetscPushErrorHandler(PetscMPIAbortErrorHandler,0);CHKERRQ(ierr)} 310e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-mpi_return_on_error",&flg1);CHKERRQ(ierr); 311e5c89e4eSSatish Balay if (flg1) { 312e5c89e4eSSatish Balay ierr = MPI_Errhandler_set(comm,MPI_ERRORS_RETURN);CHKERRQ(ierr); 313e5c89e4eSSatish Balay } 314e8fb0fc0SBarry Smith ierr = PetscOptionsHasName(PETSC_NULL,"-error_output_none",&flg1);CHKERRQ(ierr); 315e8fb0fc0SBarry Smith if (flg1) { 316e8fb0fc0SBarry Smith PetscErrorPrintf = PetscErrorPrintfNone; 317e8fb0fc0SBarry Smith } 318e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-no_signal_handler",&flg1);CHKERRQ(ierr); 319e5c89e4eSSatish Balay if (!flg1) { ierr = PetscPushSignalHandler(PetscDefaultSignalHandler,(void*)0);CHKERRQ(ierr) } 320e5c89e4eSSatish Balay 321e5c89e4eSSatish Balay /* 322e5c89e4eSSatish Balay Setup debugger information 323e5c89e4eSSatish Balay */ 324e5c89e4eSSatish Balay ierr = PetscSetDefaultDebugger();CHKERRQ(ierr); 325e5c89e4eSSatish Balay ierr = PetscOptionsGetString(PETSC_NULL,"-on_error_attach_debugger",string,64,&flg1);CHKERRQ(ierr); 326e5c89e4eSSatish Balay if (flg1) { 327e5c89e4eSSatish Balay MPI_Errhandler err_handler; 328e5c89e4eSSatish Balay 329e5c89e4eSSatish Balay ierr = PetscSetDebuggerFromString(string);CHKERRQ(ierr); 330e5c89e4eSSatish Balay ierr = MPI_Errhandler_create((MPI_Handler_function*)Petsc_MPI_DebuggerOnError,&err_handler);CHKERRQ(ierr); 331e5c89e4eSSatish Balay ierr = MPI_Errhandler_set(comm,err_handler);CHKERRQ(ierr); 332e5c89e4eSSatish Balay ierr = PetscPushErrorHandler(PetscAttachDebuggerErrorHandler,0);CHKERRQ(ierr); 333e5c89e4eSSatish Balay } 334e5c89e4eSSatish Balay ierr = PetscOptionsGetString(PETSC_NULL,"-start_in_debugger",string,64,&flg1);CHKERRQ(ierr); 335e5c89e4eSSatish Balay ierr = PetscOptionsGetString(PETSC_NULL,"-stop_for_debugger",string,64,&flg2);CHKERRQ(ierr); 336e5c89e4eSSatish Balay if (flg1 || flg2) { 337e5c89e4eSSatish Balay PetscMPIInt size; 338e5c89e4eSSatish Balay PetscInt lsize,*nodes; 339e5c89e4eSSatish Balay MPI_Errhandler err_handler; 340e5c89e4eSSatish Balay /* 341e5c89e4eSSatish Balay we have to make sure that all processors have opened 342e5c89e4eSSatish Balay connections to all other processors, otherwise once the 343e5c89e4eSSatish Balay debugger has stated it is likely to receive a SIGUSR1 344e5c89e4eSSatish Balay and kill the program. 345e5c89e4eSSatish Balay */ 346e5c89e4eSSatish Balay ierr = MPI_Comm_size(PETSC_COMM_WORLD,&size);CHKERRQ(ierr); 347e5c89e4eSSatish Balay if (size > 2) { 348e5c89e4eSSatish Balay PetscMPIInt dummy; 349e5c89e4eSSatish Balay MPI_Status status; 350e5c89e4eSSatish Balay for (i=0; i<size; i++) { 351e5c89e4eSSatish Balay if (rank != i) { 352e5c89e4eSSatish Balay ierr = MPI_Send(&dummy,1,MPI_INT,i,109,PETSC_COMM_WORLD);CHKERRQ(ierr); 353e5c89e4eSSatish Balay } 354e5c89e4eSSatish Balay } 355e5c89e4eSSatish Balay for (i=0; i<size; i++) { 356e5c89e4eSSatish Balay if (rank != i) { 357e5c89e4eSSatish Balay ierr = MPI_Recv(&dummy,1,MPI_INT,i,109,PETSC_COMM_WORLD,&status);CHKERRQ(ierr); 358e5c89e4eSSatish Balay } 359e5c89e4eSSatish Balay } 360e5c89e4eSSatish Balay } 361e5c89e4eSSatish Balay /* check if this processor node should be in debugger */ 362e5c89e4eSSatish Balay ierr = PetscMalloc(size*sizeof(PetscInt),&nodes);CHKERRQ(ierr); 363e5c89e4eSSatish Balay lsize = size; 364e5c89e4eSSatish Balay ierr = PetscOptionsGetIntArray(PETSC_NULL,"-debugger_nodes",nodes,&lsize,&flag);CHKERRQ(ierr); 365e5c89e4eSSatish Balay if (flag) { 366e5c89e4eSSatish Balay for (i=0; i<lsize; i++) { 367e5c89e4eSSatish Balay if (nodes[i] == rank) { flag = PETSC_FALSE; break; } 368e5c89e4eSSatish Balay } 369e5c89e4eSSatish Balay } 370e5c89e4eSSatish Balay if (!flag) { 371e5c89e4eSSatish Balay ierr = PetscSetDebuggerFromString(string);CHKERRQ(ierr); 372e5c89e4eSSatish Balay ierr = PetscPushErrorHandler(PetscAbortErrorHandler,0);CHKERRQ(ierr); 373e5c89e4eSSatish Balay if (flg1) { 374e5c89e4eSSatish Balay ierr = PetscAttachDebugger();CHKERRQ(ierr); 375e5c89e4eSSatish Balay } else { 376e5c89e4eSSatish Balay ierr = PetscStopForDebugger();CHKERRQ(ierr); 377e5c89e4eSSatish Balay } 378e5c89e4eSSatish Balay ierr = MPI_Errhandler_create((MPI_Handler_function*)Petsc_MPI_AbortOnError,&err_handler);CHKERRQ(ierr); 379e5c89e4eSSatish Balay ierr = MPI_Errhandler_set(comm,err_handler);CHKERRQ(ierr); 380e5c89e4eSSatish Balay } 381e5c89e4eSSatish Balay ierr = PetscFree(nodes);CHKERRQ(ierr); 382e5c89e4eSSatish Balay } 383e5c89e4eSSatish Balay 384e5c89e4eSSatish Balay ierr = PetscOptionsGetString(PETSC_NULL,"-on_error_emacs",emacsmachinename,128,&flg1);CHKERRQ(ierr); 385e5c89e4eSSatish Balay if (flg1 && !rank) {ierr = PetscPushErrorHandler(PetscEmacsClientErrorHandler,emacsmachinename);CHKERRQ(ierr)} 386e5c89e4eSSatish Balay 3876dc8fec2Sbcordonn ierr=PetscOptionsHasName(PETSC_NULL,"-z", &flgz); CHKERRQ(ierr); 3886dc8fec2Sbcordonn if(flgz){ 3896dc8fec2Sbcordonn extern int PETSC_SOCKFD; 3906dc8fec2Sbcordonn extern int PETSC_LISTENFD; 3916dc8fec2Sbcordonn extern int PETSC_LISTEN_CHECK; 392*f1384234SBarry Smith char hostname[256]; 3936dc8fec2Sbcordonn int remoteport = 9999; 3946dc8fec2Sbcordonn int listenport = 9998; 395*f1384234SBarry Smith PetscGetHostName(hostname,256); 3966dc8fec2Sbcordonn PETSC_SOCKFD = PetscOpenSocket(hostname, remoteport); 3976dc8fec2Sbcordonn PETSC_LISTEN_CHECK = 1; 3986dc8fec2Sbcordonn PETSC_LISTENFD = PetscSocketListen(hostname, listenport); 3996dc8fec2Sbcordonn PETSC_STDOUT = fdopen(PETSC_SOCKFD, "w"); 4006dc8fec2Sbcordonn fprintf(PETSC_STDOUT, "<<<start>>>"); 4016dc8fec2Sbcordonn } 4026dc8fec2Sbcordonn 403e5c89e4eSSatish Balay /* 404e5c89e4eSSatish Balay Setup profiling and logging 405e5c89e4eSSatish Balay */ 4066cf91177SBarry Smith #if defined (PETSC_USE_INFO) 4076cf91177SBarry Smith ierr = PetscOptionsHasName(PETSC_NULL,"-info",&flg1);CHKERRQ(ierr); 408e5c89e4eSSatish Balay if (flg1) { 409e5c89e4eSSatish Balay char logname[PETSC_MAX_PATH_LEN]; logname[0] = 0; 4106cf91177SBarry Smith ierr = PetscOptionsGetString(PETSC_NULL,"-info",logname,250,&flg1);CHKERRQ(ierr); 411e5c89e4eSSatish Balay if (logname[0]) { 4126cf91177SBarry Smith PetscInfoAllow(PETSC_TRUE,logname); 413e5c89e4eSSatish Balay } else { 4146cf91177SBarry Smith PetscInfoAllow(PETSC_TRUE,PETSC_NULL); 415e5c89e4eSSatish Balay } 416e5c89e4eSSatish Balay } 417865f6aa8SSatish Balay #endif 418865f6aa8SSatish Balay #if defined(PETSC_USE_LOG) 419865f6aa8SSatish Balay mname[0] = 0; 420865f6aa8SSatish Balay ierr = PetscOptionsGetString(PETSC_NULL,"-log_history",mname,PETSC_MAX_PATH_LEN,&flg1);CHKERRQ(ierr); 421865f6aa8SSatish Balay if (flg1) { 422865f6aa8SSatish Balay if (mname[0]) { 423865f6aa8SSatish Balay ierr = PetscLogOpenHistoryFile(mname,&petsc_history);CHKERRQ(ierr); 424865f6aa8SSatish Balay } else { 425865f6aa8SSatish Balay ierr = PetscLogOpenHistoryFile(0,&petsc_history);CHKERRQ(ierr); 426865f6aa8SSatish Balay } 427865f6aa8SSatish Balay } 428e5c89e4eSSatish Balay #if defined(PETSC_HAVE_MPE) 429e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-log_mpe",&flg1);CHKERRQ(ierr); 430e5c89e4eSSatish Balay if (flg1) PetscLogMPEBegin(); 431e5c89e4eSSatish Balay #endif 432e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-log_all",&flg1);CHKERRQ(ierr); 433e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-log",&flg2);CHKERRQ(ierr); 434e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-log_summary",&flg3);CHKERRQ(ierr); 435e5c89e4eSSatish Balay if (flg1) { ierr = PetscLogAllBegin();CHKERRQ(ierr); } 436e5c89e4eSSatish Balay else if (flg2 || flg3) { ierr = PetscLogBegin();CHKERRQ(ierr);} 437e5c89e4eSSatish Balay 438e5c89e4eSSatish Balay ierr = PetscOptionsGetString(PETSC_NULL,"-log_trace",mname,250,&flg1);CHKERRQ(ierr); 439e5c89e4eSSatish Balay if (flg1) { 440e5c89e4eSSatish Balay char name[PETSC_MAX_PATH_LEN],fname[PETSC_MAX_PATH_LEN]; 441e5c89e4eSSatish Balay FILE *file; 442e5c89e4eSSatish Balay if (mname[0]) { 443e5c89e4eSSatish Balay sprintf(name,"%s.%d",mname,rank); 444e5c89e4eSSatish Balay ierr = PetscFixFilename(name,fname);CHKERRQ(ierr); 445e5c89e4eSSatish Balay file = fopen(fname,"w"); 446e5c89e4eSSatish Balay if (!file) { 447e5c89e4eSSatish Balay SETERRQ1(PETSC_ERR_FILE_OPEN,"Unable to open trace file: %s",fname); 448e5c89e4eSSatish Balay } 449e5c89e4eSSatish Balay } else { 450da9f1d6bSBarry Smith file = PETSC_STDOUT; 451e5c89e4eSSatish Balay } 452e5c89e4eSSatish Balay ierr = PetscLogTraceBegin(file);CHKERRQ(ierr); 453e5c89e4eSSatish Balay } 454e5c89e4eSSatish Balay #endif 455e5c89e4eSSatish Balay 456e5c89e4eSSatish Balay /* 457e5c89e4eSSatish Balay Setup building of stack frames for all function calls 458e5c89e4eSSatish Balay */ 459e5c89e4eSSatish Balay #if defined(PETSC_USE_DEBUG) 460e5c89e4eSSatish Balay ierr = PetscStackCreate();CHKERRQ(ierr); 461e5c89e4eSSatish Balay #endif 462e5c89e4eSSatish Balay 46353acd3b1SBarry Smith ierr = PetscOptionsHasName(PETSC_NULL,"-options_gui",&PetscOptionsPublish);CHKERRQ(ierr); 464e5c89e4eSSatish Balay 465e5c89e4eSSatish Balay /* 466e5c89e4eSSatish Balay Print basic help message 467e5c89e4eSSatish Balay */ 468e5c89e4eSSatish Balay ierr = PetscOptionsHasName(PETSC_NULL,"-help",&flg1);CHKERRQ(ierr); 469e5c89e4eSSatish Balay if (flg1) { 470e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"Options for all PETSc programs:\n");CHKERRQ(ierr); 471e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -on_error_abort: cause an abort when an error is");CHKERRQ(ierr); 472e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," detected. Useful \n only when run in the debugger\n");CHKERRQ(ierr); 473e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -on_error_attach_debugger [gdb,dbx,xxgdb,ups,noxterm]\n");CHKERRQ(ierr); 474e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," start the debugger in new xterm\n");CHKERRQ(ierr); 475e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," unless noxterm is given\n");CHKERRQ(ierr); 476e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -start_in_debugger [gdb,dbx,xxgdb,ups,noxterm]\n");CHKERRQ(ierr); 477e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," start all processes in the debugger\n");CHKERRQ(ierr); 478e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -on_error_emacs <machinename>\n");CHKERRQ(ierr); 479e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," emacs jumps to error file\n");CHKERRQ(ierr); 480e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -debugger_nodes [n1,n2,..] Nodes to start in debugger\n");CHKERRQ(ierr); 481e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -debugger_pause [m] : delay (in seconds) to attach debugger\n");CHKERRQ(ierr); 482e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -stop_for_debugger : prints message on how to attach debugger manually\n");CHKERRQ(ierr); 483e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," waits the delay for you to attach\n");CHKERRQ(ierr); 484e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -display display: Location where graphics and debuggers are displayed\n");CHKERRQ(ierr); 485e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -no_signal_handler: do not trap error signals\n");CHKERRQ(ierr); 486e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -mpi_return_on_error: MPI returns error code, rather than abort on internal error\n");CHKERRQ(ierr); 487e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -fp_trap: stop on floating point exceptions\n");CHKERRQ(ierr); 488e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," note on IBM RS6000 this slows run greatly\n");CHKERRQ(ierr); 489e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -malloc_dump <optional filename>: dump list of unfreed memory at conclusion\n");CHKERRQ(ierr); 490e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -malloc: use our error checking malloc\n");CHKERRQ(ierr); 491e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -malloc no: don't use error checking malloc\n");CHKERRQ(ierr); 492e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -mallocinfo: prints total memory usage\n");CHKERRQ(ierr); 493e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -malloc_debug: enables extended checking for memory corruption\n");CHKERRQ(ierr); 494e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -options_table: dump list of options inputted\n");CHKERRQ(ierr); 495e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -options_left: dump list of unused options\n");CHKERRQ(ierr); 496e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -options_left no: don't dump list of unused options\n");CHKERRQ(ierr); 497e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -tmp tmpdir: alternative /tmp directory\n");CHKERRQ(ierr); 498e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -shared_tmp: tmp directory is shared by all processors\n");CHKERRQ(ierr); 499a8c7a070SBarry Smith ierr = (*PetscHelpPrintf)(comm," -not_shared_tmp: each processor has separate tmp directory\n");CHKERRQ(ierr); 500e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -memory_info: print memory usage at end of run\n");CHKERRQ(ierr); 501e5c89e4eSSatish Balay #if defined(PETSC_USE_LOG) 502e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -get_total_flops: total flops over all processors\n");CHKERRQ(ierr); 503e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -log[_all _summary]: logging objects and events\n");CHKERRQ(ierr); 504e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -log_trace [filename]: prints trace of all PETSc calls\n");CHKERRQ(ierr); 505e5c89e4eSSatish Balay #if defined(PETSC_HAVE_MPE) 506e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -log_mpe: Also create logfile viewable through upshot\n");CHKERRQ(ierr); 507e5c89e4eSSatish Balay #endif 5086cf91177SBarry Smith ierr = (*PetscHelpPrintf)(comm," -info <optional filename>: print informative messages about the calculations\n");CHKERRQ(ierr); 509e5c89e4eSSatish Balay #endif 510e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -v: prints PETSc version number and release date\n");CHKERRQ(ierr); 511e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -options_file <file>: reads options from file\n");CHKERRQ(ierr); 512e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm," -petsc_sleep n: sleeps n seconds before running program\n");CHKERRQ(ierr); 513e5c89e4eSSatish Balay ierr = (*PetscHelpPrintf)(comm,"-----------------------------------------------\n");CHKERRQ(ierr); 514e5c89e4eSSatish Balay } 515e5c89e4eSSatish Balay 516e5c89e4eSSatish Balay ierr = PetscOptionsGetInt(PETSC_NULL,"-petsc_sleep",&si,&flg1);CHKERRQ(ierr); 517e5c89e4eSSatish Balay if (flg1) { 518e5c89e4eSSatish Balay ierr = PetscSleep(si);CHKERRQ(ierr); 519e5c89e4eSSatish Balay } 520e5c89e4eSSatish Balay 5216cf91177SBarry Smith ierr = PetscOptionsGetString(PETSC_NULL,"-info_exclude",mname,PETSC_MAX_PATH_LEN,&flg1);CHKERRQ(ierr); 522e5c89e4eSSatish Balay ierr = PetscStrstr(mname,"null",&f);CHKERRQ(ierr); 523e5c89e4eSSatish Balay if (f) { 5246cf91177SBarry Smith ierr = PetscInfoDeactivateClass(PETSC_NULL);CHKERRQ(ierr); 525e5c89e4eSSatish Balay } 526e5c89e4eSSatish Balay 527e5c89e4eSSatish Balay PetscFunctionReturn(0); 528e5c89e4eSSatish Balay } 529e5c89e4eSSatish Balay 530e5c89e4eSSatish Balay 531