[Nek5000-users] Saving perturbation solution

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

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, <nek5000-users at lists.mcs.anl.gov> 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
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140926/856a97b0/attachment.html>