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

#! /opt/sage/sage

# @file runSLZ-05.sage

#from scipy.constants.codata import precision

def initialize_env():

    """

    Load all necessary modules.

    """

    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/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()

x = var('x')

func(x) = exp(x)

precision = 113

RRR = RealField(precision)

intervalCenter      = RRR(3/8)

icUlp               = intervalCenter.ulp()

intervalRadiusInUlp = 2^49 + 2^45

intervalRadius      = RRR(2^(-35))          

srs_run_SLZ_v05(inputFunction=func, 

                inputLowerBound         = intervalCenter - intervalRadius, 

                inputUpperBound         = intervalCenter + intervalRadius, 

                alpha                   = 6, 

                degree                  = 18, 

                precision               = precision, 

                emin                    = -16382, 

                emax                    = 16383, 

                targetHardnessToRound   = precision * 6, 

                debug                   = True)

#