xref: /petsc/include/petscerror.h (revision 55e7fe800d976e85ed2b5cd8bfdef564daa37bd9)
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 include/petsc/finclude/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 wrote options which should not be changed */
61 
62 #define PETSC_ERR_MAX_VALUE        94  /* this is always the one more than the largest error code */
63 
64 #define PetscStringizeArg(a) #a
65 #define PetscStringize(a) PetscStringizeArg(a)
66 
67 
68 /*MC
69    SETERRQ - Macro to be called when an error has been detected,
70 
71    Synopsis:
72    #include <petscsys.h>
73    PetscErrorCode SETERRQ(MPI_Comm comm,PetscErrorCode ierr,char *message)
74 
75    Collective on MPI_Comm
76 
77    Input Parameters:
78 +  comm - A communicator, use PETSC_COMM_SELF unless you know all ranks of another communicator will detect the error
79 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
80 -  message - error message
81 
82   Level: beginner
83 
84    Notes:
85     Once the error handler is called the calling function is then returned from with the given error code.
86 
87     See SETERRQ1(), SETERRQ2(), SETERRQ3() for versions that take arguments
88 
89     Experienced users can set the error handler with PetscPushErrorHandler().
90 
91    Concepts: error^setting condition
92 
93 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3()
94 M*/
95 #define SETERRQ(comm,ierr,s) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s)
96 
97 /*MC
98    SETERRMPI - Macro to be called when an error has been detected within an MPI callback function
99 
100    Synopsis:
101    #include <petscsys.h>
102    PetscErrorCode SETERRMPI(MPI_Comm comm,PetscErrorCode ierr,char *message)
103 
104    Collective on MPI_Comm
105 
106    Input Parameters:
107 +  comm - A communicator, use PETSC_COMM_SELF unless you know all ranks of another communicator will detect the error
108 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
109 -  message - error message
110 
111   Level: developer
112 
113    Notes:
114     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
115     which is registered with MPI_Add_error_code() when PETSc is initialized.
116 
117    Concepts: error^setting condition
118 
119 .seealso: SETERRQ(), CHKERRQ(), CHKERRMPI(), PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3()
120 M*/
121 #define SETERRMPI(comm,ierr,s) return (PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s),PETSC_MPI_ERROR_CODE)
122 
123 /*MC
124    SETERRQ1 - Macro that is called when an error has been detected,
125 
126    Synopsis:
127    #include <petscsys.h>
128    PetscErrorCode SETERRQ1(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg)
129 
130    Collective on MPI_Comm
131 
132    Input Parameters:
133 +  comm - A communicator, so that the error can be collective
134 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
135 .  message - error message in the printf format
136 -  arg - argument (for example an integer, string or double)
137 
138   Level: beginner
139 
140    Notes:
141     Once the error handler is called the calling function is then returned from with the given error code.
142 
143    Experienced users can set the error handler with PetscPushErrorHandler().
144 
145    Concepts: error^setting condition
146 
147 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ(), SETERRQ2(), SETERRQ3()
148 M*/
149 #define SETERRQ1(comm,ierr,s,a1) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1)
150 
151 /*MC
152    SETERRQ2 - Macro that is called when an error has been detected,
153 
154    Synopsis:
155    #include <petscsys.h>
156    PetscErrorCode SETERRQ2(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2)
157 
158    Collective on MPI_Comm
159 
160    Input Parameters:
161 +  comm - A communicator, so that the error can be collective
162 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
163 .  message - error message in the printf format
164 .  arg1 - argument (for example an integer, string or double)
165 -  arg2 - argument (for example an integer, string or double)
166 
167   Level: beginner
168 
169    Notes:
170     Once the error handler is called the calling function is then returned from with the given error code.
171 
172    Experienced users can set the error handler with PetscPushErrorHandler().
173 
174    Concepts: error^setting condition
175 
176 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ3()
177 M*/
178 #define SETERRQ2(comm,ierr,s,a1,a2) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2)
179 
180 /*MC
181    SETERRQ3 - Macro that is called when an error has been detected,
182 
183    Synopsis:
184    #include <petscsys.h>
185    PetscErrorCode SETERRQ3(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
186 
187    Collective on MPI_Comm
188 
189    Input Parameters:
190 +  comm - A communicator, so that the error can be collective
191 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
192 .  message - error message in the printf format
193 .  arg1 - argument (for example an integer, string or double)
194 .  arg2 - argument (for example an integer, string or double)
195 -  arg3 - argument (for example an integer, string or double)
196 
197   Level: beginner
198 
199    Notes:
200     Once the error handler is called the calling function is then returned from with the given error code.
201 
202     There are also versions for 4, 5, 6 and 7 arguments.
203 
204    Experienced users can set the error handler with PetscPushErrorHandler().
205 
206    Concepts: error^setting condition
207 
208 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
209 M*/
210 #define SETERRQ3(comm,ierr,s,a1,a2,a3) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3)
211 
212 /*MC
213    SETERRQ4 - Macro that is called when an error has been detected,
214 
215    Synopsis:
216    #include <petscsys.h>
217    PetscErrorCode SETERRQ4(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
218 
219    Collective on MPI_Comm
220 
221    Input Parameters:
222 +  comm - A communicator, so that the error can be collective
223 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
224 .  message - error message in the printf format
225 .  arg1 - argument (for example an integer, string or double)
226 .  arg2 - argument (for example an integer, string or double)
227 .  arg3 - argument (for example an integer, string or double)
228 -  arg4 - argument (for example an integer, string or double)
229 
230   Level: beginner
231 
232    Notes:
233     Once the error handler is called the calling function is then returned from with the given error code.
234 
235     There are also versions for 4, 5, 6 and 7 arguments.
236 
237    Experienced users can set the error handler with PetscPushErrorHandler().
238 
239    Concepts: error^setting condition
240 
241 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
242 M*/
243 #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)
244 
245 /*MC
246    SETERRQ5 - Macro that is called when an error has been detected,
247 
248    Synopsis:
249    #include <petscsys.h>
250    PetscErrorCode SETERRQ5(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
251 
252    Collective on MPI_COmm
253 
254    Input Parameters:
255 +  comm - A communicator, so that the error can be collective
256 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
257 .  message - error message in the printf format
258 .  arg1 - argument (for example an integer, string or double)
259 .  arg2 - argument (for example an integer, string or double)
260 .  arg3 - argument (for example an integer, string or double)
261 .  arg4 - argument (for example an integer, string or double)
262 -  arg5 - argument (for example an integer, string or double)
263 
264   Level: beginner
265 
266    Notes:
267     Once the error handler is called the calling function is then returned from with the given error code.
268 
269     There are also versions for 4, 5, 6 and 7 arguments.
270 
271    Experienced users can set the error handler with PetscPushErrorHandler().
272 
273    Concepts: error^setting condition
274 
275 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
276 M*/
277 #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)
278 
279 /*MC
280    SETERRQ6 - Macro that is called when an error has been detected,
281 
282    Synopsis:
283    #include <petscsys.h>
284    PetscErrorCode SETERRQ6(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
285 
286    Collective on MPI_Comm
287 
288    Input Parameters:
289 +  comm - A communicator, so that the error can be collective
290 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
291 .  message - error message in the printf format
292 .  arg1 - argument (for example an integer, string or double)
293 .  arg2 - argument (for example an integer, string or double)
294 .  arg3 - argument (for example an integer, string or double)
295 .  arg4 - argument (for example an integer, string or double)
296 .  arg5 - argument (for example an integer, string or double)
297 -  arg6 - argument (for example an integer, string or double)
298 
299   Level: beginner
300 
301    Notes:
302     Once the error handler is called the calling function is then returned from with the given error code.
303 
304     There are also versions for 4, 5, 6 and 7 arguments.
305 
306    Experienced users can set the error handler with PetscPushErrorHandler().
307 
308    Concepts: error^setting condition
309 
310 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
311 M*/
312 #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)
313 
314 /*MC
315    SETERRQ7 - Macro that is called when an error has been detected,
316 
317    Synopsis:
318    #include <petscsys.h>
319    PetscErrorCode SETERRQ7(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
320 
321    Collective on MPI_Comm
322 
323    Input Parameters:
324 +  comm - A communicator, so that the error can be collective
325 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
326 .  message - error message in the printf format
327 .  arg1 - argument (for example an integer, string or double)
328 .  arg2 - argument (for example an integer, string or double)
329 .  arg3 - argument (for example an integer, string or double)
330 .  arg4 - argument (for example an integer, string or double)
331 .  arg5 - argument (for example an integer, string or double)
332 .  arg6 - argument (for example an integer, string or double)
333 -  arg7 - argument (for example an integer, string or double)
334 
335   Level: beginner
336 
337    Notes:
338     Once the error handler is called the calling function is then returned from with the given error code.
339 
340     There are also versions for 4, 5, 6 and 7 arguments.
341 
342    Experienced users can set the error handler with PetscPushErrorHandler().
343 
344    Concepts: error^setting condition
345 
346 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
347 M*/
348 #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)
349 
350 /*MC
351    SETERRQ8 - Macro that is called when an error has been detected,
352 
353    Synopsis:
354    #include <petscsys.h>
355    PetscErrorCode SETERRQ8(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
356 
357    Collective on MPI_Comm
358 
359    Input Parameters:
360 +  comm - A communicator, so that the error can be collective
361 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
362 .  message - error message in the printf format
363 .  arg1 - argument (for example an integer, string or double)
364 .  arg2 - argument (for example an integer, string or double)
365 .  arg3 - argument (for example an integer, string or double)
366 .  arg4 - argument (for example an integer, string or double)
367 .  arg5 - argument (for example an integer, string or double)
368 .  arg6 - argument (for example an integer, string or double)
369 .  arg7 - argument (for example an integer, string or double)
370 -  arg8 - argument (for example an integer, string or double)
371 
372   Level: beginner
373 
374    Notes:
375     Once the error handler is called the calling function is then returned from with the given error code.
376 
377     There are also versions for 4, 5, 6 and 7 arguments.
378 
379    Experienced users can set the error handler with PetscPushErrorHandler().
380 
381    Concepts: error^setting condition
382 
383 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
384 M*/
385 #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)
386 
387 /*MC
388    SETERRQ9 - Macro that is called when an error has been detected,
389 
390    Synopsis:
391    #include <petscsys.h>
392    PetscErrorCode SETERRQ9(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
393 
394    Collective on MPI_Comm
395 
396    Input Parameters:
397 +  comm - A communicator, so that the error can be collective
398 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
399 .  message - error message in the printf format
400 .  arg1 - argument (for example an integer, string or double)
401 .  arg2 - argument (for example an integer, string or double)
402 .  arg3 - argument (for example an integer, string or double)
403 .  arg4 - argument (for example an integer, string or double)
404 .  arg5 - argument (for example an integer, string or double)
405 .  arg6 - argument (for example an integer, string or double)
406 .  arg7 - argument (for example an integer, string or double)
407 .  arg8 - argument (for example an integer, string or double)
408 -  arg9 - argument (for example an integer, string or double)
409 
410   Level: beginner
411 
412    Notes:
413     Once the error handler is called the calling function is then returned from with the given error code.
414 
415     There are also versions for 0 to 9 arguments.
416 
417    Experienced users can set the error handler with PetscPushErrorHandler().
418 
419    Concepts: error^setting condition
420 
421 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
422 M*/
423 #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)
424 
425 /*MC
426    SETERRABORT - Macro that can be called when an error has been detected,
427 
428    Synopsis:
429    #include <petscsys.h>
430    PetscErrorCode SETERRABORT(MPI_Comm comm,PetscErrorCode ierr,char *message)
431 
432    Collective on MPI_Comm
433 
434    Input Parameters:
435 +  comm - A communicator, so that the error can be collective
436 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
437 -  message - error message in the printf format
438 
439   Level: beginner
440 
441    Notes:
442     This function just calls MPI_Abort().
443 
444    Concepts: error^setting condition
445 
446 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
447 M*/
448 #define SETERRABORT(comm,ierr,s) do {PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s);MPI_Abort(comm,ierr);} while (0)
449 
450 /*MC
451    CHKERRQ - Checks error code, if non-zero it calls the error handler and then returns
452 
453    Synopsis:
454    #include <petscsys.h>
455    PetscErrorCode CHKERRQ(PetscErrorCode ierr)
456 
457    Not Collective
458 
459    Input Parameters:
460 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
461 
462   Level: beginner
463 
464    Notes:
465     Once the error handler is called the calling function is then returned from with the given error code.
466 
467     Experienced users can set the error handler with PetscPushErrorHandler().
468 
469     CHKERRQ(ierr) is fundamentally a macro replacement for
470          if (ierr) return(PetscError(...,ierr,...));
471 
472     Although typical usage resembles "void CHKERRQ(PetscErrorCode)" as described above, for certain uses it is
473     highly inappropriate to use it in this manner as it invokes return(PetscErrorCode). In particular,
474     it cannot be used in functions which return(void) or any other datatype.  In these types of functions,
475     you can use CHKERRV() which returns without an error code (bad idea since the error is ignored or
476          if (ierr) {PetscError(....); return(YourReturnType);}
477     where you may pass back a NULL to indicate an error. You can also call CHKERRABORT(comm,n) to have
478     MPI_Abort() returned immediately.
479 
480    Concepts: error^setting condition
481 
482 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2()
483 M*/
484 #define CHKERRQ(ierr)          do {if (PetscUnlikely(ierr)) return PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");} while (0)
485 #define CHKERRV(ierr)          do {if (PetscUnlikely(ierr)) {ierr = PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");return;}} while(0)
486 #define CHKERRABORT(comm,ierr) do {if (PetscUnlikely(ierr)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");MPI_Abort(comm,ierr);}} while (0)
487 #define CHKERRCONTINUE(ierr)   do {if (PetscUnlikely(ierr)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");}} while (0)
488 
489 
490 /*MC
491    CHKERRMPI - Checks error code, if non-zero it calls the error handler and then returns
492 
493    Synopsis:
494    #include <petscsys.h>
495    PetscErrorCode CHKERRMPI(PetscErrorCode ierr)
496 
497    Not Collective
498 
499    Input Parameters:
500 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
501 
502   Level: developer
503 
504    Notes:
505     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
506     which is registered with MPI_Add_error_code() when PETSc is initialized.
507 
508    Concepts: error^setting condition
509 
510 .seealso: CHKERRQ(), PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2()
511 M*/
512 #define CHKERRMPI(ierr)        do {if (PetscUnlikely(ierr)) return (PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," "),PETSC_MPI_ERROR_CODE);} while (0)
513 
514 #ifdef PETSC_CLANGUAGE_CXX
515 
516 /*MC
517    CHKERRXX - Checks error code, if non-zero it calls the C++ error handler which throws an exception
518 
519    Synopsis:
520    #include <petscsys.h>
521    void CHKERRXX(PetscErrorCode ierr)
522 
523    Not Collective
524 
525    Input Parameters:
526 .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
527 
528   Level: beginner
529 
530    Notes:
531     Once the error handler throws a ??? exception.
532 
533     You can use CHKERRV() which returns without an error code (bad idea since the error is ignored)
534     or CHKERRABORT(comm,n) to have MPI_Abort() returned immediately.
535 
536    Concepts: error^setting condition
537 
538 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKERRQ(), CHKMEMQ
539 M*/
540 #define CHKERRXX(ierr)  do {if (PetscUnlikely(ierr)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_IN_CXX,0);}} while(0)
541 
542 #endif
543 
544 #define CHKERRCUDA(err)   do {if (PetscUnlikely(err)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"CUDA error %d",err);} while(0)
545 #define CHKERRCUBLAS(err) do {if (PetscUnlikely(err)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"CUBLAS error %d",err);} while(0)
546 
547 /*MC
548    CHKMEMQ - Checks the memory for corruption, calls error handler if any is detected
549 
550    Synopsis:
551    #include <petscsys.h>
552    CHKMEMQ;
553 
554    Not Collective
555 
556   Level: beginner
557 
558    Notes:
559     We highly recommend using valgrind http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind for finding memory problems. This is useful
560     on systems that do not have valgrind, but much much less useful.
561 
562     Must run with the option -malloc_debug to enable this option
563 
564     Once the error handler is called the calling function is then returned from with the given error code.
565 
566     By defaults prints location where memory that is corrupted was allocated.
567 
568     Use CHKMEMA for functions that return void
569 
570    Concepts: memory corruption
571 
572 .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(),
573           PetscMallocValidate()
574 M*/
575 #define CHKMEMQ do {PetscErrorCode _7_ierr = PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__);CHKERRQ(_7_ierr);} while(0)
576 
577 #define CHKMEMA PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__)
578 
579 /*E
580   PetscErrorType - passed to the PETSc error handling routines indicating if this is the first or a later call to the error handlers
581 
582   Level: advanced
583 
584   PETSC_ERROR_IN_CXX indicates the error was detected in C++ and an exception should be generated
585 
586   Developer Notes:
587     This is currently used to decide when to print the detailed information about the run in PetscTraceBackErrorHandler()
588 
589 .seealso: PetscError(), SETERRXX()
590 E*/
591 typedef enum {PETSC_ERROR_INITIAL=0,PETSC_ERROR_REPEAT=1,PETSC_ERROR_IN_CXX = 2} PetscErrorType;
592 
593 #if defined(__clang_analyzer__)
594 __attribute__((analyzer_noreturn))
595 #endif
596 PETSC_EXTERN PetscErrorCode PetscError(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,...);
597 
598 PETSC_EXTERN PetscErrorCode PetscErrorPrintfInitialize(void);
599 PETSC_EXTERN PetscErrorCode PetscErrorMessage(int,const char*[],char **);
600 PETSC_EXTERN PetscErrorCode PetscTraceBackErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
601 PETSC_EXTERN PetscErrorCode PetscIgnoreErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
602 PETSC_EXTERN PetscErrorCode PetscEmacsClientErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
603 PETSC_EXTERN PetscErrorCode PetscMPIAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
604 PETSC_EXTERN PetscErrorCode PetscAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
605 PETSC_EXTERN PetscErrorCode PetscAttachDebuggerErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
606 PETSC_EXTERN PetscErrorCode PetscReturnErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
607 PETSC_EXTERN PetscErrorCode PetscPushErrorHandler(PetscErrorCode (*handler)(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*),void*);
608 PETSC_EXTERN PetscErrorCode PetscPopErrorHandler(void);
609 PETSC_EXTERN PetscErrorCode PetscSignalHandlerDefault(int,void*);
610 PETSC_EXTERN PetscErrorCode PetscPushSignalHandler(PetscErrorCode (*)(int,void *),void*);
611 PETSC_EXTERN PetscErrorCode PetscPopSignalHandler(void);
612 PETSC_EXTERN PetscErrorCode PetscCheckPointerSetIntensity(PetscInt);
613 
614 /*MC
615     PetscErrorPrintf - Prints error messages.
616 
617    Synopsis:
618     #include <petscsys.h>
619      PetscErrorCode (*PetscErrorPrintf)(const char format[],...);
620 
621     Not Collective
622 
623     Input Parameters:
624 .   format - the usual printf() format string
625 
626    Options Database Keys:
627 +    -error_output_stdout - cause error messages to be printed to stdout instead of the  (default) stderr
628 -    -error_output_none - to turn off all printing of error messages (does not change the way the error is handled.)
629 
630    Notes:
631     Use
632 $     PetscErrorPrintf = PetscErrorPrintfNone; to turn off all printing of error messages (does not change the way the
633 $                        error is handled.) and
634 $     PetscErrorPrintf = PetscErrorPrintfDefault; to turn it back on or you can use your own function
635 
636           Use
637      PETSC_STDERR = FILE* obtained from a file open etc. to have stderr printed to the file.
638      PETSC_STDOUT = FILE* obtained from a file open etc. to have stdout printed to the file.
639 
640           Use
641       PetscPushErrorHandler() to provide your own error handler that determines what kind of messages to print
642 
643    Level: developer
644 
645     Fortran Note:
646     This routine is not supported in Fortran.
647 
648     Concepts: error messages^printing
649     Concepts: printing^error messages
650 
651 .seealso: PetscFPrintf(), PetscSynchronizedPrintf(), PetscHelpPrintf(), PetscPrintf(), PetscPushErrorHandler(), PetscVFPrintf(), PetscHelpPrintf()
652 M*/
653 PETSC_EXTERN PetscErrorCode (*PetscErrorPrintf)(const char[],...);
654 
655 typedef enum {PETSC_FP_TRAP_OFF=0,PETSC_FP_TRAP_ON=1} PetscFPTrap;
656 PETSC_EXTERN PetscErrorCode PetscSetFPTrap(PetscFPTrap);
657 PETSC_EXTERN PetscErrorCode PetscFPTrapPush(PetscFPTrap);
658 PETSC_EXTERN PetscErrorCode PetscFPTrapPop(void);
659 
660 /*
661       Allows the code to build a stack frame as it runs
662 */
663 
664 #define PETSCSTACKSIZE 64
665 
666 typedef struct  {
667   const char      *function[PETSCSTACKSIZE];
668   const char      *file[PETSCSTACKSIZE];
669         int       line[PETSCSTACKSIZE];
670         PetscBool petscroutine[PETSCSTACKSIZE];
671         int       currentsize;
672         int       hotdepth;
673 } PetscStack;
674 
675 PETSC_EXTERN PetscStack *petscstack;
676 
677 PetscErrorCode  PetscStackCopy(PetscStack*,PetscStack*);
678 PetscErrorCode  PetscStackPrint(PetscStack *,FILE*);
679 #if defined(PETSC_SERIALIZE_FUNCTIONS)
680 #include <petsc/private/petscfptimpl.h>
681 /*
682    Registers the current function into the global function pointer to function name table
683 
684    Have to fix this to handle errors but cannot return error since used in PETSC_VIEWER_DRAW_() etc
685 */
686 #define PetscRegister__FUNCT__() do { \
687   static PetscBool __chked = PETSC_FALSE; \
688   if (!__chked) {\
689   void *ptr; PetscDLSym(NULL,PETSC_FUNCTION_NAME,&ptr);\
690   __chked = PETSC_TRUE;\
691   }} while (0)
692 #else
693 #define PetscRegister__FUNCT__()
694 #endif
695 
696 #if defined(PETSC_USE_DEBUG)
697 PETSC_STATIC_INLINE PetscBool PetscStackActive(void)
698 {
699   return(petscstack ? PETSC_TRUE : PETSC_FALSE);
700 }
701 
702 /* Stack handling is based on the following two "NoCheck" macros.  These should only be called directly by other error
703  * handling macros.  We record the line of the call, which may or may not be the location of the definition.  But is at
704  * least more useful than "unknown" because it can distinguish multiple calls from the same function.
705  */
706 
707 #define PetscStackPushNoCheck(funct,petsc_routine,hot)                        \
708   do {                                                                        \
709     PetscStackSAWsTakeAccess();                                                \
710     if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) {         \
711       petscstack->function[petscstack->currentsize]  = funct;               \
712       petscstack->file[petscstack->currentsize]      = __FILE__;            \
713       petscstack->line[petscstack->currentsize]      = __LINE__;            \
714       petscstack->petscroutine[petscstack->currentsize] = petsc_routine;    \
715       petscstack->currentsize++;                                             \
716     }                                                                         \
717     if (petscstack) {                                                        \
718       petscstack->hotdepth += (hot || petscstack->hotdepth);                \
719     }                                                                         \
720     PetscStackSAWsGrantAccess();                                               \
721   } while (0)
722 
723 #define PetscStackPopNoCheck                                            \
724   do {                                                                  \
725     PetscStackSAWsTakeAccess();                                          \
726     if (petscstack && petscstack->currentsize > 0) {                  \
727       petscstack->currentsize--;                                       \
728       petscstack->function[petscstack->currentsize]  = 0;             \
729       petscstack->file[petscstack->currentsize]      = 0;             \
730       petscstack->line[petscstack->currentsize]      = 0;             \
731       petscstack->petscroutine[petscstack->currentsize] = PETSC_FALSE;\
732     }                                                                   \
733     if (petscstack) {                                                  \
734       petscstack->hotdepth = PetscMax(petscstack->hotdepth-1,0);      \
735     }                                                                   \
736     PetscStackSAWsGrantAccess();                                         \
737   } while (0)
738 
739 /*MC
740    PetscFunctionBegin - First executable line of each PETSc function,  used for error handling. Final
741       line of PETSc functions should be PetscFunctionReturn(0);
742 
743    Synopsis:
744    #include <petscsys.h>
745    void PetscFunctionBegin;
746 
747    Not Collective
748 
749    Usage:
750 .vb
751      int something;
752 
753      PetscFunctionBegin;
754 .ve
755 
756    Notes:
757      Use PetscFunctionBeginUser for application codes.
758 
759      Not available in Fortran
760 
761    Level: developer
762 
763 .seealso: PetscFunctionReturn(), PetscFunctionBeginHot(), PetscFunctionBeginUser()
764 
765 .keywords: traceback, error handling
766 M*/
767 #define PetscFunctionBegin do {                                        \
768     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_FALSE); \
769     PetscRegister__FUNCT__();                                          \
770   } while (0)
771 
772 /*MC
773    PetscFunctionBeginHot - Substitute for PetscFunctionBegin to be used in functions that are called in
774    performance-critical circumstances.  Use of this function allows for lighter profiling by default.
775 
776    Synopsis:
777    #include <petscsys.h>
778    void PetscFunctionBeginHot;
779 
780    Not Collective
781 
782    Usage:
783 .vb
784      int something;
785 
786      PetscFunctionBeginHot;
787 .ve
788 
789    Notes:
790      Not available in Fortran
791 
792    Level: developer
793 
794 .seealso: PetscFunctionBegin, PetscFunctionReturn()
795 
796 .keywords: traceback, error handling
797 M*/
798 #define PetscFunctionBeginHot do {                                     \
799     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_TRUE);  \
800     PetscRegister__FUNCT__();                                          \
801   } while (0)
802 
803 /*MC
804    PetscFunctionBeginUser - First executable line of user provided PETSc routine
805 
806    Synopsis:
807    #include <petscsys.h>
808    void PetscFunctionBeginUser;
809 
810    Not Collective
811 
812    Usage:
813 .vb
814      int something;
815 
816      PetscFunctionBeginUser;
817 .ve
818 
819    Notes:
820       Final line of PETSc functions should be PetscFunctionReturn(0) except for main().
821 
822       Not available in Fortran
823 
824       This is identical to PetscFunctionBegin except it labels the routine as a user
825       routine instead of as a PETSc library routine.
826 
827    Level: intermediate
828 
829 .seealso: PetscFunctionReturn(), PetscFunctionBegin, PetscFunctionBeginHot
830 
831 .keywords: traceback, error handling
832 M*/
833 #define PetscFunctionBeginUser                                          \
834   do {                                                                  \
835     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_FALSE,PETSC_FALSE); \
836     PetscRegister__FUNCT__();                                           \
837   } while (0)
838 
839 
840 #define PetscStackPush(n) \
841   do {                                                                  \
842     PetscStackPushNoCheck(n,PETSC_FALSE,PETSC_FALSE);                   \
843     CHKMEMQ;                                                            \
844   } while (0)
845 
846 #define PetscStackPop                           \
847     do {                                        \
848       CHKMEMQ;                                  \
849       PetscStackPopNoCheck;                     \
850     } while (0)
851 
852 /*MC
853    PetscFunctionReturn - Last executable line of each PETSc function
854         used for error handling. Replaces return()
855 
856    Synopsis:
857    #include <petscsys.h>
858    void PetscFunctionReturn(0);
859 
860    Not Collective
861 
862    Usage:
863 .vb
864     ....
865      PetscFunctionReturn(0);
866    }
867 .ve
868 
869    Notes:
870      Not available in Fortran
871 
872    Level: developer
873 
874 .seealso: PetscFunctionBegin()
875 
876 .keywords: traceback, error handling
877 M*/
878 #define PetscFunctionReturn(a) \
879   do {                                                                \
880     PetscStackPopNoCheck;                                             \
881     return(a);} while (0)
882 
883 #define PetscFunctionReturnVoid() \
884   do {                                                                \
885     PetscStackPopNoCheck;                                             \
886     return;} while (0)
887 
888 #else
889 
890 PETSC_STATIC_INLINE PetscBool PetscStackActive(void) {return PETSC_FALSE;}
891 #define PetscStackPushNoCheck(funct,petsc_routine,hot) do {} while (0)
892 #define PetscStackPopNoCheck                           do {} while (0)
893 #define PetscFunctionBegin
894 #define PetscFunctionBeginUser
895 #define PetscFunctionBeginHot
896 #define PetscFunctionReturn(a)    return(a)
897 #define PetscFunctionReturnVoid() return
898 #define PetscStackPop             CHKMEMQ
899 #define PetscStackPush(f)         CHKMEMQ
900 
901 #endif
902 
903 /*
904     PetscStackCall - Calls an external library routine or user function after pushing the name of the routine on the stack.
905 
906    Input Parameters:
907 +   name - string that gives the name of the function being called
908 -   routine - actual call to the routine, including ierr = and CHKERRQ(ierr);
909 
910    Note: Often one should use PetscStackCallStandard() instead. This routine is intended for external library routines that DO NOT return error codes
911 
912    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.
913 
914 
915 
916 */
917 #define PetscStackCall(name,routine) do { PetscStackPush(name);routine;PetscStackPop; } while(0)
918 
919 /*
920     PetscStackCallStandard - Calls an external library routine after pushing the name of the routine on the stack.
921 
922    Input Parameters:
923 +   func-  name of the routine
924 -   args - arguments to the routine surrounded by ()
925 
926    Notes:
927     This is intended for external package routines that return error codes. Use PetscStackCall() for those that do not.
928 
929    Developer Note: this is so that when an external packge routine results in a crash or corrupts memory, they get blamed instead of PETSc.
930 
931 */
932 #define PetscStackCallStandard(func,args) do {                        \
933     PetscStackPush(#func);ierr = func args;PetscStackPop; if (ierr) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error in %s(): error code %d",#func,(int)ierr); \
934   } while (0)
935 
936 PETSC_EXTERN PetscErrorCode PetscStackCreate(void);
937 PETSC_EXTERN PetscErrorCode PetscStackView(FILE*);
938 PETSC_EXTERN PetscErrorCode PetscStackDestroy(void);
939 
940 #endif
941