<br><br><div class="gmail_quote">On Tue, Nov 13, 2012 at 10:43 PM, Hugh Aguilar <span dir="ltr"><<a href="mailto:hughaguilar96@yahoo.com" target="_blank">hughaguilar96@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div style="font-size:12pt;font-family:times new roman,new york,times,serif"><div class="im"><div>> <span style="font-size:12pt">We are doing numerical </span><span style="font-size:12pt">integration of celestial bodies over large periods of time </span><span style="font-size:12pt">(100 years is a norm).</span></div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif"><br></div></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
I'm new to Scheme, so I may be totally wrong about this --- but, isn't a numerical program like this exactly what Scheme is *not* designed for?</div></div></div></blockquote><div><br></div><div>The Supercomputer Toolkit ( <a href="http://www.hpl.hp.com/techreports/94/HPL-94-30.html">http://www.hpl.hp.com/techreports/94/HPL-94-30.html</a> ) was</div>
<div>designed for doing numerical integration of celestial bodies over large periods of time. </div><div><br></div><div>"The Toolkit's compiler uses a novel strategy based upon partial evaluation [7, 9]. This exploits the data-independence of typical numerical</div>
<div>algorithms to generate exceptionally efficient object code from source programs that are expressed in terms of highly abstract components written in</div><div>the Scheme dialect of Lisp [14]."</div><div><br></div>
<div>"The integrator and the force law were written as high level Scheme programs. The accumulation of position was implemented in</div><div>quad precision (128 bits), and the required quad precision operators were</div>
<div>written in Scheme."</div><div><br></div><div>"In hindsight, the use of quad precision appears to have been overly conservative for this</div><div>problem"</div><div><br></div></div>-- <br>~jrm<br>