intervals 0.3.56

data/ext/crlibm/configure.ac

@@ -0,0 +1,364 @@
+dnl Process this file with autoconf to produce a configure script.
+AC_INIT
+
+crlibm_version=`cat $srcdir/VERSION`
+dnl For Automake
+AM_INIT_AUTOMAKE(crlibm, $crlibm_version)
+
+AM_CONFIG_HEADER(crlibm_config.h)
+
+
+AC_SUBST(VERSION)
+
+ISODATE=`date +%Y-%m-%d`
+AC_SUBST(ISODATE)
+
+dnl Check for system information
+AC_CANONICAL_HOST
+case $host_os in
+hpux*)
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPEOS_HPUX,1,[OS type]);;
+*cygwin* | *mingw* | *uwin* | *djgpp | *emx*)
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPEOS_CYGWIN,1,[OS type]);;
+freebsd*|netbsd*)
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPEOS_BSD,1,[OS type]);;
+esac;
+
+
+case $host_cpu in
+alpha)     
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPECPU_ALPHA,1,[Processor type]) ;;
+sparc)     
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPECPU_SPARC,1,[Processor type]) ;;
+powerpc*)  
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPECPU_POWERPC,1,[Processor type]);;
+i*86|k6|k7|pentium*|athlon*)
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPECPU_X86,1,[Processor type])
+	has_ia32_de=true;;
+x86_64|amd64)
+        AC_DEFINE_UNQUOTED(CRLIBM_TYPECPU_AMD64,1,[Processor type])
+	has_ia32_de=true;;
+ia64)
+	AC_DEFINE_UNQUOTED(CRLIBM_TYPECPU_ITANIUM,1,[Processor type])
+	has_ia64_de=true;;
+esac;
+
+
+AC_CONFIG_SRCDIR([exp_accurate.c])
+AC_PROG_MAKE_SET
+
+dnl Check for compiler
+
+AC_PROG_CC(gcc-3.3 gcc-3.2 gcc icc cc)
+
+# Try to set compiler flags to get floating point arithmetics 
+# to work
+#
+if test "${GCC}" != "yes" ; then
+    case "${build}" in
+        *-*-hpux* )
+	    CFLAGS="${CFLAGS} -Wp,-H30000"
+            ;;
+        *-dec-osf* )
+            CFLAGS="${CFLAGS} -ieee"
+            ;;
+	*-apple-darwin* )
+	    CFLAGS="${CFLAGS} -fno-common"
+	    ;;
+	*-sgi-irix* )
+	    ;;
+    esac
+else
+    case "${CC}" in
+	*icc* )
+	    # -Wall causes relocation errors with icc
+	    CFLAGS="${CFLAGS} -wd269"
+	    ;;
+        * )
+            case "${build}" in
+	        *-dec-osf* )
+	            CFLAGS="${CFLAGS} -mieee -Wall"
+	            ;;
+                *alpha*-*-linux* )
+                    CFLAGS="${CFLAGS} -mieee -Wall"
+                    ;;
+                *alpha*-*-freebsd* )
+                    CFLAGS="${CFLAGS} -mieee -Wall"
+                    ;;
+	        *-apple-darwin* )
+	            CFLAGS="${CFLAGS} -fno-common -Wall"
+	            ;;
+	        * )
+	            CFLAGS="${CFLAGS} -Wall"
+	            ;;
+            esac
+	    ;;
+    esac
+fi
+
+
+
+
+
+
+dnl send this information to automake to set up compilation flags
+dnl (add your compiler here)
+AM_CONDITIONAL(COMPILER_ICC, test x$CC = xicc)
+
+dnl Checks for programs.
+AC_PROG_INSTALL
+AC_PROG_RANLIB
+AC_PROG_AWK
+AC_PROG_LN_S
+
+dnl Checks for standard C stuff.
+AC_HEADER_STDC
+AC_FUNC_MALLOC
+AC_C_CONST
+AC_C_INLINE
+
+
+AC_CHECK_FUNCS([sqrt])
+AC_CHECK_HEADERS([fenv.h])
+AC_CHECK_HEADERS([float.h])
+
+dnl adds -lm
+AC_CHECK_LIB([m], [log])
+
+dnl Other checks
+AC_C_BIGENDIAN(,,)
+dnl defines WORDS_BIGENDIAN, or not 
+
+AC_C_LONG_DOUBLE
+dnl defines HAVE_LONG_DOUBLE, or not
+
+
+AC_CHECK_SIZEOF(int)
+dnl defines SIZEOF_INT to be the size in bytes of int
+
+
+
+
+
+dnl All the rest is related to MPFR and to Ziv's lib, useful for testing.
+
+dnl Sometimes the following function is provided by aclocal, sometimes not.
+dnl Copying it here is probably a bad way of managing this problem
+AC_DEFUN(AC_ADD_PATH_LIBS,
+[
+if ` test "$1" `
+then  
+  AC_MSG_CHECKING($2 library)
+	if  test -r "$1/lib$2.a" -o -r "$1/lib$2.so"
+	then
+	  LDFLAGS="$LDFLAGS -L$1"
+          LIBS   ="$LIBS -l$2" 	
+	else
+	   AC_MSG_ERROR($2 not found)
+	fi
+  AC_MSG_RESULT(yes)
+else
+  AC_CHECK_LIB($2, main, , AC_MSG_ERROR($2 not found))
+fi
+]
+)
+
+dnl Command-line arguments to ./configure:
+
+
+dnl A switch to disable use of double-extended optimised version
+AC_ARG_ENABLE(double-extended,        
+  [  --enable-double-extended        enable double-extended optimised versions of functions if hardware supports it  [[default=yes]]],
+  [case $enableval in
+  yes|no) ;;
+  *) AC_MSG_ERROR([bad value $enableval for --enable-double-extended, need yes or no]) ;;
+  esac],
+  [enable_double_extended=yes])
+
+dnl following line needs a little bit more work on log-de.c and double-extended.h
+dnl AM_CONDITIONAL(USE_HARDWARE_DE, test x$has_ia32_de = xtrue -a x$enable_double_extended = xyes   -o  x$has_ia64_de = xtrue -a x$CC = xicc -a x$enable_double_extended = xyes  )
+AM_CONDITIONAL(USE_HARDWARE_DE, test x$has_ia32_de = xtrue -a x$enable_double_extended = xyes   )
+
+dnl If at some point we wish to include SPARC optimization back,
+dnl  the default should take the value 10 and 25 in this case 
+dnl  for around 200 bits of precision
+
+AC_ARG_ENABLE(nbdigits,
+	[  --enable-nbdigits=N      set the number of digits in the SCS structure to N [[default=8]] ],
+	[case $enableval in
+	      *[[0-9]])  AC_DEFINE_UNQUOTED(SCS_NB_WORDS, $enableval, [Number of digits in the SCS structure]);;
+	      *) AC_MSG_ERROR([bad value $enableval for --enable-nbdigits, number expected]) ;;
+	 esac],
+	[AC_DEFINE_UNQUOTED(SCS_NB_WORDS, 8, [Number of digits in the SCS structure])]
+)
+
+AC_ARG_ENABLE(digitsize,
+	[ --enable-digitsize=N      set the size of a digit to N bits [[default=30]] ],
+	[case $enableval in
+	      *[[0-9]])  AC_DEFINE_UNQUOTED(SCS_NB_BITS, $enableval, [Size of an SCS digit]);;
+	      *) AC_MSG_ERROR([bad value $enableval for --enable-digitsize, number expected]) ;;
+	 esac],
+	[AC_DEFINE_UNQUOTED(SCS_NB_BITS, 30, [Size of an SCS digit])]
+)
+
+
+dnl  ... for  MPFR
+AC_ARG_ENABLE(mpfr,       [  --enable-mpfr           enable MPFR function [[default=no]]],
+[case $enableval in
+yes|no) ;;
+*) AC_MSG_ERROR([bad value $enableval for --enable-mpfr, need yes or no]) ;;
+esac],
+[enable_mpfr=no])
+
+
+dnl  ... for GMP
+AC_ARG_ENABLE(gmp,        [  --enable-gmp            enable GMP function [[default=no]]],
+[case $enableval in
+yes|no) ;;
+*) AC_MSG_ERROR([bad value $enableval for --enable-gmp, need yes or no]) ;;
+esac],
+[enable_gmp=no])
+
+
+dnl  ... for IBM Accurate Portable Library
+AC_ARG_ENABLE(ultim,        [  --enable-ultim        enable libultim, IBM's Accurate Portable Math Library [[default=no]]],
+[case $enableval in
+yes|no) ;;
+*) AC_MSG_ERROR([bad value $enableval for --enable-ultim, need yes or no]) ;;
+esac],
+[enable_ultim=no])
+
+
+dnl  ... for Sun's libmcr
+AC_ARG_ENABLE(mcr,        [  --enable-mcr        enable Sun's LIBMCR  [[default=no]]],
+[case $enableval in
+yes|no) ;;
+*) AC_MSG_ERROR([bad value $enableval for --enable-mcr, need yes or no]) ;;
+esac],
+[enable_mcr=no])
+
+
+dnl Now consider the result of the previous to set directories:
+
+
+dnl ... for GMP
+if test "$enable_gmp" = "yes" -o "$enable_mpfr" = "yes"; then
+
+AC_ARG_WITH(gmp_include, [ --with-gmp-include=DIR   GMP include directory ],
+ with_gmp_include=$withval)
+AC_ARG_WITH(gmp_lib,     [ --with-gmp-lib=DIR       GMP lib directory ],
+ with_gmp_lib=$withval)
+AC_ARG_WITH(gmp,         [ --with-gmp=DIR           GMP directory ],
+ with_gmp_include=$withval/include with_gmp_lib=$withval/lib)
+
+if test -d "$with_gmp_include"; then
+  CPPFLAGS="$CPPFLAGS -I$with_gmp_include"
+else
+  with_gmp_include=
+fi
+AC_ADD_PATH_LIBS($with_gmp_lib, gmp)
+
+dnl If gmp.h is found, #defines HAVE_GMP_H
+AC_CHECK_HEADERS([gmp.h])
+
+dnl If libgmp is found, adds -lgmp to LIBS
+dnl correct in theory:  AC_CHECK_LIB([gmp], [__gmp_fprintf])
+AC_CHECK_LIB([gmp], [main])
+
+dnl end of test for gmp
+fi
+
+
+
+dnl  and for MPFR 
+
+if test "$enable_mpfr" = "yes"; then
+AC_ARG_WITH(mpfr_include, [ --with-mpfr-include=DIR  MPFR include directory ], 
+with_mpfr_include=$withval)
+AC_ARG_WITH(mpfr_lib,     [ --with-mpfr-lib=DIR      MPFR lib directory ], 
+with_mpfr_lib=$withval)
+AC_ARG_WITH(mpfr,         [ --with-mpfr=DIR          MPFR directory ],
+with_mpfr_include=$withval/include with_mpfr_lib=$withval/lib)
+
+if test -d "$with_mpfr_include"; then
+  CPPFLAGS="$CPPFLAGS -I$with_mpfr_include"
+else
+  with_mpfr_include=
+fi
+AC_ADD_PATH_LIBS($with_mpfr_lib, mpfr)
+
+
+dnl If mpfr.h is found, #defines HAVE_MPFR_H
+dnl this line should be
+dnl AC_CHECK_HEADERS([mpfr.h], [], [], [#if HAVE_GMP_H #include<gmp.h> #endif])
+dnl but it doesn't work properly so we assume that gmp is installed on the system
+dnl is mpfr is !!!
+AC_CHECK_HEADERS([mpfr.h], [], [], [#include<gmp.h>])
+
+dnl If MPFR is found, adds -lmpfr to LIBS
+dnl correct in theory:  AC_CHECK_LIB([mpfr], [mpfr_init2])
+AC_CHECK_LIB([mpfr], [main])
+
+dnl end of test for mpfr
+fi
+
+
+
+dnl ... for IBM
+if test "$enable_ultim" = "yes"; then
+
+AC_ARG_WITH(ultim_include, [ --with-ultim-include=DIR   ULTIM include directory ],
+ with_ultim_include=$withval)
+AC_ARG_WITH(ultim_lib,     [ --with-ultim-lib=DIR       ULTIM lib directory ],
+ with_ultim_lib=$withval)
+AC_ARG_WITH(ultim,         [ --with-ultim=DIR           ULTIM directory ],
+ with_ultim_include=$withval/include with_ultim_lib=$withval/lib)
+
+if test -d "$with_ultim_include"; then
+  CPPFLAGS="$CPPFLAGS -I$with_ultim_include"
+else
+  with_ultim_include=
+fi
+AC_ADD_PATH_LIBS($with_ultim_lib, ultim)
+
+AC_CHECK_HEADERS([MathLib.h])
+
+AC_CHECK_LIB([ultim], [main])
+
+dnl end of test for IBM libultim
+fi
+
+
+
+dnl ... for SUN
+if test "$enable_mcr" = "yes"; then
+
+AC_ARG_WITH(mcr_include, [ --with-mcr-include=DIR   LIBMCR include directory ],
+ with_mcr_include=$withval)
+AC_ARG_WITH(mcr_lib,     [ --with-mcr-lib=DIR   LIBMCR lib directory ],
+ with_mcr_lib=$withval)
+AC_ARG_WITH(mcr,         [ --with-mcr=DIR       LIBMCR directory ],
+ with_mcr_include=$withval/include with_mcr_lib=$withval/lib)
+
+if test -d "$with_mcr_include"; then
+  CPPFLAGS="$CPPFLAGS -I$with_mcr_include"
+else
+  with_mcr_include=
+fi
+AC_ADD_PATH_LIBS($with_mcr_lib, mcr)
+
+AC_CHECK_HEADERS([libmcr.h])
+
+AC_CHECK_LIB([mcr], [main])
+
+dnl end of test for IBM libultim
+fi
+
+
+#AC_CONFIG_SUBDIRS([scs_lib])
+
+
+
+AC_OUTPUT(Makefile tests/Makefile scs_lib/Makefile)
+
+

data/ext/crlibm/crlibm.h

@@ -0,0 +1,125 @@
+/*
+ * Author  : David Defour, Catherine Daramy, Florent de Dinechin
+ * Contact : David.Defour@ens-lyon.fr, catherine_daramy@ens-lyon.fr
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU Lesser General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or 
+ * (at your option) any later version.
+ * 
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  
+ */
+
+#ifndef CRLIBM_H
+#define CRLIBM_H
+
+#if defined (__cplusplus)
+extern "C" {
+#endif
+
+
+/* An init function which sets FPU flags when needed (mostly on Intel
+   architectures with default double extended) */
+extern unsigned long long crlibm_init(void);
+
+/* An exit function which restores FPU flags when needed (mostly on Intel
+   architectures with default double extended) */
+extern  void crlibm_exit(unsigned long long);
+
+
+/* Finished functions */
+/* These functions are computed in two steps and have an average
+   execution time comparable to that of a standard libm
+*/
+
+/*  exponential  */
+extern double exp_rn(double); /* to nearest  */
+extern double exp_rd(double); /* toward -inf */ 
+extern double exp_ru(double); /* toward +inf */ 
+#define exp_rz exp_rd         /* toward zero */ 
+
+/*  logarithm  */
+extern double log_rn(double); /* to nearest  */
+extern double log_rd(double); /* toward -inf */ 
+extern double log_ru(double); /* toward +inf */ 
+extern double log_rz(double); /* toward zero */ 
+
+/*  cosine  */
+extern double cos_rn(double); /* to nearest  */
+extern double cos_rd(double); /* toward -inf */ 
+extern double cos_ru(double); /* toward +inf */ 
+extern double cos_rz(double); /* toward zero */ 
+
+/*  sine  */
+extern double sin_rn(double); /* to nearest  */
+extern double sin_rd(double); /* toward -inf */ 
+extern double sin_ru(double); /* toward +inf */ 
+extern double sin_rz(double); /* toward zero */ 
+
+/*  tangent  */
+extern double tan_rn(double); /* to nearest  */
+extern double tan_rd(double); /* toward -inf */ 
+extern double tan_ru(double); /* toward +inf */
+extern double tan_rz(double); /* toward zero */
+ 
+
+/*  cotangent  */
+extern double cotan_rn(double); /* to nearest  */
+extern double cotan_rd(double); /* toward -inf */ 
+extern double cotan_ru(double); /* toward +inf */ 
+extern double cotan_rz(double); /* toward zero */ 
+
+/*  arctangent  */
+extern double atan_rn(double); /* to nearest  */
+extern double atan_rd(double); /* toward -inf */ 
+extern double atan_ru(double); /* toward +inf */ 
+extern double atan_rz(double); /* toward zero */ 
+
+/*  hyperbolic cosine*/
+extern double cosh_rn(double); /* to nearest */
+extern double cosh_rd(double); /* toward -inf */ 
+extern double cosh_ru(double); /* toward +inf */ 
+extern double cosh_rz(double); /* toward zero */ 
+
+/*  hyperbolic sine */
+extern double sinh_rn(double); /* to nearest */
+extern double sinh_rd(double); /* toward -inf */ 
+extern double sinh_ru(double); /* toward +inf */ 
+extern double sinh_rz(double); /* toward zero */ 
+
+/*  pow  */
+extern double pow_rn(double, double); /* to nearest */
+
+
+/* Unfinished functions */
+/* These functions provide correct rounding but are very slow
+   (typically 100 times slower that the standard libm) */
+
+
+extern double exp2_rn(double); /* to nearest  */
+extern double exp2_rd(double); /* toward -inf */ 
+extern double exp2_ru(double); /* toward +inf */ 
+
+
+extern double log2_rn(double); /* to nearest  */
+extern double log2_rd(double); /* toward -inf */ 
+extern double log2_ru(double); /* toward +inf */ 
+extern double log2_rz(double); /* towards zero */ 
+
+extern double log10_rn(double); /* to nearest  */
+extern double log10_rd(double); /* toward -inf */ 
+extern double log10_ru(double); /* toward +inf */ 
+extern double log10_rz(double); /* towards zero */ 
+
+#if defined (__cplusplus)
+}
+#endif
+
+#endif /* ifdef CRLIBM_H*/

data/ext/crlibm/crlibm_config.h

@@ -0,0 +1,149 @@
+/* crlibm_config.h.  Generated by configure.  */
+/* crlibm_config.h.in.  Generated from configure.ac by autoheader.  */
+
+/* Processor type */
+/* #undef CRLIBM_TYPECPU_ALPHA */
+
+/* Processor type */
+/* #undef CRLIBM_TYPECPU_AMD64 */
+
+/* Processor type */
+/* #undef CRLIBM_TYPECPU_ITANIUM */
+
+/* Processor type */
+#define CRLIBM_TYPECPU_POWERPC 1
+
+/* Processor type */
+/* #undef CRLIBM_TYPECPU_SPARC */
+
+/* Processor type */
+/* #undef CRLIBM_TYPECPU_X86 */
+
+/* OS type */
+/* #undef CRLIBM_TYPEOS_BSD */
+
+/* OS type */
+/* #undef CRLIBM_TYPEOS_CYGWIN */
+
+/* OS type */
+/* #undef CRLIBM_TYPEOS_HPUX */
+
+/* Define to 1 if you have the <fenv.h> header file. */
+#define HAVE_FENV_H 1
+
+/* Define to 1 if you have the <float.h> header file. */
+#define HAVE_FLOAT_H 1
+
+/* Define to 1 if you have the <gmp.h> header file. */
+/* #undef HAVE_GMP_H */
+
+/* Define to 1 if you have the <inttypes.h> header file. */
+#define HAVE_INTTYPES_H 1
+
+/* Define to 1 if you have the `gmp' library (-lgmp). */
+/* #undef HAVE_LIBGMP */
+
+/* Define to 1 if you have the `m' library (-lm). */
+#define HAVE_LIBM 1
+
+/* Define to 1 if you have the `mcr' library (-lmcr). */
+/* #undef HAVE_LIBMCR */
+
+/* Define to 1 if you have the <libmcr.h> header file. */
+/* #undef HAVE_LIBMCR_H */
+
+/* Define to 1 if you have the `mpfr' library (-lmpfr). */
+/* #undef HAVE_LIBMPFR */
+
+/* Define to 1 if you have the `ultim' library (-lultim). */
+/* #undef HAVE_LIBULTIM */
+
+/* Define to 1 if long double works and has more range or precision than
+   double. */
+/* #undef HAVE_LONG_DOUBLE */
+
+/* Define to 1 if your system has a GNU libc compatible `malloc' function, and
+   to 0 otherwise. */
+#define HAVE_MALLOC 1
+
+/* Define to 1 if you have the <MathLib.h> header file. */
+/* #undef HAVE_MATHLIB_H */
+
+/* Define to 1 if you have the <memory.h> header file. */
+#define HAVE_MEMORY_H 1
+
+/* Define to 1 if you have the <mpfr.h> header file. */
+/* #undef HAVE_MPFR_H */
+
+/* Define to 1 if you have the `sqrt' function. */
+#define HAVE_SQRT 1
+
+/* Define to 1 if you have the <stdint.h> header file. */
+#define HAVE_STDINT_H 1
+
+/* Define to 1 if you have the <stdlib.h> header file. */
+#define HAVE_STDLIB_H 1
+
+/* Define to 1 if you have the <strings.h> header file. */
+#define HAVE_STRINGS_H 1
+
+/* Define to 1 if you have the <string.h> header file. */
+#define HAVE_STRING_H 1
+
+/* Define to 1 if you have the <sys/stat.h> header file. */
+#define HAVE_SYS_STAT_H 1
+
+/* Define to 1 if you have the <sys/types.h> header file. */
+#define HAVE_SYS_TYPES_H 1
+
+/* Define to 1 if you have the <unistd.h> header file. */
+#define HAVE_UNISTD_H 1
+
+/* Name of package */
+#define PACKAGE "crlibm"
+
+/* Define to the address where bug reports for this package should be sent. */
+#define PACKAGE_BUGREPORT ""
+
+/* Define to the full name of this package. */
+#define PACKAGE_NAME ""
+
+/* Define to the full name and version of this package. */
+#define PACKAGE_STRING ""
+
+/* Define to the one symbol short name of this package. */
+#define PACKAGE_TARNAME ""
+
+/* Define to the version of this package. */
+#define PACKAGE_VERSION ""
+
+/* Size of an SCS digit */
+#define SCS_NB_BITS 30
+
+/* Number of digits in the SCS structure */
+#define SCS_NB_WORDS 8
+
+/* The size of a `int', as computed by sizeof. */
+#define SIZEOF_INT 4
+
+/* Define to 1 if you have the ANSI C header files. */
+#define STDC_HEADERS 1
+
+/* Version number of package */
+#define VERSION "0.10beta"
+
+/* Define to 1 if your processor stores words with the most significant byte
+   first (like Motorola and SPARC, unlike Intel and VAX). */
+#define WORDS_BIGENDIAN 1
+
+/* Define to empty if `const' does not conform to ANSI C. */
+/* #undef const */
+
+/* Define to `__inline__' or `__inline' if that's what the C compiler
+   calls it, or to nothing if 'inline' is not supported under any name.  */
+#ifndef __cplusplus
+/* #undef inline */
+#endif
+
+/* Define to rpl_malloc if the replacement function should be used. */
+/* #undef malloc */