Skip site navigation (1)Skip section navigation (2)
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+tKpWsAqH6OQDkVJjytcFg@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+tKpWsAqH6OQDkVJjytcFg@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: <http://docs.FreeBSD.org/cgi/mid.cgi?20130906194857.GA66345>