From owner-svn-src-head@FreeBSD.ORG Mon Nov 15 00:51:23 2010 Return-Path: Delivered-To: svn-src-head@freebsd.org Received: by hub.freebsd.org (Postfix, from userid 1233) id 93E9A1065670; Mon, 15 Nov 2010 00:51:23 +0000 (UTC) Date: Mon, 15 Nov 2010 00:51:23 +0000 From: Alexander Best To: Steve Kargl Message-ID: <20101115005123.GA46134@freebsd.org> References: <201011131054.oADAsA7I045096@svn.freebsd.org> <20101113125648.GA25183@freebsd.org> <20101113164645.GB69612@troutmask.apl.washington.edu> Mime-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20101113164645.GB69612@troutmask.apl.washington.edu> Cc: svn-src-head@freebsd.org, svn-src-all@freebsd.org, src-committers@freebsd.org, Ulrich Spoerlein Subject: Re: svn commit: r215237 - head/lib/msun/src X-BeenThere: svn-src-head@freebsd.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: SVN commit messages for the src tree for head/-current List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 15 Nov 2010 00:51:23 -0000 On Sat Nov 13 10, Steve Kargl wrote: > On Sat, Nov 13, 2010 at 12:56:48PM +0000, Alexander Best wrote: > > On Sat Nov 13 10, Ulrich Spoerlein wrote: > > > Author: uqs > > > Date: Sat Nov 13 10:54:10 2010 > > > New Revision: 215237 > > > URL: http://svn.freebsd.org/changeset/base/215237 > > > > > > Log: > > > Fix bug in jn(3) and jnf(3) that led to -inf results > > > > thank you very much for fixing this long outstanding issue. > > you might want to have a look at [1], where two more issues have been reported. > > > > cheers. > > alex > > > > [1] http://mailman.oakapple.net/pipermail/numeric-interest/2010-September/thread.html > > > > Ulrich, thanks. Of the two remaining issues, I can verify > one is indeed an issue. also it seems one of the regression tests in tools/regression/lib/msun fails (assert tripped): test-lrint. cheers. alex > > laptop:kargl[204] cat mn.c > #include > #include > > static const double z = 0.; > > int > main(void) > { > double x, y; > x = -1. / z; > y = exp( 0.5 * log(x)); > printf("pow(%f, 0.5) = %le, %le\n", x, pow(x, 0.5), y); > y = exp(-0.5 * log(x)); > printf("pow(%f, -0.5) = %le, %le\n", x, pow(x, -0.5), y); > return 0; > } > > laptop:kargl[205] cc -o z mn.c -fno-builtin -lm && ./z > pow(-inf, 0.5) = inf, nan > pow(-inf, -0.5) = 0.000000e+00, nan > > -- > Steve -- a13x