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

#! /opt/sage/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

    print "Checking libraries for syntax errors terminated without error."

    print

initialize_env()

objectSo = pobyso_constant_0_sa_so()

pobyso_autoprint(objectSo)

endEllipticList = pobyso_end_elliptic_list_so_sa_so(objectSo, 5)

pobyso_autoprint(endEllipticList)

        sollya_lib_append = sollya.sollya_lib_append

        sollya_lib_build_end_elliptic_list = \

            sollya.sollya_lib_build_end_elliptic_list

        sollya_lib_end_elliptic_list= \

            sollya.sollya_lib_end_elliptic_list

        sollya_lib_prepend = sollya.sollya_lib_prepend

        sollya_lib_roundcoefficients = \

            sollya.sollya_lib_roundcoefficients

    sollya_lib_append.argstypes              = [c_ulong, c_ulong]

    sollya_lib_end_elliptic_list.argtypes    = [c_ulong, c_int]

    sollya_lib_prepend.argstypes             = [c_ulong, c_ulong]

def pobyso_build_end_elliptic_list_so_so(*args):

    """

    From argumrny Sollya objects, create a Sollya end elliptic list.

    Elements of the list are "eaten" (should not be cleared individualy, 

    are cleared when the list is cleared).

    """

    if len(args) == 0:

        ## Called with an empty list produced "error".

        return sollya_lib_build_end_elliptic_list(None)

    index = 0

    ## One can not append elements to an elliptic list, prepend only is 

    #  permitted.

    for argument in reversed(args):

        if index == 0:

            listSo = sollya_lib_build_end_elliptic_list(argument, None)

        else:

            listSo = sollya_lib_prepend(argument, listSo)

        index += 1

    return listSo

# End pobyso_build_end_elliptic_list_so_so

def pobyso_compute_precision_decay_ratio_function_sa_so():

    """

    Compute the precision decay ratio function for polynomial 

    coefficient progressive trucation.

    """

    functionText = """

    proc(deg, a, b, we, wq) 

    {

      k = we * (exp(x/a)-1) + wq * (b*x)^2 + (1-we-wq) * x;

      return k/k(d);

    };

    """

    return pobyso_parse_string_sa_so(functionText)

# End  pobyso_compute_precision_decay_ratio_function.

    return pobyso_constant_from_int_sa_so(0)

    return pobyso_constant_1_sa_so()

    return pobyso_constant_from_int_sa_so(anInt)

    return sollya_lib_constant_from_int64(long(anInt))

    return constSa.value

def pobyso_end_elliptic_list_so_sa_so(objectsListSo, intCountSa):

    """

    Create a Sollya end elliptic list made of the objectListSo[0] to

     objectsListSo[intCountSa-1] objects.

    """     

    return sollya_lib_end_elliptic_list(objectSo, int(intCountSa))

def pobyso_evaluate_so_so(funcSo, argumentSo):

    """

    Evaluates funcSo for arguemntSo through sollya_lib_evaluate().

    """

    return sollya_lib_evaluate(funcSo, argumentSo)

# End pobyso_evaluate_so_so.

        print "Maximum precision of Sollya polynomial coefficients:", expressionPrecSa

    polyCoeffsListSa = polynomialSa.coefficients()

    #for coeff in polyCoeffsListSa:

    #    print coeff.abs().n()

def pobyso_free_variable():

    """

    Ultra thin wrapper around the sollya_lib_function_build_free_variable function.

    """

    return sollya_lib_build_function_free_variable()

def pobyso_inf_so_so(intervalSo):

    """

    Very thin wrapper around sollya_lib_inf().

    """

    return sollya_lib_inf(intervalSo)

# End pobyso_inf_so_so.

def pobyso_polynomial_coefficients_progressive_truncate_so_so(polySo,

                                                              funcSo,

                                                              icSo,

                                                              intervalSo,

                                                              itpSo,

                                                              ftpSo,

                                                              maxPrecSo,

                                                              maxErrSo):

    print "Input arguments:"

    pobyso_autoprint(polySo)

    pobyso_autoprint(funcSo)

    pobyso_autoprint(icSo)

    pobyso_autoprint(intervalSo)

    pobyso_autoprint(itpSo)

    pobyso_autoprint(ftpSo)

    pobyso_autoprint(maxPrecSo)

    pobyso_autoprint(maxErrSo)

    print "________________"

    precDecayFuncSo = \

        pobyso_parse_string("3/8*(exp(x/2.2)-1) + 2/8*(2*x)^2 + (1-3/8-2/8)*x")

    pobyso_autoprint(precDecayFuncSo)

    return None

def pobyso_round_coefficients_single_so_so(polySo, precSo):

    """

    Create a rounded coefficients polynomial from polynomial argument to

    the number of bits in size argument.

    All coefficients are set to the same precision.

    """

    ## TODO: check arguments.

    endEllipListSo = pobyso_build_end_elliptic_list_so_so(precSo)

    polySo = sollya_lib_roundcoefficients(polySo, endEllipListSo, None)

    #sollya_lib_clear_obj(endEllipListSo)

    return polySo

# End pobyso_round_coefficients_single_so_so

# End pobyso_set_prec_sa_so.

# End pobyso_set_prec_so_so.

def pobyso_inf_so_so(intervalSo):

    """

    Very thin wrapper around sollya_lib_inf().

    """

    return sollya_lib_inf(intervalSo)

# End pobyso_inf_so_so.

    return resultSo