- GiNaC-list - ginac.de

GiNaC 1.1.6 available!
by Jens Vollinga 22 Jan '04

22 Jan '04

Hello! GiNaC 1.1.6 "Sho Chiku Bai" is now available! It brings to you dummy functions and several bug fixes. In detail: * Added a function option "dummy()" which means "no options". This simplifies the implementation of symbolic functions which are not to be further evaluated. * Removed a bug in the numerical evaluation of Li() that caused the system to hang for certain parameter combinations. * Fixed a bug in the calculation of hash values for indices that could lead to wrong results or bogus error messages from simplify_indexed(). * Fixed a bug in the evaluation of harmonic polylogarithms for complex arguments with positive imaginary part. So we have GiNaC-Team 3 points, cimex lectularius 0 points. It's where you expect it to be: ftp://ftpthep.physik.uni-mainz.de/pub/GiNaC/ Bye, Jens

1 0

Symbols that are implicitly dependent on other symbols
by Johan Rutgeerts 16 Jan '04

16 Jan '04

Hi, Is there some concept in Ginac of "symbols which are function of some other symbol, but without making this explicit"? I'll try to explain what I mean: I wanted to: - define some function f(q_1, q_2, ... q_m) - calculate its derivatives df/d(q_j) symbolically and evaluate for specific values of q_j. But just calculating the derivatives and substituting values for q_j with subs(q_j == value) is far too slow for my application. However, I can rewrite f: f(q_1, q_2,... ) == g(r_1, r_2,.., r_n), in witch the r_i are functions of q_1, q_2,...,q_m. And I have numeric methods to calculate the value of r_i and its derivatives d(r_i)/d(q_j) for specific values of q_j. So I would like to calculate df/d(q_j) as: df/d(q_j) == sum( dg/d(r_i) * d(r_i)/d(q_j), for all i), using symbolic differentiation for the dg/d(r_i) and using the numerically calculated values of the r_i and d(r_i)/d(q_j) in the evaluation (which should go a lot faster since g(r_i) is a fairly simple function). So for this I'd need to define g as a function of r_1 to r_n, but then the r_i must somehow know that they are function of q_1 to q_m, otherwise the symbolic derivative dg/d(q_j) is zero. So my question is: is there an easy way to do this in Ginac? What I have an in mind as a possible solution is some kind of symbol (let's call it "MySymbol" to be original:) that has awareness of its numerical value and the value of its derivatives to certain variables. This would then allow me to take the derivative and would eliminate the use of subs(), making evaluation faster. The class MySymbol would have methods like: void setValue(double value); void setDerivative(symbol var, unsigned int order, double value); so that you can write: symbol q_j; MySymbol r_i; r_i.setValue(1.0); // the numeric value of r_i is 1.0 r_i.setDerivative(q_j, 1, 3.0); // the numeric value of d(r_i)/d(q_j) is 3.0 r_i.setDerivative(q_j, 2, 2.0); // the numeric value of d^2(r_i)/d(q_j)^2 is 2.0 // all other derivatives are standard zero then e.g.: r_i.diff(q_j, 1) returns a "MySymbol" with "setValue(3.0)" and "setDerivative(q_j, 1, 2.0)", r_i.diff(q_j, 2) returns a "numeric(2.0)", r_i.diff(q_j, i) with i > 2 returns a "numeric(0)", r_i.diff(l, 1) returns a "numeric(0)" cause we never did a "setDerivative()" with l as parameter. This should allow to write things like: MySymbol y = r_i + r_i.diff(q_j, 2) + r_i.diff(l, 1); double f = ex_to<double>( y ).to_double(); cout << f << endl; // should be "3.0" and: cout << sin(r_i).diff(q_j,1) << endl; // should be "cos(1.0)*3.0" but also: cout << ( sin(r_i).diff(q_j,1) ).diff(q_j,1) << endl; which should return the value of "cos(1.0) * 2.0 - sin(1.0) * (3.0)^2" Does this seem feasible? Anyone see any big flaws, has objections, knows a better name then "MySymbol"? ;) A remark: It might also be nice if one could do things like: symbol q_j; MySymbol r_i; r_i.setValue(1.0); r_i.setDerivative(q, 1, 3.0); MySymbol y = r_i + r_i.diff(q_j, 1); double f = ex_to<double>( y ).to_double(); r_i.setValue(5.0); double g = ex_to<double>( y ).to_double(); // without reassigning y cout << f << ", " << g << endl; // prints "4.0, 8.0" but this implies coupling between x and the object returned by e.g. x.diff(...). Don't know if this is good. Many greetings, Johan Rutgeerts

1 0

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

composite function
by Volker Vogel 12 Dec '03

12 Dec '03

Hello, I'am trying to define composite functions, for example this one: / 1 for all x > 1 f(x) =| x*x for all 1 >= x >= -1 \ 1 for all x <-1 Later I want to evaluate this function with different x. I tried to do this allready using symbol x("x"); string polystr="if (abs(x)>1) {1}else{x*x}"; ex poly; poly=ex(polystr,x); cout << poly; but I do not get it working. I assume that I defined my polystr incorrect, anybody knows how to do that ? Grettings, Volker Vogel ______________________________________________________________________________ Horoskop, Comics, VIPs, Wetter, Sport und Lotto im WEB.DE Screensaver1.2 Kostenlos downloaden: http://screensaver.web.de/?mc=021110

2 1

Error message during compilation of cln 1.1.5
by Chris Dams 05 Dec '03

05 Dec '03

Hi everybody, Today I installed gcc 3.3.2, so I had to recompile several things. One of the things that I compiled was cln 1.1.5. This, however, gave the error message below. It is about something in the tests directory, so this did not stop me from building GiNaC (the "consistency of numeric types" exam of GiNaC is passed, so I suppose everything is okay), but it still is a bit irritating. Perhaps the maintainer of CLN, who is believed to be reading this mailing list, knows what could be done about this. Have a nice day, Chris Dams make[1]: Entering directory `/home/chrisd/download/cln-1.1.5/tests' gcc -g -O2 -I../include -I../include -I. -I../src -I../src/base -I../src/float -I../src/float/lfloat -I../src/integer -I../src/base/random -I../src/base/digitseq -I../src/base/digit -I../src/base -I../src/base -c exam.s -o exam.o exam.s: Assembler messages: exam.s:246: Error: file number 70 already allocated make[1]: *** [exam.o] Error 1 make[1]: Leaving directory `/home/chrisd/download/cln-1.1.5/tests' make: *** [install] Error 2

1 0

Intervals in GiNaC
by Simon Langley 01 Dec '03

01 Dec '03

Carl, I have an interval class which I use in GiNaC. It isn't as tidy as I'd like but it works in the following sense: * I'm only interested in polynomials so the interval class doesn't have other functions like sin/exp/etc just field operations and powers. * It assumes that either rational or s/w floating point arithmetic is used. There is a problem with h/w floating point in GiNaC because underflow/overflow errors may not get noted. That said because of the way GiNaC works you can mix floating point and rational arithmetic. * When used in Ginac expressions, interval expressions aren't simplified. I.e. you can write ex z = interval(-1,1) + x + interval(1, 3)*4 + interval(1.3, 3.1) but z.simplify() won't add the intervals and z.subs(x == 2.3) won't reduce the expression to an interval. I've written a procedure to evaluate a Ginac expression in intervals to an interval (again assuming you only need polynomial expressions). If you want it let me know. By the way are you working with Dan Richardson or has interest in intervals at Bath spread? Simon ~~~~~ Simon Langley c/o CEMS, UWE, BS16 1QY

1 0

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