- GiNaC-devel - ginac.de

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-47-g37d324b
by git＠ginac.de 25 Jun '08

25 Jun '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 37d324b137fcb61384203b711e91d8919535f237 (commit) from 26c31eb89fe38c4d232f90b62c8acbeaaa0454f4 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 37d324b137fcb61384203b711e91d8919535f237 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Wed Jun 25 11:28:11 2008 +0200 Removed left-over debugging line. ----------------------------------------------------------------------- Summary of changes: check/exam_clifford.cpp | 2 -- 1 files changed, 0 insertions(+), 2 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

Patches
by Jens Vollinga 25 Jun '08

25 Jun '08

Hi all, I just put the gcd patch (2008-04-24) and the clifford patch (2008-04-26) into the public repository. For the clifford patch I had to change additionally three lines in clifford.cpp (removing the static keywords). Somehow wrong information is kept there, which causes an exception if more than one clifford_unit with different dimensions is to be created. Regards, Jens

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-46-g26c31eb
by git＠ginac.de 25 Jun '08

25 Jun '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 26c31eb89fe38c4d232f90b62c8acbeaaa0454f4 (commit) via fd4fbeaa409cd807e4f2f7501d278e835e964d69 (commit) from 6850c16435de33396ccd505ff26b4d4574fa74cf (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 26c31eb89fe38c4d232f90b62c8acbeaaa0454f4 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Wed Jun 25 11:17:55 2008 +0200 Make the behaviour of class function more consistent with respect to ncmul::eval(). [V.Kisil] commit fd4fbeaa409cd807e4f2f7501d278e835e964d69 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Tue Jun 24 21:26:58 2008 +0200 Improve heur_gcd() so it can handle rational polynomials, and add a test case. Previously heur_gcd() worked only with integer polynomials, and did not check if the inputs are indeed integer polynomials. That lead to an endless loop or a miscomputed gcd. [A.Sheplyakov] ----------------------------------------------------------------------- Summary of changes: check/Makefile.am | 4 ++ check/exam_clifford.cpp | 17 ++++++++ check/heur_gcd_bug.cpp | 39 +++++++++++++++++++ ginac/clifford.cpp | 11 ++++- ginac/function.pl | 11 +++++ ginac/normal.cpp | 96 ++++++++++++++++++++++++++++++++++++++-------- 6 files changed, 158 insertions(+), 20 deletions(-) create mode 100644 check/heur_gcd_bug.cpp hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, ginac_1-4, updated. release_1-4-0-51-g9968ebf
by git＠ginac.de 25 Jun '08

25 Jun '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, ginac_1-4 has been updated via 9968ebfc60c121b8373f270d79f08d85ca37ccb5 (commit) via 65196b339fbcc7d0e2eb3036fa64e5fa47b56379 (commit) from 314595d10b4a4471bd4637a9038b18dde4780513 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 9968ebfc60c121b8373f270d79f08d85ca37ccb5 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Wed Jun 25 11:17:55 2008 +0200 Make the behaviour of class function more consistent with respect to ncmul::eval(). [V.Kisil] commit 65196b339fbcc7d0e2eb3036fa64e5fa47b56379 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Tue Jun 24 21:26:58 2008 +0200 Improve heur_gcd() so it can handle rational polynomials, and add a test case. Previously heur_gcd() worked only with integer polynomials, and did not check if the inputs are indeed integer polynomials. That lead to an endless loop or a miscomputed gcd. [A.Sheplyakov] ----------------------------------------------------------------------- Summary of changes: check/Makefile.am | 4 ++ check/exam_clifford.cpp | 17 ++++++++ check/heur_gcd_bug.cpp | 39 +++++++++++++++++++ ginac/clifford.cpp | 11 ++++- ginac/function.pl | 11 +++++ ginac/normal.cpp | 96 ++++++++++++++++++++++++++++++++++++++-------- 6 files changed, 158 insertions(+), 20 deletions(-) create mode 100644 check/heur_gcd_bug.cpp hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[PATCH] Function and ncmul
by Vladimir Kisil 26 Apr '08

26 Apr '08

Dear All, Consider the following program: #include <ginac/ginac.h> #include <fstream> using namespace std; using namespace GiNaC; int main(){ realsymbol a("a"); varidx mu(symbol("mu", "\\mu"), 1); ex e = clifford_unit(mu, diag_matrix(lst(-1))), e0=e.subs(mu==0); cout << e0*e0*exp(a*e0) << endl; // -> exp(e0*a) cout << exp(a*e0)*e0*e0 << endl; // -> exp(e0*a)*e0*e0 } The expression in the second output is not simplified for the following reason. When the method ncmul::eval() collects non-commuting terms, the function exp(e0*a) gets the type of its first noncommuting argument, i.e. e0 (this is defined by function::return_type()). So all three factors are of the same type. However, later on, when ncmul::eval() calls the method eval_ncmul() for the first term in the product, e.i. exp(e0*a), it falls back to basic::eval_ncmul() not clifford::eval_ncmul() as one may expect. I am providing below a small patch which makes the behaviour of class function more consistent in this respect. Best wishes, Vladimir -- Vladimir V. Kisil email: kisilv(a)maths.leeds.ac.uk -- www: http://maths.leeds.ac.uk/~kisilv/ Parent: 6850c16435de33396ccd505ff26b4d4574fa74cf (Fixed bug in mLi summation causing premature drop-out and made Nielsen polylog) Branch: master Follows: release_1-4-0 Precedes: Method function::eval_ncmul() is defined to match the behaviour of function::return_type() ------------------------------ ginac/function.pl ------------------------------ index e91278c..3b22143 100644 @@ -507,6 +507,7 @@ public: ex expand(unsigned options=0) const; ex eval(int level=0) const; ex evalf(int level=0) const; + ex eval_ncmul(const exvector & v) const; unsigned calchash() const; ex series(const relational & r, int order, unsigned options = 0) const; ex thiscontainer(const exvector & v) const; @@ -1053,6 +1054,16 @@ ${evalf_switch_statement} throw(std::logic_error("function::evalf(): invalid nparams")); } +/* +* This method is defined to be in line with behaviour of function::return_type() +*/ +ex function::eval_ncmul(const exvector & v) const +{ + // If this function is called then the list of arguments is non-empty + // and the first argument is non-commutative, see function::return_type() + return seq.begin()->eval_ncmul(v); +} + unsigned function::calchash() const { unsigned v = golden_ratio_hash(golden_ratio_hash((p_int)tinfo()) ^ serial);

1 0

Function and ncmul
by Vladimir V. Kisil 26 Apr '08

26 Apr '08

Dear All, Consider the following program: #include <ginac/ginac.h> #include <fstream> using namespace std; using namespace GiNaC; int main(){ realsymbol a("a"); varidx mu(symbol("mu", "\\mu"), 1); ex e = clifford_unit(mu, diag_matrix(lst(-1))), e0=e.subs(mu==0), p1=e0*exp(a*e0), p2=exp(a*e0)*e0; cout << e0*p1 << endl; // -> exp(e0*a) cout << p2*e0 << endl; // -> exp(e0*a)*e0*e0 } The expression in the second output is not simplified for the following reason. When the method ncmul::eval() collects non-commuting the terms exp(e0*a) gets the type of its first noncommuting element, i.e. e0 (this is defined by function::return_type()). So all three factors are in the same type group. However when later ncmul::eval() calls the method eval_ncmul() for the first term in the product, e.i. exp(e0*a), it falls back to basic::eval_ncmul() not clifford::eval_ncmul() as one may expect. I provide below a small patch which makes the behaviour of class function more consistent in this respect. Best wishes, Vladimir -- Vladimir V. Kisil email: kisilv(a)maths.leeds.ac.uk -- www: http://maths.leeds.ac.uk/~kisilv/ Parent: 6850c16435de33396ccd505ff26b4d4574fa74cf (Fixed bug in mLi summation causing premature drop-out and made Nielsen polylog) Branch: master Follows: release_1-4-0 Precedes: Method function::eval_ncmul() is defined to match the behaviour of function::return_type() ------------------------------ ginac/function.pl ------------------------------ index e91278c..3b22143 100644 @@ -507,6 +507,7 @@ public: ex expand(unsigned options=0) const; ex eval(int level=0) const; ex evalf(int level=0) const; + ex eval_ncmul(const exvector & v) const; unsigned calchash() const; ex series(const relational & r, int order, unsigned options = 0) const; ex thiscontainer(const exvector & v) const; @@ -1053,6 +1054,16 @@ ${evalf_switch_statement} throw(std::logic_error("function::evalf(): invalid nparams")); } +/* +* This method is defined to be in line with behaviour of function::return_type() +*/ +ex function::eval_ncmul(const exvector & v) const +{ + // If this function is called then the list of arguments is non-empty + // and the first argument is non-commutative, see function::return_type() + return seq.begin()->eval_ncmul(v); +} + unsigned function::calchash() const { unsigned v = golden_ratio_hash(golden_ratio_hash((p_int)tinfo()) ^ serial);

1 0

gcd
by Jens Vollinga 24 Apr '08

24 Apr '08

Hi, there must be some bug in gcd_heur(): "gcd(5/4+n, 25/8+n^2+15/4*n)" causes an infinite loop. Any ideas? Regards, Jens

2 2

GiNaC release 1.4.3
by Jens Vollinga 04 Apr '08

04 Apr '08

Hello again everybody, in order to fulfill the increased demand for GiNaC releases we proudly present: GiNaC 1.4.3 The changes are: * Fixed bug in numerical evaluation of multiple polylogarithms and alternating multiple zeta values introduced in version 1.4.2. * Nielsen polylog now invalidates its lookup tables in case the precision (Digits) has been changed. * Added new checks for recent bugs in the numerical evaluation of Li and zeta. As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

1 0

[SCM] GiNaC -- a C++ library for symbolic computations tag, release_1-4-3, created. release_1-4-0-49-g314595d
by git＠ginac.de 04 Apr '08

04 Apr '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The tag, release_1-4-3 has been created at 314595d10b4a4471bd4637a9038b18dde4780513 (commit) - Log ----------------------------------------------------------------- commit 314595d10b4a4471bd4637a9038b18dde4780513 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Fri Apr 4 15:00:01 2008 +0200 Preparing for release 1.4.3. ----------------------------------------------------------------------- hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, ginac_1-4, updated. release_1-4-0-49-g314595d
by git＠ginac.de 04 Apr '08

04 Apr '08

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, ginac_1-4 has been updated via 314595d10b4a4471bd4637a9038b18dde4780513 (commit) from 37cce07bdc6d2508f026e66190f6a982dce2e2b0 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 314595d10b4a4471bd4637a9038b18dde4780513 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Fri Apr 4 15:00:01 2008 +0200 Preparing for release 1.4.3. ----------------------------------------------------------------------- Summary of changes: NEWS | 7 +++++++ configure.ac | 6 +++--- 2 files changed, 10 insertions(+), 3 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0