[Nek5000-users] Saving perturbation solution

[nek5000-users at lists.mcs.anl.gov]

Hello neks

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

       ifxyo = .false.

       if(istep.eq.0)then

          ifxyo = .true.

          call out_pert()

       else if(mod(istep,iostep).eq.0)then

          call out_pert()

       endif


then I get the crash

         1 U-PRES gmres:     25  3.8311E-12  5.0000E-12  1.6258E-08
3.7787E-03  1.1255E-02

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

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

 2.000000E-03     1  1.78124472E-03  1.78124472E-03  1.00000713E+00
1.00000686E+00  2.63413513E-07  3.83894325E-02 egn


Program received signal SIGFPE: Floating-point exception - erroneous
arithmetic operation.


Backtrace for this error:

#0  0x10fa6770d

#1  0x10fa67c1b

#2  0x7fff872745a9

#3  0x1042ad8fe

#4  0x1042aeafb

#5  0x1044e705c


Thanks

praveen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140909/20ac2f50/attachment.html>