/include/petscerror.h

https://bitbucket.org/petsc/petsc · C Header · 910 lines · 211 code · 60 blank · 639 comment · 34 complexity · ed7538d315b52b8a5232f25c6424cbfe MD5 · raw file 

    

  1. /*
    2. 

  2.     Contains all error handling interfaces for PETSc.
    3. 

  3. */
    4. 

  4. #if !defined(PETSCERROR_H)
    5. 

  5. #define PETSCERROR_H
    6. 

  6. 

  7. /*
    8. 

  8.      These are the generic error codes. These error codes are used
    9. 

  9.      many different places in the PETSc source code. The string versions are
    10. 

  10.      at src/sys/error/err.c any changes here must also be made there
    11. 

  11.      These are also define in include/petsc/finclude/petscerror.h any CHANGES here
    12. 

  12.      must be also made there.
    13. 

    14. 

  14. */
    15. 

  15. #define PETSC_ERR_MIN_VALUE        54   /* should always be one less then the smallest value */
    16. 

  16. 

  17. #define PETSC_ERR_MEM              55   /* unable to allocate requested memory */
    18. 

  18. #define PETSC_ERR_SUP              56   /* no support for requested operation */
    19. 

  19. #define PETSC_ERR_SUP_SYS          57   /* no support for requested operation on this computer system */
    20. 

  20. #define PETSC_ERR_ORDER            58   /* operation done in wrong order */
    21. 

  21. #define PETSC_ERR_SIG              59   /* signal received */
    22. 

  22. #define PETSC_ERR_FP               72   /* floating point exception */
    23. 

  23. #define PETSC_ERR_COR              74   /* corrupted PETSc object */
    24. 

  24. #define PETSC_ERR_LIB              76   /* error in library called by PETSc */
    25. 

  25. #define PETSC_ERR_PLIB             77   /* PETSc library generated inconsistent data */
    26. 

  26. #define PETSC_ERR_MEMC             78   /* memory corruption */
    27. 

  27. #define PETSC_ERR_CONV_FAILED      82   /* iterative method (KSP or SNES) failed */
    28. 

  28. #define PETSC_ERR_USER             83   /* user has not provided needed function */
    29. 

  29. #define PETSC_ERR_SYS              88   /* error in system call */
    30. 

  30. #define PETSC_ERR_POINTER          70   /* pointer does not point to valid address */
    31. 

  31. #define PETSC_ERR_MPI_LIB_INCOMP   87   /* MPI library at runtime is not compatible with MPI user compiled with */
    32. 

  32. 

  33. #define PETSC_ERR_ARG_SIZ          60   /* nonconforming object sizes used in operation */
    34. 

  34. #define PETSC_ERR_ARG_IDN          61   /* two arguments not allowed to be the same */
    35. 

  35. #define PETSC_ERR_ARG_WRONG        62   /* wrong argument (but object probably ok) */
    36. 

  36. #define PETSC_ERR_ARG_CORRUPT      64   /* null or corrupted PETSc object as argument */
    37. 

  37. #define PETSC_ERR_ARG_OUTOFRANGE   63   /* input argument, out of range */
    38. 

  38. #define PETSC_ERR_ARG_BADPTR       68   /* invalid pointer argument */
    39. 

  39. #define PETSC_ERR_ARG_NOTSAMETYPE  69   /* two args must be same object type */
    40. 

  40. #define PETSC_ERR_ARG_NOTSAMECOMM  80   /* two args must be same communicators */
    41. 

  41. #define PETSC_ERR_ARG_WRONGSTATE   73   /* object in argument is in wrong state, e.g. unassembled mat */
    42. 

  42. #define PETSC_ERR_ARG_TYPENOTSET   89   /* the type of the object has not yet been set */
    43. 

  43. #define PETSC_ERR_ARG_INCOMP       75   /* two arguments are incompatible */
    44. 

  44. #define PETSC_ERR_ARG_NULL         85   /* argument is null that should not be */
    45. 

  45. #define PETSC_ERR_ARG_UNKNOWN_TYPE 86   /* type name doesn't match any registered type */
    46. 

  46. 

  47. #define PETSC_ERR_FILE_OPEN        65   /* unable to open file */
    48. 

  48. #define PETSC_ERR_FILE_READ        66   /* unable to read from file */
    49. 

  49. #define PETSC_ERR_FILE_WRITE       67   /* unable to write to file */
    50. 

  50. #define PETSC_ERR_FILE_UNEXPECTED  79   /* unexpected data in file */
    51. 

  51. 

  52. #define PETSC_ERR_MAT_LU_ZRPVT     71   /* detected a zero pivot during LU factorization */
    53. 

  53. #define PETSC_ERR_MAT_CH_ZRPVT     81   /* detected a zero pivot during Cholesky factorization */
    54. 

  54. 

  55. #define PETSC_ERR_INT_OVERFLOW     84
    56. 

  56. 

  57. #define PETSC_ERR_FLOP_COUNT       90
    58. 

  58. #define PETSC_ERR_NOT_CONVERGED    91  /* solver did not converge */
    59. 

  59. #define PETSC_ERR_MISSING_FACTOR   92  /* MatGetFactor() failed */
    60. 

  60. #define PETSC_ERR_OPT_OVERWRITE    93  /* attempted to over wrote options which should not be changed */
    61. 

  61. 

  62. #define PETSC_ERR_MAX_VALUE        94  /* this is always the one more than the largest error code */
    63. 

  63. 

  64. #define PetscStringizeArg(a) #a
    65. 

  65. #define PetscStringize(a) PetscStringizeArg(a)
    66. 

  66. 

  67. 

  68. /*MC
    69. 

  69.    SETERRQ - Macro to be called when an error has been detected,
    70. 

    71. 

  71.    Synopsis:
    72. 

  72.    #include <petscsys.h>
    73. 

  73.    PetscErrorCode SETERRQ(MPI_Comm comm,PetscErrorCode ierr,char *message)
    74. 

    75. 

  75.    Collective
    76. 

    77. 

  77.    Input Parameters:
    78. 

  78. +  comm - A communicator, use PETSC_COMM_SELF unless you know all ranks of another communicator will detect the error
    79. 

  79. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    80. 

  80. -  message - error message
    81. 

    82. 

  82.   Level: beginner
    83. 

    84. 

  84.    Notes:
    85. 

  85.     Once the error handler is called the calling function is then returned from with the given error code.
    86. 

    87. 

  87.     See SETERRQ1(), SETERRQ2(), SETERRQ3() for versions that take arguments
    88. 

    89. 

  89.     Experienced users can set the error handler with PetscPushErrorHandler().
    90. 

    91. 

  91. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3()
    92. 

  92. M*/
    93. 

  93. #define SETERRQ(comm,ierr,s) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s)
    94. 

  94. 

  95. /*MC
    96. 

  96.    SETERRMPI - Macro to be called when an error has been detected within an MPI callback function
    97. 

    98. 

  98.    Synopsis:
    99. 

  99.    #include <petscsys.h>
    100. 

  100.    PetscErrorCode SETERRMPI(MPI_Comm comm,PetscErrorCode ierr,char *message)
    101. 

    102. 

  102.    Collective
    103. 

    104. 

  104.    Input Parameters:
    105. 

  105. +  comm - A communicator, use PETSC_COMM_SELF unless you know all ranks of another communicator will detect the error
    106. 

  106. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    107. 

  107. -  message - error message
    108. 

    109. 

  109.   Level: developer
    110. 

    111. 

  111.    Notes:
    112. 

  112.     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
    113. 

  113.     which is registered with MPI_Add_error_code() when PETSc is initialized.
    114. 

    115. 

  115. .seealso: SETERRQ(), CHKERRQ(), CHKERRMPI(), PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3()
    116. 

  116. M*/
    117. 

  117. #define SETERRMPI(comm,ierr,s) return (PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s),PETSC_MPI_ERROR_CODE)
    118. 

  118. 

  119. /*MC
    120. 

  120.    SETERRQ1 - Macro that is called when an error has been detected,
    121. 

    122. 

  122.    Synopsis:
    123. 

  123.    #include <petscsys.h>
    124. 

  124.    PetscErrorCode SETERRQ1(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg)
    125. 

    126. 

  126.    Collective
    127. 

    128. 

  128.    Input Parameters:
    129. 

  129. +  comm - A communicator, so that the error can be collective
    130. 

  130. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    131. 

  131. .  message - error message in the printf format
    132. 

  132. -  arg - argument (for example an integer, string or double)
    133. 

    134. 

  134.   Level: beginner
    135. 

    136. 

  136.    Notes:
    137. 

  137.     Once the error handler is called the calling function is then returned from with the given error code.
    138. 

    139. 

  139.    Experienced users can set the error handler with PetscPushErrorHandler().
    140. 

    141. 

  141. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ(), SETERRQ2(), SETERRQ3()
    142. 

  142. M*/
    143. 

  143. #define SETERRQ1(comm,ierr,s,a1) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1)
    144. 

  144. 

  145. /*MC
    146. 

  146.    SETERRQ2 - Macro that is called when an error has been detected,
    147. 

    148. 

  148.    Synopsis:
    149. 

  149.    #include <petscsys.h>
    150. 

  150.    PetscErrorCode SETERRQ2(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2)
    151. 

    152. 

  152.    Collective
    153. 

    154. 

  154.    Input Parameters:
    155. 

  155. +  comm - A communicator, so that the error can be collective
    156. 

  156. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    157. 

  157. .  message - error message in the printf format
    158. 

  158. .  arg1 - argument (for example an integer, string or double)
    159. 

  159. -  arg2 - argument (for example an integer, string or double)
    160. 

    161. 

  161.   Level: beginner
    162. 

    163. 

  163.    Notes:
    164. 

  164.     Once the error handler is called the calling function is then returned from with the given error code.
    165. 

    166. 

  166.    Experienced users can set the error handler with PetscPushErrorHandler().
    167. 

    168. 

  168. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ3()
    169. 

  169. M*/
    170. 

  170. #define SETERRQ2(comm,ierr,s,a1,a2) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2)
    171. 

  171. 

  172. /*MC
    173. 

  173.    SETERRQ3 - Macro that is called when an error has been detected,
    174. 

    175. 

  175.    Synopsis:
    176. 

  176.    #include <petscsys.h>
    177. 

  177.    PetscErrorCode SETERRQ3(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    178. 

    179. 

  179.    Collective
    180. 

    181. 

  181.    Input Parameters:
    182. 

  182. +  comm - A communicator, so that the error can be collective
    183. 

  183. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    184. 

  184. .  message - error message in the printf format
    185. 

  185. .  arg1 - argument (for example an integer, string or double)
    186. 

  186. .  arg2 - argument (for example an integer, string or double)
    187. 

  187. -  arg3 - argument (for example an integer, string or double)
    188. 

    189. 

  189.   Level: beginner
    190. 

    191. 

  191.    Notes:
    192. 

  192.     Once the error handler is called the calling function is then returned from with the given error code.
    193. 

    194. 

  194.     There are also versions for 4, 5, 6 and 7 arguments.
    195. 

    196. 

  196.    Experienced users can set the error handler with PetscPushErrorHandler().
    197. 

    198. 

  198. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    199. 

  199. M*/
    200. 

  200. #define SETERRQ3(comm,ierr,s,a1,a2,a3) return PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s,a1,a2,a3)
    201. 

  201. 

  202. /*MC
    203. 

  203.    SETERRQ4 - Macro that is called when an error has been detected,
    204. 

    205. 

  205.    Synopsis:
    206. 

  206.    #include <petscsys.h>
    207. 

  207.    PetscErrorCode SETERRQ4(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    208. 

    209. 

  209.    Collective
    210. 

    211. 

  211.    Input Parameters:
    212. 

  212. +  comm - A communicator, so that the error can be collective
    213. 

  213. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    214. 

  214. .  message - error message in the printf format
    215. 

  215. .  arg1 - argument (for example an integer, string or double)
    216. 

  216. .  arg2 - argument (for example an integer, string or double)
    217. 

  217. .  arg3 - argument (for example an integer, string or double)
    218. 

  218. -  arg4 - argument (for example an integer, string or double)
    219. 

    220. 

  220.   Level: beginner
    221. 

    222. 

  222.    Notes:
    223. 

  223.     Once the error handler is called the calling function is then returned from with the given error code.
    224. 

    225. 

  225.     There are also versions for 4, 5, 6 and 7 arguments.
    226. 

    227. 

  227.    Experienced users can set the error handler with PetscPushErrorHandler().
    228. 

    229. 

  229. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    230. 

  230. M*/
    231. 

  231. #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)
    232. 

  232. 

  233. /*MC
    234. 

  234.    SETERRQ5 - Macro that is called when an error has been detected,
    235. 

    236. 

  236.    Synopsis:
    237. 

  237.    #include <petscsys.h>
    238. 

  238.    PetscErrorCode SETERRQ5(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    239. 

    240. 

  240.    Collective
    241. 

    242. 

  242.    Input Parameters:
    243. 

  243. +  comm - A communicator, so that the error can be collective
    244. 

  244. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    245. 

  245. .  message - error message in the printf format
    246. 

  246. .  arg1 - argument (for example an integer, string or double)
    247. 

  247. .  arg2 - argument (for example an integer, string or double)
    248. 

  248. .  arg3 - argument (for example an integer, string or double)
    249. 

  249. .  arg4 - argument (for example an integer, string or double)
    250. 

  250. -  arg5 - argument (for example an integer, string or double)
    251. 

    252. 

  252.   Level: beginner
    253. 

    254. 

  254.    Notes:
    255. 

  255.     Once the error handler is called the calling function is then returned from with the given error code.
    256. 

    257. 

  257.     There are also versions for 4, 5, 6 and 7 arguments.
    258. 

    259. 

  259.    Experienced users can set the error handler with PetscPushErrorHandler().
    260. 

    261. 

  261. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    262. 

  262. M*/
    263. 

  263. #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)
    264. 

  264. 

  265. /*MC
    266. 

  266.    SETERRQ6 - Macro that is called when an error has been detected,
    267. 

    268. 

  268.    Synopsis:
    269. 

  269.    #include <petscsys.h>
    270. 

  270.    PetscErrorCode SETERRQ6(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    271. 

    272. 

  272.    Collective
    273. 

    274. 

  274.    Input Parameters:
    275. 

  275. +  comm - A communicator, so that the error can be collective
    276. 

  276. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    277. 

  277. .  message - error message in the printf format
    278. 

  278. .  arg1 - argument (for example an integer, string or double)
    279. 

  279. .  arg2 - argument (for example an integer, string or double)
    280. 

  280. .  arg3 - argument (for example an integer, string or double)
    281. 

  281. .  arg4 - argument (for example an integer, string or double)
    282. 

  282. .  arg5 - argument (for example an integer, string or double)
    283. 

  283. -  arg6 - argument (for example an integer, string or double)
    284. 

    285. 

  285.   Level: beginner
    286. 

    287. 

  287.    Notes:
    288. 

  288.     Once the error handler is called the calling function is then returned from with the given error code.
    289. 

    290. 

  290.     There are also versions for 4, 5, 6 and 7 arguments.
    291. 

    292. 

  292.    Experienced users can set the error handler with PetscPushErrorHandler().
    293. 

    294. 

  294. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    295. 

  295. M*/
    296. 

  296. #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)
    297. 

  297. 

  298. /*MC
    299. 

  299.    SETERRQ7 - Macro that is called when an error has been detected,
    300. 

    301. 

  301.    Synopsis:
    302. 

  302.    #include <petscsys.h>
    303. 

  303.    PetscErrorCode SETERRQ7(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    304. 

    305. 

  305.    Collective
    306. 

    307. 

  307.    Input Parameters:
    308. 

  308. +  comm - A communicator, so that the error can be collective
    309. 

  309. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    310. 

  310. .  message - error message in the printf format
    311. 

  311. .  arg1 - argument (for example an integer, string or double)
    312. 

  312. .  arg2 - argument (for example an integer, string or double)
    313. 

  313. .  arg3 - argument (for example an integer, string or double)
    314. 

  314. .  arg4 - argument (for example an integer, string or double)
    315. 

  315. .  arg5 - argument (for example an integer, string or double)
    316. 

  316. .  arg6 - argument (for example an integer, string or double)
    317. 

  317. -  arg7 - argument (for example an integer, string or double)
    318. 

    319. 

  319.   Level: beginner
    320. 

    321. 

  321.    Notes:
    322. 

  322.     Once the error handler is called the calling function is then returned from with the given error code.
    323. 

    324. 

  324.     There are also versions for 4, 5, 6 and 7 arguments.
    325. 

    326. 

  326.    Experienced users can set the error handler with PetscPushErrorHandler().
    327. 

    328. 

  328. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    329. 

  329. M*/
    330. 

  330. #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)
    331. 

  331. 

  332. /*MC
    333. 

  333.    SETERRQ8 - Macro that is called when an error has been detected,
    334. 

    335. 

  335.    Synopsis:
    336. 

  336.    #include <petscsys.h>
    337. 

  337.    PetscErrorCode SETERRQ8(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    338. 

    339. 

  339.    Collective
    340. 

    341. 

  341.    Input Parameters:
    342. 

  342. +  comm - A communicator, so that the error can be collective
    343. 

  343. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    344. 

  344. .  message - error message in the printf format
    345. 

  345. .  arg1 - argument (for example an integer, string or double)
    346. 

  346. .  arg2 - argument (for example an integer, string or double)
    347. 

  347. .  arg3 - argument (for example an integer, string or double)
    348. 

  348. .  arg4 - argument (for example an integer, string or double)
    349. 

  349. .  arg5 - argument (for example an integer, string or double)
    350. 

  350. .  arg6 - argument (for example an integer, string or double)
    351. 

  351. .  arg7 - argument (for example an integer, string or double)
    352. 

  352. -  arg8 - argument (for example an integer, string or double)
    353. 

    354. 

  354.   Level: beginner
    355. 

    356. 

  356.    Notes:
    357. 

  357.     Once the error handler is called the calling function is then returned from with the given error code.
    358. 

    359. 

  359.     There are also versions for 4, 5, 6 and 7 arguments.
    360. 

    361. 

  361.    Experienced users can set the error handler with PetscPushErrorHandler().
    362. 

    363. 

  363. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    364. 

  364. M*/
    365. 

  365. #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)
    366. 

  366. 

  367. /*MC
    368. 

  368.    SETERRQ9 - Macro that is called when an error has been detected,
    369. 

    370. 

  370.    Synopsis:
    371. 

  371.    #include <petscsys.h>
    372. 

  372.    PetscErrorCode SETERRQ9(MPI_Comm comm,PetscErrorCode ierr,char *formatmessage,arg1,arg2,arg3)
    373. 

    374. 

  374.    Collective
    375. 

    376. 

  376.    Input Parameters:
    377. 

  377. +  comm - A communicator, so that the error can be collective
    378. 

  378. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    379. 

  379. .  message - error message in the printf format
    380. 

  380. .  arg1 - argument (for example an integer, string or double)
    381. 

  381. .  arg2 - argument (for example an integer, string or double)
    382. 

  382. .  arg3 - argument (for example an integer, string or double)
    383. 

  383. .  arg4 - argument (for example an integer, string or double)
    384. 

  384. .  arg5 - argument (for example an integer, string or double)
    385. 

  385. .  arg6 - argument (for example an integer, string or double)
    386. 

  386. .  arg7 - argument (for example an integer, string or double)
    387. 

  387. .  arg8 - argument (for example an integer, string or double)
    388. 

  388. -  arg9 - argument (for example an integer, string or double)
    389. 

    390. 

  390.   Level: beginner
    391. 

    392. 

  392.    Notes:
    393. 

  393.     Once the error handler is called the calling function is then returned from with the given error code.
    394. 

    395. 

  395.     There are also versions for 0 to 9 arguments.
    396. 

    397. 

  397.    Experienced users can set the error handler with PetscPushErrorHandler().
    398. 

    399. 

  399. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    400. 

  400. M*/
    401. 

  401. #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)
    402. 

  402. 

  403. /*MC
    404. 

  404.    SETERRABORT - Macro that can be called when an error has been detected,
    405. 

    406. 

  406.    Synopsis:
    407. 

  407.    #include <petscsys.h>
    408. 

  408.    PetscErrorCode SETERRABORT(MPI_Comm comm,PetscErrorCode ierr,char *message)
    409. 

    410. 

  410.    Collective
    411. 

    412. 

  412.    Input Parameters:
    413. 

  413. +  comm - A communicator, so that the error can be collective
    414. 

  414. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    415. 

  415. -  message - error message in the printf format
    416. 

    417. 

  417.   Level: beginner
    418. 

    419. 

  419.    Notes:
    420. 

  420.     This function just calls MPI_Abort().
    421. 

    422. 

  422. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), CHKERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2()
    423. 

  423. M*/
    424. 

  424. #define SETERRABORT(comm,ierr,s) do {PetscError(comm,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_INITIAL,s);MPI_Abort(comm,ierr);} while (0)
    425. 

  425. 

  426. /*MC
    427. 

  427.    CHKERRQ - Checks error code, if non-zero it calls the error handler and then returns
    428. 

    429. 

  429.    Synopsis:
    430. 

  430.    #include <petscsys.h>
    431. 

  431.    PetscErrorCode CHKERRQ(PetscErrorCode ierr)
    432. 

    433. 

  433.    Not Collective
    434. 

    435. 

  435.    Input Parameters:
    436. 

  436. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    437. 

    438. 

  438.   Level: beginner
    439. 

    440. 

  440.    Notes:
    441. 

  441.     Once the error handler is called the calling function is then returned from with the given error code.
    442. 

    443. 

  443.     Experienced users can set the error handler with PetscPushErrorHandler().
    444. 

    445. 

  445.     CHKERRQ(ierr) is fundamentally a macro replacement for
    446. 

  446.          if (ierr) return(PetscError(...,ierr,...));
    447. 

    448. 

  448.     Although typical usage resembles "void CHKERRQ(PetscErrorCode)" as described above, for certain uses it is
    449. 

  449.     highly inappropriate to use it in this manner as it invokes return(PetscErrorCode). In particular,
    450. 

  450.     it cannot be used in functions which return(void) or any other datatype.  In these types of functions,
    451. 

  451.     you can use CHKERRV() which returns without an error code (bad idea since the error is ignored or
    452. 

  452.          if (ierr) {PetscError(....); return(YourReturnType);}
    453. 

  453.     where you may pass back a NULL to indicate an error. You can also call CHKERRABORT(comm,n) to have
    454. 

  454.     MPI_Abort() returned immediately.
    455. 

    456. 

  456. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2()
    457. 

  457. M*/
    458. 

  458. #define CHKERRQ(ierr)          do {if (PetscUnlikely(ierr)) return PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");} while (0)
    459. 

  459. #define CHKERRV(ierr)          do {if (PetscUnlikely(ierr)) {ierr = PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");return;}} while(0)
    460. 

  460. #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)
    461. 

  461. #define CHKERRCONTINUE(ierr)   do {if (PetscUnlikely(ierr)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_REPEAT," ");}} while (0)
    462. 

  462. 

  463. 

  464. /*MC
    465. 

  465.    CHKERRMPI - Checks error code, if non-zero it calls the error handler and then returns
    466. 

    467. 

  467.    Synopsis:
    468. 

  468.    #include <petscsys.h>
    469. 

  469.    PetscErrorCode CHKERRMPI(PetscErrorCode ierr)
    470. 

    471. 

  471.    Not Collective
    472. 

    473. 

  473.    Input Parameters:
    474. 

  474. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    475. 

    476. 

  476.   Level: developer
    477. 

    478. 

  478.    Notes:
    479. 

  479.     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
    480. 

  480.     which is registered with MPI_Add_error_code() when PETSc is initialized.
    481. 

    482. 

  482. .seealso: CHKERRQ(), PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ2()
    483. 

  483. M*/
    484. 

  484. #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)
    485. 

  485. 

  486. #ifdef PETSC_CLANGUAGE_CXX
    487. 

  487. 

  488. /*MC
    489. 

  489.    CHKERRXX - Checks error code, if non-zero it calls the C++ error handler which throws an exception
    490. 

    491. 

  491.    Synopsis:
    492. 

  492.    #include <petscsys.h>
    493. 

  493.    void CHKERRXX(PetscErrorCode ierr)
    494. 

    495. 

  495.    Not Collective
    496. 

    497. 

  497.    Input Parameters:
    498. 

  498. .  ierr - nonzero error code, see the list of standard error codes in include/petscerror.h
    499. 

    500. 

  500.   Level: beginner
    501. 

    502. 

  502.    Notes:
    503. 

  503.     Once the error handler throws a ??? exception.
    504. 

    505. 

  505.     You can use CHKERRV() which returns without an error code (bad idea since the error is ignored)
    506. 

  506.     or CHKERRABORT(comm,n) to have MPI_Abort() returned immediately.
    507. 

    508. 

  508. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKERRQ(), CHKMEMQ
    509. 

  509. M*/
    510. 

  510. #define CHKERRXX(ierr)  do {if (PetscUnlikely(ierr)) {PetscError(PETSC_COMM_SELF,__LINE__,PETSC_FUNCTION_NAME,__FILE__,ierr,PETSC_ERROR_IN_CXX,0);}} while(0)
    511. 

  511. 

  512. #endif
    513. 

  513. 

  514. #if defined(PETSC_HAVE_CUDA)
    515. 

  515. #define CHKERRCUDA(err)   do {if (PetscUnlikely(err)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"CUDA error %d",err);} while(0)
    516. 

  516. #define CHKERRCUBLAS(err) do {if (PetscUnlikely(err)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"CUBLAS error %d",err);} while(0)
    517. 

  517. #define CHKERRCUSOLVER(err) do {if (PetscUnlikely(err)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"CUSOLVER error %d",err);} while(0)
    518. 

  518. #endif
    519. 

  519. 

  520. /*MC
    521. 

  521.    CHKMEMQ - Checks the memory for corruption, calls error handler if any is detected
    522. 

    523. 

  523.    Synopsis:
    524. 

  524.    #include <petscsys.h>
    525. 

  525.    CHKMEMQ;
    526. 

    527. 

  527.    Not Collective
    528. 

    529. 

  529.   Level: beginner
    530. 

    531. 

  531.    Notes:
    532. 

  532.     We highly recommend using valgrind https://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind for finding memory problems. This is useful
    533. 

  533.     on systems that do not have valgrind, but much much less useful.
    534. 

    535. 

  535.     Must run with the option -malloc_debug (-malloc_test in debug mode; or if PetscMallocSetDebug() called) to enable this option
    536. 

    537. 

  537.     Once the error handler is called the calling function is then returned from with the given error code.
    538. 

    539. 

  539.     By defaults prints location where memory that is corrupted was allocated.
    540. 

    541. 

  541.     Use CHKMEMA for functions that return void
    542. 

    543. 

  543. .seealso: PetscTraceBackErrorHandler(), PetscPushErrorHandler(), PetscError(), SETERRQ(), CHKMEMQ, SETERRQ1(), SETERRQ2(), SETERRQ3(),
    544. 

  544.           PetscMallocValidate()
    545. 

  545. M*/
    546. 

  546. #define CHKMEMQ do {PetscErrorCode _7_ierr = PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__);CHKERRQ(_7_ierr);} while(0)
    547. 

  547. 

  548. #define CHKMEMA PetscMallocValidate(__LINE__,PETSC_FUNCTION_NAME,__FILE__)
    549. 

  549. 

  550. /*E
    551. 

  551.   PetscErrorType - passed to the PETSc error handling routines indicating if this is the first or a later call to the error handlers
    552. 

    553. 

  553.   Level: advanced
    554. 

    555. 

  555.   PETSC_ERROR_IN_CXX indicates the error was detected in C++ and an exception should be generated
    556. 

    557. 

  557.   Developer Notes:
    558. 

  558.     This is currently used to decide when to print the detailed information about the run in PetscTraceBackErrorHandler()
    559. 

    560. 

  560. .seealso: PetscError(), SETERRXX()
    561. 

  561. E*/
    562. 

  562. typedef enum {PETSC_ERROR_INITIAL=0,PETSC_ERROR_REPEAT=1,PETSC_ERROR_IN_CXX = 2} PetscErrorType;
    563. 

  563. 

  564. #if defined(__clang_analyzer__)
    565. 

  565. __attribute__((analyzer_noreturn))
    566. 

  566. #endif
    567. 

  567. PETSC_EXTERN PetscErrorCode PetscError(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,...);
    568. 

  568. 

  569. PETSC_EXTERN PetscErrorCode PetscErrorPrintfInitialize(void);
    570. 

  570. PETSC_EXTERN PetscErrorCode PetscErrorMessage(int,const char*[],char **);
    571. 

  571. PETSC_EXTERN PetscErrorCode PetscTraceBackErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    572. 

  572. PETSC_EXTERN PetscErrorCode PetscIgnoreErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    573. 

  573. PETSC_EXTERN PetscErrorCode PetscEmacsClientErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    574. 

  574. PETSC_EXTERN PetscErrorCode PetscMPIAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    575. 

  575. PETSC_EXTERN PetscErrorCode PetscAbortErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    576. 

  576. PETSC_EXTERN PetscErrorCode PetscAttachDebuggerErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    577. 

  577. PETSC_EXTERN PetscErrorCode PetscReturnErrorHandler(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*);
    578. 

  578. PETSC_EXTERN PetscErrorCode PetscPushErrorHandler(PetscErrorCode (*handler)(MPI_Comm,int,const char*,const char*,PetscErrorCode,PetscErrorType,const char*,void*),void*);
    579. 

  579. PETSC_EXTERN PetscErrorCode PetscPopErrorHandler(void);
    580. 

  580. PETSC_EXTERN PetscErrorCode PetscSignalHandlerDefault(int,void*);
    581. 

  581. PETSC_EXTERN PetscErrorCode PetscPushSignalHandler(PetscErrorCode (*)(int,void *),void*);
    582. 

  582. PETSC_EXTERN PetscErrorCode PetscPopSignalHandler(void);
    583. 

  583. PETSC_EXTERN PetscErrorCode PetscCheckPointerSetIntensity(PetscInt);
    584. 

  584. PETSC_EXTERN void PetscSignalSegvCheckPointer(void);
    585. 

  585. 

  586. /*MC
    587. 

  587.     PetscErrorPrintf - Prints error messages.
    588. 

    589. 

  589.    Synopsis:
    590. 

  590.     #include <petscsys.h>
    591. 

  591.      PetscErrorCode (*PetscErrorPrintf)(const char format[],...);
    592. 

    593. 

  593.     Not Collective
    594. 

    595. 

  595.     Input Parameters:
    596. 

  596. .   format - the usual printf() format string
    597. 

    598. 

  598.    Options Database Keys:
    599. 

  599. +    -error_output_stdout - cause error messages to be printed to stdout instead of the  (default) stderr
    600. 

  600. -    -error_output_none - to turn off all printing of error messages (does not change the way the error is handled.)
    601. 

    602. 

  602.    Notes:
    603. 

  603.     Use
    604. 

  604. $     PetscErrorPrintf = PetscErrorPrintfNone; to turn off all printing of error messages (does not change the way the
    605. 

  605. $                        error is handled.) and
    606. 

  606. $     PetscErrorPrintf = PetscErrorPrintfDefault; to turn it back on or you can use your own function
    607. 

    608. 

  608.           Use
    609. 

  609.      PETSC_STDERR = FILE* obtained from a file open etc. to have stderr printed to the file.
    610. 

  610.      PETSC_STDOUT = FILE* obtained from a file open etc. to have stdout printed to the file.
    611. 

    612. 

  612.           Use
    613. 

  613.       PetscPushErrorHandler() to provide your own error handler that determines what kind of messages to print
    614. 

    615. 

  615.    Level: developer
    616. 

    617. 

  617.     Fortran Note:
    618. 

  618.     This routine is not supported in Fortran.
    619. 

    620. 

    621. 

  621. .seealso: PetscFPrintf(), PetscSynchronizedPrintf(), PetscHelpPrintf(), PetscPrintf(), PetscPushErrorHandler(), PetscVFPrintf(), PetscHelpPrintf()
    622. 

  622. M*/
    623. 

  623. PETSC_EXTERN PetscErrorCode (*PetscErrorPrintf)(const char[],...);
    624. 

  624. 

  625. typedef enum {PETSC_FP_TRAP_OFF=0,PETSC_FP_TRAP_ON=1} PetscFPTrap;
    626. 

  626. PETSC_EXTERN PetscErrorCode PetscSetFPTrap(PetscFPTrap);
    627. 

  627. PETSC_EXTERN PetscErrorCode PetscFPTrapPush(PetscFPTrap);
    628. 

  628. PETSC_EXTERN PetscErrorCode PetscFPTrapPop(void);
    629. 

  629. 

  630. /*
    631. 

  631.       Allows the code to build a stack frame as it runs
    632. 

  632. */
    633. 

  633. 

  634. #define PETSCSTACKSIZE 64
    635. 

  635. 

  636. typedef struct  {
    637. 

  637.   const char      *function[PETSCSTACKSIZE];
    638. 

  638.   const char      *file[PETSCSTACKSIZE];
    639. 

  639.         int       line[PETSCSTACKSIZE];
    640. 

  640.         PetscBool petscroutine[PETSCSTACKSIZE];
    641. 

  641.         int       currentsize;
    642. 

  642.         int       hotdepth;
    643. 

  643. } PetscStack;
    644. 

  644. 

  645. PETSC_EXTERN PetscStack *petscstack;
    646. 

  646. 

  647. PetscErrorCode  PetscStackCopy(PetscStack*,PetscStack*);
    648. 

  648. PetscErrorCode  PetscStackPrint(PetscStack *,FILE*);
    649. 

  649. #if defined(PETSC_SERIALIZE_FUNCTIONS)
    650. 

  650. #include <petsc/private/petscfptimpl.h>
    651. 

  651. /*
    652. 

  652.    Registers the current function into the global function pointer to function name table
    653. 

    654. 

  654.    Have to fix this to handle errors but cannot return error since used in PETSC_VIEWER_DRAW_() etc
    655. 

  655. */
    656. 

  656. #define PetscRegister__FUNCT__() do { \
    657. 

  657.   static PetscBool __chked = PETSC_FALSE; \
    658. 

  658.   if (!__chked) {\
    659. 

  659.   void *ptr; PetscDLSym(NULL,PETSC_FUNCTION_NAME,&ptr);\
    660. 

  660.   __chked = PETSC_TRUE;\
    661. 

  661.   }} while (0)
    662. 

  662. #else
    663. 

  663. #define PetscRegister__FUNCT__()
    664. 

  664. #endif
    665. 

  665. 

  666. #if defined(PETSC_USE_DEBUG)
    667. 

  667. PETSC_STATIC_INLINE PetscBool PetscStackActive(void)
    668. 

  668. {
    669. 

  669.   return(petscstack ? PETSC_TRUE : PETSC_FALSE);
    670. 

  670. }
    671. 

  671. 

  672. /* Stack handling is based on the following two "NoCheck" macros.  These should only be called directly by other error
    673. 

  673.  * handling macros.  We record the line of the call, which may or may not be the location of the definition.  But is at
    674. 

  674.  * least more useful than "unknown" because it can distinguish multiple calls from the same function.
    675. 

  675.  */
    676. 

  676. 

  677. #define PetscStackPushNoCheck(funct,petsc_routine,hot)                        \
    678. 

  678.   do {                                                                        \
    679. 

  679.     PetscStackSAWsTakeAccess();                                                \
    680. 

  680.     if (petscstack && (petscstack->currentsize < PETSCSTACKSIZE)) {         \
    681. 

  681.       petscstack->function[petscstack->currentsize]  = funct;               \
    682. 

  682.       petscstack->file[petscstack->currentsize]      = __FILE__;            \
    683. 

  683.       petscstack->line[petscstack->currentsize]      = __LINE__;            \
    684. 

  684.       petscstack->petscroutine[petscstack->currentsize] = petsc_routine;    \
    685. 

  685.       petscstack->currentsize++;                                             \
    686. 

  686.     }                                                                         \
    687. 

  687.     if (petscstack) {                                                        \
    688. 

  688.       petscstack->hotdepth += (hot || petscstack->hotdepth);                \
    689. 

  689.     }                                                                         \
    690. 

  690.     PetscStackSAWsGrantAccess();                                               \
    691. 

  691.   } while (0)
    692. 

  692. 

  693. #define PetscStackPopNoCheck                                            \
    694. 

  694.   do {                                                                  \
    695. 

  695.     PetscStackSAWsTakeAccess();                                          \
    696. 

  696.     if (petscstack && petscstack->currentsize > 0) {                  \
    697. 

  697.       petscstack->currentsize--;                                       \
    698. 

  698.       petscstack->function[petscstack->currentsize]  = 0;             \
    699. 

  699.       petscstack->file[petscstack->currentsize]      = 0;             \
    700. 

  700.       petscstack->line[petscstack->currentsize]      = 0;             \
    701. 

  701.       petscstack->petscroutine[petscstack->currentsize] = PETSC_FALSE;\
    702. 

  702.     }                                                                   \
    703. 

  703.     if (petscstack) {                                                  \
    704. 

  704.       petscstack->hotdepth = PetscMax(petscstack->hotdepth-1,0);      \
    705. 

  705.     }                                                                   \
    706. 

  706.     PetscStackSAWsGrantAccess();                                         \
    707. 

  707.   } while (0)
    708. 

  708. 

  709. /*MC
    710. 

  710.    PetscFunctionBegin - First executable line of each PETSc function,  used for error handling. Final
    711. 

  711.       line of PETSc functions should be PetscFunctionReturn(0);
    712. 

    713. 

  713.    Synopsis:
    714. 

  714.    #include <petscsys.h>
    715. 

  715.    void PetscFunctionBegin;
    716. 

    717. 

  717.    Not Collective
    718. 

    719. 

  719.    Usage:
    720. 

  720. .vb
    721. 

  721.      int something;
    722. 

    723. 

  723.      PetscFunctionBegin;
    724. 

  724. .ve
    725. 

    726. 

  726.    Notes:
    727. 

  727.      Use PetscFunctionBeginUser for application codes.
    728. 

    729. 

  729.      Not available in Fortran
    730. 

    731. 

  731.    Level: developer
    732. 

    733. 

  733. .seealso: PetscFunctionReturn(), PetscFunctionBeginHot(), PetscFunctionBeginUser()
    734. 

    735. 

  735. M*/
    736. 

  736. #define PetscFunctionBegin do {                                        \
    737. 

  737.     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_FALSE); \
    738. 

  738.     PetscRegister__FUNCT__();                                          \
    739. 

  739.   } while (0)
    740. 

  740. 

  741. /*MC
    742. 

  742.    PetscFunctionBeginHot - Substitute for PetscFunctionBegin to be used in functions that are called in
    743. 

  743.    performance-critical circumstances.  Use of this function allows for lighter profiling by default.
    744. 

    745. 

  745.    Synopsis:
    746. 

  746.    #include <petscsys.h>
    747. 

  747.    void PetscFunctionBeginHot;
    748. 

    749. 

  749.    Not Collective
    750. 

    751. 

  751.    Usage:
    752. 

  752. .vb
    753. 

  753.      int something;
    754. 

    755. 

  755.      PetscFunctionBeginHot;
    756. 

  756. .ve
    757. 

    758. 

  758.    Notes:
    759. 

  759.      Not available in Fortran
    760. 

    761. 

  761.    Level: developer
    762. 

    763. 

  763. .seealso: PetscFunctionBegin, PetscFunctionReturn()
    764. 

    765. 

  765. M*/
    766. 

  766. #define PetscFunctionBeginHot do {                                     \
    767. 

  767.     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_TRUE,PETSC_TRUE);  \
    768. 

  768.     PetscRegister__FUNCT__();                                          \
    769. 

  769.   } while (0)
    770. 

  770. 

  771. /*MC
    772. 

  772.    PetscFunctionBeginUser - First executable line of user provided PETSc routine
    773. 

    774. 

  774.    Synopsis:
    775. 

  775.    #include <petscsys.h>
    776. 

  776.    void PetscFunctionBeginUser;
    777. 

    778. 

  778.    Not Collective
    779. 

    780. 

  780.    Usage:
    781. 

  781. .vb
    782. 

  782.      int something;
    783. 

    784. 

  784.      PetscFunctionBeginUser;
    785. 

  785. .ve
    786. 

    787. 

  787.    Notes:
    788. 

  788.       Final line of PETSc functions should be PetscFunctionReturn(0) except for main().
    789. 

    790. 

  790.       Not available in Fortran
    791. 

    792. 

  792.       This is identical to PetscFunctionBegin except it labels the routine as a user
    793. 

  793.       routine instead of as a PETSc library routine.
    794. 

    795. 

  795.    Level: intermediate
    796. 

    797. 

  797. .seealso: PetscFunctionReturn(), PetscFunctionBegin, PetscFunctionBeginHot
    798. 

    799. 

  799. M*/
    800. 

  800. #define PetscFunctionBeginUser                                          \
    801. 

  801.   do {                                                                  \
    802. 

  802.     PetscStackPushNoCheck(PETSC_FUNCTION_NAME,PETSC_FALSE,PETSC_FALSE); \
    803. 

  803.     PetscRegister__FUNCT__();                                           \
    804. 

  804.   } while (0)
    805. 

  805. 

  806. 

  807. #define PetscStackPush(n) \
    808. 

  808.   do {                                                                  \
    809. 

  809.     PetscStackPushNoCheck(n,PETSC_FALSE,PETSC_FALSE);                   \
    810. 

  810.     CHKMEMQ;                                                            \
    811. 

  811.   } while (0)
    812. 

  812. 

  813. #define PetscStackPop                           \
    814. 

  814.     do {                                        \
    815. 

  815.       CHKMEMQ;                                  \
    816. 

  816.       PetscStackPopNoCheck;                     \
    817. 

  817.     } while (0)
    818. 

  818. 

  819. /*MC
    820. 

  820.    PetscFunctionReturn - Last executable line of each PETSc function
    821. 

  821.         used for error handling. Replaces return()
    822. 

    823. 

  823.    Synopsis:
    824. 

  824.    #include <petscsys.h>
    825. 

  825.    void PetscFunctionReturn(0);
    826. 

    827. 

  827.    Not Collective
    828. 

    829. 

  829.    Usage:
    830. 

  830. .vb
    831. 

  831.     ....
    832. 

  832.      PetscFunctionReturn(0);
    833. 

  833.    }
    834. 

  834. .ve
    835. 

    836. 

  836.    Notes:
    837. 

  837.      Not available in Fortran
    838. 

    839. 

  839.    Level: developer
    840. 

    841. 

  841. .seealso: PetscFunctionBegin()
    842. 

    843. 

  843. M*/
    844. 

  844. #define PetscFunctionReturn(a) \
    845. 

  845.   do {                                                                \
    846. 

  846.     PetscStackPopNoCheck;                                             \
    847. 

  847.     return(a);} while (0)
    848. 

  848. 

  849. #define PetscFunctionReturnVoid() \
    850. 

  850.   do {                                                                \
    851. 

  851.     PetscStackPopNoCheck;                                             \
    852. 

  852.     return;} while (0)
    853. 

  853. 

  854. #else
    855. 

  855. 

  856. PETSC_STATIC_INLINE PetscBool PetscStackActive(void) {return PETSC_FALSE;}
    857. 

  857. #define PetscStackPushNoCheck(funct,petsc_routine,hot) do {} while (0)
    858. 

  858. #define PetscStackPopNoCheck                           do {} while (0)
    859. 

  859. #define PetscFunctionBegin
    860. 

  860. #define PetscFunctionBeginUser
    861. 

  861. #define PetscFunctionBeginHot
    862. 

  862. #define PetscFunctionReturn(a)    return(a)
    863. 

  863. #define PetscFunctionReturnVoid() return
    864. 

  864. #define PetscStackPop             CHKMEMQ
    865. 

  865. #define PetscStackPush(f)         CHKMEMQ
    866. 

  866. 

  867. #endif
    868. 

  868. 

  869. /*
    870. 

  870.     PetscStackCall - Calls an external library routine or user function after pushing the name of the routine on the stack.
    871. 

    872. 

  872.    Input Parameters:
    873. 

  873. +   name - string that gives the name of the function being called
    874. 

  874. -   routine - actual call to the routine, including ierr = and CHKERRQ(ierr);
    875. 

    876. 

  876.    Note: Often one should use PetscStackCallStandard() instead. This routine is intended for external library routines that DO NOT return error codes
    877. 

    878. 

  878.    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.
    879. 

    880. 

    881. 

    882. 

  882. */
    883. 

  883. #define PetscStackCall(name,routine) do { PetscStackPush(name);routine;PetscStackPop; } while(0)
    884. 

  884. 

  885. /*
    886. 

  886.     PetscStackCallStandard - Calls an external library routine after pushing the name of the routine on the stack.
    887. 

    888. 

  888.    Input Parameters:
    889. 

  889. +   func-  name of the routine
    890. 

  890. -   args - arguments to the routine surrounded by ()
    891. 

    892. 

  892.    Notes:
    893. 

  893.     This is intended for external package routines that return error codes. Use PetscStackCall() for those that do not.
    894. 

    895. 

  895.    Developer Note: this is so that when an external packge routine results in a crash or corrupts memory, they get blamed instead of PETSc.
    896. 

    897. 

  897. */
    898. 

  898. #define PetscStackCallStandard(func,args) do {                                                            \
    899. 

  899.     PetscErrorCode __ierr;                                                                                \
    900. 

  900.     PetscStackPush(#func);                                                                                \
    901. 

  901.     __ierr = func args;                                                                                   \
    902. 

  902.     PetscStackPop;                                                                                        \
    903. 

  903.     if (__ierr) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error in %s(): error code %d",#func,(int)__ierr); \
    904. 

  904.   } while (0)
    905. 

  905. 

  906. PETSC_EXTERN PetscErrorCode PetscStackCreate(void);
    907. 

  907. PETSC_EXTERN PetscErrorCode PetscStackView(FILE*);
    908. 

  908. PETSC_EXTERN PetscErrorCode PetscStackDestroy(void);
    909. 

  909. 

  910. #endif
    911.