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