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

#! /opt/sage/sage

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/pobyso.py")

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

    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")

def run_SLZ(function,

            lowerBound,

            upperBound,

            alpha,

            degree,

            precision,

            targetHardnessToRound,

            debug = False):

    if debug:

        print "Function                :", function

        print "Lower bound             :", lowerBound

        print "Upper bounds            :", upperBound

        print "Alpha                   :", alpha

        print "Degree                  :", degree

        print "Precision               :", precision

        print "Target hardness-to-round:", targetHardnessToRound

print "Running SLZ..."

#initialize_env()

run_SLZ(exp(x), 1/4, 1/2, 2, 2, 53, 100, True)