pseries::evalf() and pseries::eval()
Dear All, GiNaC tutorial in section "2.2 What it can do for you" says that Ginsh shall produce: > series(tgamma(x),x==0,3); x^(-1)-Euler+(1/12*Pi^2+1/2*Euler^2)*x+ (-1/3*zeta(3)-1/12*Pi^2*Euler-1/6*Euler^3)*x^2+Order(x^3) > evalf(%); x^(-1)-0.5772156649015328606+(0.9890559953279725555)*x -(0.90747907608088628905)*x^2+Order(x^3) However, with the current GiNaC versions the second output will be identical to the first. I am attaching a small patch which shall restore the expected output. Best wishes, Vladimir -- Vladimir V. Kisil http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Transformations http://goo.gl/EaG2Vu Software: Geometry of cycles http://moebinv.sourceforge.net/ Jupyter (Colab): https://github.com/vvkisil/MoebInv-notebooks Jupyter (CodeOcean): https://codeocean.com/capsule/7952650/tree
Dear Vladimir, On 29.09.19 22:35, Vladimir V. Kisil wrote:
GiNaC tutorial in section "2.2 What it can do for you" says that Ginsh shall produce:
> series(tgamma(x),x==0,3); x^(-1)-Euler+(1/12*Pi^2+1/2*Euler^2)*x+ (-1/3*zeta(3)-1/12*Pi^2*Euler-1/6*Euler^3)*x^2+Order(x^3) > evalf(%); x^(-1)-0.5772156649015328606+(0.9890559953279725555)*x -(0.90747907608088628905)*x^2+Order(x^3)
However, with the current GiNaC versions the second output will be identical to the first. I am attaching a small patch which shall restore the expected output.
Thanks for reporting this! It is a regression introduced in f8c2455fbb. I assume you just added the eval() method to it.rest in pseries::eval() because of symmetry? There shouldn't be any need to do so since 1.7.0. All my best, -richy. -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>
On Tue, 1 Oct 2019 23:55:00 +0200, "Richard B. Kreckel" <kreckel@in.terlu.de> said: RK> I assume you just added the eval() method to it.rest in RK> pseries::eval() because of symmetry? There shouldn't be any need RK> to do so since 1.7.0.
Yes, I did it formally... -- Vladimir V. Kisil http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Transformations http://goo.gl/EaG2Vu Software: Geometry of cycles http://moebinv.sourceforge.net/ Jupyter (Colab): https://github.com/vvkisil/MoebInv-notebooks Jupyter (CodeOcean): https://codeocean.com/capsule/7952650/tree
participants (2)
-
Richard B. Kreckel
-
Vladimir V. Kisil