How is it possible???
double *bds = reader->GetOutput()->GetBounds();
printf("%lf %lf %lf %lf\n", bds[0], bds[1], bds[2], bds[3]);
-> 0.00 69.00 0.00 69.00
printf("point1[0]=%f point1[1]=%f point1[2]=%f\n", point1[0], point1[1],
point1[2]);
-> point1[0]=58.196274 point1[1
I suggest you attach a debugger to the code and step through. Then
you'll be able to answer your question.
Since the velocity field of your data may be high, it is entirely
possible that the computed next position is outside of the bounds.
Perhaps reducing the step size will help. Try even Rk
Hi Fred:
As JB suggested, please try using RK45, which employs an adaptive step
size that is insensitive to the velocity magnitude and hence is capable of
capturing accurate shape of the curve (flow line). Using RK4 (still adopting
a fixed step size) does not guarantee a sufficient / desired a
Hello,
My main concern is why does the integrator return 0 instead of OUT_OF_DOMAIN?
Regards
Date: Mon, 17 Aug 2009 10:45:31 -0400
Subject: Re: [Paraview] RungeKutta4 implementation issue
From: zhanping@kitware.com
To: biddi...@cscs.ch
CC: stan1...@hotmail.fr; paraview@paraview.org
Hi Fred
saved - the courant condition is
no longer satisfied etc etc).
does that help?
JB
Regards
Date: Mon, 17 Aug 2009 10:45:31 -0400
Subject: Re: [Paraview] RungeKutta4 implementation issue
From: zhanping@kitware.com
To: biddi...@cscs.ch
CC: stan1...@hotmail.fr; paraview@paraview.org
Hi
Huh? I get over 100 emails from that list every week. There is usually more
activity on the VTK mailing list than the ParaView mailing list (according to
markmail).
-Ken
On 8/16/09 4:03 PM, "Fred Fred" wrote:
PS: This message should probably be sent to the VTK mailing list but it seems
th
...@skippingmouse.co.uk
To: stan1...@hotmail.fr
CC: zhanping@kitware.com; biddi...@cscs.ch; paraview@paraview.org
Subject: Re: [Paraview] RungeKutta4 implementation issue
My main concern is why does the integrator return 0
instead of OUT_OF_DOMAIN?
Because the (current position inside