1 /* 2 Contains all error handling interfaces for PETSc. 3 */ 4 #if !defined(__PETSCERROR_H) 5 #define __PETSCERROR_H 6 7 /* 8 Defines the function where the compiled source is located; used 9 in printing error messages. This is defined here in case the user 10 does not declare it. 11 */ 12 #if !defined(__FUNCT__) 13 #define __FUNCT__ "User provided function" 14 #endif 15 16 /* 17 These are the generic error codes. These error codes are used 18 many different places in the PETSc source code. The string versions are 19 at src/sys/error/err.c any changes here must also be made there 20 These are also define in include/petsc/finclude/petscerror.h any CHANGES here 21 must be also made there. 22 23 */ 24 #define PETSC_ERR_MIN_VALUE 54 /* should always be one less then the smallest value */ 25 26 #define PETSC_ERR_MEM 55 /* unable to allocate requested memory */ 27 #define PETSC_ERR_SUP 56 /* no support for requested operation */ 28 #define PETSC_ERR_SUP_SYS 57 /* no support for requested operation on this computer system */ 29 #define PETSC_ERR_ORDER 58 /* operation done in wrong order */ 30 #define PETSC_ERR_SIG 59 /* signal received */ 31 #define PETSC_ERR_FP 72 /* floating point exception */ 32 #define PETSC_ERR_COR 74 /* corrupted PETSc object */ 33 #define PETSC_ERR_LIB 76 /* error in library called by PETSc */ 34 #define PETSC_ERR_PLIB 77 /* PETSc library generated inconsistent data */ 35 #define PETSC_ERR_MEMC 78 /* memory corruption */ 36 #define PETSC_ERR_CONV_FAILED 82 /* iterative method (KSP or SNES) failed */ 37 #define PETSC_ERR_USER 83 /* user has not provided needed function */ 38 #define PETSC_ERR_SYS 88 /* error in system call */ 39 #define PETSC_ERR_POINTER 70 /* pointer does not point to valid address */ 40 41 #define PETSC_ERR_ARG_SIZ 60 /* nonconforming object sizes used in operation */ 42 #define PETSC_ERR_ARG_IDN 61 /* two arguments not allowed to be the same */ 43 #define PETSC_ERR_ARG_WRONG 62 /* wrong argument (but object probably ok) */ 44 #define PETSC_ERR_ARG_CORRUPT 64 /* null or corrupted PETSc object as argument */ 45 #define PETSC_ERR_ARG_OUTOFRANGE 63 /* input argument, out of range */ 46 #define PETSC_ERR_ARG_BADPTR 68 /* invalid pointer argument */ 47 #define PETSC_ERR_ARG_NOTSAMETYPE 69 /* two args must be same object type */ 48 #define PETSC_ERR_ARG_NOTSAMECOMM 80 /* two args must be same communicators */ 49 #define PETSC_ERR_ARG_WRONGSTATE 73 /* object in argument is in wrong state, e.g. unassembled mat */ 50 #define PETSC_ERR_ARG_TYPENOTSET 89 /* the type of the object has not yet been set */ 51 #define PETSC_ERR_ARG_INCOMP 75 /* two arguments are incompatible */ 52 #define PETSC_ERR_ARG_NULL 85 /* argument is null that should not be */ 53 #define PETSC_ERR_ARG_UNKNOWN_TYPE 86 /* type name doesn't match any registered type */ 54 55 #define PETSC_ERR_FILE_OPEN 65 /* unable to open file */ 56 #define PETSC_ERR_FILE_READ 66 /* unable to read from file */ 57 #define PETSC_ERR_FILE_WRITE 67 /* unable to write to file */ 58 #define PETSC_ERR_FILE_UNEXPECTED 79 /* unexpected data in file */ 59 60 #define PETSC_ERR_MAT_LU_ZRPVT 71 /* detected a zero pivot during LU factorization */ 61 #define PETSC_ERR_MAT_CH_ZRPVT 81 /* detected a zero pivot during Cholesky factorization */ 62 63 #define PETSC_ERR_INT_OVERFLOW 84 64 65 #define PETSC_ERR_FLOP_COUNT 90 66 #define PETSC_ERR_NOT_CONVERGED 91 /* solver did not converge */ 67 #define PETSC_ERR_MISSING_FACTOR 92 /* MatGetFactor() failed */ 68 #define PETSC_ERR_MAX_VALUE 93 /* this is always the one more than the largest error code */ 69 70 #define PetscStringizeArg(a) #a 71 #define PetscStringize(a) PetscStringizeArg(a) 72 73 #if defined(PETSC_USE_ERRORCHECKING) 74 75 /*MC 76 SETERRQ - Macro to be called when an error has been detected, 77 78 Synopsis: 79 #include <petscsys.h> 80 PetscErrorCode SETERRQ(MPI_Comm comm,PetscErrorCode errorcode,char *message) 81 82 Not Collective 83 84 Input Parameters: 85 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 86 - message - error message 87 88 Level: beginner 89 90 Notes: 91 Once the error handler is called the calling function is then returned from with the given error code. 92 93 See SETERRQ1(), SETERRQ2(), SETERRQ3() for versions that take arguments 94 95 In Fortran MPI_Abort() is always called 96 97 Experienced users can set the error handler with PetscPushErrorHandler(). 98 99 Concepts: error^setting condition 100 101 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3() 102 M*/ 103 #define SETERRQ(comm,n,s) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s) 104 105 /*MC 106 SETERRQ1 - Macro that is called when an error has been detected, 107 108 Synopsis: 109 #include <petscsys.h> 110 PetscErrorCode SETERRQ1(MPI_Comm comm,PetscErrorCode errorcode,char *formatmessage,arg) 111 112 Not Collective 113 114 Input Parameters: 115 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 116 . message - error message in the printf format 117 - arg - argument (for example an integer, string or double) 118 119 Level: beginner 120 121 Notes: 122 Once the error handler is called the calling function is then returned from with the given error code. 123 124 Experienced users can set the error handler with PetscPushErrorHandler(). 125 126 Concepts: error^setting condition 127 128 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ(), SETERRQ2(), SETERRQ3() 129 M*/ 130 #define SETERRQ1(comm,n,s,a1) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1) 131 132 /*MC 133 SETERRQ2 - Macro that is called when an error has been detected, 134 135 Synopsis: 136 #include <petscsys.h> 137 PetscErrorCode SETERRQ2(PetscErrorCode errorcode,char *formatmessage,arg1,arg2) 138 139 Not Collective 140 141 Input Parameters: 142 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 143 . message - error message in the printf format 144 . arg1 - argument (for example an integer, string or double) 145 - arg2 - argument (for example an integer, string or double) 146 147 Level: beginner 148 149 Notes: 150 Once the error handler is called the calling function is then returned from with the given error code. 151 152 Experienced users can set the error handler with PetscPushErrorHandler(). 153 154 Concepts: error^setting condition 155 156 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ3() 157 M*/ 158 #define SETERRQ2(comm,n,s,a1,a2) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2) 159 160 /*MC 161 SETERRQ3 - Macro that is called when an error has been detected, 162 163 Synopsis: 164 #include <petscsys.h> 165 PetscErrorCode SETERRQ3(PetscErrorCode errorcode,char *formatmessage,arg1,arg2,arg3) 166 167 Not Collective 168 169 Input Parameters: 170 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 171 . message - error message in the printf format 172 . arg1 - argument (for example an integer, string or double) 173 . arg2 - argument (for example an integer, string or double) 174 - arg3 - argument (for example an integer, string or double) 175 176 Level: beginner 177 178 Notes: 179 Once the error handler is called the calling function is then returned from with the given error code. 180 181 There are also versions for 4, 5, 6 and 7 arguments. 182 183 Experienced users can set the error handler with PetscPushErrorHandler(). 184 185 Concepts: error^setting condition 186 187 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2() 188 M*/ 189 #define SETERRQ3(comm,n,s,a1,a2,a3) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2,a3) 190 191 #define SETERRQ4(comm,n,s,a1,a2,a3,a4) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4) 192 #define SETERRQ5(comm,n,s,a1,a2,a3,a4,a5) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5) 193 #define SETERRQ6(comm,n,s,a1,a2,a3,a4,a5,a6) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6) 194 #define SETERRQ7(comm,n,s,a1,a2,a3,a4,a5,a6,a7) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6,a7) 195 #define SETERRQ8(comm,n,s,a1,a2,a3,a4,a5,a6,a7,a8) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6,a7,a8) 196 #define SETERRABORT(comm,n,s) do {PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_INITIAL,s);MPI_Abort(comm,n);} while (0) 197 198 /*MC 199 CHKERRQ - Checks error code, if non-zero it calls the error handler and then returns 200 201 Synopsis: 202 #include <petscsys.h> 203 PetscErrorCode CHKERRQ(PetscErrorCode errorcode) 204 205 Not Collective 206 207 Input Parameters: 208 . errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 209 210 Level: beginner 211 212 Notes: 213 Once the error handler is called the calling function is then returned from with the given error code. 214 215 Experienced users can set the error handler with PetscPushErrorHandler(). 216 217 CHKERRQ(n) is fundamentally a macro replacement for 218 if (n) return(PetscError(...,n,...)); 219 220 Although typical usage resembles "void CHKERRQ(PetscErrorCode)" as described above, for certain uses it is 221 highly inappropriate to use it in this manner as it invokes return(PetscErrorCode). In particular, 222 it cannot be used in functions which return(void) or any other datatype. In these types of functions, 223 you can use CHKERRV() which returns without an error code (bad idea since the error is ignored or 224 if (n) {PetscError(....); return(YourReturnType);} 225 where you may pass back a NULL to indicate an error. You can also call CHKERRABORT(comm,n) to have 226 MPI_Abort() returned immediately. 227 228 In Fortran MPI_Abort() is always called 229 230 Concepts: error^setting condition 231 232 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2() 233 M*/ 234 #define CHKERRQ(n) do {if (PetscUnlikely(n)) return PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_REPEAT," ");} while (0) 235 236 #define CHKERRV(n) do {if (PetscUnlikely(n)) {n = PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_REPEAT," ");return;}} while(0) 237 #define CHKERRABORT(comm,n) do {if (PetscUnlikely(n)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_REPEAT," ");MPI_Abort(comm,n);}} while (0) 238 #define CHKERRCONTINUE(n) do {if (PetscUnlikely(n)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_REPEAT," ");}} while (0) 239 240 #ifdef PETSC_CLANGUAGE_CXX 241 242 /*MC 243 CHKERRXX - Checks error code, if non-zero it calls the C++ error handler which throws an exception 244 245 Synopsis: 246 #include <petscsys.h> 247 void CHKERRXX(PetscErrorCode errorcode) 248 249 Not Collective 250 251 Input Parameters: 252 . errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 253 254 Level: beginner 255 256 Notes: 257 Once the error handler throws a ??? exception. 258 259 You can use CHKERRV() which returns without an error code (bad idea since the error is ignored) 260 or CHKERRABORT(comm,n) to have MPI_Abort() returned immediately. 261 262 Concepts: error^setting condition 263 264 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKERRQ(), CHKMEMQ 265 M*/ 266 #define CHKERRXX(n) do {if (PetscUnlikely(n)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,n,PETSC_ERROR_IN_CXX,0);}} while(0) 267 268 #endif 269 270 /*MC 271 CHKMEMQ - Checks the memory for corruption, calls error handler if any is detected 272 273 Synopsis: 274 #include <petscsys.h> 275 CHKMEMQ; 276 277 Not Collective 278 279 Level: beginner 280 281 Notes: 282 We highly recommend using valgrind http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind for finding memory problems. This is useful 283 on systems that do not have valgrind, but much much less useful. 284 285 Must run with the option -malloc_debug to enable this option 286 287 Once the error handler is called the calling function is then returned from with the given error code. 288 289 By defaults prints location where memory that is corrupted was allocated. 290 291 Use CHKMEMA for functions that return void 292 293 Concepts: memory corruption 294 295 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 296 PetscMallocValidate() 297 M*/ 298 #define CHKMEMQ do {PetscErrorCode _7_ierr = PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__);CHKERRQ(_7_ierr);} while(0) 299 300 #define CHKMEMA PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__) 301 302 #else /* PETSC_USE_ERRORCHECKING */ 303 304 /* 305 These are defined to be empty for when error checking is turned off, with ./configure --with-errorchecking=0 306 */ 307 308 #define SETERRQ(c,n,s) 309 #define SETERRQ1(c,n,s,a1) 310 #define SETERRQ2(c,n,s,a1,a2) 311 #define SETERRQ3(c,n,s,a1,a2,a3) 312 #define SETERRQ4(c,n,s,a1,a2,a3,a4) 313 #define SETERRQ5(c,n,s,a1,a2,a3,a4,a5) 314 #define SETERRQ6(c,n,s,a1,a2,a3,a4,a5,a6) 315 #define SETERRQ7(c,n,s,a1,a2,a3,a4,a5,a6,a7) 316 #define SETERRQ8(c,n,s,a1,a2,a3,a4,a5,a6,a7,a8) 317 #define SETERRABORT(comm,n,s) 318 319 #define CHKERRQ(n) ; 320 #define CHKERRABORT(comm,n) ; 321 #define CHKERRCONTINUE(n) ; 322 #define CHKMEMQ ; 323 324 #ifdef PETSC_CLANGUAGE_CXX 325 #define CHKERRXX(n) ; 326 #endif 327 328 #endif /* PETSC_USE_ERRORCHECKING */ 329 330 /*E 331 PetscErrorType - passed to the PETSc error handling routines indicating if this is the first or a later call to the error handlers 332 333 Level: advanced 334 335 PETSC_ERROR_IN_CXX indicates the error was detected in C++ and an exception should be generated 336 337 Developer Notes: This is currently used to decide when to print the detailed information about the run in PetscTraceBackErrorHandler() 338 339 .seealso: PetscError(), SETERRXX() 340 E*/ 341 typedef enum {PETSC_ERROR_INITIAL=0,PETSC_ERROR_REPEAT=1,PETSC_ERROR_IN_CXX = 2} PetscErrorType; 342 343 #if defined(__clang_analyzer__) 344 __attribute__((analyzer_noreturn)) 345 #endif 346 PETSC_EXTERN PetscErrorCode PetscError(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,...); 347 348 PETSC_EXTERN PetscErrorCode PetscErrorPrintfInitialize(void); 349 PETSC_EXTERN PetscErrorCode PetscErrorMessage(int,const char*[],char **); 350 PETSC_EXTERN PetscErrorCode PetscTraceBackErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 351 PETSC_EXTERN PetscErrorCode PetscIgnoreErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 352 PETSC_EXTERN PetscErrorCode PetscEmacsClientErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 353 PETSC_EXTERN PetscErrorCode PetscMPIAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 354 PETSC_EXTERN PetscErrorCode PetscAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 355 PETSC_EXTERN PetscErrorCode PetscAttachDebuggerErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 356 PETSC_EXTERN PetscErrorCode PetscReturnErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*); 357 PETSC_EXTERN PetscErrorCode PetscPushErrorHandler(PetscErrorCode (*handler)(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*),void*); 358 PETSC_EXTERN PetscErrorCode PetscPopErrorHandler(void); 359 PETSC_EXTERN PetscErrorCode PetscSignalHandlerDefault(int,void*); 360 PETSC_EXTERN PetscErrorCode PetscPushSignalHandler(PetscErrorCode (*)(int,void *),void*); 361 PETSC_EXTERN PetscErrorCode PetscPopSignalHandler(void); 362 PETSC_EXTERN PetscErrorCode PetscCheckPointerSetIntensity(PetscInt); 363 364 /*MC 365 PetscErrorPrintf - Prints error messages. 366 367 Synopsis: 368 #include <petscsys.h> 369 PetscErrorCode (*PetscErrorPrintf)(const char format[],...); 370 371 Not Collective 372 373 Input Parameters: 374 . format - the usual printf() format string 375 376 Options Database Keys: 377 + -error_output_stdout - cause error messages to be printed to stdout instead of the (default) stderr 378 - -error_output_none - to turn off all printing of error messages (does not change the way the error is handled.) 379 380 Notes: Use 381 $ PetscErrorPrintf = PetscErrorPrintfNone; to turn off all printing of error messages (does not change the way the 382 $ error is handled.) and 383 $ PetscErrorPrintf = PetscErrorPrintfDefault; to turn it back on or you can use your own function 384 385 Use 386 PETSC_STDERR = FILE* obtained from a file open etc. to have stderr printed to the file. 387 PETSC_STDOUT = FILE* obtained from a file open etc. to have stdout printed to the file. 388 389 Use 390 PetscPushErrorHandler() to provide your own error handler that determines what kind of messages to print 391 392 Level: developer 393 394 Fortran Note: 395 This routine is not supported in Fortran. 396 397 Concepts: error messages^printing 398 Concepts: printing^error messages 399 400 .seealso: PetscFPrintf(), PetscSynchronizedPrintf(), PetscHelpPrintf(), PetscPrintf(), PetscErrorHandlerPush(), PetscVFPrintf(), PetscHelpPrintf() 401 M*/ 402 PETSC_EXTERN PetscErrorCode (*PetscErrorPrintf)(const char[],...); 403 404 typedef enum {PETSC_FP_TRAP_OFF=0,PETSC_FP_TRAP_ON=1} PetscFPTrap; 405 PETSC_EXTERN PetscErrorCode PetscSetFPTrap(PetscFPTrap); 406 PETSC_EXTERN PetscErrorCode PetscFPTrapPush(PetscFPTrap); 407 PETSC_EXTERN PetscErrorCode PetscFPTrapPop(void); 408 409 /* 410 Allows the code to build a stack frame as it runs 411 */ 412 413 #define PETSCSTACKSIZE 64 414 415 typedef struct { 416 const char *function[PETSCSTACKSIZE]; 417 const char *file[PETSCSTACKSIZE]; 418 int line[PETSCSTACKSIZE]; 419 PetscBool petscroutine[PETSCSTACKSIZE]; 420 int currentsize; 421 int hotdepth; 422 } PetscStack; 423 424 PETSC_EXTERN PetscStack *petscstack; 425 426 PetscErrorCode PetscStackCopy(PetscStack*,PetscStack*); 427 PetscErrorCode PetscStackPrint(PetscStack *,FILE*); 428 #if defined(PETSC_USE_DEBUG) 429 PETSC_STATIC_INLINE PetscBool PetscStackActive(void) 430 { 431 return(petscstack ? PETSC_TRUE : PETSC_FALSE); 432 } 433 434 /* Stack handling is based on the following two "NoCheck" macros. These should only be called directly by other error 435 * handling macros. We record the line of the call, which may or may not be the location of the definition. But is at 436 * least more useful than "unknown" because it can distinguish multiple calls from the same function. 437 */ 438 439 #define PetscStackPushNoCheck(funct,petsc_routine,hot) \ 440 do { \ 441 PetscStackSAWsTakeAccess(); \ 442 if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) { \ 443 petscstack->function[petscstack->currentsize] = funct; \ 444 petscstack->file[petscstack->currentsize] = __FILE__; \ 445 petscstack->line[petscstack->currentsize] = __LINE__; \ 446 petscstack->petscroutine[petscstack->currentsize] = petsc_routine; \ 447 petscstack->currentsize++; \ 448 } \ 449 if (petscstack) { \ 450 petscstack->hotdepth += (hot || petscstack->hotdepth); \ 451 } \ 452 PetscStackSAWsGrantAccess(); \ 453 } while (0) 454 455 #define PetscStackPopNoCheck \ 456 do { \ 457 PetscStackSAWsTakeAccess(); \ 458 if (petscstack && petscstack->currentsize > 0) { \ 459 petscstack->currentsize--; \ 460 petscstack->function[petscstack->currentsize] = 0; \ 461 petscstack->file[petscstack->currentsize] = 0; \ 462 petscstack->line[petscstack->currentsize] = 0; \ 463 petscstack->petscroutine[petscstack->currentsize] = PETSC_FALSE;\ 464 } \ 465 if (petscstack) { \ 466 petscstack->hotdepth = PetscMax(petscstack->hotdepth-1,0); \ 467 } \ 468 PetscStackSAWsGrantAccess(); \ 469 } while (0) 470 471 /*MC 472 PetscFunctionBegin - First executable line of each PETSc function, used for error handling. Final 473 line of PETSc functions should be PetscFunctionReturn(0); 474 475 Synopsis: 476 #include <petscsys.h> 477 void PetscFunctionBegin; 478 479 Not Collective 480 481 Usage: 482 .vb 483 int something; 484 485 PetscFunctionBegin; 486 .ve 487 488 Notes: 489 Use PetscFunctionBeginUser for application codes. 490 491 Not available in Fortran 492 493 Level: developer 494 495 .seealso: PetscFunctionReturn(), PetscFunctionBeginHot(), PetscFunctionBeginUser() 496 497 .keywords: traceback, error handling 498 M*/ 499 #define PetscFunctionBegin do { \ 500 PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_FALSE); \ 501 PetscCheck__FUNCT__(); \ 502 PetscRegister__FUNCT__(); \ 503 } while (0) 504 505 /*MC 506 PetscFunctionBeginHot - Substitute for PetscFunctionBegin to be used in functions that are called in 507 performance-critical circumstances. Use of this function allows for lighter profiling by default. 508 509 Synopsis: 510 #include <petscsys.h> 511 void PetscFunctionBeginHot; 512 513 Not Collective 514 515 Usage: 516 .vb 517 int something; 518 519 PetscFunctionBeginHot; 520 .ve 521 522 Notes: 523 Not available in Fortran 524 525 Level: developer 526 527 .seealso: PetscFunctionBegin, PetscFunctionReturn() 528 529 .keywords: traceback, error handling 530 M*/ 531 #define PetscFunctionBeginHot do { \ 532 PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_TRUE); \ 533 PetscCheck__FUNCT__(); \ 534 PetscRegister__FUNCT__(); \ 535 } while (0) 536 537 /*MC 538 PetscFunctionBeginUser - First executable line of user provided PETSc routine 539 540 Synopsis: 541 #include <petscsys.h> 542 void PetscFunctionBeginUser; 543 544 Not Collective 545 546 Usage: 547 .vb 548 int something; 549 550 PetscFunctionBegin; 551 .ve 552 553 Notes: 554 Final line of PETSc functions should be PetscFunctionReturn(0) except for main(). 555 556 Not available in Fortran 557 558 Level: intermediate 559 560 .seealso: PetscFunctionReturn(), PetscFunctionBegin, PetscFunctionBeginHot 561 562 .keywords: traceback, error handling 563 M*/ 564 #define PetscFunctionBeginUser \ 565 do { \ 566 PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_FALSE,PETSC_FALSE); \ 567 PetscCheck__FUNCT__(); \ 568 PetscRegister__FUNCT__(); \ 569 } while (0) 570 571 572 #if defined(PETSC_SERIALIZE_FUNCTIONS) 573 #include <petsc/private/petscfptimpl.h> 574 /* 575 Registers the current function into the global function pointer to function name table 576 577 Have to fix this to handle errors but cannot return error since used in PETSC_VIEWER_DRAW_() etc 578 */ 579 #define PetscRegister__FUNCT__() do { \ 580 static PetscBool __chked = PETSC_FALSE; \ 581 if (!__chked) {\ 582 void *ptr; PetscDLSym(NULL,__FUNCT__,&ptr);\ 583 __chked = PETSC_TRUE;\ 584 }} while (0) 585 #else 586 #define PetscRegister__FUNCT__() 587 #endif 588 589 #define PetscCheck__FUNCT__() do { PetscBool _sc1,_sc2; \ 590 PetscStrcmpNoError(PETSC_FUNCTION_NAME,__FUNCT__,&_sc1);\ 591 PetscStrcmpNoError(__FUNCT__,"User provided function",&_sc2);\ 592 if (!_sc1 && !_sc2) { \ 593 printf("%s:%d: __FUNCT__=\"%s\" does not agree with %s=\"%s\"\n",__FILE__,__LINE__,__FUNCT__,PetscStringize(PETSC_FUNCTION_NAME),PETSC_FUNCTION_NAME); \ 594 } \ 595 } while (0) 596 597 #define PetscStackPush(n) \ 598 do { \ 599 PetscStackPushNoCheck(n,PETSC_FALSE,PETSC_FALSE); \ 600 CHKMEMQ; \ 601 } while (0) 602 603 #define PetscStackPop \ 604 do { \ 605 CHKMEMQ; \ 606 PetscStackPopNoCheck; \ 607 } while (0) 608 609 /*MC 610 PetscFunctionReturn - Last executable line of each PETSc function 611 used for error handling. Replaces return() 612 613 Synopsis: 614 #include <petscsys.h> 615 void PetscFunctionReturn(0); 616 617 Not Collective 618 619 Usage: 620 .vb 621 .... 622 PetscFunctionReturn(0); 623 } 624 .ve 625 626 Notes: 627 Not available in Fortran 628 629 Level: developer 630 631 .seealso: PetscFunctionBegin() 632 633 .keywords: traceback, error handling 634 M*/ 635 #define PetscFunctionReturn(a) \ 636 do { \ 637 PetscStackPopNoCheck; \ 638 return(a);} while (0) 639 640 #define PetscFunctionReturnVoid() \ 641 do { \ 642 PetscStackPopNoCheck; \ 643 return;} while (0) 644 645 #else 646 647 PETSC_STATIC_INLINE PetscBool PetscStackActive(void) {return PETSC_FALSE;} 648 #define PetscStackPushNoCheck(funct,petsc_routine,hot) do {} while (0) 649 #define PetscStackPopNoCheck do {} while (0) 650 #define PetscFunctionBegin 651 #define PetscFunctionBeginUser 652 #define PetscFunctionBeginHot 653 #define PetscFunctionReturn(a) return(a) 654 #define PetscFunctionReturnVoid() return 655 #define PetscStackPop CHKMEMQ 656 #define PetscStackPush(f) CHKMEMQ 657 658 #endif 659 660 /* 661 PetscStackCall - Calls an external library routine or user function after pushing the name of the routine on the stack. 662 663 Input Parameters: 664 + name - string that gives the name of the function being called 665 - routine - actual call to the routine, including ierr = and CHKERRQ(ierr); 666 667 Note: Often one should use PetscStackCallStandard() instead. This routine is intended for external library routines that DO NOT return error codes 668 669 Developer Note: this is so that when a user or external library routine results in a crash or corrupts memory, they get blamed instead of PETSc. 670 671 672 673 */ 674 #define PetscStackCall(name,routine) do { PetscStackPush(name);routine;PetscStackPop; } while(0) 675 676 /* 677 PetscStackCallStandard - Calls an external library routine after pushing the name of the routine on the stack. 678 679 Input Parameters: 680 + func- name of the routine 681 - args - arguments to the routine surrounded by () 682 683 Notes: This is intended for external package routines that return error codes. Use PetscStackCall() for those that do not. 684 685 Developer Note: this is so that when an external packge routine results in a crash or corrupts memory, they get blamed instead of PETSc. 686 687 */ 688 #define PetscStackCallStandard(func,args) do { \ 689 PetscStackPush(#func);ierr = func args;PetscStackPop; if (ierr) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error in %s()",#func); \ 690 } while (0) 691 692 PETSC_EXTERN PetscErrorCode PetscStackCreate(void); 693 PETSC_EXTERN PetscErrorCode PetscStackView(FILE*); 694 PETSC_EXTERN PetscErrorCode PetscStackDestroy(void); 695 696 #endif 697