xref: /petsc/include/petscerror.h (revision cd2571c5d4122701d9c8a8814fb939e6fb00b2ad)
1 /*
2     Contains all error handling interfaces for PETSc.
3 */
4 #if !defined(PETSCERROR_H)
5 #define PETSCERROR_H
6 
7 /*
8      These are the generic error codes. These error codes are used
9      many different places in the PETSc source code. The string versions are
10      at src/sys/error/err.c any changes here must also be made there
11      These are also define in src/sys/f90-mod/petscerror.h any CHANGES here
12      must be also made there.
13 
14 */
15 #define PETSC_ERR_MIN_VALUE        54   /* should always be one less then the smallest value */
16 
17 #define PETSC_ERR_MEM              55   /* unable to allocate requested memory */
18 #define PETSC_ERR_SUP              56   /* no support for requested operation */
19 #define PETSC_ERR_SUP_SYS          57   /* no support for requested operation on this computer system */
20 #define PETSC_ERR_ORDER            58   /* operation done in wrong order */
21 #define PETSC_ERR_SIG              59   /* signal received */
22 #define PETSC_ERR_FP               72   /* floating point exception */
23 #define PETSC_ERR_COR              74   /* corrupted PETSc object */
24 #define PETSC_ERR_LIB              76   /* error in library called by PETSc */
25 #define PETSC_ERR_PLIB             77   /* PETSc library generated inconsistent data */
26 #define PETSC_ERR_MEMC             78   /* memory corruption */
27 #define PETSC_ERR_CONV_FAILED      82   /* iterative method (KSP or SNES) failed */
28 #define PETSC_ERR_USER             83   /* user has not provided needed function */
29 #define PETSC_ERR_SYS              88   /* error in system call */
30 #define PETSC_ERR_POINTER          70   /* pointer does not point to valid address */
31 #define PETSC_ERR_MPI_LIB_INCOMP   87   /* MPI library at runtime is not compatible with MPI user compiled with */
32 
33 #define PETSC_ERR_ARG_SIZ          60   /* nonconforming object sizes used in operation */
34 #define PETSC_ERR_ARG_IDN          61   /* two arguments not allowed to be the same */
35 #define PETSC_ERR_ARG_WRONG        62   /* wrong argument (but object probably ok) */
36 #define PETSC_ERR_ARG_CORRUPT      64   /* null or corrupted PETSc object as argument */
37 #define PETSC_ERR_ARG_OUTOFRANGE   63   /* input argument, out of range */
38 #define PETSC_ERR_ARG_BADPTR       68   /* invalid pointer argument */
39 #define PETSC_ERR_ARG_NOTSAMETYPE  69   /* two args must be same object type */
40 #define PETSC_ERR_ARG_NOTSAMECOMM  80   /* two args must be same communicators */
41 #define PETSC_ERR_ARG_WRONGSTATE   73   /* object in argument is in wrong state, e.g. unassembled mat */
42 #define PETSC_ERR_ARG_TYPENOTSET   89   /* the type of the object has not yet been set */
43 #define PETSC_ERR_ARG_INCOMP       75   /* two arguments are incompatible */
44 #define PETSC_ERR_ARG_NULL         85   /* argument is null that should not be */
45 #define PETSC_ERR_ARG_UNKNOWN_TYPE 86   /* type name doesn't match any registered type */
46 
47 #define PETSC_ERR_FILE_OPEN        65   /* unable to open file */
48 #define PETSC_ERR_FILE_READ        66   /* unable to read from file */
49 #define PETSC_ERR_FILE_WRITE       67   /* unable to write to file */
50 #define PETSC_ERR_FILE_UNEXPECTED  79   /* unexpected data in file */
51 
52 #define PETSC_ERR_MAT_LU_ZRPVT     71   /* detected a zero pivot during LU factorization */
53 #define PETSC_ERR_MAT_CH_ZRPVT     81   /* detected a zero pivot during Cholesky factorization */
54 
55 #define PETSC_ERR_INT_OVERFLOW     84
56 
57 #define PETSC_ERR_FLOP_COUNT       90
58 #define PETSC_ERR_NOT_CONVERGED    91  /* solver did not converge */
59 #define PETSC_ERR_MISSING_FACTOR   92  /* MatGetFactor() failed */
60 #define PETSC_ERR_OPT_OVERWRITE    93  /* attempted to over write options which should not be changed */
61 #define PETSC_ERR_WRONG_MPI_SIZE   94  /* example/application run with number of MPI ranks it does not support */
62 #define PETSC_ERR_USER_INPUT       95  /* missing or incorrect user input */
63 #define PETSC_ERR_GPU_RESOURCE     96  /* unable to load a GPU resource, for example cuBLAS */
64 #define PETSC_ERR_GPU              97  /* An error from a GPU call, this may be due to lack of resources on the GPU or a true error in the call */
65 #define PETSC_ERR_MPI              98  /* general MPI error */
66 #define PETSC_ERR_MAX_VALUE        99  /* this is always the one more than the largest error code */
67 
68 #define PetscStringizeArg(a) #a
69 #define PetscStringize(a) PetscStringizeArg(a)
70 
71 
72 /*MC
73    SETERRQ - Macro to be called when an error has been detected,
74 
75    Synopsis:
76    #include <petscsys.h>
77    PetscErrorCode SETERRQ(MPI_Comm comm,PetscErrorCode ierr,char *message)
78 
79    Collective
80 
81    Input Parameters:
82 +  comm - A communicator, use PETSC_COMM_SELF unless you know all ranks of another communicator will detect the error
83 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
84 -  message - error message
85 
86   Level: beginner
87 
88    Notes:
89     Once the error handler is called the calling function is then returned from with the given error code.
90 
91     See SETERRQ1(), SETERRQ2(), SETERRQ3() for versions that take arguments
92 
93     Experienced users can set the error handler with PetscPushErrorHandler().
94 
95    Fortran Notes:
96       SETERRQ() may be called from Fortran subroutines but SETERRA() must be called from the
97       Fortran main program.
98 
99 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(), CHKERRMPI()
100 M*/
101 #define SETERRQ(comm,ierr,s) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s)
102 
103 /*
104     Returned from PETSc functions that are called from MPI, such as related to attributes
105       Do not confuse PETSC_MPI_ERROR_CODE and PETSC_ERR_MPI, the first is registered with MPI and returned to MPI as
106       an error code, the latter is a regular PETSc error code passed within PETSc code indicating an error was detected in an MPI call.
107 */
108 PETSC_EXTERN PetscMPIInt PETSC_MPI_ERROR_CLASS;
109 PETSC_EXTERN PetscMPIInt PETSC_MPI_ERROR_CODE;
110 
111 
112 /*MC
113    SETERRMPI - Macro to be called when an error has been detected within an MPI callback function
114 
115    Synopsis:
116    #include <petscsys.h>
117    PetscErrorCode SETERRMPI(MPI_Comm comm,PetscErrorCode ierr,char *message)
118 
119    Collective
120 
121    Input Parameters:
122 +  comm - A communicator, use PETSC_COMM_SELF unless you know all ranks of another communicator will detect the error
123 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
124 -  message - error message
125 
126   Level: developer
127 
128    Notes:
129     This macro is FOR USE IN MPI CALLBACK FUNCTIONS ONLY, such as those passed to MPI_Comm_create_keyval(). It always returns the error code PETSC_MPI_ERROR_CODE
130     which is registered with MPI_Add_error_code() when PETSc is initialized.
131 
132 .seealso: SETERRQ(), CHKERRQ(), CHKERRMPI(), PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3()
133 M*/
134 #define SETERRMPI(comm,ierr,s) return (PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s),PETSC_MPI_ERROR_CODE)
135 
136 /*MC
137    SETERRQ1 - Macro that is called when an error has been detected,
138 
139    Synopsis:
140    #include <petscsys.h>
141    PetscErrorCode SETERRQ1(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg)
142 
143    Collective
144 
145    Input Parameters:
146 +  comm - A communicator, so that the error can be collective
147 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
148 .  message - error message in the printf format
149 -  arg - argument (for example an integer, string or double)
150 
151   Level: beginner
152 
153    Notes:
154     Once the error handler is called the calling function is then returned from with the given error code.
155 
156    Experienced users can set the error handler with PetscPushErrorHandler().
157 
158 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ(), SETERRQ2(), SETERRQ3()
159 M*/
160 #define SETERRQ1(comm,ierr,s,a1) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1)
161 
162 /*MC
163    SETERRQ2 - Macro that is called when an error has been detected,
164 
165    Synopsis:
166    #include <petscsys.h>
167    PetscErrorCode SETERRQ2(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2)
168 
169    Collective
170 
171    Input Parameters:
172 +  comm - A communicator, so that the error can be collective
173 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
174 .  message - error message in the printf format
175 .  arg1 - argument (for example an integer, string or double)
176 -  arg2 - argument (for example an integer, string or double)
177 
178   Level: beginner
179 
180    Notes:
181     Once the error handler is called the calling function is then returned from with the given error code.
182 
183    Experienced users can set the error handler with PetscPushErrorHandler().
184 
185 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ3()
186 M*/
187 #define SETERRQ2(comm,ierr,s,a1,a2) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2)
188 
189 /*MC
190    SETERRQ3 - Macro that is called when an error has been detected,
191 
192    Synopsis:
193    #include <petscsys.h>
194    PetscErrorCode SETERRQ3(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
195 
196    Collective
197 
198    Input Parameters:
199 +  comm - A communicator, so that the error can be collective
200 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
201 .  message - error message in the printf format
202 .  arg1 - argument (for example an integer, string or double)
203 .  arg2 - argument (for example an integer, string or double)
204 -  arg3 - argument (for example an integer, string or double)
205 
206   Level: beginner
207 
208    Notes:
209     Once the error handler is called the calling function is then returned from with the given error code.
210 
211     There are also versions for 4, 5, 6 and 7 arguments.
212 
213    Experienced users can set the error handler with PetscPushErrorHandler().
214 
215 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
216 M*/
217 #define SETERRQ3(comm,ierr,s,a1,a2,a3) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3)
218 
219 /*MC
220    SETERRQ4 - Macro that is called when an error has been detected,
221 
222    Synopsis:
223    #include <petscsys.h>
224    PetscErrorCode SETERRQ4(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3,arg4)
225 
226    Collective
227 
228    Input Parameters:
229 +  comm - A communicator, so that the error can be collective
230 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
231 .  message - error message in the printf format
232 .  arg1 - argument (for example an integer, string or double)
233 .  arg2 - argument (for example an integer, string or double)
234 .  arg3 - argument (for example an integer, string or double)
235 -  arg4 - argument (for example an integer, string or double)
236 
237   Level: beginner
238 
239    Notes:
240     Once the error handler is called the calling function is then returned from with the given error code.
241 
242     There are also versions for 4, 5, 6 and 7 arguments.
243 
244    Experienced users can set the error handler with PetscPushErrorHandler().
245 
246 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
247 M*/
248 #define SETERRQ4(comm,ierr,s,a1,a2,a3,a4) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4)
249 
250 /*MC
251    SETERRQ5 - Macro that is called when an error has been detected,
252 
253    Synopsis:
254    #include <petscsys.h>
255    PetscErrorCode SETERRQ5(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3,arg4,arg5)
256 
257    Collective
258 
259    Input Parameters:
260 +  comm - A communicator, so that the error can be collective
261 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
262 .  message - error message in the printf format
263 .  arg1 - argument (for example an integer, string or double)
264 .  arg2 - argument (for example an integer, string or double)
265 .  arg3 - argument (for example an integer, string or double)
266 .  arg4 - argument (for example an integer, string or double)
267 -  arg5 - argument (for example an integer, string or double)
268 
269   Level: beginner
270 
271    Notes:
272     Once the error handler is called the calling function is then returned from with the given error code.
273 
274     There are also versions for 4, 5, 6 and 7 arguments.
275 
276    Experienced users can set the error handler with PetscPushErrorHandler().
277 
278 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
279 M*/
280 #define SETERRQ5(comm,ierr,s,a1,a2,a3,a4,a5) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5)
281 
282 /*MC
283    SETERRQ6 - Macro that is called when an error has been detected,
284 
285    Synopsis:
286    #include <petscsys.h>
287    PetscErrorCode SETERRQ6(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3,arg4,arg5,arg6)
288 
289    Collective
290 
291    Input Parameters:
292 +  comm - A communicator, so that the error can be collective
293 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
294 .  message - error message in the printf format
295 .  arg1 - argument (for example an integer, string or double)
296 .  arg2 - argument (for example an integer, string or double)
297 .  arg3 - argument (for example an integer, string or double)
298 .  arg4 - argument (for example an integer, string or double)
299 .  arg5 - argument (for example an integer, string or double)
300 -  arg6 - argument (for example an integer, string or double)
301 
302   Level: beginner
303 
304    Notes:
305     Once the error handler is called the calling function is then returned from with the given error code.
306 
307     There are also versions for 4, 5, 6 and 7 arguments.
308 
309    Experienced users can set the error handler with PetscPushErrorHandler().
310 
311 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
312 M*/
313 #define SETERRQ6(comm,ierr,s,a1,a2,a3,a4,a5,a6) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6)
314 
315 /*MC
316    SETERRQ7 - Macro that is called when an error has been detected,
317 
318    Synopsis:
319    #include <petscsys.h>
320    PetscErrorCode SETERRQ7(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3,arg4,arg5,arg6,arg7)
321 
322    Collective
323 
324    Input Parameters:
325 +  comm - A communicator, so that the error can be collective
326 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
327 .  message - error message in the printf format
328 .  arg1 - argument (for example an integer, string or double)
329 .  arg2 - argument (for example an integer, string or double)
330 .  arg3 - argument (for example an integer, string or double)
331 .  arg4 - argument (for example an integer, string or double)
332 .  arg5 - argument (for example an integer, string or double)
333 .  arg6 - argument (for example an integer, string or double)
334 -  arg7 - argument (for example an integer, string or double)
335 
336   Level: beginner
337 
338    Notes:
339     Once the error handler is called the calling function is then returned from with the given error code.
340 
341     There are also versions for 4, 5, 6 and 7 arguments.
342 
343    Experienced users can set the error handler with PetscPushErrorHandler().
344 
345 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
346 M*/
347 #define SETERRQ7(comm,ierr,s,a1,a2,a3,a4,a5,a6,a7) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6,a7)
348 
349 /*MC
350    SETERRQ8 - Macro that is called when an error has been detected,
351 
352    Synopsis:
353    #include <petscsys.h>
354    PetscErrorCode SETERRQ8(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8)
355 
356    Collective
357 
358    Input Parameters:
359 +  comm - A communicator, so that the error can be collective
360 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
361 .  message - error message in the printf format
362 .  arg1 - argument (for example an integer, string or double)
363 .  arg2 - argument (for example an integer, string or double)
364 .  arg3 - argument (for example an integer, string or double)
365 .  arg4 - argument (for example an integer, string or double)
366 .  arg5 - argument (for example an integer, string or double)
367 .  arg6 - argument (for example an integer, string or double)
368 .  arg7 - argument (for example an integer, string or double)
369 -  arg8 - argument (for example an integer, string or double)
370 
371   Level: beginner
372 
373    Notes:
374     Once the error handler is called the calling function is then returned from with the given error code.
375 
376     There are also versions for 4, 5, 6 and 7 arguments.
377 
378    Experienced users can set the error handler with PetscPushErrorHandler().
379 
380 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
381 M*/
382 #define SETERRQ8(comm,ierr,s,a1,a2,a3,a4,a5,a6,a7,a8) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6,a7,a8)
383 
384 /*MC
385    SETERRQ9 - Macro that is called when an error has been detected,
386 
387    Synopsis:
388    #include <petscsys.h>
389    PetscErrorCode SETERRQ9(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8,arg9)
390 
391    Collective
392 
393    Input Parameters:
394 +  comm - A communicator, so that the error can be collective
395 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
396 .  message - error message in the printf format
397 .  arg1 - argument (for example an integer, string or double)
398 .  arg2 - argument (for example an integer, string or double)
399 .  arg3 - argument (for example an integer, string or double)
400 .  arg4 - argument (for example an integer, string or double)
401 .  arg5 - argument (for example an integer, string or double)
402 .  arg6 - argument (for example an integer, string or double)
403 .  arg7 - argument (for example an integer, string or double)
404 .  arg8 - argument (for example an integer, string or double)
405 -  arg9 - argument (for example an integer, string or double)
406 
407   Level: beginner
408 
409    Notes:
410     Once the error handler is called the calling function is then returned from with the given error code.
411 
412     There are also versions for 0 to 9 arguments.
413 
414    Experienced users can set the error handler with PetscPushErrorHandler().
415 
416 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
417 M*/
418 #define SETERRQ9(comm,ierr,s,a1,a2,a3,a4,a5,a6,a7,a8,a9) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3,a4,a5,a6,a7,a8,a9)
419 
420 /*MC
421    SETERRABORT - Macro that can be called when an error has been detected,
422 
423    Synopsis:
424    #include <petscsys.h>
425    PetscErrorCode SETERRABORT(MPI_Comm comm,PetscErrorCode ierr,char *message)
426 
427    Collective
428 
429    Input Parameters:
430 +  comm - A communicator, so that the error can be collective
431 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
432 -  message - error message in the printf format
433 
434   Level: beginner
435 
436    Notes:
437     This function just calls MPI_Abort().
438 
439 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
440 M*/
441 #define SETERRABORT(comm,ierr,s) do {PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s);MPI_Abort(comm,ierr);} while (0)
442 
443 /*MC
444    CHKERRQ - Checks error code returned from PETSc function, if non-zero it calls the error handler and then returns. Use CHKERRMPI() for checking errors from MPI calls
445 
446    Synopsis:
447    #include <petscsys.h>
448    PetscErrorCode CHKERRQ(PetscErrorCode ierr)
449 
450    Not Collective
451 
452    Input Parameters:
453 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
454 
455   Level: beginner
456 
457    Notes:
458     Once the error handler is called the calling function is then returned from with the given error code.
459 
460     Experienced users can set the error handler with PetscPushErrorHandler().
461 
462     CHKERRQ(ierr) is fundamentally a macro replacement for
463          if (ierr) return(PetscError(...,ierr,...));
464 
465     Although typical usage resembles "void CHKERRQ(PetscErrorCode)" as described above, for certain uses it is
466     highly inappropriate to use it in this manner as it invokes return(PetscErrorCode). In particular,
467     it cannot be used in functions which return(void) or any other datatype.  In these types of functions,
468     you can use CHKERRV() which returns without an error code (bad idea since the error is ignored or
469          if (ierr) {PetscError(....); return(YourReturnType);}
470     where you may pass back a NULL to indicate an error. You can also call CHKERRABORT(comm,n) to have
471     MPI_Abort() returned immediately.
472 
473    Fortran Notes:
474       CHKERRQ() may be called from Fortran subroutines but CHKERRA() must be called from the
475       Fortran main program.
476 
477 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2()
478 M*/
479 #define CHKERRQ(ierr)          do {PetscErrorCode ierr__ = (ierr); if (PetscUnlikely(ierr__)) return PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr__,PETSC_ERROR_REPEAT," ");} while (0)
480 #define CHKERRV(ierr)          do {PetscErrorCode ierr__ = (ierr); if (PetscUnlikely(ierr__)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr__,PETSC_ERROR_REPEAT," ");return;}} while (0)
481 
482 /*MC
483    CHKERRABORT - Checks error code returned from PETSc function. If non-zero it aborts immediately.
484 
485    Synopsis:
486    #include <petscsys.h>
487    PetscErrorCode CHKERRABORT(MPI_Comm comm,PetscErrorCode ierr)
488 
489    Not Collective
490 
491    Input Parameters:
492 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
493 
494   Level: intermediate
495 
496 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2(), SETERRABORT(), CHKERRMPI()
497 M*/
498 #define CHKERRABORT(comm,ierr) do {PetscErrorCode ierr__ = (ierr); if (PetscUnlikely(ierr__)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr__,PETSC_ERROR_REPEAT," ");MPI_Abort(comm,ierr);}} while (0)
499 #define CHKERRCONTINUE(ierr)   do {PetscErrorCode ierr__ = (ierr); if (PetscUnlikely(ierr__)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr__,PETSC_ERROR_REPEAT," ");}} while (0)
500 
501 PETSC_EXTERN PetscErrorCode PetscAbortFindSourceFile_Private(const char*,PetscInt*);
502 PETSC_EXTERN PetscBool petscwaitonerrorflg,petscindebugger;
503 
504 /*MC
505    PETSCABORT - Call MPI_Abort with an informative error code
506 
507    Synopsis:
508    #include <petscsys.h>
509    PETSCABORT(MPI_Comm comm, PetscErrorCode ierr)
510 
511    Collective
512 
513    Input Parameters:
514 +  comm - A communicator, so that the error can be collective
515 -  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
516 
517    Level: advanced
518 
519    Notes:
520    We pass MPI_Abort() an error code of format XX_YYYY_ZZZ, where XX, YYYY are an index and line number of the file
521    where PETSCABORT is called, respectively. ZZZ is the PETSc error code.
522 
523    If XX is zero, this means that the call was made in the routine main().
524    If XX is one, that means 1) the file is not in PETSc (it may be in users code); OR 2) the file is in PETSc but PetscAbortSourceFiles[]
525      is out of date. PETSc developers have to update it.
526    Otherwise, look up the value of XX in the table PetscAbortSourceFiles[] in src/sys/error/err.c to map XX back to the source file where the PETSCABORT() was called.
527 
528    If the option -start_in_debugger was used then this calls abort() to stop the program in the debugger.
529 
530 M*/
531 #define PETSCABORT(comm,ierr)  \
532    do {                                                               \
533       PetscInt       idx = 0;                                         \
534       PetscMPIInt    errcode;                                         \
535       PetscAbortFindSourceFile_Private(__FILE__,&idx);                \
536       errcode = (PetscMPIInt)(0*idx*10000000 + 0*__LINE__*1000 + ierr);   \
537       if (petscwaitonerrorflg) PetscSleep(1000);                      \
538       if (petscindebugger) abort();                                   \
539       else MPI_Abort(comm,errcode);                                   \
540    } while (0)
541 
542 /*MC
543    CHKERRMPI - Checks error code returned from MPI calls, if non-zero it calls the error handler and then returns
544 
545    Synopsis:
546    #include <petscsys.h>
547    PetscErrorCode CHKERRMPI(PetscErrorCode ierr)
548 
549    Not Collective
550 
551    Input Parameters:
552 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
553 
554   Level: intermediate
555 
556    Notes:
557     Always returns the error code PETSC_ERR_MPI; the MPI error code and string are embedded in the string error message
558 
559 .seealso: CHKERRQ(), PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2(), SETERRMPI(), SETERRABORT(), CHKERRABORT()
560 M*/
561 #define CHKERRMPI(ierr) \
562 do { \
563   PetscErrorCode _7_errorcode = (ierr); \
564   if (PetscUnlikely(_7_errorcode)) { \
565     char _7_errorstring[MPI_MAX_ERROR_STRING]; \
566     PetscMPIInt _7_resultlen; \
567     MPI_Error_string(_7_errorcode,(char*)_7_errorstring,&_7_resultlen); (void)_7_resultlen; \
568     SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_MPI,"MPI error %d %s",(int)_7_errorcode,_7_errorstring); \
569   } \
570 } while (0)
571 
572 #ifdef PETSC_CLANGUAGE_CXX
573 
574 /*MC
575    CHKERRXX - Checks error code, if non-zero it calls the C++ error handler which throws an exception
576 
577    Synopsis:
578    #include <petscsys.h>
579    void CHKERRXX(PetscErrorCode ierr)
580 
581    Not Collective
582 
583    Input Parameters:
584 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
585 
586   Level: beginner
587 
588    Notes:
589     Once the error handler throws a ??? exception.
590 
591     You can use CHKERRV() which returns without an error code (bad idea since the error is ignored)
592     or CHKERRABORT(comm,n) to have MPI_Abort() returned immediately.
593 
594 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKERRQ(), CHKMEMQ
595 M*/
596 #define CHKERRXX(ierr)  do {if (PetscUnlikely(ierr)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_IN_CXX,0);}} while (0)
597 
598 #endif
599 
600 #if defined(PETSC_HAVE_CUDA)
601 #include <cusolverDn.h>
602 #include <cusolverSp.h>
603 PETSC_EXTERN const char* PetscCUSolverGetErrorName(cusolverStatus_t);
604 #define CHKERRCUSOLVER(stat)                       \
605   do {                 \
606    if (PetscUnlikely(stat)) { \
607      const char *name = PetscCUSolverGetErrorName(stat);                     \
608      if ((stat == CUSOLVER_STATUS_NOT_INITIALIZED) || (stat == CUSOLVER_STATUS_ALLOC_FAILED) || (stat == CUSOLVER_STATUS_INTERNAL_ERROR)) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_GPU_RESOURCE,"cuSolver error %d (%s). This indicates the GPU has run out resources",(int)stat,name); \
609      else SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_GPU,"cuSolver error %d (%s)",(int)stat,name); \
610    } \
611 } while (0)
612 #endif
613 
614 /* TODO: SEK:  Need to figure out the hipsolver issues */
615 #if defined(PETSC_HAVE_HIP)
616 #define CHKERRHIPSOLVER(err) do {if (PetscUnlikely(err)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"HIPSOLVER error %d",err);} while (0)
617 #endif
618 /*MC
619    CHKMEMQ - Checks the memory for corruption, calls error handler if any is detected
620 
621    Synopsis:
622    #include <petscsys.h>
623    CHKMEMQ;
624 
625    Not Collective
626 
627   Level: beginner
628 
629    Notes:
630     We highly recommend using valgrind https://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind for finding memory problems. This is useful
631     on systems that do not have valgrind, but much much less useful.
632 
633     Must run with the option -malloc_debug (-malloc_test in debug mode; or if PetscMallocSetDebug() called) to enable this option
634 
635     Once the error handler is called the calling function is then returned from with the given error code.
636 
637     By defaults prints location where memory that is corrupted was allocated.
638 
639     Use CHKMEMA for functions that return void
640 
641 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(),
642           PetscMallocValidate()
643 M*/
644 #define CHKMEMQ do {PetscErrorCode _7_ierr = PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__);CHKERRQ(_7_ierr);} while (0)
645 
646 #define CHKMEMA PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__)
647 
648 /*E
649   PetscErrorType - passed to the PETSc error handling routines indicating if this is the first or a later call to the error handlers
650 
651   Level: advanced
652 
653   PETSC_ERROR_IN_CXX indicates the error was detected in C++ and an exception should be generated
654 
655   Developer Notes:
656     This is currently used to decide when to print the detailed information about the run in PetscTraceBackErrorHandler()
657 
658 .seealso: PetscError(), SETERRXX()
659 E*/
660 typedef enum {PETSC_ERROR_INITIAL=0,PETSC_ERROR_REPEAT=1,PETSC_ERROR_IN_CXX = 2} PetscErrorType;
661 
662 #if defined(__clang_analyzer__)
663 __attribute__((analyzer_noreturn))
664 #endif
665 PETSC_EXTERN PetscErrorCode PetscError(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,...);
666 
667 PETSC_EXTERN PetscErrorCode PetscErrorPrintfInitialize(void);
668 PETSC_EXTERN PetscErrorCode PetscErrorMessage(int,const char*[],char **);
669 PETSC_EXTERN PetscErrorCode PetscTraceBackErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
670 PETSC_EXTERN PetscErrorCode PetscIgnoreErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
671 PETSC_EXTERN PetscErrorCode PetscEmacsClientErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
672 PETSC_EXTERN PetscErrorCode PetscMPIAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
673 PETSC_EXTERN PetscErrorCode PetscAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
674 PETSC_EXTERN PetscErrorCode PetscAttachDebuggerErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
675 PETSC_EXTERN PetscErrorCode PetscReturnErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
676 PETSC_EXTERN PetscErrorCode PetscPushErrorHandler(PetscErrorCode (*handler)(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*),void*);
677 PETSC_EXTERN PetscErrorCode PetscPopErrorHandler(void);
678 PETSC_EXTERN PetscErrorCode PetscSignalHandlerDefault(int,void*);
679 PETSC_EXTERN PetscErrorCode PetscPushSignalHandler(PetscErrorCode (*)(int,void *),void*);
680 PETSC_EXTERN PetscErrorCode PetscPopSignalHandler(void);
681 PETSC_EXTERN PetscErrorCode PetscCheckPointerSetIntensity(PetscInt);
682 PETSC_EXTERN void PetscSignalSegvCheckPointerOrMpi(void);
683 PETSC_DEPRECATED_FUNCTION("Use PetscSignalSegvCheckPointerOrMpi() (since version 3.13)") PETSC_STATIC_INLINE void PetscSignalSegvCheckPointer(void) {PetscSignalSegvCheckPointerOrMpi();}
684 
685 /*MC
686     PetscErrorPrintf - Prints error messages.
687 
688    Synopsis:
689     #include <petscsys.h>
690      PetscErrorCode (*PetscErrorPrintf)(const char format[],...);
691 
692     Not Collective
693 
694     Input Parameters:
695 .   format - the usual printf() format string
696 
697    Options Database Keys:
698 +    -error_output_stdout - cause error messages to be printed to stdout instead of the  (default) stderr
699 -    -error_output_none - to turn off all printing of error messages (does not change the way the error is handled.)
700 
701    Notes:
702     Use
703 $     PetscErrorPrintf = PetscErrorPrintfNone; to turn off all printing of error messages (does not change the way the
704 $                        error is handled.) and
705 $     PetscErrorPrintf = PetscErrorPrintfDefault; to turn it back on or you can use your own function
706 
707           Use
708      PETSC_STDERR = FILE* obtained from a file open etc. to have stderr printed to the file.
709      PETSC_STDOUT = FILE* obtained from a file open etc. to have stdout printed to the file.
710 
711           Use
712       PetscPushErrorHandler() to provide your own error handler that determines what kind of messages to print
713 
714    Level: developer
715 
716     Fortran Note:
717     This routine is not supported in Fortran.
718 
719 
720 .seealso: PetscFPrintf(), PetscSynchronizedPrintf(), PetscHelpPrintf(), PetscPrintf(), PetscPushErrorHandler(), PetscVFPrintf(), PetscHelpPrintf()
721 M*/
722 PETSC_EXTERN PetscErrorCode (*PetscErrorPrintf)(const char[],...);
723 
724 typedef enum {PETSC_FP_TRAP_OFF=0,PETSC_FP_TRAP_ON=1} PetscFPTrap;
725 PETSC_EXTERN PetscErrorCode PetscSetFPTrap(PetscFPTrap);
726 PETSC_EXTERN PetscErrorCode PetscFPTrapPush(PetscFPTrap);
727 PETSC_EXTERN PetscErrorCode PetscFPTrapPop(void);
728 PETSC_EXTERN PetscErrorCode PetscDetermineInitialFPTrap(void);
729 
730 /*
731       Allows the code to build a stack frame as it runs
732 */
733 
734 #define PETSCSTACKSIZE 64
735 
736 typedef struct  {
737   const char      *function[PETSCSTACKSIZE];
738   const char      *file[PETSCSTACKSIZE];
739         int       line[PETSCSTACKSIZE];
740         PetscBool petscroutine[PETSCSTACKSIZE];
741         int       currentsize;
742         int       hotdepth;
743 } PetscStack;
744 
745 PETSC_EXTERN PetscStack *petscstack;
746 
747 PetscErrorCode  PetscStackCopy(PetscStack*,PetscStack*);
748 PetscErrorCode  PetscStackPrint(PetscStack *,FILE*);
749 #if defined(PETSC_SERIALIZE_FUNCTIONS)
750 #include <petsc/private/petscfptimpl.h>
751 /*
752    Registers the current function into the global function pointer to function name table
753 
754    Have to fix this to handle errors but cannot return error since used in PETSC_VIEWER_DRAW_() etc
755 */
756 #define PetscRegister__FUNCT__() do { \
757   static PetscBool __chked = PETSC_FALSE; \
758   if (!__chked) {\
759   void *ptr; PetscDLSym(NULL,PETSC_FUNCTION_NAME,&ptr);\
760   __chked = PETSC_TRUE;\
761   }} while (0)
762 #else
763 #define PetscRegister__FUNCT__()
764 #endif
765 
766 #if defined(PETSC_USE_DEBUG)
767 PETSC_STATIC_INLINE PetscBool PetscStackActive(void)
768 {
769   return(petscstack ? PETSC_TRUE : PETSC_FALSE);
770 }
771 
772 /* Stack handling is based on the following two "NoCheck" macros.  These should only be called directly by other error
773  * handling macros.  We record the line of the call, which may or may not be the location of the definition.  But is at
774  * least more useful than "unknown" because it can distinguish multiple calls from the same function.
775  */
776 
777 #define PetscStackPushNoCheck(funct,petsc_routine,hot)                        \
778   do {                                                                        \
779     PetscStackSAWsTakeAccess();                                                \
780     if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) {         \
781       petscstack->function[petscstack->currentsize]  = funct;               \
782       petscstack->file[petscstack->currentsize]      = __FILE__;            \
783       petscstack->line[petscstack->currentsize]      = __LINE__;            \
784       petscstack->petscroutine[petscstack->currentsize] = petsc_routine;    \
785       petscstack->currentsize++;                                             \
786     }                                                                         \
787     if (petscstack) {                                                        \
788       petscstack->hotdepth += (hot || petscstack->hotdepth);                \
789     }                                                                         \
790     PetscStackSAWsGrantAccess();                                               \
791   } while (0)
792 
793 #define PetscStackPopNoCheck                                            \
794   do {                                                                  \
795     PetscStackSAWsTakeAccess();                                          \
796     if (petscstack && petscstack->currentsize > 0) {                  \
797       petscstack->currentsize--;                                       \
798       petscstack->function[petscstack->currentsize]  = NULL;             \
799       petscstack->file[petscstack->currentsize]      = NULL;             \
800       petscstack->line[petscstack->currentsize]      = 0;             \
801       petscstack->petscroutine[petscstack->currentsize] = PETSC_FALSE;\
802     }                                                                   \
803     if (petscstack) {                                                  \
804       petscstack->hotdepth = PetscMax(petscstack->hotdepth-1,0);      \
805     }                                                                   \
806     PetscStackSAWsGrantAccess();                                         \
807   } while (0)
808 
809 /*MC
810    PetscFunctionBegin - First executable line of each PETSc function,  used for error handling. Final
811       line of PETSc functions should be PetscFunctionReturn(0);
812 
813    Synopsis:
814    #include <petscsys.h>
815    void PetscFunctionBegin;
816 
817    Not Collective
818 
819    Usage:
820 .vb
821      int something;
822 
823      PetscFunctionBegin;
824 .ve
825 
826    Notes:
827      Use PetscFunctionBeginUser for application codes.
828 
829      Not available in Fortran
830 
831    Level: developer
832 
833 .seealso: PetscFunctionReturn(), PetscFunctionBeginHot(), PetscFunctionBeginUser()
834 
835 M*/
836 #define PetscFunctionBegin do {                                        \
837     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_FALSE); \
838     PetscRegister__FUNCT__();                                          \
839   } while (0)
840 
841 /*MC
842    PetscFunctionBeginHot - Substitute for PetscFunctionBegin to be used in functions that are called in
843    performance-critical circumstances.  Use of this function allows for lighter profiling by default.
844 
845    Synopsis:
846    #include <petscsys.h>
847    void PetscFunctionBeginHot;
848 
849    Not Collective
850 
851    Usage:
852 .vb
853      int something;
854 
855      PetscFunctionBeginHot;
856 .ve
857 
858    Notes:
859      Not available in Fortran
860 
861    Level: developer
862 
863 .seealso: PetscFunctionBegin, PetscFunctionReturn()
864 
865 M*/
866 #define PetscFunctionBeginHot do {                                     \
867     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_TRUE);  \
868     PetscRegister__FUNCT__();                                          \
869   } while (0)
870 
871 /*MC
872    PetscFunctionBeginUser - First executable line of user provided PETSc routine
873 
874    Synopsis:
875    #include <petscsys.h>
876    void PetscFunctionBeginUser;
877 
878    Not Collective
879 
880    Usage:
881 .vb
882      int something;
883 
884      PetscFunctionBeginUser;
885 .ve
886 
887    Notes:
888       Final line of PETSc functions should be PetscFunctionReturn(0) except for main().
889 
890       Not available in Fortran
891 
892       This is identical to PetscFunctionBegin except it labels the routine as a user
893       routine instead of as a PETSc library routine.
894 
895    Level: intermediate
896 
897 .seealso: PetscFunctionReturn(), PetscFunctionBegin, PetscFunctionBeginHot
898 
899 M*/
900 #define PetscFunctionBeginUser                                          \
901   do {                                                                  \
902     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_FALSE,PETSC_FALSE); \
903     PetscRegister__FUNCT__();                                           \
904   } while (0)
905 
906 
907 #define PetscStackPush(n) \
908   do {                                                                  \
909     PetscStackPushNoCheck(n,PETSC_FALSE,PETSC_FALSE);                   \
910     CHKMEMQ;                                                            \
911   } while (0)
912 
913 #define PetscStackPop                           \
914     do {                                        \
915       CHKMEMQ;                                  \
916       PetscStackPopNoCheck;                     \
917     } while (0)
918 
919 /*MC
920    PetscFunctionReturn - Last executable line of each PETSc function
921         used for error handling. Replaces return()
922 
923    Synopsis:
924    #include <petscsys.h>
925    void PetscFunctionReturn(0);
926 
927    Not Collective
928 
929    Usage:
930 .vb
931     ....
932      PetscFunctionReturn(0);
933    }
934 .ve
935 
936    Notes:
937      Not available in Fortran
938 
939    Level: developer
940 
941 .seealso: PetscFunctionBegin()
942 
943 M*/
944 #define PetscFunctionReturn(a) \
945   do {                                                                \
946     PetscStackPopNoCheck;                                             \
947     return(a);} while (0)
948 
949 #define PetscFunctionReturnVoid() \
950   do {                                                                \
951     PetscStackPopNoCheck;                                             \
952     return;} while (0)
953 
954 #else
955 
956 PETSC_STATIC_INLINE PetscBool PetscStackActive(void) {return PETSC_FALSE;}
957 #define PetscStackPushNoCheck(funct,petsc_routine,hot) do {} while (0)
958 #define PetscStackPopNoCheck                           do {} while (0)
959 #define PetscFunctionBegin
960 #define PetscFunctionBeginUser
961 #define PetscFunctionBeginHot
962 #define PetscFunctionReturn(a)    return(a)
963 #define PetscFunctionReturnVoid() return
964 #define PetscStackPop             CHKMEMQ
965 #define PetscStackPush(f)         CHKMEMQ
966 
967 #endif
968 
969 /*
970     PetscStackCall - Calls an external library routine or user function after pushing the name of the routine on the stack.
971 
972    Input Parameters:
973 +   name - string that gives the name of the function being called
974 -   routine - actual call to the routine, including ierr = and CHKERRQ(ierr);
975 
976    Note: Often one should use PetscStackCallStandard() instead. This routine is intended for external library routines that DO NOT return error codes
977 
978    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.
979 
980 
981 
982 */
983 #define PetscStackCall(name,routine) do { PetscStackPush(name);routine;PetscStackPop; } while (0)
984 
985 /*
986     PetscStackCallStandard - Calls an external library routine after pushing the name of the routine on the stack.
987 
988    Input Parameters:
989 +   func-  name of the routine
990 -   args - arguments to the routine surrounded by ()
991 
992    Notes:
993     This is intended for external package routines that return error codes. Use PetscStackCall() for those that do not.
994 
995    Developer Note: this is so that when an external packge routine results in a crash or corrupts memory, they get blamed instead of PETSc.
996 
997 */
998 #define PetscStackCallStandard(func,args) do {                                                            \
999     PetscErrorCode __ierr;                                                                                \
1000     PetscStackPush(#func);                                                                                \
1001     __ierr = func args;                                                                                   \
1002     PetscStackPop;                                                                                        \
1003     if (__ierr) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error in %s(): error code %d",#func,(int)__ierr); \
1004   } while (0)
1005 
1006 PETSC_EXTERN PetscErrorCode PetscStackCreate(void);
1007 PETSC_EXTERN PetscErrorCode PetscStackView(FILE*);
1008 PETSC_EXTERN PetscErrorCode PetscStackDestroy(void);
1009 
1010 #endif
1011