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