Article View: comp.lang.fortran
Article #6553Re: True probabilistic rounding [was Re: Precision ...]
From: Dr Chaos
Date: Fri, 03 Oct 2003 21:20
Date: Fri, 03 Oct 2003 21:20
26 lines
2636 bytes
2636 bytes
On 2 Oct 2003 18:36:58 GMT, Nick Maclaren <nmm1@cus.cam.ac.uk> wrote: > What is probably the case is that anyone who relies on absolutely > consistent floating-point results is living dangerously, and may be > deluding himself. I know of at least half a dozen ways where the > assumption will fail, even on consistent IEEE-like architectures. > And considerably more on older ones. > > That is part of the basis for my assertion that probabilistic rounding > is good because it forces programmers to realise that floating-point > arithmetic does not give precisely predictable results. I would like true pseudoprobabilistic rounding, one with the same statisitics as probabilistic rounding, but determinsitically repeatable if need be. Gaussian rounding would be better still. :) In dynamical systems at least there almost always is (in reality) some kind of tiny probabilistic forces, e.g. thermodynamic fluctuations on a macroscopic object. Such rounding in calculations would make numerical results which are ore like what would "physically happen", in the sense that orbits that arise are likely to be ones on the "natural" physical measure. The issue of measure of results is in fact tricky and difficult in theoretical dynamical systems, but everybody agrees that there is some kind of "natural" physical measure out there.
Technical Details
Message-ID:
<slrnbnrq08.tha.mbkennelSPAMBEGONE@lyapunov.ucsd.edu>Path:
arc1-ca.pugleaf.net!archiver1.google.com!news2.google.com!newsfeed.stanford.edu!newsfeed.berkeley.edu!ucberkeley!ihnp4.ucsd.edu!news.ucsd.edu!not-for-mailReferences:
<dG4db.598668$uu5.97984@sccrnsc04> <ble09h$pn2$1@pegasus.csx.cam.ac.uk> <45022fc8.0310021010.4aaf0e5e@posting.google.com> <blhr8a$916$1@pegasus.csx.cam.ac.uk>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: 102b75,7830916812112af0 X-Google-Attributes: gid102b75,public X-Google-Thread: 1094ba,7830916812112af0 X-Google-Attributes: gid1094ba,public X-Google-ArrivalTime: 2003-10-03 14:20:16 PST Path: arc1-ca.pugleaf.net!archiver1.google.com!news2.google.com!newsfeed.stanford.edu!newsfeed.berkeley.edu!ucberkeley!ihnp4.ucsd.edu!news.ucsd.edu!not-for-mail From: Dr Chaos <mbkennelSPAMBEGONE@NOSPAMyahoo.com> Newsgroups: comp.lang.fortran,comp.arch,bit.listserv.ibm-main Subject: Re: True probabilistic rounding [was Re: Precision ...] Followup-To: comp.lang.fortran Date: Fri, 3 Oct 2003 21:20:07 +0000 (UTC) Organization: Univ of Calif San Diego Lines: 28 Message-ID: <slrnbnrq08.tha.mbkennelSPAMBEGONE@lyapunov.ucsd.edu> References: <dG4db.598668$uu5.97984@sccrnsc04> <ble09h$pn2$1@pegasus.csx.cam.ac.uk> <45022fc8.0310021010.4aaf0e5e@posting.google.com> <blhr8a$916$1@pegasus.csx.cam.ac.uk> NNTP-Posting-Host: lyapunov.ucsd.edu X-Trace: news1.ucsd.edu 1065216007 3085 132.239.222.85 (3 Oct 2003 21:20:07 GMT) X-Complaints-To: usenet@news1.ucsd.edu NNTP-Posting-Date: Fri, 3 Oct 2003 21:20:07 +0000 (UTC) User-Agent: slrn/0.9.7.4 (Linux)