<div>Have you taken a look at numpy? [1] It was written for exactly this kind of usage.<br clear="all"></div><div><br></div><div>~/santa<br></div><div><br></div><div>[1] <a href="http://numpy.scipy.org/">http://numpy.scipy.org/</a></div>


<br><br><div class="gmail_quote">On Fri, Mar 4, 2011 at 2:32 PM, Jon Herman <span dir="ltr"><<a href="mailto:jfc.herman@gmail.com">jfc.herman@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

Hello all,<br><br>I am new to the Python language and writing a Runge-Kutta-Fellberg 7(8) integrator in Python, which requires an extreme numerical precision for my particular application. Unfortunately, I can not seem to attain it.<br>


The interesting part is if I take my exact code and translate it to Matlab code (so I use the exact same process and numbers), I get a far superior precision (the one I am expecting, in fact). This leads me to think I need to call a certain command in my Python script in order to make sure no truncation errors are building up over my integration.<br>


<br>Has anyone had similar problems? Is there a difference between how Matlab and Python store numbers, and if so how do I make Python more accurate?<br><br>I know there is a lot of packages out there, but this in fact overwhelmed me a little bit and seems to prevent me from finding the answer to my question, so I'm hoping someone with more experience will be able to enlighten me!<br>


<br>Best regards,<br><font color="#888888"><br>Jon<br>
</font><br>--<br>
<a href="http://mail.python.org/mailman/listinfo/python-list" target="_blank">http://mail.python.org/mailman/listinfo/python-list</a><br>
<br></blockquote></div><br>