<div>Hi Fred:</div>
<div> </div>
<div> 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 acurracy of numerical integration (though certainly it is better than RK2) since it is stil sensitive to the velocity magnitude.</div>
<div> </div>
<div> Recently we improved the flow line integration accuracy issue. For any problems, please feel free to let us know.</div>
<div> </div>
<div> Thanks.</div>
<div> </div>
<div> -Zhanping<br></div>
<div> -- <br>Zhanping Liu, PhD<br>Kitware, Inc.<br>28 Corporate Drive<br>Clifton Park, NY 12065-8662<br>Phone: 518-371-3971 x 138<br><a href="http://www.zhanpingliu.org">http://www.zhanpingliu.org</a><br><br></div>
<div class="gmail_quote">On Mon, Aug 17, 2009 at 5:06 AM, John Biddiscombe <span dir="ltr"><<a href="mailto:biddisco@cscs.ch">biddisco@cscs.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div text="#000000" bgcolor="#ffffff">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 type="cite">
<div>
<div></div>
<div class="h5">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 href="http://www.microsoft.com/windows/windowslive/default.aspx" target="_blank">avec vos proches</a> </div></div><pre><hr width="90%" size="4">
_______________________________________________
Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a>
Follow this link to subscribe/unsubscribe:
<a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a>
</pre></blockquote><br><font color="#888888"><br><pre cols="78">--
John Biddiscombe, email:biddisco @ <a href="http://cscs.ch/" target="_blank">cscs.ch</a>
<a href="http://www.cscs.ch/" target="_blank">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></font></div><br>_______________________________________________<br>Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a><br>
<br>Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br><br>Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>Follow this link to subscribe/unsubscribe:<br><a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br><br></blockquote></div><br><br clear="all">
<div></div><br>