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

def slz_get_intervals_and_polynomials(functionSa, degreeSa, 

                                      lowerBoundSa, 

    """

    Under the assumption that:

    - functionSa is monotonic on the [lowerBoundSa, upperBoundSa] interval;

    # Set Sollya to the necessary internal precision.

    sollyaPrecChangedSa = False

    (initialSollyaPrecSo, initialSollyaPrecSa) = pobyso_get_prec_so_so_sa()

        if internalSollyaPrecSa <= 2:

            print inspect.stack()[0][3], ": precision change <=2 requested."

        pobyso_set_prec_sa_so(internalSollyaPrecSa)

        # If the error and interval are OK a the first try, just return.

        if (boundsSa.endpoints()[1] >= scaledUpperBoundSa) and \

            (errorSa <= approxAccurSa):

                pobyso_set_prec_so_so(initialSollyaPrecSo)

            sollya_lib_clear_obj(initialSollyaPrecSo)

            sollya_lib_clear_obj(functionSo)

    return((polynomialSa, polynomialChangedVarSa, intervalSa, centerSa, errorSa))

    # End slz_interval_and_polynomial_to_sage

    """

      Check if an element (argument) of the domain of function (function)

      yields a HTRN case (rounding to next) for the target precision