Skip site navigation (1)Skip section navigation (2)
Date:      Mon, 20 Aug 2012 11:29:28 -0500
From:      Stephen Montgomery-Smith <stephen@missouri.edu>
To:        freebsd-numerics@freebsd.org
Subject:   Re: Complex arg-trig functions
Message-ID:  <503265E8.3060101@missouri.edu>
In-Reply-To: <502C0CF8.8040003@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> <20120815223631.N1751@besplex.bde.org> <502C0CF8.8040003@missouri.edu>

next in thread | previous in thread | raw e-mail | index | archive | help
On 08/15/12 15:56, Stephen Montgomery-Smith wrote:
> On 08/15/2012 08:35 AM, Bruce Evans wrote:
>> On Tue, 14 Aug 2012, Stephen Montgomery-Smith wrote:
>
>
>>> It seemed to me that there is a logic behind why the the infs and nans
>>> produce the results they do.  I noticed that do_the_hard_work()
>>> already got the answers correct for the real part *rx.  Getting the
>>> imaginary part to work as well seemed to me to be the cleanest way to
>>> make it work.  (I added all the nan and inf checking after writing the
>>> rest of the code.)
>>
>> An up-front check may still be simpler, and gives more control.  In
>> csqrt*(), I needed an explicit check and special expressions to get
>> uniform behaviour.

After thinking about it some more, I decided to go with your approach. 
The additional checks were so very few, and the code becomes quite a bit 
simpler.




Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?503265E8.3060101>