Date: Wed, 15 Aug 2012 17:11:21 -0500 From: Stephen Montgomery-Smith <stephen@missouri.edu> To: freebsd-numerics@freebsd.org Subject: Re: Complex arg-trig functions Message-ID: <502C1E89.9070408@missouri.edu> In-Reply-To: <502C194D.50903@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> <502C0998.7040004@missouri.edu> <502C194D.50903@missouri.edu>
next in thread | previous in thread | raw e-mail | index | archive | help
On 08/15/2012 04:49 PM, Stephen Montgomery-Smith wrote: > On 08/15/2012 03:42 PM, Stephen Montgomery-Smith wrote: > >> >> All your solutions depend upon using (1-tiny) with the result being >> used. But what if FE_DOWNWARD is set? Then 1-tiny becomes >> 1-DBL_EPSILON. And then if the result is used, everything is off by 1 >> ulp. >> >> And >> if ((int)(1 - tiny) == 1) >> will fail. > > How about replacing > > if (huge+ax>one && huge+bx>one) .... > > with > > if ((int)(1/ax)==0 || (int)(1/bx)==0) .... > > (We know that one of ax or bx is larger than 1.) if ((int)(1/(2+ax))==0) .... (because one of ax or bx might be 0).
Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?502C1E89.9070408>