intervals 0.3.56

data/ext/crlibm/scs_lib/missing

@@ -0,0 +1,198 @@
+#! /bin/sh
+# Common stub for a few missing GNU programs while installing.
+# Copyright (C) 1996, 1997, 2001, 2002 Free Software Foundation, Inc.
+# Franc,ois Pinard <pinard@iro.umontreal.ca>, 1996.
+
+# This program is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2, 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 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.
+
+if test $# -eq 0; then
+  echo 1>&2 "Try \`$0 --help' for more information"
+  exit 1
+fi
+
+# In the cases where this matters, `missing' is being run in the
+# srcdir already.
+if test -f configure.in; then
+  configure_ac=configure.ac
+else
+  configure_ac=configure.in
+fi
+
+case "$1" in
+
+  -h|--h|--he|--hel|--help)
+    echo "\
+$0 [OPTION]... PROGRAM [ARGUMENT]...
+
+Handle \`PROGRAM [ARGUMENT]...' for when PROGRAM is missing, or return an
+error status if there is no known handling for PROGRAM.
+
+Options:
+  -h, --help      display this help and exit
+  -v, --version   output version information and exit
+
+Supported PROGRAM values:
+  aclocal      touch file \`aclocal.m4'
+  autoconf     touch file \`configure'
+  autoheader   touch file \`config.h.in'
+  automake     touch all \`Makefile.in' files
+  bison        create \`y.tab.[ch]', if possible, from existing .[ch]
+  flex         create \`lex.yy.c', if possible, from existing .c
+  lex          create \`lex.yy.c', if possible, from existing .c
+  makeinfo     touch the output file
+  yacc         create \`y.tab.[ch]', if possible, from existing .[ch]"
+    ;;
+
+  -v|--v|--ve|--ver|--vers|--versi|--versio|--version)
+    echo "missing - GNU libit 0.0"
+    ;;
+
+  -*)
+    echo 1>&2 "$0: Unknown \`$1' option"
+    echo 1>&2 "Try \`$0 --help' for more information"
+    exit 1
+    ;;
+
+  aclocal*)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified \`acinclude.m4' or \`$configure_ac'.  You might want
+         to install the \`Automake' and \`Perl' packages.  Grab them from
+         any GNU archive site."
+    touch aclocal.m4
+    ;;
+
+  autoconf)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified \`$configure_ac'.  You might want to install the
+         \`Autoconf' and \`GNU m4' packages.  Grab them from any GNU
+         archive site."
+    touch configure
+    ;;
+
+  autoheader)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified \`acconfig.h' or \`$configure_ac'.  You might want
+         to install the \`Autoconf' and \`GNU m4' packages.  Grab them
+         from any GNU archive site."
+    files=`sed -n 's/^[ ]*A[CM]_CONFIG_HEADER(\([^)]*\)).*/\1/p' $configure_ac`
+    test -z "$files" && files="config.h"
+    touch_files=
+    for f in $files; do
+      case "$f" in
+      *:*) touch_files="$touch_files "`echo "$f" |
+				       sed -e 's/^[^:]*://' -e 's/:.*//'`;;
+      *) touch_files="$touch_files $f.in";;
+      esac
+    done
+    touch $touch_files
+    ;;
+
+  automake*)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified \`Makefile.am', \`acinclude.m4' or \`$configure_ac'.
+         You might want to install the \`Automake' and \`Perl' packages.
+         Grab them from any GNU archive site."
+    find . -type f -name Makefile.am -print |
+	   sed 's/\.am$/.in/' |
+	   while read f; do touch "$f"; done
+    ;;
+
+  bison|yacc)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified a \`.y' file.  You may need the \`Bison' package
+         in order for those modifications to take effect.  You can get
+         \`Bison' from any GNU archive site."
+    rm -f y.tab.c y.tab.h
+    if [ $# -ne 1 ]; then
+        eval LASTARG="\${$#}"
+	case "$LASTARG" in
+	*.y)
+	    SRCFILE=`echo "$LASTARG" | sed 's/y$/c/'`
+	    if [ -f "$SRCFILE" ]; then
+	         cp "$SRCFILE" y.tab.c
+	    fi
+	    SRCFILE=`echo "$LASTARG" | sed 's/y$/h/'`
+	    if [ -f "$SRCFILE" ]; then
+	         cp "$SRCFILE" y.tab.h
+	    fi
+	  ;;
+	esac
+    fi
+    if [ ! -f y.tab.h ]; then
+	echo >y.tab.h
+    fi
+    if [ ! -f y.tab.c ]; then
+	echo 'main() { return 0; }' >y.tab.c
+    fi
+    ;;
+
+  lex|flex)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified a \`.l' file.  You may need the \`Flex' package
+         in order for those modifications to take effect.  You can get
+         \`Flex' from any GNU archive site."
+    rm -f lex.yy.c
+    if [ $# -ne 1 ]; then
+        eval LASTARG="\${$#}"
+	case "$LASTARG" in
+	*.l)
+	    SRCFILE=`echo "$LASTARG" | sed 's/l$/c/'`
+	    if [ -f "$SRCFILE" ]; then
+	         cp "$SRCFILE" lex.yy.c
+	    fi
+	  ;;
+	esac
+    fi
+    if [ ! -f lex.yy.c ]; then
+	echo 'main() { return 0; }' >lex.yy.c
+    fi
+    ;;
+
+  makeinfo)
+    echo 1>&2 "\
+WARNING: \`$1' is missing on your system.  You should only need it if
+         you modified a \`.texi' or \`.texinfo' file, or any other file
+         indirectly affecting the aspect of the manual.  The spurious
+         call might also be the consequence of using a buggy \`make' (AIX,
+         DU, IRIX).  You might want to install the \`Texinfo' package or
+         the \`GNU make' package.  Grab either from any GNU archive site."
+    file=`echo "$*" | sed -n 's/.*-o \([^ ]*\).*/\1/p'`
+    if test -z "$file"; then
+      file=`echo "$*" | sed 's/.* \([^ ]*\) *$/\1/'`
+      file=`sed -n '/^@setfilename/ { s/.* \([^ ]*\) *$/\1/; p; q; }' $file`
+    fi
+    touch $file
+    ;;
+
+  *)
+    echo 1>&2 "\
+WARNING: \`$1' is needed, and you do not seem to have it handy on your
+         system.  You might have modified some files without having the
+         proper tools for further handling them.  Check the \`README' file,
+         it often tells you about the needed prerequirements for installing
+         this package.  You may also peek at any GNU archive site, in case
+         some other package would contain this missing \`$1' program."
+    exit 1
+    ;;
+esac
+
+exit 0

data/ext/crlibm/scs_lib/mkinstalldirs

@@ -0,0 +1,40 @@
+#! /bin/sh
+# mkinstalldirs --- make directory hierarchy
+# Author: Noah Friedman <friedman@prep.ai.mit.edu>
+# Created: 1993-05-16
+# Public domain
+
+# $Id: mkinstalldirs,v 1.13 1999/01/05 03:18:55 bje Exp $
+
+errstatus=0
+
+for file
+do
+   set fnord `echo ":$file" | sed -ne 's/^:\//#/;s/^://;s/\// /g;s/^#/\//;p'`
+   shift
+
+   pathcomp=
+   for d
+   do
+     pathcomp="$pathcomp$d"
+     case "$pathcomp" in
+       -* ) pathcomp=./$pathcomp ;;
+     esac
+
+     if test ! -d "$pathcomp"; then
+        echo "mkdir $pathcomp"
+
+        mkdir "$pathcomp" || lasterr=$?
+
+        if test ! -d "$pathcomp"; then
+  	  errstatus=$lasterr
+        fi
+     fi
+
+     pathcomp="$pathcomp/"
+   done
+done
+
+exit $errstatus
+
+# mkinstalldirs ends here

data/ext/crlibm/scs_lib/multiplication_scs.c

@@ -0,0 +1,456 @@
+/** Functions for SCS multiplication operations 
+@file multiplication_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"
+
+#if 0 /* used to help debugging */
+void pr(char* s,double d) {
+  db_number x;
+  x.d=d;
+  printf(s);printf("   ");
+  printf("%8x%8x . 2^%d   (%8f %8x %8x)   \n",
+	 (x.i[HI]&0x000FFFFF)+0x00100000,
+	 x.i[LO],
+	 (x.i[HI]>>20)-1023,
+	 x.d,
+	 x.i[HI],
+	 x.i[LO]);
+}
+#endif
+
+
+
+/*  Compute the carry of r1, remove it from r1, and add it to r0 */
+#define SCS_CARRY_PROPAGATE(r1,r0,tmp) \
+      {tmp = r1>>SCS_NB_BITS; r0 += tmp; r1 -= (tmp<<SCS_NB_BITS);}     
+
+
+
+
+
+
+
+/************************************************************/
+/* We have unrolled the loops for SCS_NB_WORDS==4 or 8 
+
+   We just wish gcc would do it for us ! There are option switches,
+   but they don't lead to any performance improvement. When they do,
+   this part of the source code will be removed.
+
+   In the meantime, feel free to unroll for other values. */
+
+
+/***************************/
+#if (SCS_NB_WORDS==4)
+/***************************/
+void scs_mul(scs_ptr result, scs_ptr x, scs_ptr y){
+  uint64_t     val, tmp;
+  uint64_t     r0,r1,r2,r3,r4;
+  uint64_t     x0,x1,x2,x3;
+  int                    y0,y1,y2,y3;
+    
+  R_EXP = X_EXP * Y_EXP;
+  R_SGN = X_SGN * Y_SGN;
+  R_IND = X_IND + Y_IND;
+
+  /* Partial products computation */   
+  x3=X_HW[3];  y3=Y_HW[3];  x2=X_HW[2];  y2=Y_HW[2];
+  x1=X_HW[1];  y1=Y_HW[1];  x0=X_HW[0];  y0=Y_HW[0];
+
+  r4 = x3*y1 + x2*y2 + x1*y3;
+  r3 = x3*y0 + x2*y1 + x1*y2 + x0*y3;
+  r2 = x2*y0 + x1*y1 + x0*y2;
+  r1 = x1*y0 + x0*y1 ;
+  r0 = x0*y0;
+
+  val= 0;
+  /* Carry Propagate */
+  SCS_CARRY_PROPAGATE(r4,r3,tmp)
+  SCS_CARRY_PROPAGATE(r3,r2,tmp)
+  SCS_CARRY_PROPAGATE(r2,r1,tmp)
+  SCS_CARRY_PROPAGATE(r1,r0,tmp)      
+  SCS_CARRY_PROPAGATE(r0,val,tmp)      
+ 
+  if(val != 0){
+    /* shift all the digits ! */
+    R_HW[0] = val; R_HW[1] = r0; R_HW[2] = r1;  R_HW[3] = r2;
+    R_IND += 1;
+  }
+  else {
+    R_HW[0] = r0; R_HW[1] = r1; R_HW[2] = r2; R_HW[3] = r3;
+  }
+
+}
+
+
+void scs_square(scs_ptr result, scs_ptr x){
+  uint64_t  r0,r1,r2,r3,r4;
+  uint64_t  x0,x1,x2,x3;
+  uint64_t  val, tmp;
+
+
+  R_EXP = X_EXP * X_EXP;
+  R_IND = X_IND + X_IND;
+  R_SGN = 1;
+    
+  /*
+   * Calcul des PP
+   */   
+  x3=X_HW[3];  x2=X_HW[2];  x1=X_HW[1];  x0=X_HW[0];
+
+  r0 =  x0*x0;
+  r1 = (x0*x1)* 2 ;
+  r2 =  x1*x1 + (x0*x2*2);
+  r3 = (x1*x2 +  x0*x3)* 2;
+  r4 =  x2*x2 + (x1*x3)* 2;
+
+  val= 0;
+  /* Propagation des retenues */
+  SCS_CARRY_PROPAGATE(r4,r3,tmp)
+  SCS_CARRY_PROPAGATE(r3,r2,tmp)
+  SCS_CARRY_PROPAGATE(r2,r1,tmp)
+  SCS_CARRY_PROPAGATE(r1,r0,tmp)      
+  SCS_CARRY_PROPAGATE(r0,val,tmp)      
+ 
+  if(val != 0){
+    /* shift all the digits ! */
+    R_HW[0] = val; R_HW[1] = r0; R_HW[2] = r1;  R_HW[3] = r2;
+    R_IND += 1;
+  }
+  else {
+    R_HW[0] = r0; R_HW[1] = r1; R_HW[2] = r2; R_HW[3] = r3;
+  }
+  
+}
+
+
+
+
+
+
+
+/***************************/
+#elif (SCS_NB_WORDS==8)
+/***************************/
+void scs_mul(scs_ptr result, scs_ptr x, scs_ptr y){
+  uint64_t     val, tmp;
+  uint64_t     r0,r1,r2,r3,r4,r5,r6,r7,r8;
+  uint64_t     x0,x1,x2,x3,x4,x5,x6,x7;
+  int                    y0,y1,y2,y3,y4,y5,y6,y7;
+    
+  R_EXP = X_EXP * Y_EXP;
+  R_SGN = X_SGN * Y_SGN;
+  R_IND = X_IND + Y_IND;
+
+  /* Partial products computation */   
+  x7=X_HW[7];  y7=Y_HW[7];  x6=X_HW[6];  y6=Y_HW[6];
+  x5=X_HW[5];  y5=Y_HW[5];  x4=X_HW[4];  y4=Y_HW[4];
+  x3=X_HW[3];  y3=Y_HW[3];  x2=X_HW[2];  y2=Y_HW[2];
+  x1=X_HW[1];  y1=Y_HW[1];  x0=X_HW[0];  y0=Y_HW[0];
+
+  r8 = x7*y1 + x6*y2 + x5*y3 + x4*y4 + x3*y5 + x2*y6 + x1*y7;
+  r7 = x7*y0 + x6*y1 + x5*y2 + x4*y3 + x3*y4 + x2*y5 + x1*y6 + x0*y7;
+  r6 = x6*y0 + x5*y1 + x4*y2 + x3*y3 + x2*y4 + x1*y5 + x0*y6;
+  r5 = x5*y0 + x4*y1 + x3*y2 + x2*y3 + x1*y4 + x0*y5;
+  r4 = x4*y0 + x3*y1 + x2*y2 + x1*y3 + x0*y4 ;
+  r3 = x3*y0 + x2*y1 + x1*y2 + x0*y3;
+  r2 = x2*y0 + x1*y1 + x0*y2;
+  r1 = x1*y0 + x0*y1 ;
+  r0 = x0*y0 ;
+ 
+  val= 0;
+  /* Carry Propagate */
+  SCS_CARRY_PROPAGATE(r8,r7,tmp)
+  SCS_CARRY_PROPAGATE(r7,r6,tmp)
+  SCS_CARRY_PROPAGATE(r6,r5,tmp)
+  SCS_CARRY_PROPAGATE(r5,r4,tmp)
+  SCS_CARRY_PROPAGATE(r4,r3,tmp)
+  SCS_CARRY_PROPAGATE(r3,r2,tmp)
+  SCS_CARRY_PROPAGATE(r2,r1,tmp)
+  SCS_CARRY_PROPAGATE(r1,r0,tmp)      
+  SCS_CARRY_PROPAGATE(r0,val,tmp)      
+ 
+  if(val != 0){
+    /* shift all the digits ! */
+    R_HW[0] = val; R_HW[1] = r0; R_HW[2] = r1;  R_HW[3] = r2;
+    R_HW[4] = r3;  R_HW[5] = r4; R_HW[6] = r5;  R_HW[7] = r6;
+    R_IND += 1;
+  }
+  else {
+    R_HW[0] = r0; R_HW[1] = r1; R_HW[2] = r2; R_HW[3] = r3;
+    R_HW[4] = r4; R_HW[5] = r5; R_HW[6] = r6; R_HW[7] = r7;
+  }
+
+}
+
+
+void scs_square(scs_ptr result, scs_ptr x){
+  uint64_t  r0,r1,r2,r3,r4,r5,r6,r7,r8;
+  uint64_t  x0,x1,x2,x3,x4,x5,x6,x7;
+  uint64_t  val, tmp;
+
+
+  R_EXP = X_EXP * X_EXP;
+  R_IND = X_IND + X_IND;
+  R_SGN = 1;
+    
+  /*
+   * Partial products
+   */   
+  x7=X_HW[7];  x6=X_HW[6];  x5=X_HW[5];  x4=X_HW[4];
+  x3=X_HW[3];  x2=X_HW[2];  x1=X_HW[1];  x0=X_HW[0];
+
+  r0 =  x0*x0;
+  r1 = (x0*x1)* 2 ;
+  r2 =  x1*x1 + (x0*x2*2);
+  r3 = (x1*x2 +  x0*x3)* 2;
+  r4 =  x2*x2 + (x1*x3 + x0*x4)* 2;
+  r5 = (x2*x3 +  x1*x4 + x0*x5)* 2;
+  r6 =  x3*x3 + (x2*x4 + x1*x5 + x0*x6)* 2;
+  r7 = (x3*x4 +  x2*x5 + x1*x6 + x0*x7)* 2;
+  r8 =  x4*x4 + (x3*x5 + x2*x6 + x1*x7)* 2;
+
+  val= 0;
+  /* Carry propagation */
+  SCS_CARRY_PROPAGATE(r8,r7,tmp)
+  SCS_CARRY_PROPAGATE(r7,r6,tmp)
+  SCS_CARRY_PROPAGATE(r6,r5,tmp)
+  SCS_CARRY_PROPAGATE(r5,r4,tmp)
+  SCS_CARRY_PROPAGATE(r4,r3,tmp)
+  SCS_CARRY_PROPAGATE(r3,r2,tmp)
+  SCS_CARRY_PROPAGATE(r2,r1,tmp)
+  SCS_CARRY_PROPAGATE(r1,r0,tmp)      
+  SCS_CARRY_PROPAGATE(r0,val,tmp)      
+ 
+  if(val != 0){
+    /* shift all the digits ! */
+    R_HW[0] = val; R_HW[1] = r0; R_HW[2] = r1;  R_HW[3] = r2;
+    R_HW[4] = r3;  R_HW[5] = r4; R_HW[6] = r5;  R_HW[7] = r6;
+    R_IND += 1;
+  }
+  else {
+    R_HW[0] = r0; R_HW[1] = r1; R_HW[2] = r2; R_HW[3] = r3;
+    R_HW[4] = r4; R_HW[5] = r5; R_HW[6] = r6; R_HW[7] = r7;
+  }
+  
+}
+
+
+
+/***************************/
+#else
+/***************************/
+/* From there on, the normal, unrolled case */
+
+
+void scs_mul(scs_ptr result, scs_ptr x, scs_ptr y){
+  uint64_t RES[SCS_NB_WORDS+1];
+  uint64_t val, tmp;
+  int i, j;    
+
+  R_EXP = X_EXP * Y_EXP;
+  R_SGN = X_SGN * Y_SGN;
+  R_IND = X_IND + Y_IND;
+
+  for(i=0; i<=SCS_NB_WORDS; i++)
+    RES[i]=0;
+
+  /* Compute only the first half of the partial product. See the
+     unrolled code for an example of what we compute */
+
+#ifdef CRLIBM_TYPECPU_X86
+  /* This is the only place where there is assembly code to force 64-bit
+     arithmetic. Someday gcc will catch up here, too.
+  */
+  {
+    db_number t;
+    /* i=0 */
+    for(j=0; j<(SCS_NB_WORDS); j++) {
+      __asm__ volatile("mull %3" 
+		       : "=a" (t.i[LO]), "=d" (t.i[HI])
+		       : "a" (X_HW[0]) , "g" (Y_HW[j]));
+      RES[j] += t.l;
+    }
+    /* i = 1..SCS_NB_WORDS-1 */
+    for(i=1 ; i<SCS_NB_WORDS; i++){
+      for(j=0; j<(SCS_NB_WORDS-i); j++){
+	__asm__ volatile("mull %3" 
+			 : "=a" (t.i[LO]), "=d" (t.i[HI])
+			 : "a" (X_HW[i]) , "g" (Y_HW[j]));
+	RES[i+j] += t.l;
+      }
+      __asm__ volatile("mull %3" 
+		       : "=a" (t.i[LO]), "=d" (t.i[HI])
+		       : "a" (X_HW[i]) , "g" (Y_HW[j])); 
+      /* here j==SCS_NB_WORDS-i */
+      RES[SCS_NB_WORDS] += t.l;
+    }
+ }
+
+#else /* other architectures */
+
+ /* i=0 */
+ tmp = X_HW[0];
+ for(j=0; j<(SCS_NB_WORDS); j++)
+   RES[j] += tmp * Y_HW[j];
+ /* i = 1..SCS_NB_WORDS-1 */
+ for(i=1 ; i<SCS_NB_WORDS; i++){
+      tmp = X_HW[i];
+      for(j=0; j<(SCS_NB_WORDS-i); j++)
+	RES[i+j] += tmp * Y_HW[j];
+      RES[SCS_NB_WORDS] += tmp * Y_HW[j]; /* here j==SCS_NB_WORDS-i */
+  }
+#endif/* SCS_TYPECPU_X86 */
+
+  val = 0;
+
+  /* Carry propagate */
+  for(i=SCS_NB_WORDS; i>0; i--)
+    SCS_CARRY_PROPAGATE(RES[i],RES[i-1],tmp)
+  SCS_CARRY_PROPAGATE(RES[0],val,tmp)
+
+
+  /* Store the result */
+  if(val != 0){
+    /* shift all the digits ! */     
+    R_HW[0] = val;
+    for(i=1; i<SCS_NB_WORDS; i++)
+      R_HW[i] = RES[i-1];
+  
+    R_IND += 1;
+  }else {
+    for(i=0; i<SCS_NB_WORDS; i++)
+      R_HW[i] = RES[i];
+   }
+}
+
+
+
+
+
+
+
+void scs_square(scs_ptr result, scs_ptr x){
+  uint64_t RES[SCS_NB_WORDS+1];
+  uint64_t val, tmp;
+  int i, j;
+  
+
+  R_EXP = X_EXP * X_EXP;
+  R_SGN = 1;
+  R_IND = X_IND + X_IND;
+
+  /* Set to 0 intermediate register     */
+  for(i=0; i<=SCS_NB_WORDS; i++)
+    RES[i] = 0;
+
+  /* Compute all the double partial products: 2 x_i * x_j, i!=j */
+  tmp = (uint64_t)X_HW[0];
+  for(j=1; j<SCS_NB_WORDS; j++)
+    RES[j] += tmp * X_HW[j];
+  for(i=1 ; i<(SCS_NB_WORDS+1)/2; i++){
+    tmp = (uint64_t)X_HW[i];
+    for(j=i+1; j<(SCS_NB_WORDS-i); j++)
+      RES[i+j] += tmp * X_HW[j];
+    RES[SCS_NB_WORDS] += tmp * X_HW[SCS_NB_WORDS-i];
+  }
+
+  /* All these partial products are double */
+  for(i=0; i<=SCS_NB_WORDS; i++)
+    RES[i] *=2;
+
+  /* Add partial product of the form x_i^2 */
+  for(i=0, j=0; i<=SCS_NB_WORDS; i+=2, j++){
+    RES[i]  += (uint64_t)X_HW[j] * X_HW[j];
+  }  
+
+  val = 0;
+  /* Carry propagate */
+  for(i=SCS_NB_WORDS; i>0; i--)
+      SCS_CARRY_PROPAGATE(RES[i],RES[i-1],tmp)
+ 
+  SCS_CARRY_PROPAGATE(RES[0],val,tmp)
+
+
+  /* Store the result */
+  if(val != 0){
+    /* shift all the digits ! */     
+    R_HW[0] = val;
+    for(i=1; i<SCS_NB_WORDS; i++)
+      R_HW[i] = RES[i-1];
+  
+    R_IND += 1;
+  }else {
+    for(i=0; i<SCS_NB_WORDS; i++)
+      R_HW[i] = RES[i];
+   }
+
+}
+
+
+/* 
+ * #endif corresponding to the test #if (SCS_NB_WORDS==8)
+ */
+#endif
+
+
+/*
+ Multiply x by an integer val; result is returned in x.
+ */
+ void scs_mul_ui(scs_ptr x, unsigned int val_int){
+  uint64_t val, tmp, vald, rr;
+  int i;
+
+  if (val_int == 0)
+    X_EXP = 0;
+  
+  vald = val_int;
+
+  val = 0; 
+  rr  = 0;
+  for(i=(SCS_NB_WORDS-1); i>=0; i--){
+    val    += vald * X_HW[i];
+    SCS_CARRY_PROPAGATE(val, rr, tmp)
+    X_HW[i] = val;
+    val     = rr;
+    rr      = 0; 
+  }
+
+  if(val != 0){
+    /* shift all the digits ! */ 
+    for(i=(SCS_NB_WORDS-1); i>0; i--)
+      X_HW[i] = X_HW[i-1];
+
+    X_HW[0] = (unsigned int)val;
+    X_IND  += 1;
+  }
+  
+  return;
+}
+