- GiNaC-list - ginac.de

GiNaC kudos
by Stephen Montgomery-Smith 10 Aug '10

10 Aug '10

Hi Guys, I just wrote some programs that use GiNaC to multiply polynomials whose degrees can get into the 1000's, and whose coefficients are as large as 10^10000. GiNaC performed extraordinarily well, and even went a little faster than Mathematica! I am very impressed! (If you are interested, I am trying to analyze results of Luria-Delbruck experiments using Bayesian statistics.) Stephen

2 3

arithmetic error on concat two expressions
by Kraus Philipp 07 Aug '10

07 Aug '10

Hi, I'm new at GiNaC and I use the lib for create a gradient descent algorithm (symbolic calculating). I have a formula "a*x^3+b*y^4" and I subsitute the formula in the error function "0.5*(t - <formula>)^2". After the substitution GiNaC returns the expression "(0.5)*(a*x^3+b*y^4-t)^2". I use this code: GiNaC::ex m_expression; GiNaC::symtab m_exprtable; try { m_expression = l_parser( "a*x^3+b*y^4" ); m_exprtable = l_parser.get_syms(); } catch (...) {} GiNaC::ex l_full; GiNaC::symtab l_table(m_exprtable); l_table["formula"] = m_expression; GiNaC::parser l_parser(l_table); try { l_full = l_parser( "0.5*(t - formula)^2" ); // (*) } catch (...) {} std::cout << l_full << std::endl; I using the version 1.5.8 (on 1.5.7 I had got a compiler error). Is this the correct way for concating two symbolic expression? Do I have forgot anything? Thanks Phil

4 7

is ginac thread safe for this case?
by Cristobal Navarro 26 Jul '10

26 Jul '10

will ginac be thread safe for this case? 1) threads access diferent expressions which are independant 2) each thread will increase the expression its working in.. for example thread1...expression1 += x+y thread2.... expression2 += 2*x.... and so on. 2) no function evaluation is done never. 3) however, all expressions share the same global symbols ( eg, expressions made of global symbols x and y ). 4) threads access the same common matrix of expressions, but two threads will never access the same cell. i'am testing at the moment, and is working but in the long run i am not sure if problems will arise. any information, welcome. thanks in advance. Cristobal <http://www.youtube.com/neoideo>

5 22

Store matrix, then read it?
by Cristobal Navarro 07 Jul '10

07 Jul '10

Hello, i am experiencing a little trouble when reading my stored matrix. when i write and read the same matrix i stored, it works. but if i only run the program to read the matrix file, i get this error. terminate called after throwing an instance of 'std::runtime_error' what(): unknown matrix dimensions in archive This is a one-file example source code i made for you guys so you can test, its based on the same problem i have on my bigger code. try running it as it is, then comment the "write" method and run it again so it only reads the matrices, youll get the error eventually. but you will never get it if you write and read, its weird. how can this be solved? //CODE #include <fstream> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; symbol q("q"),v("v"); void read( string filename, GiNaC::matrix* qvCol, GiNaC::matrix* qvMat ); void write( string filename, GiNaC::matrix qvCol, GiNaC::matrix qvMat ); int main(int argc, char** argv){ GiNaC::matrix qvMat(2,2), qvCol(2,1), rm1, rm2; //mat1 qvCol(0,0) = q*v; qvCol(1,0) = 2*q*v; //mat2 qvMat(0,0) = q*v; qvMat(0,1) = 2*q*v; qvMat(1,0) = 3*q*v; qvMat(1,1) = 4*q*v; cout << endl << "m1: " << endl << qvCol << endl; cout << "m2: " << endl << qvMat << endl << endl; write("example.gin", qvCol, qvMat); read("example.gin", &rm1, &rm2); cout << endl << "rm1: " << endl << rm1 << endl; cout << "rm2: " << endl << rm2 << endl << endl; return 0; } void read( string filename, GiNaC::matrix* qvCol, GiNaC::matrix* qvMat ){ cout << "Main::Read File......."; fstream fr(filename.c_str(),ios::in|ios::binary); if( !fr.is_open() ){ printf("bad file\n"); return; } GiNaC::archive_node canCol, canMat; GiNaC::lst syms(q, v); fr >> canCol; fr >> canMat; qvCol->read_archive(canCol, syms); qvMat->read_archive(canMat, syms); fr.close(); cout << "OK: " << filename << endl; } void write( string filename, GiNaC::matrix qvCol, GiNaC::matrix qvMat ){ cout << "Main::Write File......"; fflush(stdout); //abrir archivo para escritura fstream fw(filename.c_str(),ios::out|ios::binary); GiNaC::archive_node anCol, anMat; qvCol.archive( anCol ); qvMat.archive( anMat ); fw << anCol; fw << anMat; fw.close(); cout << "OK: " << filename << endl; } //END

2 2

GiNaC release 1.5.8
by Jens Vollinga 06 Jul '10

06 Jul '10

Hi everybody, GiNaC 1.5.8 is out and available. The changes are: * Changed default behaviour of conjugate() for user-defined functions to avoid wrong evaluations. A user-defined conjugate_func must now be explicitly registered. * Fixed: Parser now handles abbreviations as advertized in the manual. * Fixed two bugs in the GCD code (one infinite loop, one miscalculation). * Fixed dangerous iterator use. * Fixed a memory leak in excompiler. * Fixed several syntax issues wrt to gcc 4.5 and wrt to compile option -std=cxx0x. * Fixed a wrong URL and a texinfo problem in the manual. * Added degree_vector() utility function. As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

1 0

parser and subs problem
by Felipe Bordeu 29 Jun '10

29 Jun '10

Hello, I'm trying to use the ginac parser, but doesn't work as expected. Page 81 of GiNaC 1.5.7 (an open framework ...) tutorial (29 march 2010) " It’s also possible to map input (sub)strings to arbitrary expressions: GiNaC::symbol x, y; GiNaC::symtab table; table["x"] = x+log(y)+1; GiNaC::parser reader(table); GiNaC::ex e = reader("5*x^3 - x^2"); " but I get the error terminate called after throwing an instance of 'std::invalid_argument' what(): find_or_insert_symbol: name "x" does not correspond to a symbol Abandon function find_or_insert_symbol in file ginac/parser/parse_context.cpp:49 can handle only symbols not expressions. So I tried to use only symbols and then do some subs to get the correct expression without success. GiNaC::symbol x, y; GiNaC::ex expr1 = x*y; GiNaC::ex mat1 = GiNaC::symbolic_matrix(2, 2, "K"); GiNaC::ex mat2 = GiNaC::symbolic_matrix(2, 2, "T"); GiNaC::ex expr2 = expr1.subs(GiNaC::lst(x==mat1,y==mat2)); cout << expr2 << endl; cout << expr2.evalm() << endl; and I get the error : [[K00,K01],[K10,K11]]*[[T00,T01],[T10,T11]] terminate called after throwing an instance of 'std::runtime_error' what(): matrix::mul_scalar(): non-commutative scalar Abandon do I have to use indexed objects to get the correct answer??? Thanks -- Felipe

2 1

1by1 matrix to scalar
by Felipe Bordeu 23 Jun '10

23 Jun '10

Hello, I have a little problem, I'm dealing functions made of matrix products. My question is, how can "force" ginac to convert automatically all the 1by1 matrices to scalars?. or do I have to do it manually?? Here is a simple example. void main2(){ ex m1=symbolic_matrix(1, 3, "x"); ex m2=symbolic_matrix(3, 1, "y"); cout << m1 << std::endl; cout << m2 << std::endl; ex m3 = (m1*m2).evalm(); cout << m3 << std::endl; cout << ex_to<matrix>(m3)(0,0)<< std::endl;; ex m4 = m3-ex_to<matrix>(m3)(0,0); cout << m4 << std::endl; }; the output : [[x0,x1,x2]] [[y0],[y1],[y2]] [[y1*x1+y2*x2+y0*x0]] y1*x1+y2*x2+y0*x0 -y1*x1+[[y1*x1+y2*x2+y0*x0]]-y2*x2-y0*x0 the last line "should" be zero ?-- thanks Felipe

4 5

writing symbols and expressions into binary files
by Cristobal Navarro 13 Jun '10

13 Jun '10

Hello, is it possible to write ginac symbols and expressions into a file as binary format?? i've tried this but i get segmentation fault when printing the read variables. you can try this one file example to see the problem. it writes to example.dat as binary format, then reads from it as binary format too. //CODE #include <ginac/ginac.h> #include <iostream> #include <fstream> using namespace std; using namespace GiNaC; int main(){ //WRITE if(1){ symbol a("a"), b("b"); ex myExp = a*b; cout << "Original: Symbols are '" << a << "' and '" << b << "', myExp = " << myExp << endl; fstream fw; fw.open( "example.dat", ios::out|ios::binary ); if( fw.is_open() ){ fw.write( (char*)&a, sizeof(a) ); fw.write( (char*)&b, sizeof(b) ); fw.write( (char*)&myExp, sizeof(myExp) ); fw.close(); } else{ printf("Error Opening File for writing\n"); } } //READ (SHOULD NOT DEPEND OF WRITE CODE) if(1){ symbol ca, cb; ex copyExp; fstream fr; fr.open( "example.dat", ios::in|ios::binary ); if( fr.is_open() ){ fr.read( (char*)&ca, sizeof(ca) ); fr.read( (char*)&cb, sizeof(cb) ); fr.read( (char*)&copyExp, sizeof(copyExp) ); fr.close(); cout << "Read from file: Symbols are '" << ca << "' and '" << cb << "' myExp = " << copyExp << endl; } else{ cout << "Error opening file for reading: does it exist?" << endl; } } return 1; } //CODE thanks in advance -- Cristobal <http://www.youtube.com/neoideo>

3 5

implicit differentiation
by Felipe Bordeu 02 Jun '10

02 Jun '10

Hi everyone, I created a class for GiNaC to treat implicit differentiation, here is an example: symbol x = get_symbol("x"); ----> I use a factory ex A = dersymbol("A").dependent_of(x); -----> my new class, and I tell that is dependent of x ex dAdx = A.diff(x); cout << A << endl; -----> "A(x)" cout << "dA/dx : " << dAdx ; -----> "A(x)_1x" cout << "d2A/dx2 : " << dAdx.diff(x) ; -----> "A(x)_2x" cout << "dA/dy : " << A.diff(get_symbol("y")) ; -----> "0" Inside dersymbol I store the dependents variables and the degree of the differentiation. now my problem is that I what to make the substitution of A to something else, like x^2: cout << "dA/dx with A=x^2 is : " << dAdx.subs(A==pow(x,2)) << endl; (I know that the previous line is not correct) to get something like : "2*x" Any ideas?? Thanks (first time in this mailing list) Felipe

2 1

Term ordering and compiling C++ code
by Doug 02 Jun '10

02 Jun '10

Hi, I'm running into a term ordering problem like the one mentioned in this post. http://permalink.gmane.org/gmane.comp.mathematics.ginac.general/1107 Relevant details about my application are below, for the curious. The short of it, though, is that the C++ code I am generating with ginac is so large (50MB) that the compiler can't handle it, and moreover, the values being represented by ginac symbols are array accesses that the compiler can't optimize away. So I need to create temporary values and do sub-expression elimination on the C++ code just so the compiler can have a chance. I have successfully done this using search-and-replace on the C++ code for other applications, but I guess I was just lucky that ginac printed terms in a consistent order. Is this a problem other ginac users have confronted? From the answer to the above post, it sounds like I need to write my own print function that can ensure expressions are printed in a predictable order. It looks like I can also inspect the expression tree and possibly do a proper job of common sub-expression elimination. But if there is an easier approach, I'd be very grateful for a pointer or two. Thanks, Doug Baker ------------------- Background Details ------------------- My application is related to visual odometry. I am computing the gradient and hessian for a function of 10 variables that involves converting quaternions to rotation matrices, applying those matrices to 3D points and projecting those points onto a 2D plane. I am generating code with ginac that is initially 50MB of C++ code. (And thank you for ginac being able to handle that monstrosity with such grace.) The input symbols are all array elements (e.g. symbol x("state[t+X]")) or function calls (e.g. symbol rTi11("rTi.getValue<1,1>()")). If I take the printed output from ginac as is, the C++ compiler can't properly optimize it because it doesn't know that state[t+X] is a constant or that getValue is constant. Moreover, the expressions have huge repeated expressions, and these bog the compiler down, too. I've been able to use a simple search-and-replace approach in other problems, and on this one that same approach has gotten me from 50MB to 6MB or so, but the term ordering issue is biting me now. For example, I have one sum (a+b+c+d) that occurs 232608 times. But there are 4! = 24 permutations of that sum that I'd have to replace. It's worse than that, though, because that's not the only expression like it. Some are bigger, like this: (-2.0*sK/sKSIJ+4.0*sK*sS2/pKSIJ2+-4.0*sI*sJ*sS/pKSIJ2). This has something like 72 permutations, not even counting the permutations of pKSIJ and pKSIJ2, which are themselves expressions like (a+b+c+d). And not only are there variations on this due to term ordering, similar expressions occur with different variables (e.g. with sJ instead of sK, etc.) So when all these expressions can be output in various orders, it's no longer feasible to generate all the permutations for search-and-replace. Support NPR 20 seconds at a time. www.twentysecondsatatime.org

5 18