Date: Fri, 6 Sep 2013 12:48:57 -0700 From: Steve Kargl <sgk@troutmask.apl.washington.edu> To: Mehmet Erol Sanliturk <m.e.sanliturk@gmail.com> Cc: Baptiste Daroussin <bapt@freebsd.org>, "freebsd-numerics@FreeBSD.org" <freebsd-numerics@freebsd.org>, David Chisnall <theraven@freebsd.org>, "current@freebsd.org Current" <current@freebsd.org> Subject: Re: clang+libc++ using missing powl Message-ID: <20130906194857.GA66345@troutmask.apl.washington.edu> In-Reply-To: <CAOgwaMuE96_1OagUJ-Vf7=9Fq8M5%2BtKpWsAqH6OQDkVJjytcFg@mail.gmail.com> References: <20130904223842.GB82066@ithaqua.etoilebsd.net> <80062343-53CD-4CEF-9C47-3BF614DADB64@FreeBSD.org> <20130905210953.GA58413@troutmask.apl.washington.edu> <4DF3383B-4BE7-4947-886A-AAAD25172F68@freebsd.org> <20130906155939.GA63661@troutmask.apl.washington.edu> <CAOgwaMuE96_1OagUJ-Vf7=9Fq8M5%2BtKpWsAqH6OQDkVJjytcFg@mail.gmail.com>
next in thread | previous in thread | raw e-mail | index | archive | help
On Fri, Sep 06, 2013 at 03:02:18PM -0400, Mehmet Erol Sanliturk wrote: > On Fri, Sep 6, 2013 at 11:59 AM, Steve Kargl < > sgk@troutmask.apl.washington.edu> wrote: > > > > PS: I have working erfl and erfcl for ld80 archs. I'm still > > testing and refining the code. It turns out that computing > > the needed rational approximation is fairly difficult (at least > > for me). > > > Not only for you , because function values are not very "smooth" , and > representation by an approximation is very difficult . During many weeks , > I am trying to obtain good approximations to cumulative normal distribution > and its inverse , without very much success . > > Approximation to > inverse cumulative normal distribution is much more difficult than > cumulative normal distribution . > > I am using Fortran double precision . I think , it is necessary to use an > arbitrary precision package for optimization , but I do not have any one . > If you are using gfortran version 4.7 or newer on FreeBSD, you should be able to use a quad precision type. In addition, if you are using lang/gcc46 or newer on your FreeBSD system, then you have an arbitrary precision library installed (see math/mpfr). -- Steve
Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?20130906194857.GA66345>