September 2010 - GiNaC-devel

[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

Patch 2
by Jan Rheinländer 09 Sep '10

09 Sep '10

Problem: Return type of function is "ex". MSVC complains that it cannot convert "GiNaC::power*" to "ex" Solution: Changed expression to "GiNaC::power" and removed statusflags::dynallocated

3 2

Patch 3 updated description
by Jan Rheinländer 09 Sep '10

09 Sep '10

symmetry::compare_same_type(): const-correctness fix Problem: The method in question is const, therefore one should use const_iterator to access indices. Some compilers (for instance, MSVC) are very strict about this, and won't accept (invalid) code using non-const iterator. Solution: Changed it to const_iterator

1 0

Patch 5
by Jan Rheinländer 09 Sep '10

09 Sep '10

Problem: MSVC chokes on "struct" Solution: Removed it

1 0