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

import sys

from ctypes import *

from sage.rings.real_mpfr import *

#

## Load the Sollya library.

try:

    #sollya=CDLL("/usr/local/lib/libsollya.so")

    sollya=CDLL("/usr/local/lib/libsollya.so")    

except :

    print "\nCould not find nor load the Sollya library.\n"

    sys.exit(1)

#

try:

    # Export the functions with their name in the library (so we can

    # use them without the "sollya." prefix).

        sollya_lib_add = sollya.sollya_lib_add

        sollya_lib_abs = sollya.sollya_lib_abs

        sollya_lib_absolute = sollya.sollya_lib_absolute

        sollya_lib_append = sollya.sollya_lib_append

        sollya_lib_autoprint = sollya.sollya_lib_autoprint

        sollya_lib_build_end_elliptic_list = \

            sollya.sollya_lib_build_end_elliptic_list

        sollya_lib_build_function_abs = sollya.sollya_lib_build_function_abs

        sollya_lib_build_function_add = sollya.sollya_lib_build_function_add

        sollya_lib_build_function_div = sollya.sollya_lib_build_function_div

        sollya_lib_build_function_exp = sollya.sollya_lib_build_function_exp

        sollya_lib_build_function_free_variable = \

            sollya.sollya_lib_build_function_free_variable

        sollya_lib_build_function_mul = sollya.sollya_lib_build_function_mul

        sollya_lib_build_function_pow = sollya.sollya_lib_build_function_pow

        sollya_lib_build_function_sub = sollya.sollya_lib_build_function_sub

        sollya_lib_build_list = sollya.sollya_lib_build_list

        sollya_lib_chebyshevform = sollya.sollya_lib_chebyshevform

        sollya_lib_clear_obj = sollya.sollya_lib_clear_obj

        sollya_lib_clear_object_list = sollya.sollya_lib_clear_object_list

        sollya_lib_close = sollya.sollya_lib_close

        sollya_lib_coeff = sollya.sollya_lib_coeff

        sollya_lib_constant = sollya.sollya_lib_constant

        sollya_lib_constant_from_int = sollya.sollya_lib_constant_from_int

        sollya_lib_constant_from_int64 = sollya.sollya_lib_constant_from_int64

        sollya_lib_constant_from_mpq = sollya.sollya_lib_constant_from_mpq 

        sollya_lib_constant_from_mpz = sollya.sollya_lib_constant_from_mpz 

        sollya_lib_constant_from_uint64 = \

            sollya.sollya_lib_constant_from_uint64

        sollya_lib_copy_obj = sollya.sollya_lib_copy_obj

        sollya_lib_cos = sollya.sollya_lib_cos

        sollya_lib_degree = sollya.sollya_lib_degree

        sollya_lib_diff = sollya.sollya_lib_diff

        sollya_lib_dirtyfindzeros = \

            sollya.sollya_lib_dirtyfindzeros

        sollya_lib_dirtyinfnorm = sollya.sollya_lib_dirtyinfnorm

        sollya_lib_end_elliptic_list= \

            sollya.sollya_lib_end_elliptic_list

        sollya_lib_error=sollya.sollya_lib_error

        sollya_lib_evaluate = sollya.sollya_lib_evaluate

        sollya_lib_findzeros = sollya.sollya_lib_findzeros

        sollya_lib_free = sollya.sollya_lib_free

        sollya_lib_get_canonical = sollya.sollya_lib_get_canonical

        sollya_lib_get_constant = sollya.sollya_lib_get_constant

        sollya_lib_get_constant_as_int = \

            sollya.sollya_lib_get_constant_as_int

        sollya_lib_get_constant_as_int64 = \

            sollya.sollya_lib_get_constant_as_int64

        sollya_lib_get_constant_as_uint64 = \

            sollya.sollya_lib_get_constant_as_uint64

        sollya_lib_get_free_variable_name = \

        sollya.sollya_lib_get_free_variable_name

        sollya_lib_get_function_arity = sollya.sollya_lib_get_function_arity

        sollya_lib_get_head_function = sollya.sollya_lib_get_head_function

        sollya_lib_get_interval_from_range = \

            sollya.sollya_lib_get_interval_from_range

        sollya_lib_get_list_elements = sollya.sollya_lib_get_list_elements

        sollya_lib_get_object_list_head = \

            sollya.sollya_lib_get_object_list_head

        sollya_lib_get_object_list_tail = \

            sollya.sollya_lib_get_object_list_tail

        sollya_lib_get_prec = sollya.sollya_lib_get_prec

        sollya_lib_get_prec_of_constant = sollya.sollya_lib_get_prec_of_constant

        sollya_lib_get_prec_of_range = sollya.sollya_lib_get_prec_of_range

        sollya_lib_get_subfunctions = sollya.sollya_lib_get_subfunctions

        sollya_lib_guessdegree = sollya.sollya_lib_guessdegree

        sollya_lib_head = sollya.sollya_lib_head

        sollya_lib_inf = sollya.sollya_lib_inf

        sollya_lib_infnorm = sollya.sollya_lib_infnorm

        sollya_lib_init= sollya.sollya_lib_init

        sollya_lib_is_absolute = sollya.sollya_lib_is_absolute

        sollya_lib_is_empty_object_list = \

            sollya.sollya_lib_is_empty_object_list

        sollya_lib_is_off = sollya.sollya_lib_is_off

        sollya_lib_is_on = sollya.sollya_lib_is_on

        sollya_lib_obj_is_range = sollya.sollya_lib_obj_is_range

        sollya_lib_max = sollya.sollya_lib_max

        sollya_lib_mid = sollya.sollya_lib_mid

        sollya_lib_name_free_variable = \

            sollya.sollya_lib_name_free_variable

        sollya_lib_obj_is_function = sollya.sollya_lib_obj_is_function

        sollya_lib_obj_is_list = sollya.sollya_lib_obj_is_list

        sollya_lib_obj_is_error = sollya.sollya_lib_obj_is_error

        sollya_lib_obj_is_range = sollya.sollya_lib_obj_is_range

        sollya_lib_obj_is_structure = sollya.sollya_lib_obj_is_structure

        sollya_lib_off = sollya.sollya_lib_off

        sollya_lib_on = sollya.sollya_lib_on

        sollya_lib_max = sollya.sollya_lib_max

        sollya_lib_parse_string = sollya.sollya_lib_parse_string

        sollya_lib_precision = sollya.sollya_lib_precision

        sollya_lib_prepend = sollya.sollya_lib_prepend

        sollya_lib_range = sollya.sollya_lib_range

        sollya_lib_range_from_bounds = sollya.sollya_lib_range_from_bounds

        sollya_lib_range_from_interval = \

            sollya.sollya_lib_range_from_interval

        sollya_lib_relative = sollya.sollya_lib_relative

        sollya_lib_remez = sollya.sollya_lib_remez

        sollya_lib_roundcoefficients = \

            sollya.sollya_lib_roundcoefficients

        sollya_lib_set_canonical = sollya.sollya_lib_set_canonical

        sollya_lib_set_prec = sollya.sollya_lib_set_prec

        sollya_lib_sub = sollya.sollya_lib_sub

        sollya_lib_sup = sollya.sollya_lib_sup

        sollya_lib_supnorm = sollya.sollya_lib_supnorm

        sollya_lib_tail = sollya.sollya_lib_tail

        sollya_lib_taylor = sollya.sollya_lib_taylor

        sollya_lib_taylorform = sollya.sollya_lib_taylorform

except :

    print "\nCould not rename one of the functions.\n"

    sys.exit(1)

# Set the return type of several functions (those that have a return

# type different from c_int or void). We consider the all the functions

# returning a sollya_obj_t as returning an int.

try:

    sollya_lib_get_free_variable_name.restype = c_char_p

except :

    sys.stderr.write("\nOne of the Python-Sollya return type setting command \n")

    sys.stderr.write("has failed.\n\n")

    sys.exit(1)

#

# Set the argument type of several functions. Default arguments are int.

#

try:

    sollya_lib_add.argtypes                  = [c_ulong, c_ulong]

    sollya_lib_append.argstypes              = [c_ulong, c_ulong]

    #sollya_lib_autoprint.argtypes            = [c_ulong]

    sollya_lib_build_function_add.argtypes   = [c_int, c_int]

    sollya_lib_constant.argtypes             = [c_ulong]

    sollya_lib_constant_from_int64.argtypes  = [c_long]

    sollya_lib_constant_from_mpq.argtypes    = [c_ulong]

    sollya_lib_constant.restype              = c_ulong

    sollya_lib_dirtyinfnorm.argstypes        = [c_ulong, c_ulong]

    sollya_lib_end_elliptic_list.argtypes    = [c_ulong, c_int]

    sollya_lib_get_constant.argtypes         = [c_ulong, c_ulong]

    sollya_lib_get_constant_as_int.argtypes  = [POINTER(c_int), c_int]

    sollya_lib_get_constant_as_int64.argtypes= [POINTER(c_long), c_int]

    sollya_lib_get_function_arity.argtypes   = [POINTER(c_int), c_int]

    sollya_lib_get_head_function.argtypes    = [POINTER(c_int), c_ulong]

    sollya_lib_get_interval_from_range.argtypes = [c_ulong, c_ulong]

    sollya_lib_get_list_elements.argtypes    = [POINTER(POINTER(c_longlong)), \

                                                POINTER(c_int),\

                                                POINTER(c_int), c_int]

    sollya_lib_get_prec_of_constant.argtypes = [POINTER(c_int), c_ulong]

    sollya_lib_get_prec_of_range.argtypes    = [POINTER(c_int), c_int]

    sollya_lib_get_subfunctions.argtypes     = [c_int, POINTER(c_int), \

            POINTER(c_int), POINTER(c_int), POINTER(c_int), POINTER(c_int), \

            POINTER(c_int), POINTER(c_int), POINTER(c_int), POINTER(c_int), \

            POINTER(c_int), POINTER(c_int)]

    sollya_lib_guessdegree.argtypes          = [c_ulong, c_ulong, c_ulong, \

            c_ulong, c_ulong]

    sollya_lib_guessdegree.restype           = c_ulong        

    sollya_lib_name_free_variable.argtypes   = [POINTER(c_char)]

    sollya_lib_parse_string.restype          = c_long

    sollya_lib_prepend.argstypes             = [c_ulong, c_ulong]

    sollya_lib_range.argtypes                = [c_ulong, c_ulong]

    sollya_lib_range.restype                 = c_ulong

    sollya_lib_range_from_bounds.argtypes    = [c_ulong, c_ulong]

    sollya_lib_range_from_bounds.restype     = c_ulong

    sollya_lib_range_from_interval.argtypes  = [c_ulong]

    sollya_lib_range_from_interval.restype   = c_ulong

except :

    sys.stderr.write("\nOne of the Python-Sollya return type setting command \n")

    sys.stderr.write("has failed.\n\n")

    sys.exit(1)

#

# Give it a try!

#

sollya

sollya.sollya_lib_init(None)

sys.stderr.write("\nSuperficial Sollya library check...\n")

#print sollya.sollya_lib_get_free_variable_name()

sollyaExp = sollya_lib_parse_string("exp(x)")

#retc = sollya_lib_autoprint(sollyaExp, None)

arg = sollya_lib_constant_from_int(int(1))

#print "Type of sollya_obj_t: ",type(arg), "Value: ", arg

res = sollya_lib_evaluate(sollyaExp, arg)

#retc = sollya_lib_autoprint(res, None)

retc = sollya_lib_get_prec(None)

a = c_int(0)

sollya_lib_get_constant_as_int(byref(a), retc)

b = a.value

#print "Precision : ", a.value

#print "Address of a.value ", addressof(a)

#sollya_lib_close(None)

sys.stderr.write("\t...Sollya library check done.\n")

#

#

#

#b = RealNumber()

#print cast(retc, POINTER(c_int)).contents

#retc = cast(retc, POINTER(c_int)).contents

#arg = sollya_lib_constant_from_int(int(retc))

#retc = sollya_lib_autoprint(arg, None)

#sollya_lib_set_prec(c_int(100))