1 /* 2 Contains all error handling code for PETSc. 3 */ 4 #if !defined(__PETSCERROR_H) 5 #define __PETSCERROR_H 6 #include "petsc.h" 7 PETSC_EXTERN_CXX_BEGIN 8 9 /* 10 Defines the directory where the compiled source is located; used 11 in printing error messages. Each makefile has an entry 12 LOCDIR = thedirectory 13 and bmake/common_variables includes in CCPPFLAGS -D__SDIR__='"${LOCDIR}"' 14 which is a flag passed to the C/C++ compilers. 15 */ 16 #if !defined(__SDIR__) 17 #define __SDIR__ "unknowndirectory/" 18 #endif 19 20 /* 21 Defines the function where the compiled source is located; used 22 in printing error messages. 23 */ 24 #if !defined(__FUNCT__) 25 #define __FUNCT__ "User provided function" 26 #endif 27 28 /* 29 These are the generic error codes. These error codes are used 30 many different places in the PETSc source code. The string versions are 31 at src/sys/src/error/err.c any changes here must also be made there 32 33 */ 34 #define PETSC_ERR_MEM 55 /* unable to allocate requested memory */ 35 #define PETSC_ERR_MEM_MALLOC_0 85 /* cannot malloc zero size */ 36 #define PETSC_ERR_SUP 56 /* no support for requested operation */ 37 #define PETSC_ERR_SUP_SYS 57 /* no support for requested operation on this computer system */ 38 #define PETSC_ERR_ORDER 58 /* operation done in wrong order */ 39 #define PETSC_ERR_SIG 59 /* signal received */ 40 #define PETSC_ERR_FP 72 /* floating point exception */ 41 #define PETSC_ERR_COR 74 /* corrupted PETSc object */ 42 #define PETSC_ERR_LIB 76 /* error in library called by PETSc */ 43 #define PETSC_ERR_PLIB 77 /* PETSc library generated inconsistent data */ 44 #define PETSC_ERR_MEMC 78 /* memory corruption */ 45 #define PETSC_ERR_CONV_FAILED 82 /* iterative method (KSP or SNES) failed */ 46 #define PETSC_ERR_USER 83 /* user has not provided needed function */ 47 48 #define PETSC_ERR_ARG_SIZ 60 /* nonconforming object sizes used in operation */ 49 #define PETSC_ERR_ARG_IDN 61 /* two arguments not allowed to be the same */ 50 #define PETSC_ERR_ARG_WRONG 62 /* wrong argument (but object probably ok) */ 51 #define PETSC_ERR_ARG_CORRUPT 64 /* null or corrupted PETSc object as argument */ 52 #define PETSC_ERR_ARG_OUTOFRANGE 63 /* input argument, out of range */ 53 #define PETSC_ERR_ARG_BADPTR 68 /* invalid pointer argument */ 54 #define PETSC_ERR_ARG_NOTSAMETYPE 69 /* two args must be same object type */ 55 #define PETSC_ERR_ARG_NOTSAMECOMM 80 /* two args must be same communicators */ 56 #define PETSC_ERR_ARG_WRONGSTATE 73 /* object in argument is in wrong state, e.g. unassembled mat */ 57 #define PETSC_ERR_ARG_INCOMP 75 /* two arguments are incompatible */ 58 #define PETSC_ERR_ARG_NULL 85 /* argument is null that should not be */ 59 #define PETSC_ERR_ARG_UNKNOWN_TYPE 86 /* type name doesn't match any registered type */ 60 #define PETSC_ERR_ARG_DOMAIN 87 /* argument is not in domain of function */ 61 62 #define PETSC_ERR_FILE_OPEN 65 /* unable to open file */ 63 #define PETSC_ERR_FILE_READ 66 /* unable to read from file */ 64 #define PETSC_ERR_FILE_WRITE 67 /* unable to write to file */ 65 #define PETSC_ERR_FILE_UNEXPECTED 79 /* unexpected data in file */ 66 67 #define PETSC_ERR_MAT_LU_ZRPVT 71 /* detected a zero pivot during LU factorization */ 68 #define PETSC_ERR_MAT_CH_ZRPVT 81 /* detected a zero pivot during Cholesky factorization */ 69 70 #if defined(PETSC_USE_ERRORCHECKING) 71 72 /*MC 73 SETERRQ - Macro that is called when an error has been detected, 74 75 Not Collective 76 77 Synopsis: 78 void SETERRQ(PetscErrorCode errorcode,char *message) 79 80 81 Input Parameters: 82 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 83 - message - error message 84 85 Level: beginner 86 87 Notes: 88 Once the error handler is called the calling function is then returned from with the given error code. 89 90 See SETERRQ1(), SETERRQ2(), SETERRQ3() for versions that take arguments 91 92 93 Experienced users can set the error handler with PetscPushErrorHandler(). 94 95 Concepts: error^setting condition 96 97 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3() 98 M*/ 99 #define SETERRQ(n,s) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s);} 100 101 /*MC 102 SETERRQ1 - Macro that is called when an error has been detected, 103 104 Not Collective 105 106 Synopsis: 107 void SETERRQ1(PetscErrorCode errorcode,char *formatmessage,arg) 108 109 110 Input Parameters: 111 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 112 . message - error message in the printf format 113 - arg - argument (for example an integer, string or double) 114 115 Level: beginner 116 117 Notes: 118 Once the error handler is called the calling function is then returned from with the given error code. 119 120 Experienced users can set the error handler with PetscPushErrorHandler(). 121 122 Concepts: error^setting condition 123 124 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ(), SETERRQ2(), SETERRQ3() 125 M*/ 126 #define SETERRQ1(n,s,a1) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1);} 127 128 /*MC 129 SETERRQ2 - Macro that is called when an error has been detected, 130 131 Not Collective 132 133 Synopsis: 134 void SETERRQ2(PetscErrorCode errorcode,char *formatmessage,arg1,arg2) 135 136 137 Input Parameters: 138 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 139 . message - error message in the printf format 140 . arg1 - argument (for example an integer, string or double) 141 - arg2 - argument (for example an integer, string or double) 142 143 Level: beginner 144 145 Notes: 146 Once the error handler is called the calling function is then returned from with the given error code. 147 148 Experienced users can set the error handler with PetscPushErrorHandler(). 149 150 Concepts: error^setting condition 151 152 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ3() 153 M*/ 154 #define SETERRQ2(n,s,a1,a2) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1,a2);} 155 156 /*MC 157 SETERRQ3 - Macro that is called when an error has been detected, 158 159 Not Collective 160 161 Synopsis: 162 void SETERRQ3(PetscErrorCode errorcode,char *formatmessage,arg1,arg2,arg3) 163 164 165 Input Parameters: 166 + errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 167 . message - error message in the printf format 168 . arg1 - argument (for example an integer, string or double) 169 . arg2 - argument (for example an integer, string or double) 170 - arg3 - argument (for example an integer, string or double) 171 172 Level: beginner 173 174 Notes: 175 Once the error handler is called the calling function is then returned from with the given error code. 176 177 Experienced users can set the error handler with PetscPushErrorHandler(). 178 179 Concepts: error^setting condition 180 181 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2() 182 M*/ 183 #define SETERRQ3(n,s,a1,a2,a3) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1,a2,a3);} 184 185 #define SETERRQ4(n,s,a1,a2,a3,a4) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1,a2,a3,a4);} 186 #define SETERRQ5(n,s,a1,a2,a3,a4,a5) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1,a2,a3,a4,a5);} 187 #define SETERRQ6(n,s,a1,a2,a3,a4,a5,a6) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1,a2,a3,a4,a5,a6);} 188 #define SETERRQ7(n,s,a1,a2,a3,a4,a5,a6,a7) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s,a1,a2,a3,a4,a5,a6,a7);} 189 #define SETERRABORT(comm,n,s) {PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,1,s);MPI_Abort(comm,n);} 190 191 /*MC 192 CHKERRQ - Checks error code, if non-zero it calls the error handler and then returns 193 194 Not Collective 195 196 Synopsis: 197 void CHKERRQ(PetscErrorCode errorcode) 198 199 200 Input Parameters: 201 . errorcode - nonzero error code, see the list of standard error codes in include/petscerror.h 202 203 Level: beginner 204 205 Notes: 206 Once the error handler is called the calling function is then returned from with the given error code. 207 208 Experienced users can set the error handler with PetscPushErrorHandler(). 209 210 CHKERRQ(n) is fundamentally a macro replacement for 211 if (n) return(PetscError(...,n,...)); 212 213 Although typical usage resembles "void CHKERRQ(PetscErrorCode)" as described above, for certain uses it is 214 highly inappropriate to use it in this manner as it invokes return(PetscErrorCode). In particular, 215 it cannot be used in functions which return(void) or any other datatype. In these types of functions, 216 a more appropriate construct for using PETSc Error Handling would be 217 if (n) {PetscError(....); return(YourReturnType);} 218 219 Concepts: error^setting condition 220 221 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2() 222 M*/ 223 #define CHKERRQ(n) if (n) {return PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,0," ");} 224 225 #define CHKERRABORT(comm,n) if (n) {PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,0," ");MPI_Abort(comm,n);} 226 #define CHKERRCONTINUE(n) if (n) {PetscError(__LINE__,__FUNCT__,__FILE__,__SDIR__,n,0," ");} 227 228 /*MC 229 CHKMEMQ - Checks the memory for corruption, calls error handler if any is detected 230 231 Not Collective 232 233 Synopsis: 234 CHKMEMQ; 235 236 Level: beginner 237 238 Notes: 239 Must run with the option -malloc_debug to enable this option 240 241 Once the error handler is called the calling function is then returned from with the given error code. 242 243 By defaults prints location where memory that is corrupted was allocated. 244 245 Concepts: memory corruption 246 247 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 248 PetscMallocValidate() 249 M*/ 250 #define CHKMEMQ {PetscErrorCode _7_ierr = PetscMallocValidate(__LINE__,__FUNCT__,__FILE__,__SDIR__);CHKERRQ(_7_ierr);} 251 252 #if !defined(PETSC_SKIP_UNDERSCORE_CHKERR) 253 extern PetscErrorCode __gierr; 254 #define _ __gierr = 255 #define ___ CHKERRQ(__gierr); 256 #endif 257 258 #define PETSC_EXCEPTIONS_MAX 256 259 extern PetscErrorCode PetscErrorUncatchable[PETSC_EXCEPTIONS_MAX]; 260 extern PetscInt PetscErrorUncatchableCount; 261 extern PetscErrorCode PetscExceptions[PETSC_EXCEPTIONS_MAX]; 262 extern PetscInt PetscExceptionsCount; 263 264 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscExceptionPush(PetscErrorCode); 265 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscErrorSetCatchable(PetscErrorCode,PetscTruth); 266 EXTERN void PETSC_DLLEXPORT PetscExceptionPop(PetscErrorCode); 267 268 /*MC 269 PetscExceptionCaught - Indicates if exception zierr was caught. 270 271 Not Collective 272 273 Synopsis: 274 PetscTruth PetscExceptionCaught(PetscErrorCode xierr,PetscErrorCode zierr); 275 276 Input Parameters: 277 + xierr - error code returned from PetscExceptionTry1() 278 - zierr - error code you want it to be 279 280 Level: advanced 281 282 Notes: 283 PETSc must not be configured using the option --with-errorchecking=0 for this to work 284 285 Concepts: exceptions, exception hanlding 286 287 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 288 CHKERRQ(), PetscExceptionTry1(), PetscExceptionValue() 289 M*/ 290 static inline PetscTruth PetscExceptionCaught(PetscErrorCode xierr,PetscErrorCode zierr) { 291 PetscInt i; 292 if (xierr != zierr) return PETSC_FALSE; 293 for (i=0; i<PetscErrorUncatchableCount; i++) { 294 if (PetscErrorUncatchable[i] == zierr) { 295 return PETSC_FALSE; 296 } 297 } 298 return PETSC_TRUE; 299 } 300 301 /*MC 302 PetscExceptionValue - Indicates if the error code is one that is currently being tried 303 304 Not Collective 305 306 Synopsis: 307 PetscTruth PetscExceptionValue(PetscErrorCode xierr); 308 309 Input Parameters: 310 . xierr - error code 311 312 Level: developer 313 314 Notes: 315 PETSc must not be configured using the option --with-errorchecking=0 for this to work 316 317 Concepts: exceptions, exception hanlding 318 319 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 320 CHKERRQ(), PetscExceptionTry1(), PetscExceptionCaught() 321 M*/ 322 static inline PetscTruth PetscExceptionValue(PetscErrorCode zierr) { 323 PetscInt i; 324 for (i=0; i<PetscExceptionsCount; i++) { 325 if (PetscExceptions[i] == zierr) { 326 return PETSC_TRUE; 327 } 328 } 329 return PETSC_FALSE; 330 } 331 332 /*MC 333 PetscExceptionTry1 - Runs the routine, causing a particular error code to be treated as an exception, 334 rather than an error. That is if that error code is treated the program returns to this level, 335 but does not call the error handlers 336 337 Not Collective 338 339 Synopsis: 340 PetscExceptionTry1(PetscErrorCode routine(....),PetscErrorCode); 341 342 Level: advanced 343 344 Notes: 345 PETSc must not be configured using the option --with-errorchecking=0 for this to work 346 347 Note: In general, the outer most try on an exception is the one that will be caught (that is trys down in 348 PETSc code will not usually handle an exception that was issued above). See SNESSolve() for an example 349 of how the local try is ignored if a higher (in the stack) one is also in effect. 350 351 Concepts: exceptions, exception hanlding 352 353 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 354 CHKERRQ(), PetscExceptionCaught() 355 M*/ 356 extern PetscErrorCode PetscExceptionTmp; 357 #define PetscExceptionTry1(a,b) (PetscExceptionTmp = PetscExceptionPush(b)) ? PetscExceptionTmp : (PetscExceptionTmp = a , PetscExceptionPop(b),PetscExceptionTmp) 358 359 #else 360 #define SETERRQ(n,s) ; 361 #define SETERRQ1(n,s,a1) ; 362 #define SETERRQ2(n,s,a1,a2) ; 363 #define SETERRQ3(n,s,a1,a2,a3) ; 364 #define SETERRQ4(n,s,a1,a2,a3,a4) ; 365 #define SETERRQ5(n,s,a1,a2,a3,a4,a5) ; 366 #define SETERRQ6(n,s,a1,a2,a3,a4,a5,a6) ; 367 #define SETERRABORT(comm,n,s) ; 368 369 #define CHKERRQ(n) ; 370 #define CHKERRABORT(comm,n) ; 371 #define CHKERRCONTINUE(n) ; 372 373 #define CHKMEMQ ; 374 375 #if !defined(PETSC_SKIP_UNDERSCORE_CHKERR) 376 #define _ 377 #define ___ 378 #endif 379 380 #define PetscExceptionTry1(a,b) a 381 #endif 382 383 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscErrorPrintfInitialize(void); 384 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscErrorMessage(int,const char*[],char **); 385 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscTraceBackErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 386 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscIgnoreErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 387 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscEmacsClientErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 388 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStopErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 389 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscAbortErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 390 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscAttachDebuggerErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 391 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscError(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,...) PETSC_PRINTF_FORMAT_CHECK(7,8); 392 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPushErrorHandler(PetscErrorCode (*handler)(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*),void*); 393 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPopErrorHandler(void); 394 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscDefaultSignalHandler(int,void*); 395 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPushSignalHandler(PetscErrorCode (*)(int,void *),void*); 396 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPopSignalHandler(void); 397 398 typedef enum {PETSC_FP_TRAP_OFF=0,PETSC_FP_TRAP_ON=1} PetscFPTrap; 399 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscSetFPTrap(PetscFPTrap); 400 401 /* 402 Allows the code to build a stack frame as it runs 403 */ 404 #if defined(PETSC_USE_DEBUG) 405 406 #define PETSCSTACKSIZE 15 407 408 typedef struct { 409 const char *function[PETSCSTACKSIZE]; 410 const char *file[PETSCSTACKSIZE]; 411 const char *directory[PETSCSTACKSIZE]; 412 int line[PETSCSTACKSIZE]; 413 int currentsize; 414 } PetscStack; 415 416 extern PETSC_DLLEXPORT PetscStack *petscstack; 417 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackCopy(PetscStack*,PetscStack*); 418 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackPrint(PetscStack*,FILE* fp); 419 420 #define PetscStackActive (petscstack != 0) 421 422 423 /*MC 424 PetscFunctionBegin - First executable line of each PETSc function 425 used for error handling. 426 427 Synopsis: 428 void PetscFunctionBegin; 429 430 Usage: 431 .vb 432 int something; 433 434 PetscFunctionBegin; 435 .ve 436 437 Notes: 438 Not available in Fortran 439 440 Level: developer 441 442 .seealso: PetscFunctionReturn() 443 444 .keywords: traceback, error handling 445 M*/ 446 #define PetscFunctionBegin \ 447 {\ 448 if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) { \ 449 petscstack->function[petscstack->currentsize] = __FUNCT__; \ 450 petscstack->file[petscstack->currentsize] = __FILE__; \ 451 petscstack->directory[petscstack->currentsize] = __SDIR__; \ 452 petscstack->line[petscstack->currentsize] = __LINE__; \ 453 petscstack->currentsize++; \ 454 }} 455 456 #define PetscStackPush(n) \ 457 {if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) { \ 458 petscstack->function[petscstack->currentsize] = n; \ 459 petscstack->file[petscstack->currentsize] = "unknown"; \ 460 petscstack->directory[petscstack->currentsize] = "unknown"; \ 461 petscstack->line[petscstack->currentsize] = 0; \ 462 petscstack->currentsize++; \ 463 }} 464 465 #define PetscStackPop \ 466 {if (petscstack && petscstack->currentsize > 0) { \ 467 petscstack->currentsize--; \ 468 petscstack->function[petscstack->currentsize] = 0; \ 469 petscstack->file[petscstack->currentsize] = 0; \ 470 petscstack->directory[petscstack->currentsize] = 0; \ 471 petscstack->line[petscstack->currentsize] = 0; \ 472 }}; 473 474 /*MC 475 PetscFunctionReturn - Last executable line of each PETSc function 476 used for error handling. Replaces return() 477 478 Synopsis: 479 void PetscFunctionReturn(0); 480 481 Usage: 482 .vb 483 .... 484 PetscFunctionReturn(0); 485 } 486 .ve 487 488 Notes: 489 Not available in Fortran 490 491 Level: developer 492 493 .seealso: PetscFunctionBegin() 494 495 .keywords: traceback, error handling 496 M*/ 497 #define PetscFunctionReturn(a) \ 498 {\ 499 PetscStackPop; \ 500 return(a);} 501 502 #define PetscFunctionReturnVoid() \ 503 {\ 504 PetscStackPop; \ 505 return;} 506 507 508 #else 509 510 #define PetscFunctionBegin 511 #define PetscFunctionReturn(a) return(a) 512 #define PetscFunctionReturnVoid() return 513 #define PetscStackPop 514 #define PetscStackPush(f) 515 #define PetscStackActive 0 516 517 #endif 518 519 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackCreate(void); 520 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackView(PetscViewer); 521 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackDestroy(void); 522 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackPublish(void); 523 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackDepublish(void); 524 525 526 PETSC_EXTERN_CXX_END 527 #endif 528