cpow and clog

Steve Kargl sgk at troutmask.apl.washington.edu
Mon Nov 6 20:41:29 UTC 2017


On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote:
> Hello,
> 
> I would like to have some feedback on my attempt to import OpenBSD
> code for cpow and clog:
> 
> https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341
> https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693
> 
> What happened to the alternative implementation mentioned in the thread below?

bde has an implementation of clog[fl].  He may someday 
commit it.  I don't know if anyone ever worked on cpow[fl].
I stopped working on powl and tgammal when I returned my
commit bit due to differences with "higher-ranking" committers.

> And what had stopped the developers from just reusing the Net-i
> or OpenBSD code?

How have you tested the NetBSD and/or OpenBSD code?  What is the
quality?  Have the long double clogl and cpowl been tested on both
ld80 and ld128 hardware?  See FreeBSD's lib/msun/src/math_private.h
for a discussion of possible issues of using I from complex.h in this
code.

> Far as I see, FreeBSD's msun source resembles those other libm's a
> lot anyway?

Well, of course, the various libm's look alike.  They all started
life from Sun Microsystems fdlibm code.  You can get the last
release from Netlib.

http://www.netlib.org/fdlibm/index.html

-- 
Steve


More information about the freebsd-numerics mailing list