Use of C99 extra long double math functions after r236148

Steve Kargl sgk at troutmask.apl.washington.edu
Sun Aug 12 22:56:06 UTC 2012


On Tue, Jul 17, 2012 at 02:21:25PM +1000, Peter Jeremy wrote:
> On 2012-Jul-16 21:01:18 -0700, Steve Kargl <sgk at troutmask.apl.washington.edu> wrote:
> >On Mon, Jul 16, 2012 at 10:40:25PM -0500, Stephen Montgomery-Smith wrote:
> >> 
> >> I came up with pseudo code that looks a bit like this.
> ...
> >Stop.  Please see msun/src/math_private.h. You cannot
> >use I in any expression.
> 
> Note the "pseudo code" reference.  I agree that the actual code has to
> jump through hoops to avoid compiler issues with complex arithmetic
> but doing that for pseudocode just obfuscates it.
> 
> >Also, consult n1256.pdf for x,y = +-0, +-inf, nan. 
> >There are specific requirements that must be met.
> 
> Again, handling the special cases listed in G.6 is all just
> boilerplate code that we can take as assumed for pseudocode.  IMO, it
> would be nice if we could come up with a formal, compact (one/line per
> rule) representation of G.6 that could be pre-processed into wrapper
> code that handles all the 0/Inf/NaN special-cases and then calls
> FOO_finite() which has the hand-written code to handle "normal" cases.
> 

As someone who spent 10+ years getting sqrtl(), cbrtl(), ccosh(), sinl(),
cosl(), tanl(), etc into FreeBSD, I respectfully disagree with your
take that it is just boilerplate.  Getting this stuff right is much harder
than I think some people understand.  Oh well, I'll back to lurking and
working on things I need.
 
-- 
Steve


More information about the freebsd-numerics mailing list