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