[src/trunk]: src fenv support for sparc64 from Stathis Kamperis

[christos <christos%NetBSD.org@localhost>]

details:   https://anonhg.NetBSD.org/src/rev/0ff71046578b
branches:  trunk
changeset: 761517:0ff71046578b
user:      christos <christos%NetBSD.org@localhost>
date:      Mon Jan 31 00:19:33 2011 +0000

description:
fenv support for sparc64 from Stathis Kamperis

diffstat:

 distrib/sets/lists/comp/md.sparc64 |    3 +-
 include/fenv.h                     |    4 +-
 lib/libm/Makefile                  |    6 +-
 lib/libm/arch/sparc64/fenv.c       |  329 +++++++++++++++++++++++++++++++++++++
 sys/arch/sparc64/include/Makefile  |    4 +-
 sys/arch/sparc64/include/fenv.h    |   80 ++++++++
 6 files changed, 419 insertions(+), 7 deletions(-)

diffs (truncated from 497 to 300 lines):

diff -r e784614b4256 -r 0ff71046578b distrib/sets/lists/comp/md.sparc64
--- a/distrib/sets/lists/comp/md.sparc64        Mon Jan 31 00:10:06 2011 +0000
+++ b/distrib/sets/lists/comp/md.sparc64        Mon Jan 31 00:19:33 2011 +0000
@@ -1,4 +1,4 @@
-# $NetBSD: md.sparc64,v 1.79 2011/01/18 21:07:51 pooka Exp $
+# $NetBSD: md.sparc64,v 1.80 2011/01/31 00:19:33 christos Exp $
 ./usr/include/ieeefp.h                         comp-c-include
 ./usr/include/sparc                            comp-c-include
 ./usr/include/sparc/_G_config.h                        comp-obsolete           obsolete
@@ -85,6 +85,7 @@
 ./usr/include/sparc64/endian_machdep.h         comp-c-include
 ./usr/include/sparc64/fbio.h                   comp-obsolete           obsolete
 ./usr/include/sparc64/fbvar.h                  comp-obsolete           obsolete
+./usr/include/sparc64/fenv.h                   comp-c-include
 ./usr/include/sparc64/float.h                  comp-c-include
 ./usr/include/sparc64/frame.h                  comp-c-include
 ./usr/include/sparc64/fsr.h                    comp-c-include
diff -r e784614b4256 -r 0ff71046578b include/fenv.h
--- a/include/fenv.h    Mon Jan 31 00:10:06 2011 +0000
+++ b/include/fenv.h    Mon Jan 31 00:19:33 2011 +0000
@@ -1,4 +1,4 @@
-/*     $NetBSD: fenv.h,v 1.2 2010/08/02 17:49:04 joerg Exp $   */
+/*     $NetBSD: fenv.h,v 1.3 2011/01/31 00:19:33 christos Exp $        */
 /*
  * Copyright (c) 2010 The NetBSD Foundation, Inc.
  * All rights reserved.
@@ -25,7 +25,7 @@
  * POSSIBILITY OF SUCH DAMAGE.
  */
 
-#if !defined(__i386__) && !defined(__amd64__)
+#if !defined(__i386__) && !defined(__amd64__) && !defined(__sparc__)
 #error "fenv.h is currently not supported for this architecture"
 #endif
 
diff -r e784614b4256 -r 0ff71046578b lib/libm/Makefile
--- a/lib/libm/Makefile Mon Jan 31 00:10:06 2011 +0000
+++ b/lib/libm/Makefile Mon Jan 31 00:19:33 2011 +0000
@@ -1,4 +1,4 @@
-#  $NetBSD: Makefile,v 1.101 2011/01/12 23:03:56 joerg Exp $
+#  $NetBSD: Makefile,v 1.102 2011/01/31 00:19:33 christos Exp $
 #
 #  @(#)Makefile 5.1beta 93/09/24
 #
@@ -50,6 +50,8 @@
 .if (${MACHINE_ARCH} == "alpha")
 .PATH: ${.CURDIR}/arch/alpha
 ARCH_SRCS = s_copysign.S s_copysignf.S lrint.S
+.elif (${MACHINE_ARCH} == "sparc64")
+.PATH: ${.CURDIR}/arch/sparc64
 .elif ((${MACHINE_ARCH} == "i386") || (${MACHINE_ARCH} == "x86_64"))
 
 .if (${MACHINE_ARCH} == "x86_64")
@@ -174,7 +176,7 @@
 .if (${MACHINE_ARCH} == "vax")
 SRCS=  ${NOIEEE_SRCS} ${NOIEEE_ARCH}
 .else
-SRCS=  ${COMMON_SRCS}
+SRCS=  ${COMMON_SRCS} fenv.c
 .endif
 
 .ifdef ARCH_ADDS
diff -r e784614b4256 -r 0ff71046578b lib/libm/arch/sparc64/fenv.c
--- /dev/null   Thu Jan 01 00:00:00 1970 +0000
+++ b/lib/libm/arch/sparc64/fenv.c      Mon Jan 31 00:19:33 2011 +0000
@@ -0,0 +1,329 @@
+/*     $NetBSD: fenv.c,v 1.1 2011/01/31 00:19:33 christos Exp $        */
+
+/*-
+ * Copyright (c) 2004-2005 David Schultz <das%FreeBSD.ORG@localhost>
+ * All rights reserved.
+ *
+ * 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.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR 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 AUTHOR 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
+ */
+#include <sys/cdefs.h>
+__RCSID("$NetBSD: fenv.c,v 1.1 2011/01/31 00:19:33 christos Exp $");
+
+#include <assert.h>
+#include <fenv.h>
+
+/* Load floating-point state register (all 64bits) */
+#define        __ldxfsr(__r)   __asm__ __volatile__            \
+       ("ldx %0, %%fsr" : : "m" (__r))
+
+/* Save floating-point state register (all 64bits) */
+#define        __stxfsr(__r)   __asm__ __volatile__            \
+       ("stx %%fsr, %0" : "=m" (*(__r)))
+
+/*
+ * The feclearexcept() function clears the supported floating-point exceptions
+ * represented by `excepts'.
+ */
+int
+feclearexcept(int excepts)
+{
+       fexcept_t r;
+       int ex;
+
+       _DIAGASSERT((excepts & ~FE_ALL_EXCEPT) == 0);
+
+       ex = excepts & FE_ALL_EXCEPT;
+
+       __stxfsr(&r);
+       r &= ~ex;
+       __ldxfsr(r);
+
+       /* Success */
+       return 0;
+}
+
+/*
+ * The fegetexceptflag() function stores an implementation-defined
+ * representation of the states of the floating-point status flags indicated
+ * by the argument excepts in the object pointed to by the argument flagp.
+ */
+int
+fegetexceptflag(fexcept_t *flagp, int excepts)
+{
+       fexcept_t r;
+       int ex;
+
+       _DIAGASSERT(flagp != NULL);
+       _DIAGASSERT((excepts & ~_FE_ALL_EXCEPT) == 0);
+
+       ex = excepts & FE_ALL_EXCEPT;
+
+       __stxfsr(&r);
+       *flagp = r & ex;
+
+       /* Success */
+       return 0;
+}
+
+
+/*
+ * This function sets the floating-point status flags indicated by the argument
+ * `excepts' to the states stored in the object pointed to by `flagp'. It does
+ * NOT raise any floating-point exceptions, but only sets the state of the flags.
+ */
+int
+fesetexceptflag(const fexcept_t *flagp, int excepts)
+{
+       fexcept_t r;
+       int ex;
+
+       _DIAGASSERT(flagp != NULL);
+       _DIAGASSERT((excepts & ~FE_ALL_EXCEPT) == 0);
+
+       ex = excepts & FE_ALL_EXCEPT;
+
+       __stxfsr(&r);
+       r &= ~ex;
+       r |= *flagp & ex;
+       __ldxfsr(r);
+
+       /* Success */
+       return 0;
+}
+
+/*
+ * The feraiseexcept() function raises the supported floating-point exceptions
+ * represented by the argument `excepts'.
+ *
+ * The order in which these floating-point exceptions are raised is unspecified
+ * (by the standard).
+ */
+int
+feraiseexcept(int excepts)
+{
+       volatile double d;
+       int ex;
+
+       _DIAGASSERT((excepts & ~FE_ALL_EXCEPT) == 0);
+
+       ex = excepts & FE_ALL_EXCEPT;
+
+       /*
+        * With a compiler that supports the FENV_ACCESS pragma properly, simple
+        * expressions like '0.0 / 0.0' should be sufficient to generate traps.
+        * Unfortunately, we need to bring a volatile variable into the equation
+        * to prevent incorrect optimizations.
+        */
+       if (ex & FE_INVALID) {
+               d = 0.0;
+               d = 0.0 / d;
+       }
+       if (ex & FE_DIVBYZERO) {
+               d = 0.0;
+               d = 1.0 / d;
+       }
+       if (ex & FE_OVERFLOW) {
+               d = 0x1.ffp1023;
+               d *= 2.0;
+       }
+       if (ex & FE_UNDERFLOW) {
+               d = 0x1p-1022;
+               d /= 0x1p1023;
+       }
+       if (ex & FE_INEXACT) {
+               d = 0x1p-1022;
+               d += 1.0;
+       }
+
+       /* Success */
+       return 0;
+}
+
+/*
+ * The fetestexcept() function determines which of a specified subset of the
+ * floating-point exception flags are currently set. The `excepts' argument
+ * specifies the floating-point status flags to be queried.
+ */
+int
+fetestexcept(int excepts)
+{
+       fexcept_t r;
+
+       _DIAGASSERT((excepts & ~FE_ALL_EXCEPT) == 0);
+
+       __stxfsr(&r);
+
+       return r & (excepts & FE_ALL_EXCEPT);
+}
+
+/*
+ * The fegetround() function gets the current rounding direction.
+ */
+int
+fegetround(void)
+{
+       fenv_t r;
+
+       __stxfsr(&r);
+
+       return (r >> _ROUND_SHIFT) & _ROUND_MASK;
+}
+
+/*
+ * The fesetround() function establishes the rounding direction represented by
+ * its argument `round'. If the argument is not equal to the value of a rounding
+ * direction macro, the rounding direction is not changed.
+ */
+int
+fesetround(int round)
+{
+       fenv_t r;
+
+       _DIAGASSERT((round & ~_ROUND_MASK) == 0);
+       if (round & ~_ROUND_MASK)
+               return -1;
+
+       __stxfsr(&r);
+       r &= ~(_ROUND_MASK << _ROUND_SHIFT);
+       r |= round << _ROUND_SHIFT;
+       __ldxfsr(r);
+
+       /* Success */
+       return 0;
+}
+
+/*
+ * The fegetenv() function attempts to store the current floating-point
+ * environment in the object pointed to by envp.
+ */
+int
+fegetenv(fenv_t *envp)
+{
+       _DIAGASSERT(envp != NULL);
+
+       __stxfsr(envp);
+
+       /* Success */
+       return 0;
+}
+
+
+/*
+ * The feholdexcept() function saves the current floating-point environment
+ * in the object pointed to by envp, clears the floating-point status flags, and
+ * then installs a non-stop (continue on floating-point exceptions) mode, if
+ * available, for all floating-point exceptions.
+ */
+int