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.

    """

    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/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 *

#

x = var('x')

func(x) = exp(x)

precision = 53

RRR = RealField(precision)

intervalCenter      = RRR(3/2)

icUlp               = intervalCenter.ulp()

intervalRadiusInUlp = 2^8

intervalRadius      = RRR(2^(-35))          

srs_run_SLZ_v06(inputFunction           = func, 

                inputLowerBound         = intervalCenter - intervalRadiusInUlp * icUlp, 

                inputUpperBound         = intervalCenter + intervalRadiusInUlp * icUlp, 

                alpha                   = 2, 

                degree                  = 2, 

                precision               = precision, 

                emin                    = -1022, 

                emax                    = 1023, 

                targetHardnessToRound   = precision + 10, 

                debug                   = True)

#