[<prev] [next>] [thread-next>] [day] [month] [year] [list]
Message-ID: <21ca5c9-b1e-71b5-87b-a37f81f691ab@esi.com.au>
Date: Mon, 12 Aug 2024 00:01:01 +1000 (AEST)
From: Damian McGuckin <damianm@....com.au>
To: MUSL <musl@...ts.openwall.com>
Subject: catan(z)
In this routine, there are 2 lines of code
t = 0.5 * atan2(2.0 * x, a);
w = _redupi(t);
The first computes atan2() which returns a number in the range [-pi,+pi]
which means that t is a number in the range [-pi/2,+pi/2].
As far as I understand, the routine _redupi(t) accepts a argument and
reduces it into the range [-pi, +pi]. Am I mistaken?
If I am not wrong, and 't' is already in the range that _redupi(t) will
return, then _redupi() just returns its argument (within a rounding
error). So it seems like a no-op to me.
Can anybody comment?
The routine is originally from Moshier's Cephes library. I have Moshier's
book but I cannot find any hints there.
Thanks - Damian
Powered by blists - more mailing lists
Confused about mailing lists and their use?
Read about mailing lists on Wikipedia
and check out these
guidelines on proper formatting of your messages.
