/pobysoPythonSage/src/pobyso.py - Diff - Pobyso - Forge du Centre Blaise Pascal

Révision 215 pobysoPythonSage/src/pobyso.py

         return rangeSo
     # End pobyso_bounds_to_range_sa_so
     def pobyso_build_end_elliptic_list_so_so(*args):
         """
         From argumrny Sollya objects, create a Sollya end elliptic list.
         Elements of the list are "eaten" (should not be cleared individualy,
         are cleared when the list is cleared).
         """
         if len(args) == 0:
             ## Called with an empty list produced "error".
             return sollya_lib_build_end_elliptic_list(None)
         index = 0
         ## One can not append elements to an elliptic list, prepend only is
         #  permitted.
         for argument in reversed(args):
             if index == 0:
                 listSo = sollya_lib_build_end_elliptic_list(argument, None)
             else:
                 listSo = sollya_lib_prepend(argument, listSo)
             index += 1
         return listSo
     # End pobyso_build_end_elliptic_list_so_so
     def pobyso_build_function_sub_so_so(exp1So, exp2So):
         return(sollya_lib_build_function_sub(exp1So, exp2So))
-...
         return(imageIntervalSa)
     # End pobyso_compute_pos_function_abs_val_bounds_sa_sa
     def pobyso_compute_precision_decay_ratio_function_sa_so():
         """
         Compute the precision decay ratio function for polynomial
         coefficient progressive trucation.
         """
         functionText = """
         proc(deg, a, b, we, wq)
+        {
           k = we * (exp(x/a)-1) + wq * (b*x)^2 + (1-we-wq) * x;
           return k/k(d);
         };
         """
         return pobyso_parse_string_sa_so(functionText)
     # End  pobyso_compute_precision_decay_ratio_function.
     def pobyso_constant(rnArg):
         """ Legacy function. See pobyso_constant_sa_so. """
         return(pobyso_constant_sa_so(rnArg))
-...
         """
         Obvious.
         """
         return(pobyso_constant_from_int_sa_so(0))
         return pobyso_constant_from_int_sa_so(0)
     def pobyso_constant_1():
         """
         Obvious.
         Legacy function. See pobyso_constant_so_so.
         """
         return(pobyso_constant_1_sa_so())
         return pobyso_constant_1_sa_so()
     def pobyso_constant_1_sa_so():
         """
-...
     def pobyso_constant_from_int(anInt):
         """ Legacy function. See pobyso_constant_from_int_sa_so. """
         return(pobyso_constant_from_int_sa_so(anInt))
         return pobyso_constant_from_int_sa_so(anInt)
     def pobyso_constant_from_int_sa_so(anInt):
         """
         Get a Sollya constant from a Sage int.
         """
         return(sollya_lib_constant_from_int64(long(anInt)))
         return sollya_lib_constant_from_int64(long(anInt))
     def pobyso_constant_from_int_so_sa(constSo):
         """
-...
         """
         constSa = c_long(0)
         sollya_lib_get_constant_as_int64(byref(constSa), constSo)
         return(constSa.value)
         return constSa.value
     # End pobyso_constant_from_int_so_sa
     def pobyso_constant_from_mpq_sa_so(rationalSa):
-...
         currentSollyaPrecSa = pobyso_get_prec_so_sa()
         return pobyso_constant_sa_so(rnArgSa, currentSollyaPrecSa)
     # End pobyso_constant_sollya_prec_sa_so
     def pobyso_end_elliptic_list_so_sa_so(objectsListSo, intCountSa):
         """
         Create a Sollya end elliptic list made of the objectListSo[0] to
          objectsListSo[intCountSa-1] objects.
         """
         return sollya_lib_end_elliptic_list(objectSo, int(intCountSa))
     def pobyso_error_so():
         return sollya_lib_error(None)
     # End pobyso_error().
     def pobyso_evaluate_so_so(funcSo, argumentSo):
         """
         Evaluates funcSo for arguemntSo through sollya_lib_evaluate().
         """
         return sollya_lib_evaluate(funcSo, argumentSo)
     # End pobyso_evaluate_so_so.
     def pobyso_float_poly_sa_so(polySa, precSa = None):
         """
         Create a Sollya polynomial from a Sage RealField polynomial.
-...
         """
         if realFieldSa is None:
             expressionPrecSa = pobyso_get_max_prec_of_exp_so_sa(polySo)
             print "Maximum precision of Sollya polynomial coefficients:", expressionPrecSa
             realFieldSa      = RealField(expressionPrecSa)
         #print "Sollya expression before...",
         #pobyso_autoprint(polySo)
-...
         #print polyRingSa
         # Do not use the polynomial(expressionSa, ring=polyRingSa) form!
         polynomialSa = polyRingSa(expressionSa)
         polyCoeffsListSa = polynomialSa.coefficients()
         #for coeff in polyCoeffsListSa:
         #    print coeff.abs().n()
         return polynomialSa
     # End pobyso_float_poly_so_sa
     def pobyso_free_variable():
         """
         Ultra thin wrapper around the sollya_lib_function_build_free_variable function.
         """
         return sollya_lib_build_function_free_variable()
     def pobyso_function_type_as_string(funcType):
         """ Legacy function. See pobyso_function_type_as_string_so_sa. """
-...
                 return int(degreeIntervalSa.upper())
         # End pobyso_guess_degree_so_sa
     def pobyso_inf_so_so(intervalSo):
         """
         Very thin wrapper around sollya_lib_inf().
         """
         return sollya_lib_inf(intervalSo)
     # End pobyso_inf_so_so.
     def pobyso_infnorm_so_so(func, interval, file = None, intervalList = None):
         print "Do not use this function. User pobyso_supnorm_so_so instead."
         return None
-...
         sollya_lib_clear_obj(precisionSo)
         return precisionSa
     # End pobyso_precision_so_sa
     def pobyso_polynomial_coefficients_progressive_truncate_so_so(polySo,
                                                                   funcSo,
                                                                   icSo,
                                                                   intervalSo,
                                                                   itpSo,
                                                                   ftpSo,
                                                                   maxPrecSo,
                                                                   maxErrSo):
         print "Input arguments:"
         pobyso_autoprint(polySo)
         pobyso_autoprint(funcSo)
         pobyso_autoprint(icSo)
         pobyso_autoprint(intervalSo)
         pobyso_autoprint(itpSo)
         pobyso_autoprint(ftpSo)
         pobyso_autoprint(maxPrecSo)
         pobyso_autoprint(maxErrSo)
         print "________________"
         precDecayFuncSo = \
             pobyso_parse_string("3/8*(exp(x/2.2)-1) + 2/8*(2*x)^2 + (1-3/8-2/8)*x")
         pobyso_autoprint(precDecayFuncSo)
         return None
     def pobyso_range(rnLowerBound, rnUpperBound):
         """ Legacy function. See pobyso_range_sa_so. """
         return pobyso_range_sa_so(rnLowerBound, rnUpperBound)
-...
         return(sollya_lib_remez(funcSo, degreeSo, rangeSo, weightSo, qualitySo, None))
     # End pobyso_remez_canonical_so_so.
     def pobyso_round_coefficients_single_so_so(polySo, precSo):
         """
         Create a rounded coefficients polynomial from polynomial argument to
         the number of bits in size argument.
         All coefficients are set to the same precision.
         """
         ## TODO: check arguments.
         endEllipListSo = pobyso_build_end_elliptic_list_so_so(precSo)
         polySo = sollya_lib_roundcoefficients(polySo, endEllipListSo, None)
         #sollya_lib_clear_obj(endEllipListSo)
         return polySo
     # End pobyso_round_coefficients_single_so_so
     def pobyso_set_canonical_off():
         sollya_lib_set_canonical(sollya_lib_off())
-...
         a = c_int(p)
         precSo = c_void_p(sollya_lib_constant_from_int(a))
         sollya_lib_set_prec(precSo, None)
     # End pobyso_set_prec_sa_so.
     def pobyso_set_prec_so_so(newPrecSo):
         sollya_lib_set_prec(newPrecSo, None)
     # End pobyso_set_prec_so_so.
     def pobyso_inf_so_so(intervalSo):
         """
         Very thin wrapper around sollya_lib_inf().
         """
         return sollya_lib_inf(intervalSo)
     # End pobyso_inf_so_so.
     def pobyso_supnorm_so_so(polySo, funcSo, intervalSo, errorTypeSo = None,\
                              accuracySo = None):
         """
-...
                                              intervalSo, errorTypeSo, None)
         if createdErrorType:
             sollya_lib_clear_obj(errorTypeSo)
         return(resultSo)
         return resultSo
     def pobyso_univar_polynomial_print_reverse(polySa):

Formats disponibles : Unified diff

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

Révision 215 pobysoPythonSage/src/pobyso.py