November 2003 - GiNaC-list

differential operators
by Kai Ludwig 29 Feb '04

29 Feb '04

Hi All, I'm new to GiNac. Is there a way to execute symbolic differentiation of expressions that includes differential operators ? e.g. suppose an expression ex F = u * dx*u When differentiate that expression F.diff(u) it should result in the expression dx*u - u*dx Regards Kai -- http://echempp.sourceforge.net Kai Ludwig Universität Tübingen Institut für Organische Chemie Auf der Morgenstelle 18 D-72076 Tübingen Tel.: 07071/29-73049 Mail: kai.ludwig(a)uni-tuebingen.de

4 10

complex conjugation
by Andrius Kurtinaitis 12 Dec '03

12 Dec '03

Hello, it appears that the GiNaC expressions do not support the complex conjugation... :-( I need it, but I do not know yet how to implement it. Has anybody tried to do this? Clearly conjugation is a function ;-) Ok, i'll try to read the tutorial. Let's see if I get this working... Best regards Andrius Kurtinaitis

3 6

automatic differentiation and change expression ordering
by Andrius Kurtinaitis 20 Nov '03

20 Nov '03

Hello, I have two questions: 1. is it possible to prevent automatic differentiation of some expression? For example, when I do (with ginac 1.1.4): DECLARE_FUNCTION_1P( intensity ) static ex intensity_eval( const ex & r ) { return intensity(r).hold(); } // here i tried to hold it :-) REGISTER_FUNCTION( intensity, eval_func(intensity_eval) ) int main() { symbol r("r"); ex der1 = mul( diff( intensity(r), r ), r ).hold(); ex deriv = der1.diff(r); cout << deriv << endl; return 0; } I get the following output: D[0,0](intensity)(r)*r+D[0](intensity)(r) But I would like to get something like this: D[0]( D[0](intensity)(r)*r )(r) Is it possible with ginac? 2. What is the simplest way to change the order in multiplication so that the differential operators print after all other terms: r*D[0,0](intensity)(r) instead of D[0,0](intensity)(r)*r And the third one ;-) Is there maybe some work done to output fderivative in latex style? Thanks in advance Andrius

2 1

CLN + GiNaC on MinGW
by Jonathan Brandmeyer 05 Nov '03

05 Nov '03

I am attempting to build CLN (and later, GiNaC) on my Windows XP machine with Mingw. I am using g++ 3.2.3 and the latest binary utilities. I have a build of GMP 4.2.1 for use as a backend for CLN (installed under c:\msys\1.0\local\). CLN was configured with CPPFLAGS=-I/usr/local/include LIBS=/usr/local/lib/libgmp-3.dll CXXFLAGS="-O2 -finline-functions -fomit-frame-pointer -march=i686 -finline-limit=1000" ./configure --with-gmp. The build progressed and after a couple hours the link failed for the shared library. It looks a lot like what was described on ginac-devel when a fellow had problems on Alpha (thread begins here: http://www.ginac.de/lists/ginac-devel/msg00091.html) except that I am using a relatively modern toolchain. The linker errors are of the form: cl_t_c2.lo(.text+0x4e):cl_t_c2.cc: undefined reference to `_GLOBAL__I_cl_module__cl_prin_globals__firstglobalfun' Any thoughts? -Jonathan Brandmeyer P.S. I am using msys to drive the build, but I do not want to use Cygwin.

3 8

GiNaC 1.1.5 now available
by Christian Bauer 05 Nov '03

05 Nov '03

Hello! "Release Early, Release Often" The GiNaC Group is pleased to announce the immediate availability of GiNaC 1.1.5 which contains lots of improvements to the nestedsums functions: - Harmonic polylogarithms now numerically evaluate for arbitrary arguments (parameter must still be positive integers). - The zeta function now can also be given a lst as a parameter in which case it becomes a multiple zeta value. The use of mZeta is deprecated. - The order of parameters for the multiple polylogarithm has been corrected. - Documentation for the nested sums functions zeta, harmonic polylog, multiple polylog, etc. has been added. Run, don't walk, to everyone's favorite FTP server: ftp://ftpthep.physik.uni-mainz.de/pub/GiNaC/ Have a nice day, Christian -- / Physics is an algorithm \/ http://www.uni-mainz.de/~bauec002/

1 0

CLN + GiNaC on Apple OSX
by Jonathan Brandmeyer 04 Nov '03

04 Nov '03

In addition to building on Windows/MinGW, I am also trying to get this going on Mac OSX. I have had problems with Apple's early releases of GCC, so I am using a local build of GCC-3.3.1. I believe that the binary utilities are mostly GNU with Apple's linker. (as reports 1.38 as its version). I am building CLN against GMP 4.1.2. ( I said 4.2.1 earlier, but I really installed 4.1.2). CLN was configured with CC=gcc-3.3 CXX=g++-3.3 CXXFLAGS="-O2 -finline-functions" CFFLAGS="-Dunix -D__unix" (-D__unix was needed for a random number constructor) ./configure --with-gmp The build fails with cascading syntax errors in modinteger/cl_MI.cc and several other files. Thank you for any insight you may have on this problem, -Jonathan Brandmeyer g++-3.3 -v: Reading specs from /usr/local/lib/gcc-lib/powerpc-apple-darwin6.8/3.3.1/specs Configured with: ../gcc-3.3.1/configure --enable-languages=c,c++ --enable-threads=posix --program-suffix=-3.3 Thread model: posix gcc version 3.3.1 The build fails in the modinteger directory with this error: make SUBDIR=modinteger /bin/sh ../libtool --mode=compile g++-3.3 -O2 -finline-functions -Dunix -D__unix -I../include -I../include -I./integer/hash -I./integer -I./base/hash -I./base/digitseq -I./base/digit -Ibase -I./base -c ./modinteger/cl_MI.cc rm -f .libs/cl_MI.lo g++-3.3 -O2 -finline-functions -Dunix -D__unix -I../include -I../include -I./integer/hash -I./integer -I./base/hash -I./base/digitseq -I./base/digit -Ibase -I./base -c ./modinteger/cl_MI.cc -fno-common -DPIC -o .libs/cl_MI.lo In file included from modinteger/cl_MI.cc:80: modinteger/cl_MI_pow2.h:16: error: parse error before numeric constant modinteger/cl_MI_pow2.h: In function `const cln::cl_I cln::pow2_reduce_modulo(...)': modinteger/cl_MI_pow2.h:19: error: `x' undeclared (first use this function) modinteger/cl_MI_pow2.h:19: error: (Each undeclared identifier is reported only once for each function it appears in.) modinteger/cl_MI_pow2.h: In function `const cln::_cl_MI cln::pow2_canonhom(cln::cl_heap_modint_ring*, const cln::cl_I&)': modinteger/cl_MI_pow2.h:24: warning: cannot pass objects of non-POD type `const class cln::cl_I' through `...'; call will abort at runtime modinteger/cl_MI_pow2.h: At global scope: modinteger/cl_MI_pow2.h:27: error: parse error before numeric constant modinteger/cl_MI_pow2.h: In function `const cln::_cl_MI cln::pow2_plus(...)': modinteger/cl_MI_pow2.h:30: error: `y' undeclared (first use this function) modinteger/cl_MI_pow2.h: At global scope: modinteger/cl_MI_pow2.h:34: error: parse error before numeric constant modinteger/cl_MI_pow2.h:41: error: parse error before numeric constant modinteger/cl_MI_pow2.h:48: error: parse error before numeric constant modinteger/cl_MI_pow2.h:54: error: parse error before numeric constant modinteger/cl_MI_pow2.h:61: error: parse error before numeric constant modinteger/cl_MI_pow2.h:80: error: parse error before numeric constant modinteger/cl_MI_pow2.h:105: error: parse error before numeric constant modinteger/cl_MI_pow2.h:126: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:126: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:126: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:138: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*)' modinteger/cl_MI_pow2.h:138: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:138: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:138: error: invalid conversion from `const cln::cl_MI_x (*)(...)' to `const cln::cl_MI_x (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:138: error: invalid conversion from `const cln::cl_MI_x (*)(...)' to `const cln::cl_MI_x (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_pow2.h:138: error: invalid conversion from `const cln::cl_I (*)(...)' to `const cln::cl_I (*)(cln::cl_heap_modint_ring*, const cln::cl_I&)' In file included from modinteger/cl_MI.cc:81: modinteger/cl_MI_pow2m1.h:16: error: parse error before numeric constant modinteger/cl_MI_pow2m1.h: In function `const cln::_cl_MI cln::pow2m1_canonhom(cln::cl_heap_modint_ring*, const cln::cl_I&)': modinteger/cl_MI_pow2m1.h:46: warning: cannot pass objects of non-POD type ` const class cln::cl_I' through `...'; call will abort at runtime modinteger/cl_MI_pow2m1.h: At global scope: modinteger/cl_MI_pow2m1.h:49: error: parse error before numeric constant modinteger/cl_MI_pow2m1.h:58: error: parse error before numeric constant modinteger/cl_MI_pow2m1.h:79: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_pow2m1.h:79: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_pow2m1.h:79: error: invalid conversion from `const cln::cl_I (*)(...)' to `const cln::cl_I (*)(cln::cl_heap_modint_ring*, const cln::cl_I&)' In file included from modinteger/cl_MI.cc:82: modinteger/cl_MI_pow2p1.h:16: error: parse error before numeric constant modinteger/cl_MI_pow2p1.h: In function `const cln::_cl_MI cln::pow2p1_canonhom(cln::cl_heap_modint_ring*, const cln::cl_I&)': modinteger/cl_MI_pow2p1.h:53: warning: cannot pass objects of non-POD type ` const class cln::cl_I' through `...'; call will abort at runtime modinteger/cl_MI_pow2p1.h: At global scope: modinteger/cl_MI_pow2p1.h:56: error: parse error before numeric constant modinteger/cl_MI_pow2p1.h:67: error: parse error before numeric constant modinteger/cl_MI_pow2p1.h:90: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_pow2p1.h:90: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_pow2p1.h:90: error: invalid conversion from `const cln::cl_I (*)(...)' to `const cln::cl_I (*)(cln::cl_heap_modint_ring*, const cln::cl_I&)' In file included from modinteger/cl_MI.cc:83: modinteger/cl_MI_montgom.h:87: error: parse error before numeric constant modinteger/cl_MI_montgom.h:93: error: parse error before numeric constant modinteger/cl_MI_montgom.h:99: error: parse error before numeric constant modinteger/cl_MI_montgom.h:106: error: parse error before numeric constant modinteger/cl_MI_montgom.h:112: error: parse error before numeric constant modinteger/cl_MI_montgom.h:118: error: parse error before numeric constant modinteger/cl_MI_montgom.h:132: error: parse error before numeric constant modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*)' modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::cl_I&)' modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::_cl_MI (*)(...)' to `const cln::_cl_MI (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::cl_MI_x (*)(...)' to `const cln::cl_MI_x (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::cl_MI_x (*)(...)' to `const cln::cl_MI_x (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&, const cln::_cl_MI&)' modinteger/cl_MI_montgom.h:158: error: invalid conversion from `const cln::cl_I (*)(...)' to `const cln::cl_I (*)(cln::cl_heap_modint_ring*, const cln::_cl_MI&)' make[2]: *** [cl_MI.lo] Error 1

2 3