Thanks very much for the quick fix, thanks again!
Feng Feng,
On 13.04.19 07:06, Feng Feng wrote:
I am sorry for that I didn't make it clear.
the test program shows that s1 output is Order(ep^(-1)), there is no other term in the output, so s1 output indicates that there is no ep^(-2) in expr, I am not sure I interpret the s1 output correctly or not?
I'm sure you interpret it correctly. There is indeed a bug in GiNaC.
A reduced example of the same problem is this:
symbol x("x"), ep("?"); ex expr = x * (tgamma(ep) / ep - sin(ep) / ep);
cout << expr.series(ep, -2) << endl; // good: Order(?^(-2))
cout << expr.series(ep, -1) << endl; // D'oh: Order(?^(-1))
cout << expr.series(ep, 0) << endl; // good: (x)*?^(-2)+(-Euler*x)*?^(-1)+Order(1)
It seems it was broken by commit db81420a61 in 2004.
Best, -richy.
------------------------------
Subject: Digest Footer
_______________________________________________ GiNaC-list mailing list GiNaC-list@ginac.de https://www.cebix.net/mailman/listinfo/ginac-list
------------------------------
End of GiNaC-list Digest, Vol 140, Issue 8 ******************************************