From owner-freebsd-bugs@FreeBSD.ORG Wed Jan 20 22:20:09 2010 Return-Path: Delivered-To: freebsd-bugs@hub.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:4f8:fff6::34]) by hub.freebsd.org (Postfix) with ESMTP id 3D329106566B for ; Wed, 20 Jan 2010 22:20:09 +0000 (UTC) (envelope-from gnats@FreeBSD.org) Received: from freefall.freebsd.org (freefall.freebsd.org [IPv6:2001:4f8:fff6::28]) by mx1.freebsd.org (Postfix) with ESMTP id 6F6618FC17 for ; Wed, 20 Jan 2010 22:20:06 +0000 (UTC) Received: from freefall.freebsd.org (localhost [127.0.0.1]) by freefall.freebsd.org (8.14.3/8.14.3) with ESMTP id o0KMK65v077067 for ; Wed, 20 Jan 2010 22:20:06 GMT (envelope-from gnats@freefall.freebsd.org) Received: (from gnats@localhost) by freefall.freebsd.org (8.14.3/8.14.3/Submit) id o0KMK6hg077066; Wed, 20 Jan 2010 22:20:06 GMT (envelope-from gnats) Resent-Date: Wed, 20 Jan 2010 22:20:06 GMT Resent-Message-Id: <201001202220.o0KMK6hg077066@freefall.freebsd.org> Resent-From: FreeBSD-gnats-submit@FreeBSD.org (GNATS Filer) Resent-To: freebsd-bugs@FreeBSD.org Resent-Reply-To: FreeBSD-gnats-submit@FreeBSD.org, Paul Guyot Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:4f8:fff6::34]) by hub.freebsd.org (Postfix) with ESMTP id B8BA01065692 for ; Wed, 20 Jan 2010 22:13:25 +0000 (UTC) (envelope-from nobody@FreeBSD.org) Received: from www.freebsd.org (www.freebsd.org [IPv6:2001:4f8:fff6::21]) by mx1.freebsd.org (Postfix) with ESMTP id A850F8FC19 for ; Wed, 20 Jan 2010 22:13:25 +0000 (UTC) Received: from www.freebsd.org (localhost [127.0.0.1]) by www.freebsd.org (8.14.3/8.14.3) with ESMTP id o0KMDPpG030954 for ; Wed, 20 Jan 2010 22:13:25 GMT (envelope-from nobody@www.freebsd.org) Received: (from nobody@localhost) by www.freebsd.org (8.14.3/8.14.3/Submit) id o0KMDPY5030953; Wed, 20 Jan 2010 22:13:25 GMT (envelope-from nobody) Message-Id: <201001202213.o0KMDPY5030953@www.freebsd.org> Date: Wed, 20 Jan 2010 22:13:25 GMT From: Paul Guyot To: freebsd-gnats-submit@FreeBSD.org X-Send-Pr-Version: www-3.1 Cc: Subject: kern/143037: strtod doesn't always round to the nearest X-BeenThere: freebsd-bugs@freebsd.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Bug reports List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 20 Jan 2010 22:20:09 -0000 >Number: 143037 >Category: kern >Synopsis: strtod doesn't always round to the nearest >Confidential: no >Severity: non-critical >Priority: low >Responsible: freebsd-bugs >State: open >Quarter: >Keywords: >Date-Required: >Class: sw-bug >Submitter-Id: current-users >Arrival-Date: Wed Jan 20 22:20:06 UTC 2010 >Closed-Date: >Last-Modified: >Originator: Paul Guyot >Release: 7.1 >Organization: >Environment: FreeBSD papaya.local 7.1-RELEASE-p4 FreeBSD 7.1-RELEASE-p4 #0: Fri Apr 17 05:48:23 UTC 2009 root@papaya.local:/usr/obj/usr/src/sys/GENERIC amd64 >Description: Calling strtod with "2.225073858507201e-308" yields 2^-1022 (0,16,0,0,0,0,0,0), while it should yield 2^-1022 - 2^-1074 (0,15,255,255,255,255,255,255), which is the nearest representable value. gcc on the very same FreeBSD box rounds this to the proper value (see sample code below). >How-To-Repeat: #include #include #include int main() { if (fegetround() == FE_TONEAREST) { printf("rounding to nearest\n"); } else { printf("not rounding to nearest\n"); } const char* fp_str = "2.225073858507201e-308"; printf("string = %s\n", fp_str); double fp_gcc = 2.225073858507201e-308; printf("gcc = %.20g\n", fp_gcc); double parsed = strtod(fp_str, NULL); printf("parsed = %.20g\n", parsed); double fp = 0.0; char* fp_p = (char*) &fp; fp_p[0] = 255; fp_p[1] = 255; fp_p[2] = 255; fp_p[3] = 255; fp_p[4] = 255; fp_p[5] = 255; fp_p[6] = 15; fp_p[7] = 0; printf("smaller = %.20g\n", fp); if (fp == parsed) { printf("actually rounded to nearest (and towards zero)\n"); } fp_p[0] = 0; fp_p[1] = 0; fp_p[2] = 0; fp_p[3] = 0; fp_p[4] = 0; fp_p[5] = 0; fp_p[6] = 16; fp_p[7] = 0; printf("larger = %.20g\n", fp); if (fp == parsed) { printf("actually rounded towards +infinity\n"); } return 0; } >Fix: >Release-Note: >Audit-Trail: >Unformatted: