- GiNaC-list - ginac.de

GiNaC Match Documentation
by Mark Round 15 Aug '13

15 Aug '13

I have been encountering problems with the GiNaC function match. Basically I can not understand the documentation. I wanted to detail my story in the hope that perhaps my problems are solved and in the future someone else can avoid my pitfalls. The documentation informs me that I need an associative array. I can not find how to construct such an array elsewhere in the documentation so I assume that I am supposed to use exhashmap. That was a point of confusion for me. Then I designed the following code snipet, symbol x("x"), y("y"); ex tmp=x+y; exhashmap <ex> array; tmp.match(wild(1)+y); which apparently works fine. Now I wish to fill 'array' with the replacement table. So I modified to, tmp.match(wild(1)+y,array); which returned the following compiler messages, momentum_integration.cpp: In function ‘GiNaC::ex bring_to_parameterisable_form(GiNaC::ex)’: momentum_integration.cpp:13:33: error: no matching function for call to ‘GiNaC::ex::match(const GiNaC::ex, GiNaC::exhashmap<GiNaC::ex>&)’ momentum_integration.cpp:13:33: note: candidates are: In file included from /usr/include/ginac/ginac.h:30:0, from momentum_integration.h:14, from momentum_integration.cpp:1: /usr/include/ginac/ex.h:150:7: note: bool GiNaC::ex::match(const GiNaC::ex&) const /usr/include/ginac/ex.h:150:7: note: candidate expects 1 argument, 2 provided /usr/include/ginac/ex.h:151:7: note: bool GiNaC::ex::match(const GiNaC::ex&, GiNaC::exmap&) const /usr/include/ginac/ex.h:151:7: note: no known conversion for argument 2 from ‘GiNaC::exhashmap<GiNaC::ex>’ to ‘GiNaC::exmap& {aka std::map<GiNaC::ex, GiNaC::ex, GiNaC::ex_is_less>&}’ make: *** [momentum_integration.o] Error 1 I attempted to use the command exmap but I can not find anything on that in the documentation. I tried several permutations around, exmap <ex,ex> array; which made the compiler return, momentum_integration.cpp: In function ‘GiNaC::ex bring_to_parameterisable_form(GiNaC::ex)’: momentum_integration.cpp:11:8: error: ‘GiNaC::exmap {aka std::map<GiNaC::ex, GiNaC::ex, GiNaC::ex_is_less>}’ is not a template make: *** [momentum_integration.o] Error 1 At which point it is unclear to me what else I should try. The code was ran on a 64-bit Linux Mint-Debian distribution with ginac-1.6.2.

2 1

ginac 1.6.2-2 on mingw-32.
by testhandle 07 Aug '13

07 Aug '13

(most of this text was written using google translate. english is not good at me. so there's a possibility that the meaning is not transmitted.) there may be unnecessary operation and defects. however, it is forced by this operation, but it have confirmed the operation of the tutorial examples. * step 0, dependencies. gmp 5.0.1 cln (gmp) 1.3.3 gcc 4.7.2 ld 2.23.1 m4 1.4.16 autoconf 2.68 autoreconf 2.68 automake 1.11 libtool 2.4 python 2.7.3 bison 2.4.2 flex 2.5.35 * step 1. $ git clone git://www.ginac.de/ginac.git * step 2, modify the 'ginac/config/config.rpath'. line 67: replace mingw* | pw32* | os2*) -> pw32* | os2*) line 57: insert mingw*) wl='Wl,' ;; line 135: replace cygwin* | mingw* | pw32*) -> cygwin* | pw32*) line 185: replace cygwin* | mingw* | pw32*) -> cygwin* | pw32*) line 199: replace linux*) -> linux* | mingw*) * step 3. $ cd ginac $ autoreconf -i * step 4, modify the 'ginac/configure'. line 15539: delete * step 5, modify the 'ginsh/ginsh_lexer.lpp'. line 37: replace #include "ginsh_parser.hpp" -> #include "ginsh_parser.h" * step 6. $ export CLN_LIBS="-L/usr/local/lib -lcln -lgmp" $ export CLN_CFLAGS="-I/usr/local/include" $ ./configure $ make $ make install * test 1. // hello.cpp ---------------- #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main() { symbol x("x"), y("y"); ex poly; for (int i=0; i<3; ++i) poly += factorial(i+16)*pow(x,i)*pow(y,2-i); cout << poly << endl; return 0; } $ g++ -c -O3 -I/usr/local/include -L/usr/local/bin hello.cpp $ g++ -o hello hello.o -lginac -lcln -lgmp &> result.txt $ ./hello 355687428096000*x*y+20922789888000*y^2+6402373705728000*x^2 * test 2. // hello2.cpp ---------------- #include <iostream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; ex HermitePoly(const symbol & x, int n) { ex HKer=exp(-pow(x, 2)); // uses the identity H_n(x) == (-1)^n exp(x^2) (d/dx)^n exp(-x^2) return normal(pow(-1, n) * diff(HKer, x, n) / HKer); } int main() { symbol z("z"); for (int i=0; i<6; ++i) cout << "H_" << i << "(z) == " << HermitePoly(z,i) << endl; return 0; } $ g++ -c -O3 -I/usr/local/include -L/usr/local/bin hello2.cpp $ g++ -o hello2 hello2.o -lginac -lcln -lgmp &> result.txt $ ./hello2 H_0(z) == 1 H_1(z) == 2*z H_2(z) == -2+4*z^2 H_3(z) == -12*z+8*z^3 H_4(z) == 12+16*z^4-48*z^2 H_5(z) == 120*z+32*z^5-160*z^3

1 0

Numeric and Constants
by David Narvaez 13 Jul '13

13 Jul '13

Hi, While implementing a GiNaC based backend for my software I decided to replace doubles by GiNaC::exs. That seemed to work fine until I actually needed the final native value of the expression, so I changed my mind to use GiNaC::numerics instead, which works fine except for constants, since I'm unable to do things like GiNaC::numeric angle = 2 * GiNaC::Pi; is there a quick way to pass from a constant to a numeric? (or alternatively a way to turn expressions into doubles and I should go back to GiNaC::ex). Thanks in advance. David E. Narvaez

1 0

variable isolation
by Andrea Monzali 12 Jun '13

12 Jun '13

Good evening everybody, I'm new to ginac, i read the tutorial but i've a question, if I have an equation expression, something like symbol x("x"), y("y"); ex eq = 2*x + 3*x*y - 5*y == 0 is there a method to isolate a single variable, i.e. cout << isolate(eq, x) << endl; // x == 5*y*(2 + 3*y)^(-1) cout << isolate(eq, y) << endl; // y == 2*x*(5 - 3*x)^(-1) thanks for all replies.. -am

2 1

Fwd: find series roots with Ginac
by Jey Nadarasa 22 Apr '13

22 Apr '13

Hello everyone, I would like to know if it is possible to find series roots with Ginac ? That's all for the moment, Advance thanks for your answer, Best Regards, JN

1 0

[(Possible) Bug Report] g++-4.8 (slightly) breaks GiNaC::symbol
by Song Gao 25 Mar '13

25 Mar '13

Hey, I had a weird problem using GiNaC with g++-4.8. It seems that GiNac::symbol will be deleted twice if it's not dynamically declared. In other words, some resources in it is free'd before the automatic destructor is called when the scope finishes. The minimum code to re-produce the issue is: #include <ginac/ginac.h> int main() { GiNaC::symbol x("x"); return 0; } $ g++-4.8 -g -o main main0.cc -lginac -lcln && ./main main(69018) malloc: *** error for object 0x7fff73275570: pointer being freed was not allocated It seems this issue only exists in g++-4.8. g++-4.7 is fine. I briefly looked at symbol.h but did not manage to locate the problem. I was hoping someone here more familiar with the code base could fix it or provide some direction for me. Thanks. -- Song Gao http://song.gao.io/

2 1

GiNaC Multithreading
by Florian 09 Mar '13

09 Mar '13

Hi, I have read the issue of multithreading from 2010 and I wonder whether there has been attempt to make GiNaC thread-safe. I think that it shouldn't be a big deal to make the shared data like the refcounter atomic (e.g. using std::atomic< T >) or to the synchronize it's manipulation (e.g. using std::mutex). Best regards, Florian -- Florian Corzilius RWTH Aachen University Computer Science Department, Informatik 2 D-52056 Aachen, Germany http://www-i2.informatik.rwth-aachen.de/i2/corzilius/ phone: +49 241 80-21244

1 0

cmake msvc build
by Philipp Kraus 25 Feb '13

25 Feb '13

Hello, I love GiNaC and I can build it with MSYS very well. My question is for a feature request is: Can a cmake build added, so that I can build GiNaC native under MSVC? I think this can be a nice solution. Thanks Phil

3 3

subs does not work
by Zhao Li 21 Feb '13

21 Feb '13

Here I am trying to make a substitution. For testing, the function "has" works as shown below. But the "subs" just does not implement the substitution. cout << has(temp, JF1(wild(1), GA(wild(5)), wild(2))*JF1(wild(3), GA(wild(5)), wild(4)) , has_options::algebraic ) << endl; cout << temp.subs( JF1(wild(1),GA(wild(5)),wild(2))*JF1(wild(3),GA(wild(5)),wild(4)) == 2, subs_options::algebraic ) << endl; The functions "JF1" and "GA" are both declared using DECLARE_FUNCTION and commutative. DECLARE_FUNCTION_1P( GA ); REGISTER_FUNCTION( GA, dummy() ); Best regards, Zhao Theory Division Institute of High Energy Physics Chinese Academy of Sciences

2 1

Creating a number of FUNCTION_1P
by Dale Lukas Peterson 20 Feb '13

20 Feb '13

I am trying to declare N GiNaC 1 parameter functions (DECLARE_FUNCTION_1P), each of which represents a different function of the same single parameter (in my case, this parameter represents time). My primary motivation for doing this is to be able to easily create a bunch of symbol-like objects that all depend on time so that I can trust that the chain rule will work properly when differentiating with respect to this single parameter. I think this has come up before, but I haven't seen a way to declare a vector or list of these functions without N explicit statements to DECLARE_FUNCTION_1P. Is there a simple way to do this in GiNaC? Thanks, Luke

1 1