Article View: comp.lang.fortran
Article #6531Re: Precision of floating point multiply
From: Glen Herrmannsfeldt
Date: Mon, 29 Sep 2003 07:41
Date: Mon, 29 Sep 2003 07:41
36 lines
2840 bytes
2840 bytes
"Robert Myers" <rmyers@rustuck.com> wrote in message news:iu6fnvch0g09q9f335tt89ogfcj5jm2ovh@4ax.com... > On Sat, 27 Sep 2003 00:17:13 GMT, "Glen Herrmannsfeldt" > <gah@ugcs.caltech.edu> wrote: > > <snip> > > > >PL/I has the MULTIPLY function where one can specify the desired result > >precision. At least one compiler I know of expands both operands before > >doing a double precision multiply when double precision is desired. > > > >If someone wanted to comment about other hardware or compilers in use in the > >1960's and 1970's, that might also be useful. > > > > When I took numerical methods as a graduate student (mid-seventies), I > knew only PL/I and 360 Assembler. Everybody was running on IBM 360 > (or whatever the model was at the time). My classmates, who all used > Fortran IV (while I was coding in PL/I), invariably got different > answers (in the least significant digits) doing the same problem on > the same machine. Since I was not majoring in computer science, > neither the professor nor I had the competence or the interest to > explore the details, although I might take a different attitude today. Since my post on Fortran G and H, I have also run the same test on PL/I (F). That might be the compiler you used. For C=C+A*B; single precision A and B, and double precision C, it does the same thing Fortran did. It would be interesting to see the results from newer Fortran and PL/I compilers. -- glen
Technical Details
Message-ID:
<kmRdb.626671$YN5.452649@sccrnsc01>Path:
arc1-ca.pugleaf.net!news1.google.com!newsfeed.stanford.edu!logbridge.uoregon.edu!arclight.uoregon.edu!wn13feed!wn11feed!worldnet.att.net!204.127.198.203!attbi_feed3!attbi.com!sccrnsc01.POSTED!not-for-mailReferences:
<dG4db.598668$uu5.97984@sccrnsc04> <iu6fnvch0g09q9f335tt89ogfcj5jm2ovh@4ax.com>Show Full Headers
Complete Article Headers
X-Google-Language: ENGLISH,ASCII-7-bit X-Google-Thread: fd4af,7830916812112af0 X-Google-Attributes: gidfd4af,public X-Google-Thread: fd705,7830916812112af0 X-Google-Attributes: gidfd705,public X-Google-Thread: 1094ba,7830916812112af0 X-Google-Attributes: gid1094ba,public X-Google-Thread: 102b75,7830916812112af0 X-Google-Attributes: gid102b75,public X-Google-ArrivalTime: 2003-09-29 00:41:05 PST Path: arc1-ca.pugleaf.net!news1.google.com!newsfeed.stanford.edu!logbridge.uoregon.edu!arclight.uoregon.edu!wn13feed!wn11feed!worldnet.att.net!204.127.198.203!attbi_feed3!attbi.com!sccrnsc01.POSTED!not-for-mail From: "Glen Herrmannsfeldt" <gah@ugcs.caltech.edu> Newsgroups: comp.lang.fortran,comp.arch,bit.listserv.ibm-main,alt.folklore.computers References: <dG4db.598668$uu5.97984@sccrnsc04> <iu6fnvch0g09q9f335tt89ogfcj5jm2ovh@4ax.com> Subject: Re: Precision of floating point multiply Lines: 38 X-Priority: 3 X-MSMail-Priority: Normal X-Newsreader: Microsoft Outlook Express 5.00.2919.6700 X-MimeOLE: Produced By Microsoft MimeOLE V5.00.2919.6700 Message-ID: <kmRdb.626671$YN5.452649@sccrnsc01> NNTP-Posting-Host: 12.207.204.17 X-Complaints-To: abuse@comcast.net X-Trace: sccrnsc01 1064821264 12.207.204.17 (Mon, 29 Sep 2003 07:41:04 GMT) NNTP-Posting-Date: Mon, 29 Sep 2003 07:41:04 GMT Organization: Comcast Online Date: Mon, 29 Sep 2003 07:41:04 GMT