10/15/2010, 04:03 PM
(This post was last modified: 10/15/2010, 09:47 PM by nuninho1980.)
lastest your code for faster.
I tried bases 1.5, 1.68, 1.7 and 100 and worked (not work on your older code). it's better. but I get errors after try init(1.6) and init(200).
errors here:
base 1.6
base 200
I tried base 1.67 but I get 3 bits precision (last loop (from 2 until 5) is less bits precision than previous loop)
may you fix this code.
I tried bases 1.5, 1.68, 1.7 and 100 and worked (not work on your older code). it's better. but I get errors after try init(1.6) and init(200).
errors here:
base 1.6
PHP Code:
(15:51) gp > init(1.6);loop(1)
base 1.60000000000000000000000
fixed point 1.77824328112897128908095 + 1.46912237320727818425047*I
Pseudo Period 8.97674947562306210989089 + 1.04999719021153010066853*I
3 strm(s) out of 50 sexp(z) generates 36 Riemann samples
2 rtrm(s) out of 18 riemaprx(z) generates 12 sexp samples
*** exp: exponent (expo) overflow
base 200
PHP Code:
(15:46) gp > init(200)
base 200
fixed point -0.167639843463924048976342 + 0.375686229934324598314735*I
Pseudo Period 2.73715856720609094587638 + 1.07145668424439506771849*I
%1 = 1
(15:46) gp > loop(1)
3 strm(s) out of 50 sexp(z) generates 36 Riemann samples
3 rtrm(s) out of 18 riemaprx(z) generates 12 sexp samples
sexp(-0.5)= 0.00071918820319672353153203228317623
0.3364955007 Riemann/sexp binary precision bits I=0.1305261922*I
1.758198575 recenter/renorm 0
(15:46) gp > loop(4)
1 strm(s) out of 12 sexp(z) generates 36 Riemann samples
2 rtrm(s) out of 18 riemaprx(z) generates 12 sexp samples
*** exp: exponent (expo) overflow
I tried base 1.67 but I get 3 bits precision (last loop (from 2 until 5) is less bits precision than previous loop)
may you fix this code.