- GiNaC-list - ginac.de

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

Re: GiNaC error-message
by Chris Dams 11 Oct '06

11 Oct '06

Dear Martijn, On Tue, 10 Oct 2006, Martijn Frijlink wrote: > I am using GiNaC 1.3.3 and while running the executable of the attached > file blockrec.cc I got the following message : > > $ ./blockrec 11 > Discrete frequencies 0 ... 10 will be used > > 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! The problem is the line CosC[i] = ex_to<numeric>(Lsigma*a_m).to_double(); The expression Lsigma*a_m is not a numeric but e.g. something like 6*sqrt(3)*Pi^(-1). To evaluate this into a numeric you have to use evalf. Best wishes, Chris

1 0

memory allocation, matrices
by Marko Riedel 11 Sep '06

11 Sep '06

Hi folks, I'm going to ask a very basic question, namely if you can point me to a tutorial where GiNaC's memory managment is described. I read the tutorial from your website but I could not find any info on this. I am familiar with the retain/release mechanism of used by Cocoa and GNUstep, and the use of autorelease pools. If you could compare these concepts to their GiNaC equivalents, that would be of great help. E.g. in Objective C I can create a tree of objects and either free them myself or put them into an autorelease pool that will be emptied by the runloop. How would I do this in GiNaC? Suppose I create a matrix of rational functions, a single entry of which I'm interested in retaining. I would like to free the matrix and retain that one function. Can you provide a code snippet that does this? I would also like to know what the limitations of the linear equation solver are. How many variables can it handle? What are the time and space complexities of the solver in terms of the number of variables and the number of equations? Thanks in advance for any help you may be able to give. Best regards, Marko Riedel +-------------------------------------------------------------+ | Marko Riedel, EDV Neue Arbeit gGmbH, mriedel(a)lsi.upc.es | | http://www.geocities.com/markoriedelde/index.html | +-------------------------------------------------------------+

3 8

CXXABI_1.3.1
by Marko Riedel 07 Sep '06

07 Sep '06

Hi folks, I downloaded, compiled and installed GiNaC using the following version of gcc: gcc -v Using built-in specs. Target: i686-pc-linux-gnu Configured with: ../configure Thread model: posix gcc version 4.0.2 I now appear to have a working setup. E.g. I can run "ginsh": ginsh - GiNaC Interactive Shell (ginac V1.3.5) __, _______ Copyright (C) 1999-2006 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. <-------------' For details type `warranty;'. When I compile and run my test program, however, I get the following error message: /usr/lib/libstdc++.so.6: version `CXXABI_1.3.1' not found \ (required by /usr/local/lib/libginac-1.3.so.2) Am I perhaps missing some environment variables or compiler switches here? I have the following files: /usr/include/c++/3.3/cxxabi.h /usr/local/include/c++/4.0.2/i686-pc-linux-gnu/bits/cxxabi_tweaks.h /usr/local/include/c++/4.0.2/cxxabi.h I compiled my program with the following command: c++ -g complete41.cpp -o complete41 -I/usr/local/include -lcln -lginac Any help would be appreciated. Best regards, Marko Riedel -- +------------------------------------------------------------+ | Marko Riedel, EDV Neue Arbeit gGmbH, mriedel(a)neuearbeit.de | | http://www.geocities.com/markoriedelde/index.html | +------------------------------------------------------------+

2 2

compiling cln 1.1.13
by Marko Riedel 07 Sep '06

07 Sep '06

Hi folks, I am trying to compile cln on Debian 3.1 (Kernel 2.6.12). This is my version of gcc: gcc -v Using built-in specs. Target: i686-pc-linux-gnu Configured with: ../configure Thread model: posix gcc version 4.0.2 This is the output from "make". I get the impression that I have to set an environment variable to specify endianess. cd src && make all make[1]: Entering directory `/usr/local/src/cln-1.1.13/src' make SUBDIR=base make[2]: Entering directory `/usr/local/src/cln-1.1.13/src' /bin/sh ../libtool --mode=compile g++ -O -I../include -I../include \ -I./base -c ./base/cl_N_err_d0.cc mkdir .libs g++ -O -I../include -I../include -I./base -c ./base/cl_N_err_d0.cc \ -fPIC -DPIC -o .libs/cl_N_err_d0.o In file included from ./base/cl_N_err_d0.cc:4: ./base/cl_sysdep.h:30:4: error: #error "Bogus CL_CPU_BIG_ENDIAN_P!" In file included from ./base/cl_macros.h:6, from ./base/cl_sysdep.h:56, from ./base/cl_N_err_d0.cc:4: ../include/cln/types.h:18:6: error: #error "No 8 bit integer type?" ../include/cln/types.h:25:6: error: #error "No 16 bit integer type?" ../include/cln/types.h:35:6: error: #error "No 32 bit integer type?" ../include/cln/types.h:57: error: 'sint8' does not name a type ../include/cln/types.h:58: error: 'uint8' does not name a type ../include/cln/types.h:60: error: 'sint16' does not name a type ../include/cln/types.h:61: error: 'uint16' does not name a type ../include/cln/types.h:63: error: 'sint32' does not name a type ../include/cln/types.h:64: error: 'uint32' does not name a type ../include/cln/types.h:128: error: 'sint32' does not name a type ../include/cln/types.h:129: error: 'uint32' does not name a type ../include/cln/number.h:194: error: expected constructor, destructor, \ or type conversion before 'CL_DEFINE_LONG_ASSIGNMENT_OPERATORS' make[2]: *** [cl_N_err_d0.lo] Error 1 make[2]: Leaving directory `/usr/local/src/cln-1.1.13/src' make[1]: *** [base.target] Error 2 make[1]: Leaving directory `/usr/local/src/cln-1.1.13/src' make: *** [all] Error 2 Can you help me out here? Thanks! -- +------------------------------------------------------------+ | Marko Riedel, EDV Neue Arbeit gGmbH, mriedel(a)neuearbeit.de | | http://www.geocities.com/markoriedelde/index.html | +------------------------------------------------------------+

2 2

libstdc++5 vs 6
by Marko Riedel 06 Sep '06

06 Sep '06

Hi folks, I installed the GiNaC Debian package on my machine. I now have ||/ Name Version Description +++-==============-==========-============== ii libginac-dev 1.3.0-2 The GiNaC framework (development files) I get the following warning when I link my test program: /usr/bin/ld: warning: libstdc++.so.5, \ needed by /usr/bin/../lib/libginac.so, \ may conflict with libstdc++.so.6 The program segfaults in GiNaC::matrix::operator= (). How do I tell the linker to link against libstdc++.so.5 only, thereby avoiding conflicts? This is my version of GCC: gcc -v Using built-in specs. Target: i686-pc-linux-gnu Configured with: ../configure Thread model: posix gcc version 4.0.2 Thanks! -- +------------------------------------------------------------+ | Marko Riedel, EDV Neue Arbeit gGmbH, mriedel(a)neuearbeit.de | | http://www.geocities.com/markoriedelde/index.html | +------------------------------------------------------------+

2 1