- GiNaC-devel - ginac.de

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-229-g9506030
by git＠ginac.de 15 Sep '10

15 Sep '10

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 9506030d2059ef868d5da8dcbe714c9c5b50ff50 (commit) from 8cce492692eb74d636a50f420fd79877e5d4981a (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 9506030d2059ef868d5da8dcbe714c9c5b50ff50 Author: Richard Kreckel <kreckel(a)ginac.de> Date: Wed Sep 15 09:11:57 2010 +0200 Be more careful about final top-level substitution. Substituting x==log(x) in exp(x) erroneously returned log(x) because of a final subst(x==log(x)) after having eval'ed exp(log(x)) -> x. This final substitution is wrong in the general case. On the other hand, the intent is to syntactically substitute functions of a given kind etc. This patch suppresses the final top-level substitution unless the intermediate result is a container. Thanks to Burcin Erocal for reporting this bug (originally described by Kees van Schaik on sage-support(a)googlegroops.com) ----------------------------------------------------------------------- Summary of changes: check/exam_misc.cpp | 9 +++++++++ doc/tutorial/ginac.texi | 2 +- ginac/container.h | 24 ++++++++++++++++++++---- 3 files changed, 30 insertions(+), 5 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, ginac_1-5, updated. release_1-4-0-249-g90ad10b
by git＠ginac.de 15 Sep '10

15 Sep '10

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, ginac_1-5 has been updated via 90ad10b58d02365a407b2d84d8b93e50030feaa5 (commit) from db47e259edabbed2b8eb2fa1e4ed93588c26e304 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 90ad10b58d02365a407b2d84d8b93e50030feaa5 Author: Richard Kreckel <kreckel(a)ginac.de> Date: Wed Sep 15 09:11:57 2010 +0200 Be more careful about final top-level substitution. Substituting x==log(x) in exp(x) erroneously returned log(x) because of a final subst(x==log(x)) after having eval'ed exp(log(x)) -> x. This final substitution is wrong in the general case. On the other hand, the intent is to syntactically substitute functions of a given kind etc. This patch suppresses the final top-level substitution unless the intermediate result is a container. Thanks to Burcin Erocal for reporting this bug (originally described by Kees van Schaik on sage-support(a)googlegroops.com) ----------------------------------------------------------------------- Summary of changes: check/exam_misc.cpp | 9 +++++++++ doc/tutorial/ginac.texi | 2 +- ginac/container.h | 24 ++++++++++++++++++++---- 3 files changed, 30 insertions(+), 5 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

Patch 6
by Jan private 14 Sep '10

14 Sep '10

Problem: MSVC creates different symbols for "extern const ex _ex0" if it is declared at top-level (inside GiNaC namespace) and if it is declared inside the body of a function (also inside GiNaC namespace). When linking, unresolved external symbols occur. Solution: Moved the declaration out of the body of the function Note: Instead of #include utils.h the same effect would probably be achieved by namespace GiNaC { extern const ex _ex0; } at the beginning of the file.

2 3

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-228-g8cce492
by git＠ginac.de 13 Sep '10

13 Sep '10

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 8cce492692eb74d636a50f420fd79877e5d4981a (commit) from 5f896fa7f59bbce727e4bba23df9c4bbdbb55c29 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 8cce492692eb74d636a50f420fd79877e5d4981a Author: Richard Kreckel <kreckel(a)ginac.de> Date: Mon Sep 13 22:39:05 2010 +0200 Fix autoconf warning about AC_INIT use. Apparently, autoconf doesn't like the angle brackets in the bug-report argument any more. It comlains "not a literal: <ginac-list(a)ginac.de>". So let's remove it. And while at it, also provide tarname and url arguments. ----------------------------------------------------------------------- Summary of changes: configure.ac | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, ginac_1-5, updated. release_1-4-0-248-gdb47e25
by git＠ginac.de 13 Sep '10

13 Sep '10

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, ginac_1-5 has been updated via db47e259edabbed2b8eb2fa1e4ed93588c26e304 (commit) from cc94094751459129e1392a93dfa0264adac789a5 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit db47e259edabbed2b8eb2fa1e4ed93588c26e304 Author: Richard Kreckel <kreckel(a)ginac.de> Date: Mon Sep 13 22:39:05 2010 +0200 Fix autoconf warning about AC_INIT use. Apparently, autoconf doesn't like the angle brackets in the bug-report argument any more. It comlains "not a literal: <ginac-list(a)ginac.de>". So let's remove it. And while at it, also provide tarname and url arguments. ----------------------------------------------------------------------- Summary of changes: configure.ac | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

Patch7
by Jan private 11 Sep '10

11 Sep '10

Hello Alexei, here is another hairy one, which is a compiler bug in my opinion: Problem: MSVC does not include the exprseq::info() method, giving unresolved symbols when linking Solution (rather, workaround): Add an unarchiver :-) Note: Someone else suggested adding a dummy method in exprseq.cpp basic* dummy_func(void) { return new exprseq(); } but I thought that maybe an unarchiver is at least remotely useful... Have fun! Jan

2 1

Patch 8
by Jan private 11 Sep '10

11 Sep '10

Hello Alexei, this is the last issue of GiNaC. All 58 tests passed sucessfully. CLN has a lot of issues left, though :-( Problem: MSVC cannot handle string constants above 16K Solution: Split them up into pieces Best regards, Jan

1 0

Patch 4
by Jan private 10 Sep '10

10 Sep '10

Hello Alexei, OK, so using compiler.h the patch reduces to a fourliner in that file because compiler.h is being included in all the concerned files (parse_binop_rhs.cpp, parser.cpp) already anyway. The sizeof() problem miraculously disappeared... Problem: 1. MSVC does not have __func__, but has __FUNCTION__ 2. MSVC does not have __alignof__, but has __alignof Solution: Added #defines to compiler.h

1 0

proposed integration algorithm
by Jose Antonio Garcia Peiro 10 Sep '10

10 Sep '10

Hi I am very interested in GiNaC, but I also need symbolic integration, and GiNaC only compute basic integrals. I wrote a simple recursive routine that applies a basic integration rules, matching the input function. It is a good way to implement integration algoritm? Thanks int(c,x) = c*x int(x,x) = x^2/2 int(f+g,x) = int(f,x) + int(g,x) int(c*f,x) = c*int(f,x) int( f^c*1,x) = f^(c+1)/(c+1) or log(f) int( f^c*f',x) = f^(c+1)/(c+1) or log(f) some examples: cout << "int(1+x,x) = " << myinteg( 1 + x, x ) << endl; // 1/2*x^2+x cout << "int(3*x^2+x+6,x) = " << myinteg( 3*pow(x,2)+x+6, x ) << endl; // 1/2*x^2+x^3+6*x cout << "int((3*x+5)^2*3,x) = " << myinteg( pow(3*x+5,2)*3, x ) << endl; // 1/3*(5+3*x)^3 cout << "int((3*x+5)^(-1)*3,x) = " << myinteg( pow(3*x+5,-1)*3, x ) << endl; // log(5+3*x) cout << "int((3*x+5)^b*3,x) = " << myinteg( pow(3*x+5,b)*3, x ) << endl; // (1+b)^(-1)*(5+3*x)^(1+b) the code (with debug messages): ex myinteg( const ex &f, const symbol &x ){ exmap repls; cout << endl << "f(x) = " << f << endl; //int(c,x) c*x if( !has( f, x ) ){ cout << "match int(c,x)" << endl; return f*x; } //int(x,x) x^2/2 if( f.is_equal(x) ){ cout << "match int(x,x)" << endl; return pow(x,2)/2; } //int(f+g,x) int(f,x)+int(g,x) if( f.match( wild(0)+wild(1), repls ) ){ cout << "match int(f+g,x)" << endl; cout << "wild(0)+wild(1) = " << repls << endl; ex f, g; for( exmap::const_iterator i=repls.begin(); i!=repls.end(); ++i){ if(i->first.is_equal(wild(0)) ){ f = i->second; }else if(i->first.is_equal(wild(1))){ g = i->second; } } cout << "f = " << f << endl; cout << "g = " << g << endl; return myinteg( f, x ) + myinteg( g, x ); } //int(c*f,x) c*int(f,x) if( f.match( wild(0)*wild(1), repls ) ){ cout << "match int(c*f,x)" << endl; cout << "wild(0)*wild(1) = " << repls << endl; ex c, f; for( exmap::const_iterator i=repls.begin(); i!=repls.end(); ++i){ if(i->first.is_equal(wild(0)) ){ c = i->second; }else if(i->first.is_equal(wild(1))){ f = i->second; } } if( !has( c, x ) && has( f, x )){ return c*myinteg( f, x ); }if( has( c, x ) && !has( f, x )){ return f*myinteg( c, x ); } } //int( pow( f, c ), x ) pow( f, c+1 ) * pow( c+1, -1 ) if c != -1, or ln( f ) if c == -1 if( f.match( pow(wild(0),wild(1)), repls ) ){ cout << "match int( pow( f, c ), x )" << endl; cout << "pow(wild(0),wild(1)) = " << repls << endl; ex f, c; for( exmap::const_iterator i=repls.begin(); i!=repls.end(); ++i){ if(i->first.is_equal(wild(0)) ){ f = i->second; }else if(i->first.is_equal(wild(1))){ c = i->second; } } cout << "f = " << f << endl; cout << "c = " << c << endl; if( has( f, x ) && !has( c, x ) ){ ex dft = f.diff(x,1); ex df = 1; cout << "df/dft = " << df/dft << endl; if( !has(df/dft,x) ){ cout << "!has(df/dft,x)" << endl; if( c.is_equal(-1) ){ return log( f ) * (df/dft); }else{ cout << "df = " << df << endl; cout << "dft = " << dft << endl; cout << "(df/dft) = " << (df/dft) << endl; cout << "(1/dft) = " << (1/dft) << endl; cout << "pow(dft,-1) = " << pow(dft,-1) << endl; cout << "f = " << f << endl; cout << "c = " << c << endl; cout << "(pow( f, c+1 )/( c+1 )) = " << (pow( f, c+1 )/( c+1 )) << endl; cout << "pow( f, c+1 )*pow( c+1, -1) = " << pow( f, c+1 )*pow( c+1, -1) << endl; cout << "(pow( f, c+1 )/( c+1 ))*(df/dft) = " << (pow( f, c+1 )/( c+1 ))*(df/dft) << endl; cout << "pow( f, c+1 )*pow( c+1, -1) = " << pow( f, c+1 )*pow( c+1, -1)*df*pow( dft, -1) << endl; return (pow( f, c+1 )/( c+1 )) * (df/dft); } } } } //int( pow( f, c ) * f', x ) pow( f, c+1 ) * pow( c+1, -1 ) if c != -1, or ln( f ) if c == -1 if( f.match( wild(0)*pow(wild(1),wild(2)), repls ) ){ cout << "match int( pow( f, c ) * f', x )" << endl; cout << "wild(0)*pow(wild(1),wild(2)) = " << repls << endl; ex df, f, c; for( exmap::const_iterator i=repls.begin(); i!=repls.end(); ++i){ if(i->first.is_equal(wild(0)) ){ df = i->second; }else if(i->first.is_equal(wild(1))){ f = i->second; }else if(i->first.is_equal(wild(2))){ c = i->second; } } cout << "f' = " << df << endl; cout << "f = " << f << endl; cout << "c = " << c << endl; if( has( f, x ) && !has( c, x ) ){ ex dft = f.diff(x,1); cout << "df/dft = " << df/dft << endl; if( !has(df/dft,x) ){ cout << "!has(df/dft,x)" << endl; if( c.is_equal(-1) ){ return log( f ) * (df/dft); }else{ return pow( f, c+1 )*pow( c+1, -1) * (df/dft); } } } } return 0; }

1 0

Patch 4
by Jan Rheinländer 09 Sep '10

09 Sep '10

MSVC unknown functions/symbols Problem: 1. MSVC does not know __func__ 2. MSVC does not like sizeof(void*) 3. MSVC does not have __alignof__ Solution: 1. a) Put AC_CPP_FUNC in acinclude.m4 b) Added check to configure.ac c) Changed parser/debug.h 2. a) Added AC_CHECK_SIZEOF(void*) to configure.ac b) Changed polynomial/primes_factory.h 3. a) Added AC_CHECK_ALIGNOF(void*) to configure.ac b) Changed polynomial/primes_factory.h

3 4