[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] RK8pd nan?

From: Zine Smith
Subject: [Help-gsl] RK8pd nan?
Date: Wed, 13 Apr 2005 08:54:41 -0700 (PDT)

I have struggled with receiving a nan inside rk8 that
cause my integration to bomb. It occurs during the
first call to 'evolve'.  Using GDB, I've found the
error to occur during the K4 step on the first element
of the state (dim=6).  Here is the output from GDB for
the elements in the K4 i=0 step.

(gdb) printf  "%g + %g * (%g * %g + %g * %g + %g *
1.37814e+07 + 0.1 * (0.3125 * 1 + -1.17188 * 0.982509
+ 1.17188 * 0.973763)

Given these numbers I would not expect a nan, but when
I display the ytmp[0] element, I get the following in

(gdb) printf "%g\n",ytmp[0]

Any help would be greatly appreciated.  I will be more
than happy to provide any code or information that is
needed to help with this issue.

Zine Smith

Do you Yahoo!? 
Yahoo! Small Business - Try our new resources site!

reply via email to

[Prev in Thread] Current Thread [Next in Thread]