<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Abhinav -- it is still done today.   For Intel's MKL we must have a team of programmers that specialize in writing math at the lowest levels.  DEC, CDC, Cray, IBM did the same thing back in the day.   Check out:  <a href="https://software.intel.com/en-us/mkl">Intel Math Kernel Library (<i>a.k.a.</i> MKL)</a>.  </div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 19, 2019 at 2:34 PM Abhinav Rajagopalan <<a href="mailto:abhinavrajagopalan@gmail.com">abhinavrajagopalan@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-family:"courier new",monospace;font-size:small"><div class="gmail_default" style="color:rgb(0,0,0)">Forgive me for both hijacking this thread, and to address my amateurish gnawing concern, but how was it be possible to write differential/integral equations at an assembly/machine level at the time, especially in machines such as the PDP-7 and such which had IIRC just 16 instructions and operated on the basis of mere words, especially the floating point math being done. Surmising from some personal experience that writing mathematical programs is hard even now, although there exist certain functional paradigms, and specialised environments such as MATLAB or Mathematica. The complexity seems to remain the same if not more now, due to the vast oodles of data to handle stemming from the nature of the world.</div><div class="gmail_default" style="color:rgb(0,0,0)"><br></div><div class="gmail_default" style="color:rgb(0,0,0)">Were they loaded as just words as any other instruction or were there separate coprocessors that did the number crunching? I'm guessing Fortran-ish kind of implementations were done, but the hardware level computation itself I just can't process.</div><div class="gmail_default" style="color:rgb(0,0,0)"><br></div><div class="gmail_default" style="color:rgb(0,0,0)">It just blows my mind now thinking backwards in terms of those monster machines being loaded with trails of paper tape instructions to play Space Travel. Being born in the late 90's doesn't help me too.</div><div class="gmail_default" style="color:rgb(0,0,0)"><br></div><div class="gmail_default" style="color:rgb(0,0,0)">Also, on a related note, don't know if you've watched the <a href="https://youtu.be/EY6q5dv_B-o" target="_blank">interview</a> of Ken done by Brian at the Vintage Comptuer Federation 2019, there might be a few surprises lurking around the middle of that when they discuss pipes and grep. </div><div class="gmail_default" style="color:rgb(0,0,0)"><br></div><div class="gmail_default" style="color:rgb(0,0,0)">Thank you!</div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 19, 2019 at 8:11 PM Doug McIlroy <<a href="mailto:doug@cs.dartmouth.edu" target="_blank">doug@cs.dartmouth.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">I was about to add a footnote to history about<br>
how the broad interests and collegiality of<br>
Bell Labs staff made Space Travel work, when<br>
I saw that Ken beat me to telling how he got<br>
help from another Turing Award winner.<br>
<br>
> while writing "space travel,"<br>
> i could not get the space ship integration<br>
> around a planet to keep from either gaining or<br>
> losing energy due to floating point errors.<br>
> i asked dick hamming if he could help. after<br>
> a couple hours, he came back with a formula.<br>
> i tried it and it worked perfectly. it was some<br>
> weird simple double integration that self<br>
> corrected for fp round off. as near as i can<br>
> ascertain, the formula was never published<br>
> and no one i have asked (including me) has<br>
> been able to recreate it.<br>
<br>
If I remember correctly, the cause of Ken's<br>
difficulty was not roundoff error. It<br>
was discretization error in the integration<br>
formula--probably f(t+dt)=f(t)+f'(t)dt.<br>
Dick saw that the formula did not conserve<br>
energy and found an alternative that did.<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><span style="font-family:"courier new",monospace"><br></span></div><div><span style="font-family:"courier new",monospace"><span style="color:rgb(0,0,0)">Abhinav Rajagopalan</span></span></div><div><span><span style="color:rgb(68,68,68)"><br></span></span></div><div><font style="background-color:rgb(255,0,0)" face="tahoma, sans-serif" color="#666666"><br></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div>