- GiNaC-list - ginac.de

Release of GiNaC 1.3.6
by Jens Vollinga 13 Dec '06

13 Dec '06

Hello everybody, GiNaC 1.3.6 is out and available. The changes are: - Better worst case behavior in heuristic gcd. - Improved collecting for sparse multivariate polynomials. - Improved collect_common_factors. - Fixed bug in numerical evaluation of multiple polylogs. - Fixed numerical integration of complex functions. - Fixed bug in tensor::replace_contr_index(). As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

1 0

Internal error: statement in file ./real/conv/cl_R_to_double.cc, line 55 has been reached!!
by Pedro Cruz 08 Dec '06

08 Dec '06

Hi, -------------------- this output -------------------- Pi=3.1415926535897932385 Internal error: statement in file ./real/conv/cl_R_to_double.cc, line 55 has been reached!! Please send the authors of the program a description how you produced this error! -------------------- has been produced by -------------------- #include <iostream> #include <sstream> #include <stdexcept> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { cout << "Pi=" << Pi.evalf() << endl; if (is_a<constant>(Pi)) { double f = ex_to<numeric>(Pi).to_double(); cout << f << endl; } return 0; } --------------------------------------------- on debian (etch testing) with libginac1.3c2a. --------------------------------------------- It was a first example i tried to do in order to understand ginac/cln concepts: <ex> <constant> and <numeric>. How to convert from Pi to C++ double type in one C++ line ? Thanks, Pedro Cruz

2 1

GiNaC as a polynomials manipulator?
by Francesco Biscani 03 Dec '06

03 Dec '06

Hello, I'd like to know if GiNaC can be used as a polynomials manipulator for celestial mechanics. It's needed to multiplicate and add polynomials by polynomials. An important requisite is to be able to extract the degree of the monomials in order to decide whether that particular monomial should be eliminated from the polynomial or not. Example: x^2*y+x^2*y^2+x^3*z If the truncation threshold is set to order 4 this polynomial gets truncated to x^2*y I imagine that a function to place polynomials in some kind of "canonical" form is needed, and then it is needed to inspect each monomial to see if it will be kept or not. Is this possible with GiNaC? Thanks and best regards, Francesco

2 1

Fwd: compute groebner basis?
by fabian＠fseoane.net 03 Dec '06

03 Dec '06

Hello! I'm new to GiNaC and my question is: Has GiNaC got a routine to compute Groebner basis? If the answer is no, then if i tried to implement one, would it be welcome ? Greetings ------------- Fabian

2 1

Problem with scalar_products : Scalar product of two vectors in index notation
by aravind b 09 Nov '06

09 Nov '06

Hi, Hope everyone's doing fine ! Please go through this code : #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { symbol i_sym("i"), j_sym("j"), e("e"), a("a"), b("b") ; idx i(i_sym, 3), j( j_sym, 3) ; //defining indices i and j scalar_products sp; sp.add( indexed(e, i), indexed(e, j), delta_tensor(i,j) ); // e.i* e.j = delta_tensor.i.j cout << (indexed(e,i)*indexed(e,j)).simplify_indexed(sp) <<endl ; return 0; } does NOT give delta.j.i instead just gives e.i * e.j - why? Basically, i want to define an orthonormal basis vector set {e.i}. If i have vector_a = indexed(a, i) * indexed(e, i) and vector_b = indexed(b, j) * indexed(e, j) , then i want adotb to reduce to a.i * b.i ( ofcoure after i apply simplify_indexed(sp) to adotb ) Thanks a lot for going through this mail.

2 3

Problem with linking ginac
by aravind b 08 Nov '06

08 Nov '06

Hello, I am using quantian *0.7.9.2 *live dvd and i am a total newbie to linux. I am really unable to use ginac. Firstly, locate ginac-config doesnt give any results. whereas locate ginac gives: *locate: warning: database `/var/cache/locate/locatedb' is more than 8 days old /usr/lib/libginac-1.3.so.2 /usr/lib/libginac-1.3.so.2.1.0 /usr/share/doc/ginac-tools /usr/share/doc/ginac-tools/AUTHORS /usr/share/doc/ginac-tools/changelog.Debian.gz /usr/share/doc/ginac-tools/changelog.gz /usr/share/doc/ginac-tools/copyright /usr/share/doc/ginac-tools/NEWS.gz /usr/share/doc/ginac-tools/README /usr/share/doc/libginac1.3c2a /usr/share/doc/libginac1.3c2a/AUTHORS /usr/share/doc/libginac1.3c2a/changelog.Debian.gz /usr/share/doc/libginac1.3c2a/changelog.gz /usr/share/doc/libginac1.3c2a/copyright /usr/share/doc/libginac1.3c2a/NEWS.gz /usr/share/doc/libginac1.3c2a/README /usr/share/menu/ginac-tools /usr/share/pixmaps/ginac.xpm /var/lib/dpkg/info/ginac-tools.list /var/lib/dpkg/info/ginac-tools.md5sums /var/lib/dpkg/info/ginac-tools.postinst /var/lib/dpkg/info/ginac-tools.postrm /var/lib/dpkg/info/libginac1.3c2.list /var/lib/dpkg/info/libginac1.3c2.postrm /var/lib/dpkg/info/libginac1.3c2a.list /var/lib/dpkg/info/libginac1.3c2a.md5sums /var/lib/dpkg/info/libginac1.3c2a.postinst /var/lib/dpkg/info/libginac1.3c2a.postrm /var/lib/dpkg/info/libginac1.3c2a.shlibs *Guys, please help me out here. How do i link ginac to my C++ code using g++? Thanks a lot for going through this mail with regards, aravind * *

3 6

Problem with linking ginac
by aravind b 05 Nov '06

05 Nov '06

Hello, I am using quantian *0.7.9.2 *live dvd and i am a total newbie to linux and to the mailing lists. I am really unable to use ginac. Firstly, locate ginac-config doesnt give any results. whereas locate ginac gives: *locate: warning: database `/var/cache/locate/locatedb' is more than 8 days old /usr/lib/libginac-1.3.so.2 /usr/lib/libginac-1.3.so.2.1.0 /usr/share/doc/ginac-tools /usr/share/doc/ginac-tools/AUTHORS /usr/share/doc/ginac-tools/changelog.Debian.gz /usr/share/doc/ginac-tools/changelog.gz /usr/share/doc/ginac-tools/copyright /usr/share/doc/ginac-tools/NEWS.gz /usr/share/doc/ginac-tools/README /usr/share/doc/libginac1.3c2a /usr/share/doc/libginac1.3c2a/AUTHORS /usr/share/doc/libginac1.3c2a/changelog.Debian.gz /usr/share/doc/libginac1.3c2a/changelog.gz /usr/share/doc/libginac1.3c2a/copyright /usr/share/doc/libginac1.3c2a/NEWS.gz /usr/share/doc/libginac1.3c2a/README /usr/share/menu/ginac-tools /usr/share/pixmaps/ginac.xpm /var/lib/dpkg/info/ginac-tools.list /var/lib/dpkg/info/ginac-tools.md5sums /var/lib/dpkg/info/ginac-tools.postinst /var/lib/dpkg/info/ginac-tools.postrm /var/lib/dpkg/info/libginac1.3c2.list /var/lib/dpkg/info/libginac1.3c2.postrm /var/lib/dpkg/info/libginac1.3c2a.list /var/lib/dpkg/info/libginac1.3c2a.md5sums /var/lib/dpkg/info/libginac1.3c2a.postinst /var/lib/dpkg/info/libginac1.3c2a.postrm /var/lib/dpkg/info/libginac1.3c2a.shlibs *Guys, please help me out here. How do i link ginac to my C++ code using g++? Thanks a lot for going through this mail with regards, aravind

1 0

amd64 - test failed
by Peper 22 Oct '06

22 Oct '06

Hello, test failed on amd64, any idea? GiNaC will now take an exam with specific input (like a pupils' exam): examining several historic failures just out of paranoia................ passed examining consistency of numeric types.....Error: caught exception syntax error, unexpected ')', expecting ']' or ',' at ) examining power laws..... passed examining consistency of symbolic functions.... passed examining consistency of nestedsums functions..................................Error: caught exception syntax error, unexpected T_EOF, expecting ']' or ',' at examining symbolic differentiationError: caught exception syntax error, unexpected T_EOF, expecting ']' or ',' at examining polynomial GCD computation......... passed examining rational function normalization..... passed examining series expansion............ passed examining symbolic matrix manipulations........ passed examining linear solve...... passed examining indexed objects........ passed examining color objects... passed examining clifford objects.......... passed examining archiving system. passed examining structure template.. passed examining hash maps............ passed examining miscellaneous other things........ passed Error: something went wrong. (3 individual failures) please check exams.out against exams.ref for more details. happy debugging! ./exams.ref exams.out differ: byte 90, line 4 FAIL: run_exams GiNaC will now run through some rather costly random consistency checks: checking consistency of numeric types.. passed checking consistency of symbolic functions.... passed checking symbolic matrix manipulations..... passed checking linear solve.......... passed PASS: run_checks GiNaC will now run through some basic timings: timing commutative expansion and substitution.... passed size: 50 100 200 400 time/s: 0.03 0.13 0.54 2.42 timing Laurent series expansion of Gamma function.... passed order: 10 15 20 25 time/s: 0.05 0.09 0.38 1.47 timing determinant of univariate symbolic Vandermonde matrices.... passed dim: 6x6 8x8 10x10 12x12 time/s: 0.006 0.05 0.5 4.61 timing determinant of polyvariate symbolic Toeplitz matricesError: caught exception syntax error, unexpected T_EOF, expecting ']' or ',' at timing hash map operations.... passed size: 10000 50000 100000 500000 insert/s: 0.02 0.09 0.19 0.86 find/s: 0.01 0.05 0.1 0.5 erase/s: 0.01 0.04 0.08 0.43 timing Lewis-Wester test A (divide factorials). passed 0.05s timing Lewis-Wester test B (sum of rational numbers). passed 0.003s timing Lewis-Wester test C (gcd of big integers). passed 0.024s timing Lewis-Wester test D (normalized sum of rational fcns). passed 0.08s timing Lewis-Wester test E (normalized sum of rational fcns). passed 0.04s timing Lewis-Wester test F (gcd of 2-var polys). passed 0.012s timing Lewis-Wester test G (gcd of 3-var polys). passed 0.38s timing Lewis-Wester test H (det of 80x80 Hilbert). passed 1.24s timing Lewis-Wester test I (invert rank 40 Hilbert). passed 0.37s timing Lewis-Wester test J (check rank 40 Hilbert). passed 0.22s timing Lewis-Wester test K (invert rank 70 Hilbert). passed 2.4s timing Lewis-Wester test L (check rank 70 Hilbert). passed 1.28s timing Lewis-Wester test M1 (26x26 sparse, det)Error: caught exception parser stack overflow at ^ timing Lewis-Wester test M2 (101x101 sparse, det) disabled timing Lewis-Wester test N (poly at rational fcns) disabled timing Lewis-Wester test O1 (three 15x15 dets)... passed 12.486s (average) timing Lewis-Wester test O2 (Resultant) disabled timing Lewis-Wester test P (det of sparse rank 101). passed 0.56s timing Lewis-Wester test P' (det of less sparse rank 101). passed 0.84s timing Lewis-Wester test Q (charpoly(P)). passed 19.01s timing Lewis-Wester test Q' (charpoly(P')). passed 37.2s timing computation of antipodes in Yukawa theory...... passed 21.85s (total) timing Fateman's polynomial expand benchmark. passed 47.85s Error: something went wrong. (2 individual failures) please check times.out against times.ref for more details. happy debugging! ./times.ref times.out differ: byte 261, line 8 FAIL: run_times ====================================== 2 of 3 tests failed Please report to <ginac-list(a)ginac.de> ====================================== make[3]: *** [check-TESTS] Error 1 make[3]: Leaving directory `/var/tmp/portage/sci-mathematics/ginac-1.3.4/work/ginac-1.3.4/check' make[2]: *** [check-am] Error 2 make[2]: Leaving directory `/var/tmp/portage/sci-mathematics/ginac-1.3.4/work/ginac-1.3.4/check' make[1]: *** [check-recursive] Error 1 make[1]: Leaving directory `/var/tmp/portage/sci-mathematics/ginac-1.3.4/work/ginac-1.3.4' make: *** [check] Error 2

2 1

Suggestion for the FAQ.
by Chris Dams 17 Oct '06

17 Oct '06

Dear all, I think the following should added to the FAQ. Q: I am getting an error message that looks more or less like Internal error: statement in file real/conv/cl_R_to_double.cc, line 55 has been reached!! Please send the authors of the program a description how you produced this error! What should I do? A: This is an error message of the CLN library that GiNaC uses to handle its numerics. The problem usually is that an object was passed to a CLN-function that should have been a number but is not. Your best bet is to check for code that looks like ex_to<numeric>(x) and to check whether x really is a number. The just given code fragment cannot be right if is_a<numeric>(x) yields false. Thanks! Chris

4 4

hash tables and sets
by Marko Riedel 12 Oct '06

12 Oct '06

Hi folks, I would like to write a set implementation for expressions, or borrow one, if it has already been written. I have several questions. 1. What C++ type is best suited to representing sets? 2. Is there a hash function implemented for GiNaC expressions? 3. Can GiNaC expressions be put into some kind of canonical form? 4. (Requirement.) The hash function should be such that e.g. a + b + c * (d + e) and b + (e + d) * c + a hash to the same value. Thanks! Marko Riedel +-------------------------------------------------------------+ | Marko Riedel, EDV Neue Arbeit gGmbH, mriedel(a)lsi.upc.edu | | http://www.geocities.com/markoriedelde/index.html | +-------------------------------------------------------------+

3 7