<div dir="ltr">Hello neks<div><br></div><div>I am running the perturbation code by modifying the os7000 example. The code runs. But if I try to save the perturbation solution like this</div><div><br></div><div><p style="margin:0px;font-size:11px;font-family:Menlo">       ifxyo <span style="color:rgb(52,189,38)">=</span> <span style="color:rgb(195,55,32)">.false.</span></p>
<p style="margin:0px;font-size:11px;font-family:Menlo">       <span style="color:rgb(52,189,38)">if</span>(istep<span style="color:rgb(52,189,38)">.eq.</span><span style="color:rgb(195,55,32)">0</span>)<span style="color:rgb(52,189,38)">then</span></p>
<p style="margin:0px;font-size:11px;font-family:Menlo">          ifxyo <span style="color:rgb(52,189,38)">=</span> <span style="color:rgb(195,55,32)">.true.</span></p>
<p style="margin:0px;font-size:11px;font-family:Menlo">          <span style="color:rgb(52,187,199)">call</span> out_pert()</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">       <span style="color:rgb(52,189,38)">else</span> <span style="color:rgb(52,189,38)">if</span>(<span style="color:rgb(52,187,199)">mod</span>(istep,iostep)<span style="color:rgb(52,189,38)">.eq.</span><span style="color:rgb(195,55,32)">0</span>)<span style="color:rgb(52,189,38)">then</span></p>
<p style="margin:0px;font-size:11px;font-family:Menlo">          <span style="color:rgb(52,187,199)">call</span> out_pert()</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">       <span style="color:rgb(52,189,38)">endif</span></p><p style="margin:0px;font-size:11px;font-family:Menlo"><span style="color:rgb(52,189,38)"><br></span></p>then I get the crash</div><div><br></div><div><p style="margin:0px;font-size:11px;font-family:Menlo">         1 U-PRES gmres:     25  3.8311E-12  5.0000E-12  1.6258E-08  3.7787E-03  1.1255E-02</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">filt amp 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0500</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">filt trn 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9500</p>
<p style="margin:0px;font-size:11px;font-family:Menlo"> 2.000000E-03     1  1.78124472E-03  1.78124472E-03  1.00000713E+00  1.00000686E+00  2.63413513E-07  3.83894325E-02 egn</p>
<p style="margin:0px;font-size:11px;font-family:Menlo;min-height:13px"><br></p>
<p style="margin:0px;font-size:11px;font-family:Menlo">Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.</p>
<p style="margin:0px;font-size:11px;font-family:Menlo;min-height:13px"><br></p>
<p style="margin:0px;font-size:11px;font-family:Menlo">Backtrace for this error:</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">#0  0x10fa6770d</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">#1  0x10fa67c1b</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">#2  0x7fff872745a9</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">#3  0x1042ad8fe</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">#4  0x1042aeafb</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">#5  0x1044e705c</p><p style="margin:0px;font-size:11px;font-family:Menlo"><br></p>Thanks<br><br>praveen</div></div>