intervals 0.3.56

data/ext/crlibm/scs_lib/division_scs.c

@@ -0,0 +1,110 @@
+/** Functions for SCS inverse and division
+
+@file division_scs.c
+
+@author Defour David David.Defour@ens-lyon.fr
+@author Florent de Dinechin Florent.de.Dinechin@ens-lyon.fr 
+ 
+This file is part of the SCS library.
+
+
+*/
+
+/*
+Copyright (C) 2002  David Defour and Florent de Dinechin
+
+    This library 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.1 of the License, or (at your option) any later version.
+
+    This library 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
+    Lesser General Public License for more details.
+
+    You should have received a copy of the GNU Lesser General Public
+    License along with this library; if not, write to the Free Software
+    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+
+ */
+#include "scs.h"
+#include "scs_private.h"
+
+
+void scs_div_2(scs_ptr num) {
+  /* small function to divide by 2 any SCS number */
+  unsigned int carry, mask, old_value;
+  int i;
+  carry = 0x00000000;
+  mask = (unsigned int)(((0x1) << SCS_NB_BITS)-1);/*we now have a mask for the used bits in a word*/
+  /* if it's a normal number, i.e. not zero nor NaN */
+  if( (num->exception).d == (double) 1) {
+    /* first, a loop to rotate all numbers to the right*/
+    for(i = 0; i < SCS_NB_WORDS; i++) {
+      old_value = num->h_word[i];
+      num->h_word[i] = (old_value & !mask) | carry | ((old_value >> 1) & mask);
+      carry = old_value & 0x00000001;/* it can be interesting to keep the last bit of each number =) */
+      carry = carry << (SCS_NB_BITS-1);
+    }
+    /* In the SCS format, the first number can't be zero, so we must handle this particular case */
+    if (num->h_word[0] == 0){
+      num->index = num->index - 1;
+      for(i = 1; i < SCS_NB_WORDS; i++) {
+	num->h_word[i-1] = num->h_word[i];
+      }
+      num->h_word[SCS_NB_WORDS-1] = 0;
+    }
+  }
+  else {
+    (num->exception).d = (num->exception).d / 2;/* zero, NaN, ... */
+  }
+}
+
+
+
+
+/*
+ * Compute 1/x with a Newton scheme
+ */
+void scs_inv(scs_ptr result, scs_ptr x){
+  scs_t tmp, res, res1, scstwo;
+  double app_x, inv;
+  
+  scs_set(tmp, x);  tmp->index = 0;  scs_get_d(&app_x, tmp);
+
+  scs_set_si(scstwo, 2);
+  /* inv is a 53-bit approximation of 1/x */
+  inv = 1/app_x;
+  
+  scs_set_d(res, inv);
+  res->index -= x->index;
+
+  /* First Newton Iteration */
+  scs_mul(res1, x, res);
+  scs_sub(res1, scstwo, res1);
+  scs_mul(res, res, res1);
+
+  /* Second Newton Iteration */
+  scs_mul(res1, x, res);
+  scs_sub(res1, scstwo, res1);
+  scs_mul(result, res, res1); 
+
+  return;
+}
+
+/*
+ * Compute result = x/y; 
+ */
+void scs_div(scs_ptr result, scs_ptr x, scs_ptr y){ 
+  scs_t res;
+ 
+  if (X_EXP != 1){
+    R_EXP = X_EXP / Y_EXP;
+    return;
+  }
+
+  scs_inv(res, y);
+  scs_mul(result, res, x);
+  return;
+}

data/ext/crlibm/scs_lib/double2scs.c

@@ -0,0 +1,174 @@
+/** Conversion of floating-point double to SCS
+@file double2scs.c
+
+@author Defour David David.Defour@ens-lyon.fr
+@author Florent de Dinechin Florent.de.Dinechin@ens-lyon.fr 
+
+This file is part of the SCS library.
+*/
+
+/*
+Copyright (C) 2002  David Defour and Florent de Dinechin
+
+    This library 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.1 of the License, or (at your option) any later version.
+
+    This library 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
+    Lesser General Public License for more details.
+
+    You should have received a copy of the GNU Lesser General Public
+    License along with this library; if not, write to the Free Software
+    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+
+ */
+#include "scs.h"
+#include "scs_private.h"
+
+
+/** Convert a double precision number in it SCS multiprecision
+  representation
+ */
+
+void scs_set_d(scs_ptr result, double x){
+  db_number nb, mantissa;
+  int exponent, exponent_remainder;
+  int ind, i;
+
+  if(x>=0){R_SGN = 1;    nb.d = x;}
+  else    {R_SGN = -1;   nb.d = -x;}
+
+  exponent = nb.i[HI] & 0x7ff00000 ;
+
+  if (exponent == 0x7ff00000)  {
+    /*
+     * x = +/- Inf, s/qNAN
+     */
+    R_EXP = x;  
+    for(i=0; i<SCS_NB_WORDS; i++)
+      R_HW[i] = 0;
+ 
+    R_IND = 0;
+    R_SGN = 1;
+  } 
+
+  else {    /* Normals,  denormals, +/- 0.  */
+
+    /* This number is not an exception */
+    R_EXP = 1;
+
+#if 1
+
+    if (exponent == 0){
+      /* x is a denormal number : bring it back to the normal range */
+      nb.d = nb.d * SCS_RADIX_TWO_DOUBLE;      /* 2^(2.SCS_NB_BITS) */
+      exponent = nb.i[HI] & 0x7ff00000 ;
+      R_IND = -2;
+    }else {
+      R_IND = 0;
+    }
+
+    exponent = exponent >> 20;  /* get the actual value */
+
+    ind = ((exponent +((100*SCS_NB_BITS)-1023))/SCS_NB_BITS) - 100 ;  
+    /* means : = (exponent -1023 + 100*SCS_NB_BITS)/SCS_NB_BITS -100 
+     The business with 100*SCS_NB_BITS is to stay within the positive
+     range for exponent_remainder between 1 and SCS_NB_BITS */
+
+    exponent_remainder = exponent - 1022 - (SCS_NB_BITS*ind);
+
+    R_IND += ind;
+
+    /* now get the mantissa and add the implicit 1 in fp. format*/
+    mantissa.l = (nb.l & ULL(000fffffffffffff)) | ULL(0010000000000000);
+
+
+    /* and spread it over the structure
+       Everything here is 64-bit arithmetic */
+    R_HW[0] = (unsigned int) (mantissa.l >> (53 - exponent_remainder) );
+
+    /* 11 = 64-53 */
+    mantissa.l =  (mantissa.l << (exponent_remainder+11));
+    R_HW[1] = (mantissa.i[HI] >> (32 - SCS_NB_BITS))& SCS_MASK_RADIX ;
+    mantissa.l =  (mantissa.l << SCS_NB_BITS);
+    R_HW[2] = (mantissa.i[HI] >> (32 - SCS_NB_BITS))& SCS_MASK_RADIX ;
+#if SCS_NB_BITS < 27
+    mantissa.l =  (mantissa.l << SCS_NB_BITS);
+    R_HW[3] = (mantissa.i[HI] >> (32 - SCS_NB_BITS))& SCS_MASK_RADIX ;
+#else
+    R_HW[3] = 0 ;
+#endif
+
+#if (SCS_NB_WORDS==8)
+      R_HW[4] = 0; R_HW[5] = 0; R_HW[6] = 0; R_HW[7] = 0;
+#else
+    for(i=4; i<SCS_NB_WORDS; i++)
+      R_HW[i] = 0;
+#endif
+
+
+
+#else /* Other algorithm as in the research report. Slower */
+    R_IND = 0;
+
+    while(nb.d>SCS_RADIX_ONE_DOUBLE) {
+	R_IND++;
+	nb.d *= SCS_RADIX_MONE_DOUBLE;
+      }
+
+    while(nb.d<1)  {
+	R_IND--;
+	nb.d *= SCS_RADIX_ONE_DOUBLE;
+      }
+
+    i=0;
+    while(nb.d != 0){
+      R_HW[i] = (unsigned int) nb.d;
+      nb.d = (nb.d - (double)R_HW[i]) * SCS_RADIX_ONE_DOUBLE;
+      i++;
+    }
+    for(; i<SCS_NB_WORDS; i++)
+      R_HW[i] = 0;
+
+#endif
+
+  } /* end if test NaN etc */
+
+  return;
+}
+
+
+/**
+  Convert an integer number in it scs multiprecision
+  representation
+ */
+void scs_set_si(scs_ptr result, int x){
+  unsigned int ux;
+  int i;
+  
+  if(x>=0){R_SGN = 1;   ux = (unsigned int)x;}
+  else    {R_SGN = -1;  ux = (unsigned int)-x;}
+  
+
+  if (ux > SCS_RADIX){
+    R_IND   = 1;
+    R_HW[0] = (ux - SCS_RADIX) >> SCS_NB_BITS;
+    R_HW[1] =  ux - (R_HW[0] << SCS_NB_BITS);
+  }else {
+    R_IND   = 0;
+    R_HW[0] = ux;
+    R_HW[1] = 0;
+  }
+
+  for(i=2; i<SCS_NB_WORDS; i++)
+    R_HW[i] = 0;
+
+  if (x != 0)  R_EXP = 1;
+  else         R_EXP = 0;
+  
+  return;
+}
+

data/ext/crlibm/scs_lib/install-sh

@@ -0,0 +1,269 @@
+#!/bin/sh
+#
+# install - install a program, script, or datafile
+#
+# This originates from X11R5 (mit/util/scripts/install.sh), which was
+# later released in X11R6 (xc/config/util/install.sh) with the
+# following copyright and license.
+#
+# Copyright (C) 1994 X Consortium
+#
+# Permission is hereby granted, free of charge, to any person obtaining a copy
+# of this software and associated documentation files (the "Software"), to
+# deal in the Software without restriction, including without limitation the
+# rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
+# sell copies of the Software, and to permit persons to whom the Software is
+# furnished to do so, subject to the following conditions:
+#
+# The above copyright notice and this permission notice shall be included in
+# all copies or substantial portions of the Software.
+#
+# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL THE
+# X CONSORTIUM BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN
+# AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNEC-
+# TION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+#
+# Except as contained in this notice, the name of the X Consortium shall not
+# be used in advertising or otherwise to promote the sale, use or other deal-
+# ings in this Software without prior written authorization from the X Consor-
+# tium.
+#
+#
+# FSF changes to this file are in the public domain.
+#
+# Calling this script install-sh is preferred over install.sh, to prevent
+# `make' implicit rules from creating a file called install from it
+# when there is no Makefile.
+#
+# This script is compatible with the BSD install script, but was written
+# from scratch.  It can only install one file at a time, a restriction
+# shared with many OS's install programs.
+
+
+# set DOITPROG to echo to test this script
+
+# Don't use :- since 4.3BSD and earlier shells don't like it.
+doit="${DOITPROG-}"
+
+
+# put in absolute paths if you don't have them in your path; or use env. vars.
+
+mvprog="${MVPROG-mv}"
+cpprog="${CPPROG-cp}"
+chmodprog="${CHMODPROG-chmod}"
+chownprog="${CHOWNPROG-chown}"
+chgrpprog="${CHGRPPROG-chgrp}"
+stripprog="${STRIPPROG-strip}"
+rmprog="${RMPROG-rm}"
+mkdirprog="${MKDIRPROG-mkdir}"
+
+transformbasename=""
+transform_arg=""
+instcmd="$mvprog"
+chmodcmd="$chmodprog 0755"
+chowncmd=""
+chgrpcmd=""
+stripcmd=""
+rmcmd="$rmprog -f"
+mvcmd="$mvprog"
+src=""
+dst=""
+dir_arg=""
+
+while [ x"$1" != x ]; do
+    case $1 in
+	-c) instcmd="$cpprog"
+	    shift
+	    continue;;
+
+	-d) dir_arg=true
+	    shift
+	    continue;;
+
+	-m) chmodcmd="$chmodprog $2"
+	    shift
+	    shift
+	    continue;;
+
+	-o) chowncmd="$chownprog $2"
+	    shift
+	    shift
+	    continue;;
+
+	-g) chgrpcmd="$chgrpprog $2"
+	    shift
+	    shift
+	    continue;;
+
+	-s) stripcmd="$stripprog"
+	    shift
+	    continue;;
+
+	-t=*) transformarg=`echo $1 | sed 's/-t=//'`
+	    shift
+	    continue;;
+
+	-b=*) transformbasename=`echo $1 | sed 's/-b=//'`
+	    shift
+	    continue;;
+
+	*)  if [ x"$src" = x ]
+	    then
+		src=$1
+	    else
+		# this colon is to work around a 386BSD /bin/sh bug
+		:
+		dst=$1
+	    fi
+	    shift
+	    continue;;
+    esac
+done
+
+if [ x"$src" = x ]
+then
+	echo "install:	no input file specified"
+	exit 1
+else
+	true
+fi
+
+if [ x"$dir_arg" != x ]; then
+	dst=$src
+	src=""
+	
+	if [ -d $dst ]; then
+		instcmd=:
+		chmodcmd=""
+	else
+		instcmd=mkdir
+	fi
+else
+
+# Waiting for this to be detected by the "$instcmd $src $dsttmp" command
+# might cause directories to be created, which would be especially bad 
+# if $src (and thus $dsttmp) contains '*'.
+
+	if [ -f $src -o -d $src ]
+	then
+		true
+	else
+		echo "install:  $src does not exist"
+		exit 1
+	fi
+	
+	if [ x"$dst" = x ]
+	then
+		echo "install:	no destination specified"
+		exit 1
+	else
+		true
+	fi
+
+# If destination is a directory, append the input filename; if your system
+# does not like double slashes in filenames, you may need to add some logic
+
+	if [ -d $dst ]
+	then
+		dst="$dst"/`basename $src`
+	else
+		true
+	fi
+fi
+
+## this sed command emulates the dirname command
+dstdir=`echo $dst | sed -e 's,[^/]*$,,;s,/$,,;s,^$,.,'`
+
+# Make sure that the destination directory exists.
+#  this part is taken from Noah Friedman's mkinstalldirs script
+
+# Skip lots of stat calls in the usual case.
+if [ ! -d "$dstdir" ]; then
+defaultIFS=' 	
+'
+IFS="${IFS-${defaultIFS}}"
+
+oIFS="${IFS}"
+# Some sh's can't handle IFS=/ for some reason.
+IFS='%'
+set - `echo ${dstdir} | sed -e 's@/@%@g' -e 's@^%@/@'`
+IFS="${oIFS}"
+
+pathcomp=''
+
+while [ $# -ne 0 ] ; do
+	pathcomp="${pathcomp}${1}"
+	shift
+
+	if [ ! -d "${pathcomp}" ] ;
+        then
+		$mkdirprog "${pathcomp}"
+	else
+		true
+	fi
+
+	pathcomp="${pathcomp}/"
+done
+fi
+
+if [ x"$dir_arg" != x ]
+then
+	$doit $instcmd $dst &&
+
+	if [ x"$chowncmd" != x ]; then $doit $chowncmd $dst; else true ; fi &&
+	if [ x"$chgrpcmd" != x ]; then $doit $chgrpcmd $dst; else true ; fi &&
+	if [ x"$stripcmd" != x ]; then $doit $stripcmd $dst; else true ; fi &&
+	if [ x"$chmodcmd" != x ]; then $doit $chmodcmd $dst; else true ; fi
+else
+
+# If we're going to rename the final executable, determine the name now.
+
+	if [ x"$transformarg" = x ] 
+	then
+		dstfile=`basename $dst`
+	else
+		dstfile=`basename $dst $transformbasename | 
+			sed $transformarg`$transformbasename
+	fi
+
+# don't allow the sed command to completely eliminate the filename
+
+	if [ x"$dstfile" = x ] 
+	then
+		dstfile=`basename $dst`
+	else
+		true
+	fi
+
+# Make a temp file name in the proper directory.
+
+	dsttmp=$dstdir/#inst.$$#
+
+# Move or copy the file name to the temp name
+
+	$doit $instcmd $src $dsttmp &&
+
+	trap "rm -f ${dsttmp}" 0 &&
+
+# and set any options; do chmod last to preserve setuid bits
+
+# If any of these fail, we abort the whole thing.  If we want to
+# ignore errors from any of these, just make sure not to ignore
+# errors from the above "$doit $instcmd $src $dsttmp" command.
+
+	if [ x"$chowncmd" != x ]; then $doit $chowncmd $dsttmp; else true;fi &&
+	if [ x"$chgrpcmd" != x ]; then $doit $chgrpcmd $dsttmp; else true;fi &&
+	if [ x"$stripcmd" != x ]; then $doit $stripcmd $dsttmp; else true;fi &&
+	if [ x"$chmodcmd" != x ]; then $doit $chmodcmd $dsttmp; else true;fi &&
+
+# Now rename the file to the real destination.
+
+	$doit $rmcmd -f $dstdir/$dstfile &&
+	$doit $mvcmd $dsttmp $dstdir/$dstfile 
+
+fi &&
+
+
+exit 0