Source-Changes-HG archive
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index][Old Index]
[src/trunk]: src/lib/libc/arch Make the C version of frexp() consitent across...
details: https://anonhg.NetBSD.org/src/rev/a16b38269cbd
branches: trunk
changeset: 475902:a16b38269cbd
user: mycroft <mycroft%NetBSD.org@localhost>
date: Sun Aug 29 18:30:14 1999 +0000
description:
Make the C version of frexp() consitent across platforms, and handle the
Inf/NaN case.
diffstat:
lib/libc/arch/alpha/gen/frexp.c | 86 ++++++++++++++++++++++++++------------
lib/libc/arch/arm32/gen/frexp.c | 7 +-
lib/libc/arch/i386/gen/frexp.c | 39 ++++++++++++----
lib/libc/arch/m68k/gen/frexp.c | 39 ++++++++++++----
lib/libc/arch/mips/gen/frexp.c | 39 ++++++++++++----
lib/libc/arch/powerpc/gen/frexp.c | 37 ++++++++++++----
lib/libc/arch/sparc64/gen/frexp.c | 4 +-
7 files changed, 175 insertions(+), 76 deletions(-)
diffs (truncated from 461 to 300 lines):
diff -r 3544cfde68fb -r a16b38269cbd lib/libc/arch/alpha/gen/frexp.c
--- a/lib/libc/arch/alpha/gen/frexp.c Sun Aug 29 18:23:26 1999 +0000
+++ b/lib/libc/arch/alpha/gen/frexp.c Sun Aug 29 18:30:14 1999 +0000
@@ -1,50 +1,82 @@
-/* $NetBSD: frexp.c,v 1.2 1999/03/10 08:15:43 mycroft Exp $ */
-
+/* $NetBSD: frexp.c,v 1.3 1999/08/29 18:30:14 mycroft Exp $ */
/*
- * Copyright (c) 1994, 1995 Carnegie-Mellon University.
- * All rights reserved.
+ * Copyright (c) 1992, 1993
+ * The Regents of the University of California. All rights reserved.
+ *
+ * This software was developed by the Computer Systems Engineering group
+ * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
+ * contributed to Berkeley.
*
- * Author: Chris G. Demetriou
- *
- * Permission to use, copy, modify and distribute this software and
- * its documentation is hereby granted, provided that both the copyright
- * notice and this permission notice appear in all copies of the
- * software, derivative works or modified versions, and any portions
- * thereof, and that both notices appear in supporting documentation.
- *
- * CARNEGIE MELLON ALLOWS FREE USE OF THIS SOFTWARE IN ITS "AS IS"
- * CONDITION. CARNEGIE MELLON DISCLAIMS ANY LIABILITY OF ANY KIND
- * FOR ANY DAMAGES WHATSOEVER RESULTING FROM THE USE OF THIS SOFTWARE.
- *
- * Carnegie Mellon requests users of this software to return to
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ * notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ * notice, this list of conditions and the following disclaimer in the
+ * documentation and/or other materials provided with the distribution.
+ * 3. All advertising materials mentioning features or use of this software
+ * must display the following acknowledgement:
+ * This product includes software developed by the University of
+ * California, Berkeley and its contributors.
+ * 4. Neither the name of the University nor the names of its contributors
+ * may be used to endorse or promote products derived from this software
+ * without specific prior written permission.
*
- * Software Distribution Coordinator or Software.Distribution%CS.CMU.EDU@localhost
- * School of Computer Science
- * Carnegie Mellon University
- * Pittsburgh PA 15213-3890
+ * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.
*
- * any improvements or extensions that they make and grant Carnegie the
- * rights to redistribute these changes.
+ * from: Header: frexp.c,v 1.1 91/07/07 04:45:01 torek Exp
*/
+#include <sys/cdefs.h>
+#if defined(LIBC_SCCS) && !defined(lint)
+#if 0
+static char sccsid[] = "@(#)frexp.c 8.1 (Berkeley) 6/4/93";
+#else
+__RCSID("$NetBSD: frexp.c,v 1.3 1999/08/29 18:30:14 mycroft Exp $");
+#endif
+#endif /* LIBC_SCCS and not lint */
+
#include <sys/types.h>
#include <machine/ieee.h>
#include <math.h>
+/*
+ * Split the given value into a fraction in the range [0.5, 1.0) and
+ * an exponent, such that frac * (2^exp) == value. If value is 0,
+ * return 0.
+ */
double
frexp(value, eptr)
double value;
int *eptr;
{
- union doub {
+ union {
double v;
struct ieee_double s;
- } u;
+ } u;
if (value) {
+ /*
+ * Fractions in [0.5..1.0) have an exponent of 2^-1.
+ * Leave Inf and NaN alone, however.
+ * WHAT ABOUT DENORMS?
+ */
u.v = value;
- *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
- u.s.dbl_exp = DBL_EXP_BIAS - 1;
+ if (u.s.dbl_exp != DBL_EXP_INFNAN) {
+ *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
+ u.s.dbl_exp = DBL_EXP_BIAS - 1;
+ }
return (u.v);
} else {
*eptr = 0;
diff -r 3544cfde68fb -r a16b38269cbd lib/libc/arch/arm32/gen/frexp.c
--- a/lib/libc/arch/arm32/gen/frexp.c Sun Aug 29 18:23:26 1999 +0000
+++ b/lib/libc/arch/arm32/gen/frexp.c Sun Aug 29 18:30:14 1999 +0000
@@ -1,5 +1,4 @@
-/* $NetBSD: frexp.c,v 1.5 1999/03/10 08:15:43 mycroft Exp $ */
-
+/* $NetBSD: frexp.c,v 1.6 1999/08/29 18:30:16 mycroft Exp $ */
/*
* Copyright (c) 1992, 1993
* The Regents of the University of California. All rights reserved.
@@ -36,7 +35,7 @@
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
* SUCH DAMAGE.
*
- * from: frexp.c,v 1.1 1993/10/07 00:21:33 cgd Exp $
+ * from: Header: frexp.c,v 1.1 91/07/07 04:45:01 torek Exp
*/
#include <sys/cdefs.h>
@@ -44,7 +43,7 @@
#if 0
static char sccsid[] = "@(#)frexp.c 8.1 (Berkeley) 6/4/93";
#else
-__RCSID("$NetBSD: frexp.c,v 1.5 1999/03/10 08:15:43 mycroft Exp $");
+__RCSID("$NetBSD: frexp.c,v 1.6 1999/08/29 18:30:16 mycroft Exp $");
#endif
#endif /* LIBC_SCCS and not lint */
diff -r 3544cfde68fb -r a16b38269cbd lib/libc/arch/i386/gen/frexp.c
--- a/lib/libc/arch/i386/gen/frexp.c Sun Aug 29 18:23:26 1999 +0000
+++ b/lib/libc/arch/i386/gen/frexp.c Sun Aug 29 18:30:14 1999 +0000
@@ -1,8 +1,11 @@
-/* $NetBSD: frexp.c,v 1.3 1999/03/10 08:14:44 mycroft Exp $ */
-
-/*-
- * Copyright (c) 1991 The Regents of the University of California.
- * All rights reserved.
+/* $NetBSD: frexp.c,v 1.4 1999/08/29 18:30:16 mycroft Exp $ */
+/*
+ * Copyright (c) 1992, 1993
+ * The Regents of the University of California. All rights reserved.
+ *
+ * This software was developed by the Computer Systems Engineering group
+ * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
+ * contributed to Berkeley.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
@@ -31,14 +34,16 @@
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
* SUCH DAMAGE.
+ *
+ * from: Header: frexp.c,v 1.1 91/07/07 04:45:01 torek Exp
*/
#include <sys/cdefs.h>
#if defined(LIBC_SCCS) && !defined(lint)
#if 0
-static char sccsid[] = "@(#)frexp.c 5.1 (Berkeley) 3/7/91";
+static char sccsid[] = "@(#)frexp.c 8.1 (Berkeley) 6/4/93";
#else
-__RCSID("$NetBSD: frexp.c,v 1.3 1999/03/10 08:14:44 mycroft Exp $");
+__RCSID("$NetBSD: frexp.c,v 1.4 1999/08/29 18:30:16 mycroft Exp $");
#endif
#endif /* LIBC_SCCS and not lint */
@@ -46,6 +51,11 @@
#include <machine/ieee.h>
#include <math.h>
+/*
+ * Split the given value into a fraction in the range [0.5, 1.0) and
+ * an exponent, such that frac * (2^exp) == value. If value is 0,
+ * return 0.
+ */
double
frexp(value, eptr)
double value;
@@ -53,13 +63,20 @@
{
union {
double v;
- struct ieee_double s;
- } u;
+ struct ieee_double s;
+ } u;
if (value) {
+ /*
+ * Fractions in [0.5..1.0) have an exponent of 2^-1.
+ * Leave Inf and NaN alone, however.
+ * WHAT ABOUT DENORMS?
+ */
u.v = value;
- *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
- u.s.dbl_exp = DBL_EXP_BIAS - 1;
+ if (u.s.dbl_exp != DBL_EXP_INFNAN) {
+ *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
+ u.s.dbl_exp = DBL_EXP_BIAS - 1;
+ }
return (u.v);
} else {
*eptr = 0;
diff -r 3544cfde68fb -r a16b38269cbd lib/libc/arch/m68k/gen/frexp.c
--- a/lib/libc/arch/m68k/gen/frexp.c Sun Aug 29 18:23:26 1999 +0000
+++ b/lib/libc/arch/m68k/gen/frexp.c Sun Aug 29 18:30:14 1999 +0000
@@ -1,8 +1,11 @@
-/* $NetBSD: frexp.c,v 1.4 1999/03/10 08:14:45 mycroft Exp $ */
-
-/*-
- * Copyright (c) 1991 The Regents of the University of California.
- * All rights reserved.
+/* $NetBSD: frexp.c,v 1.5 1999/08/29 18:30:16 mycroft Exp $ */
+/*
+ * Copyright (c) 1992, 1993
+ * The Regents of the University of California. All rights reserved.
+ *
+ * This software was developed by the Computer Systems Engineering group
+ * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
+ * contributed to Berkeley.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
@@ -31,14 +34,16 @@
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
* SUCH DAMAGE.
+ *
+ * from: Header: frexp.c,v 1.1 91/07/07 04:45:01 torek Exp
*/
#include <sys/cdefs.h>
#if defined(LIBC_SCCS) && !defined(lint)
#if 0
-static char sccsid[] = "from: @(#)frexp.c 5.1 (Berkeley) 3/6/91";
+static char sccsid[] = "@(#)frexp.c 8.1 (Berkeley) 6/4/93";
#else
-__RCSID("$NetBSD: frexp.c,v 1.4 1999/03/10 08:14:45 mycroft Exp $");
+__RCSID("$NetBSD: frexp.c,v 1.5 1999/08/29 18:30:16 mycroft Exp $");
#endif
#endif /* LIBC_SCCS and not lint */
@@ -46,6 +51,11 @@
#include <machine/ieee.h>
#include <math.h>
+/*
+ * Split the given value into a fraction in the range [0.5, 1.0) and
+ * an exponent, such that frac * (2^exp) == value. If value is 0,
+ * return 0.
+ */
double
frexp(value, eptr)
double value;
@@ -53,13 +63,20 @@
{
union {
double v;
- struct ieee_double s;
- } u;
+ struct ieee_double s;
+ } u;
if (value) {
+ /*
+ * Fractions in [0.5..1.0) have an exponent of 2^-1.
+ * Leave Inf and NaN alone, however.
+ * WHAT ABOUT DENORMS?
+ */
u.v = value;
- *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
- u.s.dbl_exp = DBL_EXP_BIAS - 1;
+ if (u.s.dbl_exp != DBL_EXP_INFNAN) {
+ *eptr = u.s.dbl_exp - (DBL_EXP_BIAS - 1);
+ u.s.dbl_exp = DBL_EXP_BIAS - 1;
+ }
return (u.v);
} else {
*eptr = 0;
diff -r 3544cfde68fb -r a16b38269cbd lib/libc/arch/mips/gen/frexp.c
--- a/lib/libc/arch/mips/gen/frexp.c Sun Aug 29 18:23:26 1999 +0000
+++ b/lib/libc/arch/mips/gen/frexp.c Sun Aug 29 18:30:14 1999 +0000
@@ -1,8 +1,11 @@
-/* $NetBSD: frexp.c,v 1.4 1999/08/29 18:23:26 mycroft Exp $ */
-
-/*-
Home |
Main Index |
Thread Index |
Old Index