From owner-freebsd-current@FreeBSD.ORG  Fri Sep  6 19:48:59 2013
Return-Path: <owner-freebsd-current@FreeBSD.ORG>
Delivered-To: current@freebsd.org
Received: from mx1.freebsd.org (mx1.freebsd.org [8.8.178.115])
 (using TLSv1 with cipher ADH-AES256-SHA (256/256 bits))
 (No client certificate requested)
 by hub.freebsd.org (Postfix) with ESMTP id 075B4582;
 Fri,  6 Sep 2013 19:48:59 +0000 (UTC)
 (envelope-from sgk@troutmask.apl.washington.edu)
Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu
 [128.95.76.21])
 (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits))
 (No client certificate requested)
 by mx1.freebsd.org (Postfix) with ESMTPS id D72942D61;
 Fri,  6 Sep 2013 19:48:58 +0000 (UTC)
Received: from troutmask.apl.washington.edu (localhost.apl.washington.edu
 [127.0.0.1])
 by troutmask.apl.washington.edu (8.14.7/8.14.7) with ESMTP id r86Jmv2x066355; 
 Fri, 6 Sep 2013 12:48:57 -0700 (PDT)
 (envelope-from sgk@troutmask.apl.washington.edu)
Received: (from sgk@localhost)
 by troutmask.apl.washington.edu (8.14.7/8.14.7/Submit) id r86JmvcM066354;
 Fri, 6 Sep 2013 12:48:57 -0700 (PDT) (envelope-from sgk)
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>
Subject: Re: clang+libc++ using missing powl
Message-ID: <20130906194857.GA66345@troutmask.apl.washington.edu>
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>
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Content-Disposition: inline
In-Reply-To: <CAOgwaMuE96_1OagUJ-Vf7=9Fq8M5+tKpWsAqH6OQDkVJjytcFg@mail.gmail.com>
User-Agent: Mutt/1.5.21 (2010-09-15)
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>
X-BeenThere: freebsd-current@freebsd.org
X-Mailman-Version: 2.1.14
Precedence: list
List-Id: Discussions about the use of FreeBSD-current
 <freebsd-current.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/options/freebsd-current>, 
 <mailto:freebsd-current-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/freebsd-current>
List-Post: <mailto:freebsd-current@freebsd.org>
List-Help: <mailto:freebsd-current-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/freebsd-current>,
 <mailto:freebsd-current-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Fri, 06 Sep 2013 19:48:59 -0000

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