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

    maxErrorSa = pobyso_get_constant_as_rn_with_rf_so_sa(maxErrorSo)

    while maxErrorSa > approxPrecSa:

        #print "++Approximation error:", maxErrorSa

        sollya_lib_clear_obj(polySo)

        sollya_lib_clear_obj(intervalCenterSo)

        shrinkFactorSa = RRR('5')/(maxErrorSa/approxPrecSa).log2().abs()

        #shrinkFactorSa = 1.5/(maxErrorSa/approxPrecSa)

        #errorRatioSa = approxPrecSa/maxErrorSa

                                   actualShrinkFactorSa

        #print "Current upper bound:", currentUpperBoundSa

        sollya_lib_clear_obj(currentRangeSo)

        if currentUpperBoundSa <= currentLowerBoundSa or \

              currentUpperBoundSa == currentLowerBoundSa.nextabove():

            sollya_lib_clear_obj(absoluteErrorTypeSo)

                                                      currentUpperBoundSa)

        # print "New interval:",

        # pobyso_autoprint(currentRangeSo)

        (polySo, intervalCenterSo, maxErrorSo) = \

            pobyso_taylor_expansion_no_change_var_so_so(functionSo, degreeSo, 

                                                        currentRangeSo, 

    - the corresponding approximation error.

    TODO: fix endless looping for some parameters sets.

    """

    # Set Sollya to the necessary internal precision.

    currentSollyaPrecSo = pobyso_get_prec_so()

    currentSollyaPrecSa = pobyso_constant_from_int_so_sa(currentSollyaPrecSo)

    if internalSollyaPrecSa > currentSollyaPrecSa:

        pobyso_set_prec_sa_so(internalSollyaPrecSa)

    #

    x = functionSa.variables()[0] # Actual variable name can be anything.

    # Scaled function: [1=,2] -> [1,2].

                                            upperBoundSa,                     \

                                            floatingPointPrecSa)

    #

    print "Approximation precision: ", RR(approxPrecSa)

    # Prepare the arguments for the Taylor expansion computation with Sollya.

    functionSo = pobyso_parse_string_sa_so(fff._assume_str())

                                        approxPrecSa, internalSollyaPrecSa)

    if polySo is None:

        print "slz_get_intervals_and_polynomials: Aborting and returning None!"

            sollya_lib_clear_obj(currentSollyaPrecSo)

        sollya_lib_clear_obj(functionSo)

        sollya_lib_clear_obj(degreeSo)

        sollya_lib_clear_obj(scaledBoundsSo)

        #print "Approximation error:", errorSa

        return resultArray

    # The following ratio is always >= 1

    currentErrorRatio = approxPrecSa / errorSa

    # Starting point for the next upper bound.

    return timing

# End test_pobyso_bounds_to_range_sa_so

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

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

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

    return timing

# End test_pobyso_get_list_elements    

        precSa = pobyso_get_prec_so_sa()

# End test_pobyso_get_prec_so_sa  

        precSa = pobyso_get_prec_of_constant_so_sa(constSo)

    sollya_lib_clear_obj(constSo)

# End test_pobyso_get_prec_of_constant_so_sa

    return timing

# End test_pobyso_taylor_expansion_no_change_var_so_so

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

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

    functionSo  = pobyso_parse_string('exp(x)')

    wrapped = test_pobyso_wrapper(test)

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

    return timing

    # pobyso_get_list_elements_so_so clears taylorFormSo.

    (taylorFormListSaSo, numElementsSa, isEndEllipticSa) = \

        pobyso_get_list_elements_so_so(taylorFormSo)

    sollya_lib_clear_obj(taylorFormSo)

    #polySo = taylorFormListSaSo[0]

    #errorRangeSo = sollya_lib_copy_obj(taylorFormListSaSo[2])

    errorRangeSo = taylorFormListSaSo[2]