Complex arg-trig functions

[Stephen Montgomery-Smith]

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.