Dear neks I am still stuck with the problem of not being able to save the perturbation solution using out_pert(). Do you have any suggestions ? Thanks praveen On Tue, Sep 9, 2014 at 6:57 PM, <[email protected]> wrote:
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
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users