Date: Tue, 14 Aug 2012 11:17:27 -0500 From: Stephen Montgomery-Smith <stephen@missouri.edu> To: freebsd-numerics@freebsd.org Subject: Re: Complex arg-trig functions Message-ID: <502A7A17.1090803@missouri.edu> In-Reply-To: <502A780B.2010106@missouri.edu> References: <5017111E.6060003@missouri.edu> <501C361D.4010807@missouri.edu> <20120804165555.X1231@besplex.bde.org> <501D51D7.1020101@missouri.edu> <20120805030609.R3101@besplex.bde.org> <501D9C36.2040207@missouri.edu> <20120805175106.X3574@besplex.bde.org> <501EC015.3000808@missouri.edu> <20120805191954.GA50379@troutmask.apl.washington.edu> <20120807205725.GA10572@server.rulingia.com> <20120809025220.N4114@besplex.bde.org> <5027F07E.9060409@missouri.edu> <20120814003614.H3692@besplex.bde.org> <50295F5C.6010800@missouri.edu> <20120814072946.S5260@besplex.bde.org> <50297CA5.5010900@missouri.edu> <50297E43.7090309@missouri.edu> <20120814201105.T934@besplex.bde.org> <502A780B.2010106@missouri.edu>
next in thread | previous in thread | raw e-mail | index | archive | help
On 08/14/2012 11:08 AM, Stephen Montgomery-Smith wrote: > On 08/14/2012 05:46 AM, Bruce Evans wrote: >> On Mon, 13 Aug 2012, Stephen Montgomery-Smith wrote: >> >>> >>> if (sqrt_huge+x>one && sqrt_huge+y>one) >> >> x and y can be DBL_MAX, giving overflow. > > Why? When x is DBL_MAX, sqrt_huge is so very much smaller than DBL_MAX > that DBL_MAX+sqrt_huge should be DBL_MAX within floating point > precision. So no overflow. I wrote a short test program. It seems to work.
Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?502A7A17.1090803>