- GiNaC-devel - ginac.de

infinite loop in simplify_indexed()
by PG CLARK 21 May '11

21 May '11

I've finally tracked down the source of my troubles in getting the little program I've been working on to run and I'm 99% convinced it's a bug in the library. Please take a look at this brief sample program I wrote that illustrates the problem. #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { symbol i_sym("i"), j_sym("j"), k_sym("k"), l_sym("l"), m_sym("m"), n_sym("n"); idx i(i_sym, 2), j(j_sym, 2), k(k_sym, 2), l(l_sym, 2), m(m_sym, 2), n(n_sym, 2); symbol alpha_sym("\\alpha"), beta_sym("\\beta"), gamma_sym("\\gamma"); idx Va(alpha_sym, 3), Vb(beta_sym, 3), Vg(gamma_sym, 3); symbol U("U"), v("v"), h("h"); exvector IJKLMNVec; idx IJKLMNIA[] ={i, j, k, l, m, n}; IJKLMNVec.assign (IJKLMNIA,IJKLMNIA+6); ex myEx = (indexed(h, sy_symm(sy_none(0, 1), sy_none(2, 3),sy_none(4, 5)), IJKLMNVec)*indexed(U, i, Va)*indexed(v, j, Va)*indexed(U, k, Vb) *indexed(v, l, Vb)*indexed(U, m, Vg)*indexed(v, n, Vg)).simplify_indexed(); cout << myEx << endl; return 0; } You'll find it never terminates. simplify_indexed() has an infinite loop that was unanticipated caused by some rare forms of indexed expressions. This program will just run and run trapped in implify_indexed() and never print anything out. If a fix is not forthcoming any ideas for a work around? regards peter clark ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

3 5

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-345-g56bbb99
by git＠ginac.de 21 May '11

21 May '11

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 56bbb99adc43dad22ee615380d4446a87b713cb8 (commit) from f320e27e9cfe287168c879af5991babffaa1e9c8 (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 56bbb99adc43dad22ee615380d4446a87b713cb8 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Sat May 21 01:30:08 2011 +0200 Limiting the costly symmetrization tests inside simplify_indexed() to cases where at least one antisymmetric or cyclic non-free index is involved. Thanks to P.G.Clark for reporting the problem. ----------------------------------------------------------------------- Summary of changes: ginac/indexed.cpp | 37 +++++++++++++++++++++++-------------- ginac/symmetry.cpp | 12 ++++++++++++ ginac/symmetry.h | 2 ++ 3 files changed, 37 insertions(+), 14 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-344-gf320e27
by git＠ginac.de 20 May '11

20 May '11

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 f320e27e9cfe287168c879af5991babffaa1e9c8 (commit) from 0c2f0f4c6d4118e817c5b9a14a9d0f3ada82e9fe (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 f320e27e9cfe287168c879af5991babffaa1e9c8 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Fri May 20 23:30:38 2011 +0200 Changed naming convention for the library. Now, it is only libginac.so without the release version included (instead of libginac-1.5.so, for example). Since the previous commit broke the ABI, we are free to implement this long planned change now. ----------------------------------------------------------------------- Summary of changes: configure.ac | 6 ------ ginac/Makefile.am | 2 +- 2 files changed, 1 insertions(+), 7 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-343-g0c2f0f4
by git＠ginac.de 20 May '11

20 May '11

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 0c2f0f4c6d4118e817c5b9a14a9d0f3ada82e9fe (commit) from 293ff6f6fe89c3e5413c8da57feedf3bd2016cd7 (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 0c2f0f4c6d4118e817c5b9a14a9d0f3ada82e9fe Author: Richard Kreckel <kreckel(a)ginac.de> Date: Fri May 20 22:20:57 2011 +0200 Fixed a bug introduced in last commit. Before, sqrt(2)*x used to be a polynomial in x. Now, it is a polynomial in x again. ----------------------------------------------------------------------- Summary of changes: check/exam_paranoia.cpp | 14 ++++++++++---- ginac/add.cpp | 10 ++++++++++ ginac/add.h | 1 + ginac/expairseq.cpp | 22 ---------------------- ginac/expairseq.h | 1 - ginac/mul.cpp | 11 +++++++++++ ginac/mul.h | 1 + 7 files changed, 33 insertions(+), 27 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-342-g293ff6f
by git＠ginac.de 20 May '11

20 May '11

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 293ff6f6fe89c3e5413c8da57feedf3bd2016cd7 (commit) from c152ba30b2938b3e1fe6927574883eff08dae88e (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 293ff6f6fe89c3e5413c8da57feedf3bd2016cd7 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Fri May 20 13:55:35 2011 +0200 Fixed a bug in is_polynomial() that caused expressions like sqrt(x*x+1)*sqrt(x+1) to be wrongly identified as polynomials in x. Thanks to Florent Hivert for reporting. ----------------------------------------------------------------------- Summary of changes: check/exam_paranoia.cpp | 13 +++++++++++++ ginac/expairseq.cpp | 19 +++++++++++++++---- 2 files changed, 28 insertions(+), 4 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-341-gc152ba3
by git＠ginac.de 20 May '11

20 May '11

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 c152ba30b2938b3e1fe6927574883eff08dae88e (commit) via 46e10827e1cb1888ab5419b16d3381d29161daa5 (commit) from ab1aedf73a5d062bf50ed1e740f38a4503d7fead (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 c152ba30b2938b3e1fe6927574883eff08dae88e Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Fri May 20 11:54:27 2011 +0200 ./configure now gives more information/warnings about missing utilities. commit 46e10827e1cb1888ab5419b16d3381d29161daa5 Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Thu May 19 10:29:14 2011 +0200 Avoiding a spurious missing ginac_parser.h error when doing a parallel build (make -j..). ----------------------------------------------------------------------- Summary of changes: acinclude.m4 | 8 +++++--- configure.ac | 21 +++++++++++++++++---- ginsh/Makefile.am | 2 +- 3 files changed, 23 insertions(+), 8 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-339-gab1aedf
by git＠ginac.de 19 May '11

19 May '11

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 ab1aedf73a5d062bf50ed1e740f38a4503d7fead (commit) from f9afb6aca6a971650dff63b11ca8c2ef18506690 (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 ab1aedf73a5d062bf50ed1e740f38a4503d7fead Author: Jens Vollinga <jensv(a)nikhef.nl> Date: Thu May 19 09:44:11 2011 +0200 The patch improves the run time at the expense of using more RAM in some situations. Please note: it doesn't improve the actual algorithm (iteration over all permutations). Thanks to Alexei Sheplyakov. ----------------------------------------------------------------------- Summary of changes: ginac/symmetry.cpp | 7 +++++-- 1 files changed, 5 insertions(+), 2 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

expairseq::conjugate
by Burcin Erocal 15 May '11

15 May '11

Hi, the following problem (in Sage syntax) was reported to the Sage issue tracker [1]: ginsh - GiNaC Interactive Shell (ginac V1.5.8) <snip> > conjugate(I*sqrt(-3)); -I*sqrt(-3) [1] http://trac.sagemath.org/sage_trac/ticket/10964 The problem seems to be the expair::conjugate() method which is called from the conjugateepvector() function in expairseq.cpp. This function is called from the pseries class as well. So, moving the latter in the expairseq class to use the virtual functions recombine_pair_to_ex() etc. does not work. Any ideas for a clean solution to this? Thank you. Burcin

3 5

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-338-gf9afb6a
by git＠ginac.de 15 May '11

15 May '11

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 f9afb6aca6a971650dff63b11ca8c2ef18506690 (commit) from e16348f2f6edb1747dca9e4c47300e0436eea617 (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 f9afb6aca6a971650dff63b11ca8c2ef18506690 Author: Richard Kreckel <kreckel(a)ginac.de> Date: Sun May 15 18:40:44 2011 +0200 Fix mul::conjugate(). Overwrite conjugate() for mul objects because the base class's conjugate function is incorrect at branch cuts when exponents are fractional (roots). This has been reported as Sage #10964. While at it, fix some funny indentation and remove unneeded tests. ----------------------------------------------------------------------- Summary of changes: check/exam_paranoia.cpp | 15 +++++++++++++++ ginac/expairseq.cpp | 4 +--- ginac/expairseq.h | 4 ++-- ginac/inifcns_trans.cpp | 2 +- ginac/mul.cpp | 32 ++++++++++++++++++++++++++++++++ ginac/mul.h | 1 + ginac/pseries.cpp | 6 ++---- 7 files changed, 54 insertions(+), 10 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

a clarification regarding the problem with simplify_indexed.
by PG CLARK 14 May '11

14 May '11

In fact the code doesn't even care about the tensor symmetry. This code also fails to terminate. #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { symbol i_sym("i"), j_sym("j"), k_sym("k"), l_sym("l"), m_sym("m"), n_sym("n"); idx i(i_sym, 2), j(j_sym, 2), k(k_sym, 2), l(l_sym, 2), m(m_sym, 2), n(n_sym, 2); symbol alpha_sym("\\alpha"), beta_sym("\\beta"), gamma_sym("\\gamma"); idx Va(alpha_sym, 3), Vb(beta_sym, 3), Vg(gamma_sym, 3); symbol U1("U1"), U2("U2"), U3("U3"), U4("U4"), U5("U5"), U6("U6"), h("h"); exvector IJKLMNVec; idx IJKLMNIA[] ={i, j, k, l, m, n}; IJKLMNVec.assign (IJKLMNIA,IJKLMNIA+6); ex myEx1 = (indexed(h, IJKLMNVec)*indexed(U1, i, Va)*indexed(U5, m, Vg) *indexed(U3, k, Vb)*indexed(U2, j, Va)*indexed(U4, l, Vb)*indexed(U6, n, Vg)); ex myEx2 = (indexed(h, IJKLMNVec)*indexed(U1, i, Va)*indexed(U1, m, Vg) *indexed(U1, k, Vb)*indexed(U1, j, Va)*indexed(U1, l, Vb)*indexed(U1, n, Vg)); ex myEx3 = (indexed(h, IJKLMNVec)*indexed(U1, i, Va)*indexed(U2, m, Vg) *indexed(U1, k, Vb)*indexed(U2, j, Va)*indexed(U1, l, Vb)*indexed(U2, n, Vg)); // swap comments to see .simplify_indexed() fails in all cases cout << myEx2.simplify_indexed() << endl; //cout << myEx3.simplify_indexed() << endl; //cout << myEx1.simplify_indexed() << endl; return 0; } Could this be an error in simplify_indexed_product? ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

1 0