Révision 227 pobysoPythonSage/src/sageSLZ/sageRunSLZ.sage
| sageRunSLZ.sage (revision 227) | ||
|---|---|---|
| 2966 | 2966 |
2- once a non null resultant is found, check for roots; |
| 2967 | 2967 |
3- constant resultant == no root. |
| 2968 | 2968 |
""" |
| 2969 |
|
|
| 2970 | 2969 |
if debug: |
| 2971 | 2970 |
print "Function :", inputFunction |
| 2972 | 2971 |
print "Lower bound :", inputLowerBound |
| ... | ... | |
| 3046 | 3045 |
internalSollyaPrec = ceil((RR('1.5') * targetHardnessToRound) / 64) * 64
|
| 3047 | 3046 |
if internalSollyaPrec < 192: |
| 3048 | 3047 |
internalSollyaPrec = 192 |
| 3048 |
pobyso_lib_init() |
|
| 3049 | 3049 |
pobyso_set_prec_sa_so(internalSollyaPrec) |
| 3050 | 3050 |
print "Sollya internal precision:", internalSollyaPrec |
| 3051 | 3051 |
targetPlusOnePrecRF = RealField(RRR.prec()+1) |
| ... | ... | |
| 3108 | 3108 |
iterCount += 1 |
| 3109 | 3109 |
print "[", lb, ",", ub, "]", ((ub - lb) / intervalUlp).log2().n(), \ |
| 3110 | 3110 |
"log2(numbers)." |
| 3111 |
#print "Debugging..." |
|
| 3111 | 3112 |
### Compute a Sollya polynomial that will honor the Taylor condition. |
| 3112 |
prceSo = slz_compute_polynomial_and_interval_01(scaledfSo,
|
|
| 3113 |
prceSo = slz_compute_polynomial_and_interval_02(scaledfSo,
|
|
| 3113 | 3114 |
degreeSo, |
| 3114 | 3115 |
lb, |
| 3115 | 3116 |
ub, |
| 3116 | 3117 |
polyApproxAccur, |
| 3117 | 3118 |
debug=True) |
| 3118 | 3119 |
if debug: |
| 3119 |
print "Sollya Taylor polynomial:", pobyso_autoprint(prceSo[0]) |
|
| 3120 |
print "Sollya interval :", pobyso_autoprint(prceSo[1]) |
|
| 3121 |
print "Sollya interval center :", pobyso_autoprint(prceSo[2]) |
|
| 3122 |
print "Sollya Taylor error :", pobyso_autoprint(prceSo[3]) |
|
| 3120 |
print "Sollya Taylor polynomial:", ; pobyso_autoprint(prceSo[0])
|
|
| 3121 |
print "Sollya interval :", ; pobyso_autoprint(prceSo[1])
|
|
| 3122 |
print "Sollya interval center :", ; pobyso_autoprint(prceSo[2])
|
|
| 3123 |
print "Sollya Taylor error :", ; pobyso_autoprint(prceSo[3])
|
|
| 3123 | 3124 |
|
| 3124 | 3125 |
### Convert back the data into Sage space. |
| 3125 | 3126 |
(floatP, floatPcv, intvl, ic, terr) = \ |
Formats disponibles : Unified diff