[GiNaC-devel] pseries::evalf() and pseries::eval()

[Richard B. Kreckel]

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/>