Source-Changes-HG archive
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index][Old Index]
[src/trunk]: src/lib/libc/arch/m68k/gen Reimplement the FSCALE version as C-w...
details: https://anonhg.NetBSD.org/src/rev/afae062413c7
branches: trunk
changeset: 475916:afae062413c7
user: mycroft <mycroft%NetBSD.org@localhost>
date: Sun Aug 29 19:41:53 1999 +0000
description:
Reimplement the FSCALE version as C-with-assembler. (Make the compiler deal
with ABI crap.)
diffstat:
lib/libc/arch/m68k/gen/Makefile.inc | 13 +-
lib/libc/arch/m68k/gen/ldexp.S | 62 --------------
lib/libc/arch/m68k/gen/ldexp.c | 152 ------------------------------------
lib/libc/arch/m68k/gen/ldexp_040.c | 68 ++++++++++++++++
lib/libc/arch/m68k/gen/ldexp_881.c | 62 ++++++++++++++
5 files changed, 136 insertions(+), 221 deletions(-)
diffs (truncated from 388 to 300 lines):
diff -r 8c7a0eb17e57 -r afae062413c7 lib/libc/arch/m68k/gen/Makefile.inc
--- a/lib/libc/arch/m68k/gen/Makefile.inc Sun Aug 29 19:41:27 1999 +0000
+++ b/lib/libc/arch/m68k/gen/Makefile.inc Sun Aug 29 19:41:53 1999 +0000
@@ -1,4 +1,4 @@
-# $NetBSD: Makefile.inc,v 1.10 1999/08/29 18:08:30 mycroft Exp $
+# $NetBSD: Makefile.inc,v 1.11 1999/08/29 19:41:53 mycroft Exp $
SRCS+= alloca.S fabs.S frexp.c infinity.c isinf.c isnan.c modf.S
SRCS+= flt_rounds.S fpgetmask.S fpgetround.S fpgetsticky.S fpsetmask.S \
@@ -14,13 +14,12 @@
SRCS+= _setjmp.S
SRCS+= sigsetjmp.S __sigsetjmp14.S
-.ifdef(M68040)
-SRCS+= ldexp.c
-.elifdef(M68060)
-SRCS+= ldexp.c
+.ifdef M68040
+SRCS+= ldexp_040.c
+.elifdef M68060
+SRCS+= ldexp_040.c
.else
-SRCS+= ldexp.S
+SRCS+= ldexp_881.c
.endif
-
LSRCS+= Lint_bswap16.c Lint_bswap32.c Lint_bswap64.c
DPSRCS+= Lint_bswap16.c Lint_bswap32.c Lint_bswap64.c
diff -r 8c7a0eb17e57 -r afae062413c7 lib/libc/arch/m68k/gen/ldexp.S
--- a/lib/libc/arch/m68k/gen/ldexp.S Sun Aug 29 19:41:27 1999 +0000
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,62 +0,0 @@
-/* $NetBSD: ldexp.S,v 1.6 1999/05/02 20:17:21 kleink Exp $ */
-
-/*-
- * Copyright (c) 1990 The Regents of the University of California.
- * All rights reserved.
- *
- * This code is derived from software contributed to Berkeley by
- * the Systems Programming Group of the University of Utah Computer
- * Science Department.
- *
- * 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.
- *
- * 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.
- */
-
-#include <machine/asm.h>
-
-#if defined(LIBC_SCCS) && !defined(lint)
-#if 0
- RCSID("from: @(#)ldexp.s 5.1 (Berkeley) 5/12/90")
-#else
- RCSID("$NetBSD: ldexp.S,v 1.6 1999/05/02 20:17:21 kleink Exp $")
-#endif
-#endif /* LIBC_SCCS and not lint */
-
-/*
- * double ldexp(val, exp)
- * returns: val * (2**exp), for integer exp
- */
-ENTRY(ldexp)
- fmoved sp@(4),fp0
- fscalel sp@(12),fp0
-#ifndef __SVR4_ABI__
- fmoved fp0,sp@-
- movel sp@+,d0
- movel sp@+,d1
-#endif
- rts
diff -r 8c7a0eb17e57 -r afae062413c7 lib/libc/arch/m68k/gen/ldexp.c
--- a/lib/libc/arch/m68k/gen/ldexp.c Sun Aug 29 19:41:27 1999 +0000
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,152 +0,0 @@
-/* $NetBSD: ldexp.c,v 1.6 1999/08/29 19:04:03 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
- * 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.
- *
- * 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.
- *
- * from: Header: ldexp.c,v 1.1 91/07/07 04:28:19 torek Exp
- */
-
-#include <sys/cdefs.h>
-#if defined(LIBC_SCCS) && !defined(lint)
-#if 0
-static const char sccsid[] = "@(#)ldexp.c 8.1 (Berkeley) 6/4/93";
-#else
-__RCSID("$NetBSD: ldexp.c,v 1.6 1999/08/29 19:04:03 mycroft Exp $");
-#endif
-#endif /* LIBC_SCCS and not lint */
-
-#include <sys/types.h>
-#include <machine/ieee.h>
-#include <errno.h>
-#include <math.h>
-
-/*
- * double ldexp(double val, int exp)
- * returns: val * (2**exp)
- */
-double
-ldexp(val, exp)
- double val;
- int exp;
-{
- register int oldexp, newexp, mulexp;
- union {
- double v;
- struct ieee_double s;
- } u, mul;
-
- /*
- * If input is zero, or no change, just return input.
- * Likewise, if input is Inf or NaN, just return it.
- */
- u.v = val;
- oldexp = u.s.dbl_exp;
- if (val == 0.0 || exp == 0 || oldexp == DBL_EXP_INFNAN)
- return (val);
-
- /*
- * Compute new exponent and check for over/under flow.
- * Underflow, unfortunately, could mean switching to denormal.
- * If result out of range, set ERANGE and return 0 if too small
- * or Inf if too big, with the same sign as the input value.
- */
- newexp = oldexp + exp;
- if (newexp >= DBL_EXP_INFNAN) {
- /* u.s.dbl_sign = val < 0.0; -- already set */
- u.s.dbl_exp = DBL_EXP_INFNAN;
- u.s.dbl_frach = u.s.dbl_fracl = 0;
- errno = ERANGE;
- return (u.v); /* Inf */
- }
- if (newexp <= 0) {
- /*
- * The output number is either a denormal or underflows
- * (see comments in machine/ieee.h).
- */
- if (newexp <= -DBL_FRACBITS) {
- errno = ERANGE;
- return (0.0);
- }
- /*
- * We are going to produce a denorm. Our `exp' argument
- * might be as small as -2097, and we cannot compute
- * 2^-2097, so we may have to do this as many as three
- * steps (not just two, as for positive `exp's below).
- */
- mul.v = 1.0;
- while (exp <= -DBL_EXP_BIAS) {
- mul.s.dbl_exp = 1;
- val *= mul.v;
- exp += DBL_EXP_BIAS - 1;
- }
- mul.s.dbl_exp = exp + DBL_EXP_BIAS;
- val *= mul.v;
- return (val);
- }
-
- /*
- * Newexp is positive.
- *
- * If oldexp is zero, we are starting with a denorm, and simply
- * adjusting the exponent will produce bogus answers. We need
- * to fix that first.
- */
- if (oldexp == 0) {
- /*
- * Multiply by 2^mulexp to make the number normalizable.
- * We cannot multiply by more than 2^1023, but `exp'
- * argument might be as large as 2046. A single
- * adjustment, however, will normalize the number even
- * for huge `exp's, and then we can use exponent
- * arithmetic just as for normal `double's.
- */
- mulexp = exp <= DBL_EXP_BIAS ? exp : DBL_EXP_BIAS;
- mul.v = 1.0;
- mul.s.dbl_exp = mulexp + DBL_EXP_BIAS;
- val *= mul.v;
- if (mulexp == exp)
- return (val);
- u.v = val;
- newexp -= mulexp;
- }
-
- /*
- * Both oldexp and newexp are positive; just replace the
- * old exponent with the new one.
- */
- u.s.dbl_exp = newexp;
- return (u.v);
-}
diff -r 8c7a0eb17e57 -r afae062413c7 lib/libc/arch/m68k/gen/ldexp_040.c
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/lib/libc/arch/m68k/gen/ldexp_040.c Sun Aug 29 19:41:53 1999 +0000
@@ -0,0 +1,68 @@
+/* $NetBSD: ldexp_040.c,v 1.1 1999/08/29 19:41:53 mycroft Exp $ */
+
+/*-
+ * Copyright (c) 1998 The NetBSD Foundation, Inc.
+ * All rights reserved.
+ *
+ * This code is derived from software contributed to The NetBSD Foundation
+ * by Charles M. Hannum.
+ *
+ * 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 NetBSD
+ * Foundation, Inc. and its contributors.
+ * 4. Neither the name of The NetBSD Foundation nor the names of its
+ * contributors may be used to endorse or promote products derived
+ * from this software without specific prior written permission.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. 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 FOUNDATION 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.
+ */
+
+#include <sys/cdefs.h>
+#if defined(LIBC_SCCS) && !defined(lint)
+__RCSID("$NetBSD: ldexp_040.c,v 1.1 1999/08/29 19:41:53 mycroft Exp $");
+#endif /* LIBC_SCCS and not lint */
+
+#include <sys/types.h>
+#include <machine/ieee.h>
+#include <math.h>
Home |
Main Index |
Thread Index |
Old Index