<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
I suggest you attach a debugger to the code and step through. Then
you'll be able to answer your question. <br>
<br>
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 Rk4.5 with adaptive
steps, this may help if the velocity gradient at the edge of the domain
is very large..<br>
<br>
JB<br>
<br>
<blockquote cite="mid:BAY109-W39D69AB1A871CCE56004B6B2000@phx.gbl"
type="cite">
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>How is it possible???<br>
<br>
double *bds = reader->GetOutput()->GetBounds();<br>
printf("%lf %lf %lf %lf\n", bds[0], bds[1], bds[2], bds[3]);<br>
<br>
-> 0.000000 69.000000 0.000000 69.000000<br>
<br>
printf("point1[0]=%f point1[1]=%f point1[2]=%f\n", point1[0],
point1[1], point1[2]);<br>
<br>
-> point1[0]=58.196274 point1[1]=68.862898 point1[2]=0.000000<br>
<br>
if (integrator->ComputeNextStep(point1, point2, 0, delT.Interval,
stepTaken, minStep, maxStep, this->MaximumError, error) != 0)<br>
break;<br>
printf("point2[0]=%f point2[1]=%f point2[2]=%f\n", point2[0],
point2[1], point2[2]);<br>
<br>
-> point2[0]=59.002164 point2[1]=69.454414 point2[2]=0.000000<br>
<br>
PS: This message should probably be sent to the VTK mailing list but it
seems that it is no longer so active...<br>
<br>
<br>
<hr>Découvrez toutes les possibilités de communication <a
moz-do-not-send="true"
href="http://www.microsoft.com/windows/windowslive/default.aspx"
target="_new">avec vos proches</a>
<pre wrap="">
<hr size="4" width="90%">
_______________________________________________
Powered by <a class="moz-txt-link-abbreviated" href="http://www.kitware.com">www.kitware.com</a>
Visit other Kitware open-source projects at <a class="moz-txt-link-freetext" href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a>
Please keep messages on-topic and check the ParaView Wiki at: <a class="moz-txt-link-freetext" href="http://paraview.org/Wiki/ParaView">http://paraview.org/Wiki/ParaView</a>
Follow this link to subscribe/unsubscribe:
<a class="moz-txt-link-freetext" href="http://www.paraview.org/mailman/listinfo/paraview">http://www.paraview.org/mailman/listinfo/paraview</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="78">--
John Biddiscombe, email:biddisco @ cscs.ch
<a class="moz-txt-link-freetext" href="http://www.cscs.ch/">http://www.cscs.ch/</a>
CSCS, Swiss National Supercomputing Centre | Tel: +41 (91) 610.82.07
Via Cantonale, 6928 Manno, Switzerland | Fax: +41 (91) 610.82.82</pre>
</body>
</html>