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

import timeit

def test_pobyso_wrapper(func, *args, **kwargs):

    def wrapped():

        return func(*args, **kwargs)

    return wrapped

def test_pobyso_bounds_to_range_sa_so(repeat=1000, number=10):

    functionName =  inspect.stack()[0][3]

    print "Running", inspect.stack()[0][3], "..."

    #pobyso_autoprint(precSo)

    #

    def test(lowerBoundSa, upperBoundSa):

    #

    wrapped = test_pobyso_wrapper(test,lowerBoundSa, upperBoundSa)

    timing = min(timeit.repeat(wrapped, repeat=repeat, number=number))

    #

    #pobyso_autoprint(precSo)

    print "\t...", functionName, "done."

    return timing

def test_pobyso_get_list_elements_so_so(repeat = 1000, number=10):

    functionName = inspect.stack()[0][3]

    print "Running", functionName, "..."

    def test():

        listSo = pobyso_parse_string_sa_so("[|0,1,2,3,4,5,6|];")

        (listSaSo, elementsNumSa, isEndEllipticSa) = \

            pobyso_get_list_elements(listSo)

        for elementSo in listSaSo:

            sollya_lib_clear_obj(elementSo)

        # Notice that we use sollya_lib_clear_obj and not

        # sollya_lib_clear_object_list.

        sollya_lib_clear_obj(listSo, None)

    #

    wrapped = test_pobyso_wrapper(test)

    timing = min(timeit.repeat(wrapped, repeat=repeat, number=number))

    print "\t...", functionName, "done."

    return timing

# End test_pobyso_get_list_elements    

def test_pobyso_taylor_expansion_no_change_var_so_so(repeat=1000, number=10):

    functionName = inspect.stack()[0][3]

    print "Running", functionName, "..."

    def test(functionSo, degreeSo, rangeSo, errorTypeSo):

    # End test

    wrapped = test_pobyso_wrapper(test, 

                                  functionSo, 

                                  degreeSo, 

                                  rangeSo, 

                                  errorTypeSo)

    timing = min(timeit.repeat(wrapped, repeat=repeat, number=number))

    print "\t...", functionName, "done."

    return timing

def test_pobyso_taylor_form_so_so_2(repeat=1000, number=10):

    functionName = inspect.stack()[0][3]

    print "Running", functionName, "..."

    functionSo  = pobyso_parse_string('exp(x)')

    degreeSo    = pobyso_constant_from_int_sa_so(20)

    rangeSo     = pobyso_bounds_to_range_sa_so(RR(1),RR(2))

    intervalCenterSo = sollya_lib_mid(rangeSo, None)

    errorTypeSo = pobyso_absolute_so_so()

    def test(functionSo, degreeSo, intervalCenterSo, errorTypeSo):

        taylorformSo = pobyso_taylorform_so_so(functionSo, 

                                               degreeSo, 

                                               intervalCenterSo, 

                                               errorTypeSo)

        sollya_lib_clear_obj(taylorformSo)

    wrapped = test_pobyso_wrapper(test, 

                                  functionSo, 

                                  degreeSo, 

                                  intervalCenterSo, 

                                  errorTypeSo)

    timing = min(timeit.repeat(wrapped, repeat=repeat, number=number))

    sollya_lib_clear_obj(functionSo,None)

    sollya_lib_clear_obj(degreeSo,None)

    sollya_lib_clear_obj(rangeSo,None)

    sollya_lib_clear_obj(intervalCenterSo,None)

    sollya_lib_clear_obj(errorTypeSo,None)

    return timing

# End   test_pobyso_taylor_so_so_2

def test_pobyso_dummy_test(iterationsNum=1, executionsNum=1):

    aString = "2+3"

    def test():

        return 2+3

    wrapped = test_pobyso_wrapper(test)

    timing = timeit.timeit(wrapped, number=executionsNum)

    return timing

def pobyso_clear_taylorform_sa_so(taylorFormSaSo):

    """

    This method is necessary to correctly clean up the memory from Taylor forms.

    These are made of a Sollya object, a Sollya object list, a Sollya object.

    For no clearly understood reason, sollya_lib_clear_object_list crashed 

    when applied to the object list.

    Here, we decompose it into Sage list of Sollya objects references and we

     clear them one by one. 

    """

    sollya_lib_clear_obj(taylorFormSaSo[0])

    (coefficientsErrorsListSaSo, numElementsSa, isEndEllipticSa) = \

        pobyso_get_list_elements_so_so(taylorFormSaSo[1])

    for element in coefficientsErrorsListSaSo:

        sollya_lib_clear_obj(element)

    sollya_lib_clear_obj(taylorFormSaSo[1])

    sollya_lib_clear_obj(taylorFormSaSo[2])

# End pobyso_clear_taylorform_sa_so 

    Get a Sage int from a Sollya int constant.

def pobyso_get_list_elements_so_so(objectListSo):

    listAsSageList = []

                                          objectListSo)

       listAsSageList.append(sollya_lib_copy_obj(listAddress[i]))

       # Clear each of the elements returned by Sollya.

       sollya_lib_clear_obj(listAddress[i])

    # Free the list itself.   

    sollya_lib_free(listAddress)

    return(listAsSageList, numElements.value, isEndElliptic.value)

    # taylorFormListSaSo is a Python list of Sollya objects references that 

    # are copies of the elements of taylorFormSo.

    # pobyso_get_list_elements_so_so clears taylorFormSo.

def pobyso_taylor_expansion_no_change_var_so_so(functionSo, degreeSo, rangeSo,

                                                errorTypeSo=None, 

    intervalCenterSo = sollya_lib_mid(rangeSo, None)

    taylorFormSo = sollya_lib_taylorform(functionSo, degreeSo,

                                         intervalCenterSo,

    sollya_lib_clear_obj(taylorFormSo)

    #polySo = taylorFormListSaSo[0]

    pobyso_clear_taylorform_sa_so(taylorFormListSaSo)

        sollya_lib_free = sollya.sollya_lib_free