Complex arg-trig functions

Stephen Montgomery-Smith stephen at missouri.edu
Fri Sep 21 01:41:06 UTC 2012


On 09/18/2012 09:15 AM, Bruce Evans wrote:

> @      if (ax == 1 && ay < DBL_EPSILON) {
> @ -        if ((int)ay==0) { /* raise inexact */
> @ +        if (1) { /* inexact will be raised by log() */
> @              /*
> @               * If ay == 0, divide-by-zero will be (correctly)
>
> I didn't re-indent this.

I have put back the old optimizations in catrig.c.  The only change I 
have made so far is that I did re-indent this.


More information about the freebsd-numerics mailing list