Consider this: higgsino:~/projects/GiNaC/ginsh$ ./ginsh ginsh - GiNaC Interactive Shell (GiNaC V0.5.2) __, _______ Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, (__) * | Germany. This is free software with ABSOLUTELY NO WARRANTY. ._) i N a C | You are welcome to redistribute it under certain conditions; <-------------' see the file COPYING for details. Type ?? for a list of help topics.
series(gamma(x),x,0,2); collect(series(gamma(x),x,0,2),x); -EulerGamma+x*(1/12*Pi^2+1/2*EulerGamma^2)+x^(-1)+Order(x^2) -EulerGamma+x*(1/12*Pi^2+1/2*EulerGamma^2)+x^2*Order(1)+x^(-1)
The other bug is with the fraction free elimination in ginac/matrix.cpp: If I try to solve the system [[ [[0,2-a]], [[5-4*a,3+2*a]] ]] * [[ [[x]], [[y]] ]] == [[ [[-2+5*a]], [[-5-2*a]] ]] I run into trouble within assert-blocks. I will have a look at the second problem. -rbk. -- Richard Kreckel <Richard.Kreckel@Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>
Hi! On Mon, 21 Feb 2000, Richard B. Kreckel wrote:
series(gamma(x),x,0,2); collect(series(gamma(x),x,0,2),x); -EulerGamma+x*(1/12*Pi^2+1/2*EulerGamma^2)+x^(-1)+Order(x^2) -EulerGamma+x*(1/12*Pi^2+1/2*EulerGamma^2)+x^2*Order(1)+x^(-1)
Fixed. Bye, Christian -- / Coding on PowerPC and proud of it \/ http://www.uni-mainz.de/~bauec002/
On Mon, 21 Feb 2000, Richard B. Kreckel wrote:
The other bug is with the fraction free elimination in ginac/matrix.cpp: If I try to solve the system [[ [[0,2-a]], [[5-4*a,3+2*a]] ]] * [[ [[x]], [[y]] ]] == [[ [[-2+5*a]], [[-5-2*a]] ]]
Looks like the reason is a problem in normal(): symbol a("a"); ex z=((-5-2*a)-((2-5*a)/(-2+a))*(3+2*a))/(5-4*a); cout << "z=" << z << endl; cout << "normal(z)=" << normal(z) << endl; z=(5-4*a)^(-1)*(-5-2*a-(3+2*a)*(2-5*a)*(-2+a)^(-1)) normal(z)=-(5-4*a)^(-1)*(-22+2*a*(-2+a)-10*a^2)*(-2+a)^(-1) -- Alexander Frink E-Mail: Alexander.Frink@Uni-Mainz.DE Institut fuer Physik Phone: +49-6131-3923391 Johannes-Gutenberg-Universitaet D-55099 Mainz, Germany
Hi! On Mon, 21 Feb 2000, Alexander Frink wrote:
z=(5-4*a)^(-1)*(-5-2*a-(3+2*a)*(2-5*a)*(-2+a)^(-1)) normal(z)=-(5-4*a)^(-1)*(-22+2*a*(-2+a)-10*a^2)*(-2+a)^(-1)
Fixed. Bye, Christian -- / Coding on PowerPC and proud of it \/ http://www.uni-mainz.de/~bauec002/
participants (3)
-
Alexander Frink
-
Christian Bauer
-
Richard B. Kreckel