- GiNaC-list - ginac.de

g++-3.3 and template specializations
by Kai Ludwig 08 Jun '04

08 Jun '04

Hi, there are some error messages, g++-3.3 gives me concerning missing template<> statements in GiNaC-1.2.1 Kai /usr/local/include/ginac/container.h:665: error: explicit specialization not preceded by `template <>' In file included from /usr/local/include/ginac/normal.h:29, from /usr/local/include/ginac/ginac.h:31, from ../genericRateLaw.hpp:32, from ../ecco.hpp:41, from eccoTest.cpp:1: /usr/local/include/ginac/lst.h:35: error: explicit specialization not preceded by `template <>' /usr/local/include/ginac/lst.h:38: error: explicit specialization not preceded by `template <>' /usr/local/include/ginac/lst.h:41: error: explicit specialization not preceded by `template <>' /usr/local/include/ginac/lst.h:44: error: explicit specialization not preceded by `template <>' In file included from /usr/local/include/ginac/ncmul.h:26, from /usr/local/include/ginac/structure.h:29, from /usr/local/include/ginac/ginac.h:42, from ../genericRateLaw.hpp:32, from ../ecco.hpp:41, from eccoTest.cpp:1: /usr/local/include/ginac/exprseq.h:35: error: explicit specialization not preceded by `template <>' -- http://echempp.sourceforge.net Kai Ludwig Institut für Organische Chemie Auf der Morgenstelle 18 72076 Tübingen Tel.: 07071/29-73049 Mail: kai.ludwig(a)uni-tuebingen.de

1 1

Re: Error while building GiNaC-1.2.1
by Kai Ludwig 08 Jun '04

08 Jun '04

Hi, seems that bison/flex (yacc/lex) have to be installed on the system, before GiNaC can be compiled. Nothing is said about that in the README/INSTALL files. The configure script checks it, but does not fail. Anyway, after providing them, g++-3.3 compiles GiNaC-1.2.1 fine. Cheers, Kai -------- Original Message -------- Subject: Error while building GiNaC-1.2.1 From: "Kai Ludwig" <kai.ludwig(a)uni-tuebingen.de> Date: Tue, June 8, 2004 9:42 am To: <ginac-devel(a)ginac.de> Hi All, First of all: You're doing a great job with GiNaC! I've tried to compile the GiNaC-1.2.1 version with g++-3.3.3-41 under a SuSE9.1 Linux system. Got the following: g++ -DHAVE_CONFIG_H -I. -I. -I.. -g -O2 -MT input_lexer.lo -MD -MP -MF .deps/input_lexer.Tpo -c input_lexer.cc -fPIC -DPIC -o .libs/input_lexer.oIn file included from /home/cbauer/projects/GiNaC-1.2/ginac/input_lexer.ll:35:input_lexer.h:43: error: conflicting types for `char ginac_yytext[]' input_lexer.cc:446: error: previous declaration as `char*ginac_yytext' input_lexer.h:43: error: storage size of `ginac_yytext' isn't known input_lexer.h:43: error: storage size of `ginac_yytext' isn't known make[2]: *** [input_lexer.lo] Error 1 make[2]: Leaving directory `/usr/local/src/ginac/GiNaC-1.2.1/ginac' make[1]: *** [all-recursive] Error 1 make[1]: Leaving directory `/usr/local/src/ginac/GiNaC-1.2.1' make: *** [all] Error 2 Cheers, Kai -- http://echempp.sourceforge.net Kai Ludwig Institut für Organische Chemie Auf der Morgenstelle 18 72076 Tübingen Tel.: 07071/29-73049 Mail: kai.ludwig(a)uni-tuebingen.de -- http://echempp.sourceforge.net Kai Ludwig Institut für Organische Chemie Auf der Morgenstelle 18 72076 Tübingen Tel.: 07071/29-73049 Mail: kai.ludwig(a)uni-tuebingen.de

1 0

[Fwd: error in cln check]
by Kai Ludwig 08 Jun '04

08 Jun '04

Hi All, I posted the message below to Richard last week. I recognized on his homepage that he could be a bit busy. Anyway, maybe it is also relevant for GiNaC users. Therefore, I decided to post it again to the ginac lists. Meanwhile, I've installed cln-1.1.5 from a rpm-package located at "http://mirror.mcs.anl.gov/suse/i386/9.1/suse/i586/" That should work. Kai -------- Original Message -------- Subject: error in cln check From: "Kai Ludwig" <kai.ludwig(a)uni-tuebingen.de> Date: Fri, June 4, 2004 4:17 pm To: <richard.kreckel(a)ginac.de> Hello Richard, I've just compiled the cln-1.1.7 with the g++-3.3.3-41 on a SuSE Linux 9.1 platform and got the following linker error: g++ -g -O2 -I../include -I../include -I. -I../src -I../src/base -I../src/float -I../src/float/lfloat -I../src/integer -I../src/base/random -I../src/base/digitseq -I../src/base/digit -I../src/base -I../src/base -S exam_I_sqrtp.cc -o exam_I_sqrtp.sgcc -g -O2 -I../include -I../include -I. -I../src -I../src/base -I../src/float -I../src/float/lfloat -I../src/integer -I../src/base/random -I../src/base/digitseq -I../src/base/digit -I../src/base -I../src/base -c exam_I_sqrtp.s -o exam_I_sqrtp.o/bin/sh ../libtool --mode=link g++ -g -O2 exam.o exam_I.o exam_RA.o exam_SF.o exam_FF.o exam_DF.o exam_LF.o exam_I_gcd.o exam_I_sqrtp.o .../src/libcln.la -lm -o exammkdir .libs g++ -g -O2 exam.o exam_I.o exam_RA.o exam_SF.o exam_FF.o exam_DF.o exam_LF.o exam_I_gcd.o exam_I_sqrtp.o -o .libs/exam .../src/.libs/libcln.so /usr/lib/libgmp.so -lm -Wl,--rpath -Wl,/usr/local/libexam.o(.text+0x42): In function `_GLOBAL__I__Z15test_elementaryv': .../include/cln/float.h:649: undefined reference to `_GLOBAL__I_cl_module__cl_ieee__firstglobalfun'exam_SF.o(.text+0x41): In function `_GLOBAL__I__Z11test_sfloatv': .../include/cln/float.h:649: undefined reference to `_GLOBAL__I_cl_module__cl_ieee__firstglobalfun'exam_FF.o(.text+0x41): In function `_GLOBAL__I__Z11test_ffloatv': .../include/cln/float.h:649: undefined reference to `_GLOBAL__I_cl_module__cl_ieee__firstglobalfun'exam_FF.o(.text+0x46): In function `_GLOBAL__I__Z11test_ffloatv': .../include/cln/ffloat.h:289: undefined reference to `_GLOBAL__I_cl_module__cl_ieee__firstglobalfun'exam_DF.o(.text+0x41): In function `_GLOBAL__I__Z11test_dfloatv': .../include/cln/float.h:649: undefined reference to `_GLOBAL__I_cl_module__cl_ieee__firstglobalfun'exam_DF.o(.text+0x46):../include/cln/dfloat.h:290:more undefined references to `_GLOBAL__I_cl_module__cl_ieee__firstglobalfun' follow../src/.libs/libcln.so: undefined reference to `_GLOBAL__I_cl_module__cl_SV_ringelt__firstglobalfun'../src/.libs/libcln.so: undefined reference to `_GLOBAL__I_cl_module__cl_GV_number__firstglobalfun'../src/.libs/libcln.so: undefined reference to `_GLOBAL__I_cl_module__cl_GV_I__firstglobalfun'../src/.libs/libcln.so: undefined reference to `_GLOBAL__I_cl_module__cl_SV_number__firstglobalfun'../src/.libs/libcln.so: undefined reference to `_GLOBAL__I_cl_module__cl_symbol__firstglobalfun' Cheers Kai -- http://echempp.sourceforge.net Kai Ludwig Institut für Organische Chemie Auf der Morgenstelle 18 72076 Tübingen Tel.: 07071/29-73049 Mail: kai.ludwig(a)uni-tuebingen.de

1 0

power rules
by econrad 22 May '04

22 May '04

Here's an annoying issue I have: I'm using GiNaC as a symbolic front end for a dynamical systems tool I'm building that does alot of numerical bifurcation analysis and such (I'm using GiNaC to generate all the symbolic info I need up front, translating this into numerical prodedures, and then using LAPACK/BLAS and some new continuation algorithms on the what I've generated). Models typically have terms like "a*k^n/(p^n+k^n)" in the RHS of the ODEs, where a, k, p are real numbers and n is typically an integer. Calculation of the jacobian then has terms like -a*k^n*n*p^(n-1)/(p^n+k^n)^2. Problem is, GiNaC uses the safe assumption and calculates -a*k^n*n*p^n*(p^n+k^n)^(-2)*p^(-1), which introduces a singularity into the derivatives that isn't really there. I understand the rational behind not automatically making p^n/p^1=p^(n-1) when one doesn't know anything about n, but nevertheless, I know that I'll only be using a real n>1 (or in any case, I'll only use n that makes this a valid reduction). How can I force this power transformation to take place... I'm sure I can do some convoluted matching and replacing to get this to work, but a user should be able to give an option to "eval" or to "normal" to get this operation performed at their own risk, yes? Have I missed something? Any help would be appreciated. Emery Emery Conrad 4089 Derring Hall Virginia Tech, Blacksburg VA 540.231.5508

2 1

GiNaC 1.2.1 now available
by Christian Bauer 23 Apr '04

23 Apr '04

Hello! GiNaC 1.2.1 is now available! The changes are: - Fixed infinite recursion in atan2_evalf() and improved atan2_eval(). - Added automatic evaluations for trigonometric functions with negative arguments (e.g. sin(-2) -> -sin(2)). - Fixed a static initialization order goof-up. - Fixed various bugs in series expansion. If you had problems using GiNaC 1.2 on certain platforms because you were getting segmentation faults all the way, this release should solve the issue. You can download GiNaC 1.2.1 as always from ftp://ftpthep.physik.uni-mainz.de/pub/GiNaC Bye, Jens and Christian -- / Physics is an algorithm \/ http://www.uni-mainz.de/~bauec002/

1 0

Extracting monomials from a polynomial
by Anish Muttreja 11 Apr '04

11 Apr '04

Hi I am using ginac in a symbolic regression program. It generates different polynomials using genetic programming, mulitplies each monomial by a weight and then fits the resulting expression using regression. I need to extract each monomial from the polynomial. Until now , I have been getting by only dealing with the case where the top level class is "add", and the program worked well enough for me to publish a paper on the use of symbolic regression in software modelling. Yes , I wanted to let the authors know, they have another citation :-). But this restricts the program. Is there a better way to do this ? Will ex::expand () convert the expression to a polynomial I should mention that the ginac expression itself is generated from a genetic evolved tree, so that I can generate all the monomials while doing this conversion, but that will leave me with the task of collecting like terms, something I would rather have done for me :-). Thanks Anish -- Anish Muttreja ~ Graduate Student Electrical Engineering ~ Princeton University ======================================================== As an adolescent I aspired to lasting fame, I craved factual certainty, and I thirsted for a meaningful vision of human life - so I became a scientist. This is like becoming an archbishop so you can meet girls. -Matt Cartmill, anthropology professor and author (1943- )

2 2

Vectors
by Olivier Parisy 09 Apr '04

09 Apr '04

Hi, I'd be interested in representing cartesian vectors in GiNaC. I suppose I could represent one using a matrix, but having a specific node for it and, as an example, for cross product or vectors sum would allow for higher-level expressions representation. I didn't see any hint for this in the current GiNaC distribution or in its todo list, but perhaps anyone already did it ? Any advice ? Regards, Olivier Parisy.

1 0

series expansion once again...
by varg＠thsun1.jinr.ru 05 Apr '04

05 Apr '04

Hello! > Series expansion now works predictably. All terms with the exponent > of the expansion variable smaller than the given order are > calculated exactly. If the series is not terminating, the Order > function is (at least) of the given order. ( GiNaC 1.1.7 changelog ) That sounds nice, but let's consider a simple program. #include <iostream> #include <stdexcept> using namespace std; #include <ginac/ginac.h> using namespace GiNaC; int main(int argc, char** argv) { try { symbol eps("epsilon"); ex part1 = tgamma(2-eps)*pow(eps, -1)*pow(tgamma(3-eps), -1); ex part2 = tgamma(2-eps)*pow(eps, 2)*pow(tgamma(3-eps), -1); ex test = part1 + part2; cout << test << " = " << test.series(eps, 1) << endl; cout << part1 << " = " << part1.series(eps, 1) << endl; cout << part2 << " = " << part2.series(eps, 1) << endl; return 0; } catch ( exception & ee) { cerr << ee.what() << endl; return 1; } } My naive expectations are: 1) first and second lines should be the same 2) third line should be Order(epsilon) But the output is: tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^2+tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^(-1) = Order(epsilon^(-1)) tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^(-1) = 1/2*epsilon^(-1)+1/4+Order(epsilon) tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^2 = Order(epsilon^(-1)) Is this a bug or a feature? Thanks, Alexei. P.S. $ ginac-config --version 1.1.7 $ g++ --version | head -n 1 g++ (GCC) 3.3.3 (Debian 20040401)

3 2

How to calculate the area of a circle?
by Tom Brinkman 29 Mar '04

29 Mar '04

How do i solve for the area of a circle? Example: symbol r("r"); ex e = pow(r,2) * Pi == 30; ex solve = lsolve(e,r); //This example does not work. An error message indicating that this is not a linear equation is returned. How do I solve for "r" in this case. //Could you also show an example of how to solve for a simple polynomial, such as: ex = pow(a,2) + a = 3 How do i solve for "a". regards, Tom Brikman __________________________________ Do you Yahoo!? Yahoo! Finance Tax Center - File online. File on time. http://taxes.yahoo.com/filing.html

2 1

troubles with ginac 1.2 on debian unstable
by Andrius Kurtinaitis 25 Mar '04

25 Mar '04

Hello, i was so happy about 1.2 release, that I downloaded the new 1.2 version of ginac, compiled and installed on my debian unstable machine, using usual configure, make, su, make install. I also have another ginac version 1.1.6-1 installed. When I build a simple proggy using ginac, I use new (1.2) ginac lib and the old cln: g++ -o charpoly charpoly.cpp /usr/local/lib/libginac.a -lcln The resulting program crashes with the segmentation fault. Even if it does nothing. Including <ginac/ginac.h> and compiling the empty program is enough to produce the segfault. The current installation contains the cln version 1.1.6-1. Any hints? I thought, it is because of version incompatibility between the cln library and the new ginac. Then i tried to download the newest cln, which, according to its site should be 1.1.7. But the download site ftp://ftpthep.physik.uni-mainz.de/pub/gnu/ does not contain this version... Or is it because of different compiler versions used to compile cln and my own ginac 1.2? Any ideas? Andrius Kurtinaitis

4 8