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

    TODO: completely rework and test.

    """

    pobyso = pobyso_name_free_variable_sa_so(funcSa.variables()[0])

    rangeSo = pobyso_range_sa_so(lowerBoundSa, upperBoundSa)

    infnormSo = pobyso_infnorm_so_so(funcSo,rangeSo)

    # Sollya return the infnorm as an interval.

    binadeTopLimitSa = 2**ceil(fMaxUpperBoundSa.log2())

    # Put up together the function to use to compute the lower bound.

    funcAuxSo = pobyso_parse_string(str(binadeTopLimitSa) +  \

    pobyso_autoprint(funcAuxSo)

    # Clear the Sollya range before a new call to infnorm and issue the call.

    sollya_lib_clear_obj(infnormSo)

    Sa_sa variant of the solly_guessdegree function.

    Return 0 if something goes wrong.

    """

    functionSo = pobyso_parse_string_sa_so(functionAsStringSa)

    if pobyso_is_error_so_sa(functionSo):

        sollya_lib_clear_obj(functionSo)

    else:

        approxErrorSo = pobyso_constant_sa_so(RR(approxErrorSa))

    if not weightSa is None:

        weightSo = pobyso_parse_string_sa_so(weightAsStringSa)

        if pobyso_is_error(weightSo):

            sollya_lib_clear_obj(functionSo)

        if len(localFuncSa.variables()) > 0:

            currentVariableNameSa = localFuncSa.variables()[0]

            sollya_lib_name_free_variable(str(currentVariableNameSa))

    # Expression test.

    elif type(func) == type(zorglub):

        # Until we are able to translate Sage expressions into Sollya 

        if len(func.variables()) > 0:

            currentVariableNameSa = func.variables()[0]

            sollya_lib_name_free_variable(str(currentVariableNameSa))

    else:

        return(None)

    if weight is None: # No weight given -> 1.

    elif parent(weight) == parent("string"): # Weight given as string: parse it.

        weightSo = sollya_lib_parse_string(func)

    elif type(weight) == type(zorglub): # Weight given as symbolice expression.

    else:

        return(None)

    degreeSo = pobyso_constant_from_int(degree)

    else: # Interval.

        pointSo = pobyso_interval_to_range_sa_so(pointSa, pointPrecSa)

    # Sollyafy the function.

    if sollya_lib_obj_is_error(functionSo):

        print "pobyso_tailorform: function string can't be parsed!"

        return None

print "sagePolynomialOperations loading..."

def spo_add_polynomial_coeffs_to_matrix_row(poly, 

                                            knownMonomials, 

    #

    print "Approximation precision: ", RR(approxPrecSa)

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

    degreeSo = pobyso_constant_from_int_sa_so(degreeSa)

    scaledBoundsSo = pobyso_bounds_to_range_sa_so(scaledLowerBoundSa, 

                                                  scaledUpperBoundSa)

    except Exception:

        print "denominators:", rationalList, \

                "does not understand the len() function."

    if listLength == 0:

    if listLength == 1:

        try:

            return([QQ(rationalList[0]).denominator()])

        except Exception:

            print "denominators:", rationalList[0], \

                    "has no \"denominator()\" member." 

    denominatorsList = []    

    for i in xrange(listLength):

        try:

        except Exception:

            print "denominators:", rationalList[i], \

                    "has no \"denominator()\" member." 

    return(denominatorsList)

# End sro_denominators

    except Exception:

        print "numerators:", rationalList, \

                "does not understand the len() function."

    if listLength == 0:

    if listLength == 1:

        try:

            return(QQ(rationalList[0]).numerator())

        except Exception:

            print "denominators:", rationalList[0], \

                    "has no \"numerator()\" member." 

    numeratorsList = []    

    for i in xrange(listLength):

        try:

        except Exception:

            print "numerators:", rationalList[i], \

                    "has no \"numerator()\" member." 

# End sro_numerators

print "\t...sageRationalOperations loaded"

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

    intervalSa = RIF(1,2)

    functionSo = pobyso_parse_string("exp(x)")

    #RealIntervalField

    pointSo    = pobyso_constant_sa_so(RR("1.5"))

    intervalSo = pobyso_interval_to_range_sa_so(intervalSa)