Laboratoire de l'Informatique et du Parallélisme » Pobyso

#! /bin/bash

#@file argsForRunSLZ-53-01proj.sh

DEGREE=3

ALPHA=3

HARDNESS_TO_ROUND="53+50"

BINADE="-31"

#

printf -- "$DEGREE \

$ALPHA \

$HARDNESS_TO_ROUND \

$BINADE\n"

#! /opt/sage/sage

# @file runSLZ-113proj.sage

#

#@par Changes from runSLZ-113.sage

# LLL reduction is not performed on the matrix itself but rather on the

# product of the matrix with a uniform random matrix.

# The reduced matrix obtained is discarded but the transformation matrix 

# obtained is used to multiply the original matrix in order to reduced it.

# If a sufficient level of reduction is obtained, we stop here. If not

# the product matrix obtained above is LLL reduced. But as it has been

# pre-reduced at the above step, reduction is supposed to be much faster.

#

# Both reductions combined should hopefully be faster than a straight single

# reduction.

#

# Run SLZ for p=113

#from scipy.constants.codata import precision

def initialize_env():

    """

    Load all necessary modules.

    """

    compiledSpyxDir = "/home/storres/recherche/arithmetique/pobysoPythonSage/compiledSpyx"

    if compiledSpyxDir not in sys.path:

        sys.path.append(compiledSpyxDir)

    if not 'mpfi' in sage.misc.cython.standard_libs:

        sage.misc.cython.standard_libs.append('mpfi')

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sollya_lib.sage")

#    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageMpfr.spyx")

#    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageGMP.spyx")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/pobyso.py")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageSLZ.sage")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageNumericalOperations.sage")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageRationalOperations.sage")

    # Matrix operations are loaded by polynomial operations.

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sagePolynomialOperations.sage")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageRunSLZ.sage")

print "Running SLZ..."

initialize_env()

from sageMpfr import *

from sageGMP  import *

import sys

from subprocess import call

#

## Main variables and parameters.

x            = var('x')

func(x)      = exp(x)

precision    = 53

emin         = -1022

emax         = 1023 

RRR          = RealField(precision)

degree       = 0

alpha        = 0

htrn         = 0

signAsString = ""

binade       = 0

debugMode    = False          

## Local functions

#

def usage():

    write = sys.stderr.write

    write("\nUsage:\n")

    write("  " + scriptName + " <degree> <alpha> <htrn> <intervalCenter>\n")

    write("               <numberOfNumbers> [debug]\n")

    write("\nArguments:\n")

    write("  degree          the degree of the polynomial (integer)\n")

    write("  alpha           alpha (integer)\n")

    write("  htrn            hardness-to-round - a number of bits (integer)\n")

    write("  binade          the binade we want to  explore\n")

    write("                  as a positive integral expression\n")

    write("  sign            the sign of the interval we want to expolore\n")

    write("                  either \"pos\" or \"neg\"\n")

    write("  debug           debug mode (\"debug\", in any case)\n\n")

    sys.exit(2)

# End usage.

#

argsCount = len(sys.argv)

scriptName = os.path.basename(__file__)

if argsCount < 6:

    usage()

for index in xrange(1,argsCount):

    if index == 1:

        degree = int(sys.argv[index])

    elif index == 2:

        alpha = int(sys.argv[index])

    elif index == 3:

        htrn = int(eval(sys.argv[index]))

    elif index == 4:

        binade = int(eval(sys.argv[index]))

    elif index == 5:

        signAsString = sys.argv[index]

    elif index == 6:

        debugMode = sys.argv[index].upper()

        debugMode = (debugMode == "DEBUG")

# Done with command line arguments collection.

#

## Debug printing

print "degree         :", degree

print "alpha          :", alpha

print "htrn           :", htrn

print "binade         :", binade

print "sign           :", signAsString

print "debug mode     :", debugMode

print

#

## Set the terminal window title.

terminalWindowTitle = ['stt', 

                        str(degree), 

                        str(alpha), 

                        str(htrn), 

                        str(binade),

                        signAsString]

call(terminalWindowTitle)

#

binadeUlp    = 2^(binade - precision + 1)

signAsString = signAsString.lower()

if signAsString == "pos": 

    lowerBound   = RRR(2^binade)

    upperBound   = RRR(2^(binade+1)) - binadeUlp

elif signAsString == "neg":

    lowerBound  = -RRR(2^(binade+1)) + binadeUlp

    upperBound  =  -RRR(2^binade)

    print lowerBound

    print upperBound

else:

    exceptionErrorMess = "\"" + signAsString + "\" is an invalid sign." \

    ". Aborting!\n"

    raise Exception(exceptionErrorMess)

srs_run_SLZ_v05(inputFunction=func, 

                     inputLowerBound         = lowerBound, 

                     inputUpperBound         = upperBound, 

                     alpha                   = alpha, 

                     degree                  = degree, 

                     precision               = precision, 

                     emin                    = emin, 

                     emax                    = emax, 

                     targetHardnessToRound   = htrn, 

                     debug                   = debugMode)

        ### If uncommented, the following statement avoids performing

        #   an actual LLL reduction. This allows for demonstrating

        #   the behavior of our pseudo-reduction alone.

        #return ()

        print "First step of reduction afforded enough vectors"

        print "Less than 2 Coppersmith condition compliant vectors after extra reduction."

    Compute a uniform RV in {-1, 1} (not zero).

    Compute a uniform RV in [-2^(n-1), 2^(n-1)-1] (zero is included).

# This file was *autogenerated* from the file ./runSLZ-53.sage

from sage.all_cmdline import *   # import sage library

_sage_const_3 = Integer(3); _sage_const_2 = Integer(2); _sage_const_1 = Integer(1); _sage_const_0 = Integer(0); _sage_const_6 = Integer(6); _sage_const_5 = Integer(5); _sage_const_4 = Integer(4); _sage_const_1023 = Integer(1023); _sage_const_1022 = Integer(1022); _sage_const_10 = Integer(10); _sage_const_53 = Integer(53)#! /opt/sage/sage

# @file runSLZ-113.sage

#

# Run SLZ for p=113

#from scipy.constants.codata import precision

def initialize_env():

    """

    Load all necessary modules.

    """

    compiledSpyxDir = "/home/storres/recherche/arithmetique/pobysoPythonSage/compiledSpyx"

    if compiledSpyxDir not in sys.path:

        sys.path.append(compiledSpyxDir)

    if not 'mpfi' in sage.misc.cython.standard_libs:

        sage.misc.cython.standard_libs.append('mpfi')

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sollya_lib.sage")

#    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageMpfr.spyx")

#    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageGMP.spyx")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/pobyso.py")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageSLZ.sage")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageNumericalOperations.sage")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageRationalOperations.sage")

    # Matrix operations are loaded by polynomial operations.

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sagePolynomialOperations.sage")

    load("/home/storres/recherche/arithmetique/pobysoPythonSage/src/sageSLZ/sageRunSLZ.sage")

print "Running SLZ..."

initialize_env()

from sageMpfr import *

from sageGMP  import *

import sys

from subprocess import call

#

## Main variables and parameters.

x         = var('x')

__tmp__=var("x"); func = symbolic_expression(exp(x)).function(x)

precision = _sage_const_53 

emin      = -_sage_const_1022 

emax      = _sage_const_1023  

RRR = RealField(precision)

degree              = _sage_const_0 

alpha               = _sage_const_0 

htrn                = _sage_const_0 

intervalCenter      = _sage_const_0 

intervalRadius      = _sage_const_0 

debugMode           = False          

## Local functions

#

def usage():

    write = sys.stderr.write

    write("\nUsage:\n")

    write("  " + scriptName + " <degree> <alpha> <htrn> <intervalCenter>\n")

    write("               <numberOfNumbers> [debug]\n")

    write("\nArguments:\n")

    write("  degree          the degree of the polynomial (integer)\n")

    write("  alpha           alpha (integer)\n")

    write("  htrn            hardness-to-round - a number of bits (integer)\n")

    write("  intervalCenter  the interval center (a floating-point number)\n")

    write("  numberOfNumbers the number of floating-point numbers in the interval\n")

    write("                  as a positive integral expression\n")

    write("  debug           debug mode (\"debug\", in any case)\n\n")

    sys.exit(_sage_const_2 )

# End usage.

#

argsCount = len(sys.argv)

scriptName = os.path.basename(__file__)

if argsCount < _sage_const_5 :

    usage()

for index in xrange(_sage_const_1 ,argsCount):

    if index == _sage_const_1 :

        degree = int(sys.argv[index])

    elif index == _sage_const_2 :

        alpha = int(sys.argv[index])

    elif index == _sage_const_3 :

        htrn = int(eval(sys.argv[index]))

    elif index == _sage_const_4 :

        try:

            intervalCenter = QQ(sage_eval(sys.argv[index]))

        except:

            intervalCenter = RRR(sys.argv[index])

        intervalCenter = RRR(intervalCenter)

    elif index == _sage_const_5 :

        ## Can be read as rational number but must end up as an integer.

        numberOfNumbers = QQ(sage_eval(sys.argv[index]))

        if numberOfNumbers != numberOfNumbers.round():

            raise Exception("Invalid number of numbers: " + sys.argv[index] + ".")

        numberOfNumbers = numberOfNumbers.round()

        ## The number must be strictly positive.

        if numberOfNumbers <= _sage_const_0 :

            raise Exception("Invalid number of numbers: " + sys.argv[index] + ".")

    elif index == _sage_const_6 :

        debugMode = sys.argv[index].upper()

        debugMode = (debugMode == "DEBUG")

# Done with command line arguments collection.

#

## Debug printing

print "degree         :", degree

print "alpha          :", alpha

print "htrn           :", htrn

print "interval center:", intervalCenter.n(prec=_sage_const_10 ).str(truncate=False)

print "num of nums    :", RR(numberOfNumbers).log2().n(prec=_sage_const_10 ).str(truncate=False)

print "debug mode     :", debugMode

print

#

## Set the terminal window title.

terminalWindowTitle = ['stt', str(degree), str(alpha), str(htrn), 

                       intervalCenter.n(prec=_sage_const_10 ).str(truncate=False), 

                       RRR(numberOfNumbers).log2().n(prec=_sage_const_10 ).str(truncate=False)]

call(terminalWindowTitle)

#

intervalCenterBinade = slz_compute_binade(intervalCenter)

intervalRadius       = \

    _sage_const_2 **intervalCenterBinade * _sage_const_2 **(-precision + _sage_const_1 ) * numberOfNumbers / _sage_const_2 

srs_run_SLZ_v05(inputFunction=func, 

                inputLowerBound         = intervalCenter - intervalRadius, 

                inputUpperBound         = intervalCenter + intervalRadius, 

                alpha                   = alpha, 

                degree                  = degree, 

                precision               = precision, 

                emin                    = emin, 

                emax                    = emax, 

                targetHardnessToRound   = htrn, 

                debug                   = debugMode)

#! /bin/bash

#@file argsForRunSLZ-53proj-01.sh

#

# These data correspond to the small example of the 

# SLZ-implementation chapter.

#

DEGREE=2

ALPHA=2

HARDNESS_TO_ROUND="53+10"

INTERVAL_CENTER=3/2

NUM_OF_NUMBERS="2^9"

#

printf -- "$DEGREE \

$ALPHA \

$HARDNESS_TO_ROUND \

$INTERVAL_CENTER \

$NUM_OF_NUMBERS\n"

#@file argsForRunSLZ-53-01.sh

# These data correspond to the small example of the 

# SLZ-implementation chapter.

#

DEGREE=2

ALPHA=2

HARDNESS_TO_ROUND="53+10"

INTERVAL_CENTER=3/2

NUM_OF_NUMBERS="2^9"

#