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

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

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

    return sollya_lib_build_function_sub(exp1So, exp2So)

def pobyso_build_list_of_ints_sa_so(*args):

    """

    Build a Sollya list from Sage integral arguments. 

    """

    if len(args) == 0:

        return pobyso_build_list_so_so()

    ## Make a Sage list of integral constants.

    intsList = []

    for intElem in args:

        intsList.append(pobyso_constant_from_int_sa_so(intElem))

    return pobyso_build_list_so_so(*intsList) 

def pobyso_build_list_so_so(*args):

    """

    Make a Sollya list out of Sollya objects passed as arguments.

    If one wants to call it with a list argument, should prepend a "*"

    before the list variable name.

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

    are cleared when the list is cleared).

    """

    if len(args) == 0:

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

        return sollya_lib_build_list(None)

    index = 0

    ## Append the Sollya elements one by one.

    for elementSo in args:

        if index == 0:

            listSo = sollya_lib_build_list(elementSo, None)

        else:

            listSo = sollya_lib_append(listSo, elementSo)

        index += 1

    return listSo

# End pobyso_build list_so_so

def pobyso_clear_obj(objSo):

    """

    Free a Sollya object's memory.

    Very thin wrapper around sollya_lib_clear_obj().

    """

    sollya_lib_clear_obj(objSo)

# End pobyso_clear_obj. 

     clear them one at a time. 

    Deprecated, use pobyso_cmp_sa_so_sa instead.

    """

    print "Deprecated, use pobyso_cmp_sa_so_sa instead."

    return pobyso_cmp_sa_so_sa(rnArgSa, cteSo)

# End  pobyso_cmp

def pobyso_cmp_sa_so_sa(rnArgSa, cteSo):

    """

    ## Compare the Sage RealNumber version of the Sollya constant with rnArg

    #  through a direct comparison of underlying MPFR numbers.

    return cmp_rn_value(rnArgSa, rnLocal)

# End pobyso_smp_sa_so_sa

def pobyso_dirty_inf_norm_so_so(funcSo, rangeSo, preSo=None):

    """

    Thin wrapper around sollya_dirtyinfnorm().

    """

    # TODO: manage the precision change.

    return sollya_lib_dirtyinfnorm(funcSo, rangeSo)

# End pobyso_dirty_inf_norm_so_so

def pobyso_remez_exponents_list_sa_so(func,     \

                                 exponentsList, \

                                 lowerBound,    \

                                 upperBound,    \

                                 weight = None, \

                                 quality = None):

    """

    All arguments are Sage/Python.

    The functions (func and weight) must be passed as expressions or strings.

    Otherwise the function fails. 

    The return value is a pointer to a Sollya function.

    lowerBound and upperBound mus be reals.

    """

    var('zorglub')    # Dummy variable name for type check only. Type of

    # zorglub is "symbolic expression".

    currentVariableNameSa = None

    # The func argument can be of different types (string, 

    # symbolic expression...)

    if parent(func) == parent("string"):

        localFuncSa = sage_eval(func,globals())

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

            currentVariableNameSa = localFuncSa.variables()[0]

            sollya_lib_name_free_variable(str(currentVariableNameSa))

            functionSo = \

              sollya_lib_parse_string(localFuncSa._assume_str().replace('_SAGE_VAR_', ''))

    # Expression test.

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

        # Until we are able to translate Sage expressions into Sollya 

        # expressions : parse the string version.

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

            currentVariableNameSa = func.variables()[0]

            sollya_lib_name_free_variable(str(currentVariableNameSa))

            functionSo = \

              sollya_lib_parse_string(func._assume_str().replace('_SAGE_VAR_', ''))

    else:

        return(None)

    ## Deal with the weight, much in the same way as with the function.

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

        weightSo = pobyso_constant_1_sa_so()

    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.

        functionSo = \

          sollya_lib_parse_string(weight._assume_str().replace('_SAGE_VAR_', '',100))

    else:

        return(None)

    rangeSo = pobyso_bounds_to_range_sa_so(lowerBound, upperBound)

    if not quality is None:

        qualitySo= pobyso_constant_sa_so(quality)

    else:

        qualitySo = None

    #

    ## Tranform the Sage list of exponents into a Sollya list.

    exponentsListSo = pobyso_build_list_of_ints_sa_so(*exponentsList)

    remezPolySo = sollya_lib_remez(functionSo, \

                                   exponentsListSo, \

                                   rangeSo, \

                                   weightSo, \

                                   qualitySo, \

                                   None)

    sollya_lib_clear_obj(functionSo)

    sollya_lib_clear_obj(exponentsListSo)

    sollya_lib_clear_obj(rangeSo)

    sollya_lib_clear_obj(weightSo)

    if not qualitySo is None:

        sollya_lib_clear_obj(qualitySo)

    return(remezPolySo)

# End pobyso_remez_exponentsList_sa_so

        # Notice the **: we are in Pythonland here!