<div dir="ltr">Dear neks<div><br></div><div>I am still stuck with the problem of not being able to save the perturbation solution using out_pert(). Do you have any suggestions ?</div><div><br></div><div>Thanks</div><div>praveen</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 9, 2014 at 6:57 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>