From owner-freebsd-standards@FreeBSD.ORG  Fri Jun  4 20:01:11 2004
Return-Path: <owner-freebsd-standards@FreeBSD.ORG>
Delivered-To: freebsd-standards@hub.freebsd.org
Received: from mx1.FreeBSD.org (mx1.freebsd.org [216.136.204.125])
	by hub.freebsd.org (Postfix) with ESMTP id 6C5CD16A4CE
	for <freebsd-standards@hub.freebsd.org>;
	Fri,  4 Jun 2004 20:01:11 -0700 (PDT)
Received: from freefall.freebsd.org (freefall.freebsd.org [216.136.204.21])
	by mx1.FreeBSD.org (Postfix) with ESMTP id 6641543D2F
	for <freebsd-standards@hub.freebsd.org>;
	Fri,  4 Jun 2004 20:01:11 -0700 (PDT)
	(envelope-from gnats@FreeBSD.org)
Received: from freefall.freebsd.org (gnats@localhost [127.0.0.1])
	i5530lfi044236	for <freebsd-standards@freefall.freebsd.org>;
	Fri, 4 Jun 2004 20:00:47 -0700 (PDT)
	(envelope-from gnats@freefall.freebsd.org)
Received: (from gnats@localhost)
	by freefall.freebsd.org (8.12.11/8.12.11/Submit) id i5530loM044235;
	Fri, 4 Jun 2004 20:00:47 -0700 (PDT)
	(envelope-from gnats)
Date: Fri, 4 Jun 2004 20:00:47 -0700 (PDT)
Message-Id: <200406050300.i5530loM044235@freefall.freebsd.org>
To: freebsd-standards@FreeBSD.org
From: David Schultz <das@FreeBSD.ORG>
Subject: Re: standards/59797: Implement C99's round[f]() math fucntions
X-BeenThere: freebsd-standards@freebsd.org
X-Mailman-Version: 2.1.1
Precedence: list
Reply-To: David Schultz <das@FreeBSD.ORG>
List-Id: Standards compliance <freebsd-standards.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/listinfo/freebsd-standards>,
	<mailto:freebsd-standards-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/freebsd-standards>
List-Post: <mailto:freebsd-standards@freebsd.org>
List-Help: <mailto:freebsd-standards-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/freebsd-standards>,
	<mailto:freebsd-standards-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Sat, 05 Jun 2004 03:01:11 -0000

The following reply was made to PR standards/59797; it has been noted by GNATS.

From: David Schultz <das@FreeBSD.ORG>
To: "Steven G. Kargl" <kargl@troutmask.apl.washington.edu>
Cc: FreeBSD-gnats-submit@FreeBSD.ORG, freebsd-standards@FreeBSD.ORG
Subject: Re: standards/59797: Implement C99's round[f]() math fucntions
Date: Fri, 4 Jun 2004 19:53:59 -0700

 Sorry, I've put this off way too long.  The good news is that I'm
 now going to do something about it.  The bad news is that I found
 a significant bug in the proposed implementation.  Namely, round()
 and roundf() often get the wrong answer for halfway cases.  In
 IEEE-754 round-to-nearest mode, numbers that are halfway between
 two representable numbers are supposed to be rounded to even.  For
 instance, 0.5 becomes 0, and 1.5 becomes 2.  I don't see an easy
 way to make this work in the present implementation without
 fiddling with the underlying bits.  Perhaps we need to implement
 it more similarly to fdlibm's rint().
 
 BTW, benchmarking shows that using the sample implementation that
 appears in the C99 standard results in a slowdown of two orders of
 magnitude over your round() implementation and four orders of
 magnitude over the x87 frndint instruction.  Just setting the
 rounding mode and calling rint() also results in a significant
 slowdown.  Thus, we definitely want something that's in the spirit
 of what you wrote, but perhaps one that operates on the bits directly.