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/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 #define CHKMEMA {PetscMallocValidate(__LINE__,__FUNCT__,__FILE__,__SDIR__);} 253 254 #if !defined(PETSC_SKIP_UNDERSCORE_CHKERR) 255 extern PetscErrorCode __gierr; 256 #define _ __gierr = 257 #define ___ CHKERRQ(__gierr); 258 #endif 259 260 #define PETSC_EXCEPTIONS_MAX 256 261 extern PetscErrorCode PetscErrorUncatchable[PETSC_EXCEPTIONS_MAX]; 262 extern PetscInt PetscErrorUncatchableCount; 263 extern PetscErrorCode PetscExceptions[PETSC_EXCEPTIONS_MAX]; 264 extern PetscInt PetscExceptionsCount; 265 266 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscExceptionPush(PetscErrorCode); 267 EXTERN void PETSC_DLLEXPORT PetscExceptionPop(PetscErrorCode); 268 269 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscErrorSetCatchable(PetscErrorCode,PetscTruth); 270 271 /*MC 272 PetscExceptionCaught - Indicates if exception zierr was caught. 273 274 Not Collective 275 276 Synopsis: 277 PetscTruth PetscExceptionCaught(PetscErrorCode xierr,PetscErrorCode zierr); 278 279 Input Parameters: 280 + xierr - error code returned from PetscExceptionTry1() 281 - zierr - error code you want it to be 282 283 Level: advanced 284 285 Notes: 286 PETSc must not be configured using the option --with-errorchecking=0 for this to work 287 288 Concepts: exceptions, exception hanlding 289 290 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 291 CHKERRQ(), PetscExceptionTry1(), PetscExceptionValue() 292 M*/ 293 PETSC_STATIC_INLINE PetscTruth PetscExceptionCaught(PetscErrorCode xierr,PetscErrorCode zierr) { 294 PetscInt i; 295 if (xierr != zierr) return PETSC_FALSE; 296 for (i=0; i<PetscErrorUncatchableCount; i++) { 297 if (PetscErrorUncatchable[i] == zierr) { 298 return PETSC_FALSE; 299 } 300 } 301 return PETSC_TRUE; 302 } 303 304 /*MC 305 PetscExceptionValue - Indicates if the error code is one that is currently being tried 306 307 Not Collective 308 309 Synopsis: 310 PetscTruth PetscExceptionValue(PetscErrorCode xierr); 311 312 Input Parameters: 313 . xierr - error code 314 315 Level: developer 316 317 Notes: 318 PETSc must not be configured using the option --with-errorchecking=0 for this to work 319 320 Concepts: exceptions, exception hanlding 321 322 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 323 CHKERRQ(), PetscExceptionTry1(), PetscExceptionCaught() 324 M*/ 325 PETSC_STATIC_INLINE PetscTruth PetscExceptionValue(PetscErrorCode zierr) { 326 PetscInt i; 327 for (i=0; i<PetscExceptionsCount; i++) { 328 if (PetscExceptions[i] == zierr) { 329 return PETSC_TRUE; 330 } 331 } 332 return PETSC_FALSE; 333 } 334 335 /*MC 336 PetscExceptionTry1 - Runs the routine, causing a particular error code to be treated as an exception, 337 rather than an error. That is if that error code is treated the program returns to this level, 338 but does not call the error handlers 339 340 Not Collective 341 342 Synopsis: 343 PetscExceptionTry1(PetscErrorCode routine(....),PetscErrorCode); 344 345 Level: advanced 346 347 Notes: 348 PETSc must not be configured using the option --with-errorchecking=0 for this to work 349 350 Note: In general, the outer most try on an exception is the one that will be caught (that is trys down in 351 PETSc code will not usually handle an exception that was issued above). See SNESSolve() for an example 352 of how the local try is ignored if a higher (in the stack) one is also in effect. 353 354 Concepts: exceptions, exception hanlding 355 356 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), 357 CHKERRQ(), PetscExceptionCaught() 358 M*/ 359 extern PetscErrorCode PetscExceptionTmp; 360 #define PetscExceptionTry1(a,b) (PetscExceptionTmp = PetscExceptionPush(b)) ? PetscExceptionTmp : (PetscExceptionTmp = a , PetscExceptionPop(b),PetscExceptionTmp) 361 362 #else 363 #define SETERRQ(n,s) ; 364 #define SETERRQ1(n,s,a1) ; 365 #define SETERRQ2(n,s,a1,a2) ; 366 #define SETERRQ3(n,s,a1,a2,a3) ; 367 #define SETERRQ4(n,s,a1,a2,a3,a4) ; 368 #define SETERRQ5(n,s,a1,a2,a3,a4,a5) ; 369 #define SETERRQ6(n,s,a1,a2,a3,a4,a5,a6) ; 370 #define SETERRABORT(comm,n,s) ; 371 372 #define CHKERRQ(n) ; 373 #define CHKERRABORT(comm,n) ; 374 #define CHKERRCONTINUE(n) ; 375 376 #define CHKMEMQ ; 377 378 #if !defined(PETSC_SKIP_UNDERSCORE_CHKERR) 379 #define _ 380 #define ___ 381 #endif 382 383 #define PetscErrorSetCatchable(a,b) 0 384 #define PetscExceptionCaught(a,b) PETSC_FALSE 385 #define PetscExceptionValue(a) PETSC_FALSE 386 #define PetscExceptionTry1(a,b) a 387 #endif 388 389 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscErrorPrintfInitialize(void); 390 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscErrorMessage(int,const char*[],char **); 391 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscTraceBackErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 392 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscIgnoreErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 393 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscEmacsClientErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 394 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStopErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 395 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscAbortErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 396 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscAttachDebuggerErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 397 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscReturnErrorHandler(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*); 398 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscError(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,...) PETSC_PRINTF_FORMAT_CHECK(7,8); 399 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPushErrorHandler(PetscErrorCode (*handler)(int,const char*,const char*,const char*,PetscErrorCode,int,const char*,void*),void*); 400 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPopErrorHandler(void); 401 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscDefaultSignalHandler(int,void*); 402 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPushSignalHandler(PetscErrorCode (*)(int,void *),void*); 403 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscPopSignalHandler(void); 404 405 typedef enum {PETSC_FP_TRAP_OFF=0,PETSC_FP_TRAP_ON=1} PetscFPTrap; 406 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscSetFPTrap(PetscFPTrap); 407 408 /* 409 Allows the code to build a stack frame as it runs 410 */ 411 #if defined(PETSC_USE_DEBUG) 412 413 #define PETSCSTACKSIZE 15 414 415 typedef struct { 416 const char *function[PETSCSTACKSIZE]; 417 const char *file[PETSCSTACKSIZE]; 418 const char *directory[PETSCSTACKSIZE]; 419 int line[PETSCSTACKSIZE]; 420 int currentsize; 421 } PetscStack; 422 423 extern PETSC_DLLEXPORT PetscStack *petscstack; 424 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackCopy(PetscStack*,PetscStack*); 425 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackPrint(PetscStack*,FILE* fp); 426 427 #define PetscStackActive (petscstack != 0) 428 429 430 /*MC 431 PetscFunctionBegin - First executable line of each PETSc function 432 used for error handling. 433 434 Synopsis: 435 void PetscFunctionBegin; 436 437 Usage: 438 .vb 439 int something; 440 441 PetscFunctionBegin; 442 .ve 443 444 Notes: 445 Not available in Fortran 446 447 Level: developer 448 449 .seealso: PetscFunctionReturn() 450 451 .keywords: traceback, error handling 452 M*/ 453 #define PetscFunctionBegin \ 454 {\ 455 if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) { \ 456 petscstack->function[petscstack->currentsize] = __FUNCT__; \ 457 petscstack->file[petscstack->currentsize] = __FILE__; \ 458 petscstack->directory[petscstack->currentsize] = __SDIR__; \ 459 petscstack->line[petscstack->currentsize] = __LINE__; \ 460 petscstack->currentsize++; \ 461 }} 462 463 #define PetscStackPush(n) \ 464 {if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) { \ 465 petscstack->function[petscstack->currentsize] = n; \ 466 petscstack->file[petscstack->currentsize] = "unknown"; \ 467 petscstack->directory[petscstack->currentsize] = "unknown"; \ 468 petscstack->line[petscstack->currentsize] = 0; \ 469 petscstack->currentsize++; \ 470 }} 471 472 #define PetscStackPop \ 473 {if (petscstack && petscstack->currentsize > 0) { \ 474 petscstack->currentsize--; \ 475 petscstack->function[petscstack->currentsize] = 0; \ 476 petscstack->file[petscstack->currentsize] = 0; \ 477 petscstack->directory[petscstack->currentsize] = 0; \ 478 petscstack->line[petscstack->currentsize] = 0; \ 479 }}; 480 481 /*MC 482 PetscFunctionReturn - Last executable line of each PETSc function 483 used for error handling. Replaces return() 484 485 Synopsis: 486 void PetscFunctionReturn(0); 487 488 Usage: 489 .vb 490 .... 491 PetscFunctionReturn(0); 492 } 493 .ve 494 495 Notes: 496 Not available in Fortran 497 498 Level: developer 499 500 .seealso: PetscFunctionBegin() 501 502 .keywords: traceback, error handling 503 M*/ 504 #define PetscFunctionReturn(a) \ 505 {\ 506 PetscStackPop; \ 507 return(a);} 508 509 #define PetscFunctionReturnVoid() \ 510 {\ 511 PetscStackPop; \ 512 return;} 513 514 515 #else 516 517 #define PetscFunctionBegin 518 #define PetscFunctionReturn(a) return(a) 519 #define PetscFunctionReturnVoid() return 520 #define PetscStackPop 521 #define PetscStackPush(f) 522 #define PetscStackActive 0 523 524 #endif 525 526 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackCreate(void); 527 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackView(PetscViewer); 528 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackDestroy(void); 529 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackPublish(void); 530 EXTERN PetscErrorCode PETSC_DLLEXPORT PetscStackDepublish(void); 531 532 533 PETSC_EXTERN_CXX_END 534 #endif 535