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 /* Linux functions CPU_SET and others don't work if sched.h is not included before 410 including pthread.h. Also, these functions are active only if either _GNU_SOURCE 411 or __USE_GNU is not set (see /usr/include/sched.h and /usr/include/features.h), hence 412 set these first. 413 */ 414 #if defined(PETSC_HAVE_PTHREADCLASSES) || defined (PETSC_HAVE_OPENMP) 415 #if defined(PETSC_HAVE_SCHED_H) 416 #ifndef _GNU_SOURCE 417 #define _GNU_SOURCE 418 #endif 419 #include <sched.h> 420 #endif 421 #include <pthread.h> 422 #endif 423 424 /* 425 This code is for managing thread local global variables. Each of Linux, Microsoft WINDOWS, OpenMP, and Apple OS X have 426 different ways to indicate this. On OS X each thread local global is accessed by using a pthread_key_t for that variable. 427 Thus we have functions for creating destroying and using the keys. Except for OS X these access functions merely directly 428 acess the thread local variable. 429 */ 430 431 #if defined(PETSC_HAVE_PTHREADCLASSES) && !defined(PETSC_PTHREAD_LOCAL) 432 typedef pthread_key_t PetscThreadKey; 433 /* Get the value associated with key */ 434 PETSC_STATIC_INLINE void* PetscThreadLocalGetValue(PetscThreadKey key) 435 { 436 return pthread_getspecific(key); 437 } 438 439 /* Set the value for key */ 440 PETSC_STATIC_INLINE void PetscThreadLocalSetValue(PetscThreadKey *key,void* value) 441 { 442 pthread_setspecific(*key,(void*)value); 443 } 444 445 /* Create pthread thread local key */ 446 PETSC_STATIC_INLINE void PetscThreadLocalRegister(PetscThreadKey *key) 447 { 448 pthread_key_create(key,NULL); 449 } 450 451 /* Delete pthread thread local key */ 452 PETSC_STATIC_INLINE void PetscThreadLocalDestroy(PetscThreadKey key) 453 { 454 pthread_key_delete(key); 455 } 456 #else 457 typedef void* PetscThreadKey; 458 PETSC_STATIC_INLINE void* PetscThreadLocalGetValue(PetscThreadKey key) 459 { 460 return key; 461 } 462 463 PETSC_STATIC_INLINE void PetscThreadLocalSetValue(PetscThreadKey *key,void* value) 464 { 465 *key = value; 466 } 467 468 PETSC_STATIC_INLINE void PetscThreadLocalRegister(PETSC_UNUSED PetscThreadKey *key) 469 { 470 } 471 472 PETSC_STATIC_INLINE void PetscThreadLocalDestroy(PETSC_UNUSED PetscThreadKey key) 473 { 474 } 475 #endif 476 477 /* 478 Allows the code to build a stack frame as it runs 479 */ 480 481 #define PETSCSTACKSIZE 64 482 483 typedef struct { 484 const char *function[PETSCSTACKSIZE]; 485 const char *file[PETSCSTACKSIZE]; 486 int line[PETSCSTACKSIZE]; 487 PetscBool petscroutine[PETSCSTACKSIZE]; 488 int currentsize; 489 int hotdepth; 490 } PetscStack; 491 492 #if defined(PETSC_HAVE_PTHREADCLASSES) 493 #if defined(PETSC_PTHREAD_LOCAL) 494 PETSC_EXTERN PETSC_PTHREAD_LOCAL PetscStack *petscstack; 495 #else 496 PETSC_EXTERN PetscThreadKey petscstack; 497 #endif 498 #elif defined(PETSC_HAVE_OPENMP) 499 PETSC_EXTERN PetscStack *petscstack; 500 #pragma omp threadprivate(petscstack) 501 #else 502 PETSC_EXTERN PetscStack *petscstack; 503 #endif 504 505 PETSC_EXTERN PetscErrorCode PetscStackCopy(PetscStack*,PetscStack*); 506 PETSC_EXTERN PetscErrorCode PetscStackPrint(PetscStack*,FILE* fp); 507 508 #if defined(PETSC_USE_DEBUG) 509 PETSC_STATIC_INLINE PetscBool PetscStackActive(void) 510 { 511 return(PetscThreadLocalGetValue(petscstack) ? PETSC_TRUE : PETSC_FALSE); 512 } 513 514 /* Stack handling is based on the following two "NoCheck" macros. These should only be called directly by other error 515 * handling macros. We record the line of the call, which may or may not be the location of the definition. But is at 516 * least more useful than "unknown" because it can distinguish multiple calls from the same function. 517 */ 518 519 #define PetscStackPushNoCheck(funct,petsc_routine,hot) \ 520 do { \ 521 PetscStack* petscstackp; \ 522 PetscStackSAWsTakeAccess(); \ 523 petscstackp = (PetscStack*)PetscThreadLocalGetValue(petscstack); \ 524 if (petscstackp && (petscstackp->currentsize < PETSCSTACKSIZE)) { \ 525 petscstackp->function[petscstackp->currentsize] = funct; \ 526 petscstackp->file[petscstackp->currentsize] = __FILE__; \ 527 petscstackp->line[petscstackp->currentsize] = __LINE__; \ 528 petscstackp->petscroutine[petscstackp->currentsize] = petsc_routine; \ 529 petscstackp->currentsize++; \ 530 } \ 531 if (petscstackp) { \ 532 petscstackp->hotdepth += (hot || petscstackp->hotdepth); \ 533 } \ 534 PetscStackSAWsGrantAccess(); \ 535 } while (0) 536 537 #define PetscStackPopNoCheck \ 538 do {PetscStack* petscstackp; \ 539 PetscStackSAWsTakeAccess(); \ 540 petscstackp = (PetscStack*)PetscThreadLocalGetValue(petscstack); \ 541 if (petscstackp && petscstackp->currentsize > 0) { \ 542 petscstackp->currentsize--; \ 543 petscstackp->function[petscstackp->currentsize] = 0; \ 544 petscstackp->file[petscstackp->currentsize] = 0; \ 545 petscstackp->line[petscstackp->currentsize] = 0; \ 546 petscstackp->petscroutine[petscstackp->currentsize] = PETSC_FALSE;\ 547 } \ 548 if (petscstackp) { \ 549 petscstackp->hotdepth = PetscMax(petscstackp->hotdepth-1,0); \ 550 } \ 551 PetscStackSAWsGrantAccess(); \ 552 } while (0) 553 554 /*MC 555 PetscFunctionBegin - First executable line of each PETSc function, used for error handling. Final 556 line of PETSc functions should be PetscFunctionReturn(0); 557 558 Synopsis: 559 #include <petscsys.h> 560 void PetscFunctionBegin; 561 562 Not Collective 563 564 Usage: 565 .vb 566 int something; 567 568 PetscFunctionBegin; 569 .ve 570 571 Notes: 572 Use PetscFunctionBeginUser for application codes. 573 574 Not available in Fortran 575 576 Level: developer 577 578 .seealso: PetscFunctionReturn(), PetscFunctionBeginHot(), PetscFunctionBeginUser() 579 580 .keywords: traceback, error handling 581 M*/ 582 #define PetscFunctionBegin do { \ 583 PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_FALSE); \ 584 PetscCheck__FUNCT__(); \ 585 PetscRegister__FUNCT__(); \ 586 } while (0) 587 588 /*MC 589 PetscFunctionBeginHot - Substitute for PetscFunctionBegin to be used in functions that are called in 590 performance-critical circumstances. Use of this function allows for lighter profiling by default. 591 592 Synopsis: 593 #include <petscsys.h> 594 void PetscFunctionBeginHot; 595 596 Not Collective 597 598 Usage: 599 .vb 600 int something; 601 602 PetscFunctionBeginHot; 603 .ve 604 605 Notes: 606 Not available in Fortran 607 608 Level: developer 609 610 .seealso: PetscFunctionBegin, PetscFunctionReturn() 611 612 .keywords: traceback, error handling 613 M*/ 614 #define PetscFunctionBeginHot do { \ 615 PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_TRUE); \ 616 PetscCheck__FUNCT__(); \ 617 PetscRegister__FUNCT__(); \ 618 } while (0) 619 620 /*MC 621 PetscFunctionBeginUser - First executable line of user provided PETSc routine 622 623 Synopsis: 624 #include <petscsys.h> 625 void PetscFunctionBeginUser; 626 627 Not Collective 628 629 Usage: 630 .vb 631 int something; 632 633 PetscFunctionBegin; 634 .ve 635 636 Notes: 637 Final line of PETSc functions should be PetscFunctionReturn(0) except for main(). 638 639 Not available in Fortran 640 641 Level: intermediate 642 643 .seealso: PetscFunctionReturn(), PetscFunctionBegin, PetscFunctionBeginHot 644 645 .keywords: traceback, error handling 646 M*/ 647 #define PetscFunctionBeginUser \ 648 do { \ 649 PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_FALSE,PETSC_FALSE); \ 650 PetscCheck__FUNCT__(); \ 651 PetscRegister__FUNCT__(); \ 652 } while (0) 653 654 655 #if defined(PETSC_SERIALIZE_FUNCTIONS) 656 #include <petsc/private/petscfptimpl.h> 657 /* 658 Registers the current function into the global function pointer to function name table 659 660 Have to fix this to handle errors but cannot return error since used in PETSC_VIEWER_DRAW_() etc 661 */ 662 #define PetscRegister__FUNCT__() do { \ 663 static PetscBool __chked = PETSC_FALSE; \ 664 if (!__chked) {\ 665 void *ptr; PetscDLSym(NULL,__FUNCT__,&ptr);\ 666 __chked = PETSC_TRUE;\ 667 }} while (0) 668 #else 669 #define PetscRegister__FUNCT__() 670 #endif 671 672 #define PetscCheck__FUNCT__() do { PetscBool _sc1,_sc2; \ 673 PetscStrcmpNoError(PETSC_FUNCTION_NAME,__FUNCT__,&_sc1);\ 674 PetscStrcmpNoError(__FUNCT__,"User provided function",&_sc2);\ 675 if (!_sc1 && !_sc2) { \ 676 printf("%s:%d: __FUNCT__=\"%s\" does not agree with %s=\"%s\"\n",__FILE__,__LINE__,__FUNCT__,PetscStringize(PETSC_FUNCTION_NAME),PETSC_FUNCTION_NAME); \ 677 } \ 678 } while (0) 679 680 #define PetscStackPush(n) \ 681 do { \ 682 PetscStackPushNoCheck(n,PETSC_FALSE,PETSC_FALSE); \ 683 CHKMEMQ; \ 684 } while (0) 685 686 #define PetscStackPop \ 687 do { \ 688 CHKMEMQ; \ 689 PetscStackPopNoCheck; \ 690 } while (0) 691 692 /*MC 693 PetscFunctionReturn - Last executable line of each PETSc function 694 used for error handling. Replaces return() 695 696 Synopsis: 697 #include <petscsys.h> 698 void PetscFunctionReturn(0); 699 700 Not Collective 701 702 Usage: 703 .vb 704 .... 705 PetscFunctionReturn(0); 706 } 707 .ve 708 709 Notes: 710 Not available in Fortran 711 712 Level: developer 713 714 .seealso: PetscFunctionBegin() 715 716 .keywords: traceback, error handling 717 M*/ 718 #define PetscFunctionReturn(a) \ 719 do { \ 720 PetscStackPopNoCheck; \ 721 return(a);} while (0) 722 723 #define PetscFunctionReturnVoid() \ 724 do { \ 725 PetscStackPopNoCheck; \ 726 return;} while (0) 727 728 #else 729 730 PETSC_STATIC_INLINE PetscBool PetscStackActive(void) {return PETSC_FALSE;} 731 #define PetscStackPushNoCheck(funct,petsc_routine,hot) do {} while (0) 732 #define PetscStackPopNoCheck do {} while (0) 733 #define PetscFunctionBegin 734 #define PetscFunctionBeginUser 735 #define PetscFunctionBeginHot 736 #define PetscFunctionReturn(a) return(a) 737 #define PetscFunctionReturnVoid() return 738 #define PetscStackPop CHKMEMQ 739 #define PetscStackPush(f) CHKMEMQ 740 741 #endif 742 743 /* 744 PetscStackCall - Calls an external library routine or user function after pushing the name of the routine on the stack. 745 746 Input Parameters: 747 + name - string that gives the name of the function being called 748 - routine - actual call to the routine, including ierr = and CHKERRQ(ierr); 749 750 Note: Often one should use PetscStackCallStandard() instead. This routine is intended for external library routines that DO NOT return error codes 751 752 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. 753 754 755 756 */ 757 #define PetscStackCall(name,routine) do { PetscStackPush(name);routine;PetscStackPop; } while(0) 758 759 /* 760 PetscStackCallStandard - Calls an external library routine after pushing the name of the routine on the stack. 761 762 Input Parameters: 763 + func- name of the routine 764 - args - arguments to the routine surrounded by () 765 766 Notes: This is intended for external package routines that return error codes. Use PetscStackCall() for those that do not. 767 768 Developer Note: this is so that when an external packge routine results in a crash or corrupts memory, they get blamed instead of PETSc. 769 770 */ 771 #define PetscStackCallStandard(func,args) do { \ 772 PetscStackPush(#func);ierr = func args;PetscStackPop; if (ierr) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error in %s()",#func); \ 773 } while (0) 774 775 PETSC_EXTERN PetscErrorCode PetscStackCreate(void); 776 PETSC_EXTERN PetscErrorCode PetscStackView(FILE*); 777 PETSC_EXTERN PetscErrorCode PetscStackDestroy(void); 778 779 #endif 780