<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Aptos;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:Monaco;
panose-1:0 0 0 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:12.0pt;
font-family:"Aptos",sans-serif;
mso-ligatures:standardcontextual;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Aptos",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#467886" vlink="#96607D" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">After updating to PETSc 3.23 (from a quite old version, ~3.8), I found that my old way of exception handling not working any more, on my Mac.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I would like to learn the proper way to handle exceptions in PETSc code.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Here is my pseudo code:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">================================================================<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">PetscErrorCode SNESFormFunction(SNES, Vec, Vec, void*);<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">double my_function();<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">================================================================<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">double my_function()<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">{<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> if (all_good)<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">return 1;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> else<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> {<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">throw 199;<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">return 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> }<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">}<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">PetscErrorCode SNESFormFunction(SNES snes, Vec u, Vec r, void* ctx)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">{<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> double my_value = my_function();<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> // compute residuals<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">}<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">int main(int argc, char **argv)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">{<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> Initialize_PETSc();<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> double dt = 1, dt_min = 0.001;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> while (dt > dt_min)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> {<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">try<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">{<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco"> SNESSolve(AppCtx.snes, NULL, AppCtx.u);<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">}<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">catch (int err)<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">{<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco"> dt *= 0.5;<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:19.5pt"><span style="font-size:8.0pt;font-family:Monaco">}<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> }<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco"> PetscFinalize();<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">}<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:Monaco">================================================================<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">This piece of logic used to work well, but now giving me the following error:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Current time (the starting time of this time step) = 0.. <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> NL step = 0, SNES Function norm = 8.60984E+03<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;background:yellow;mso-highlight:yellow">libc++abi: terminating due to uncaught exception of type int<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;background:yellow;mso-highlight:yellow">Abort trap: 6</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Q1: why this exception catch logic not working any more?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Q2: is there any good example of PETSc exception handling I can follow?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Best,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-Ling<o:p></o:p></span></p>
</div>
</body>
</html>