January 2001 - GiNaC-devel

polynomial arithmetic over ginac
by Parisse Bernard 23 May '01

23 May '01

I have now a little bit improved the package of my polynomial and rational fractions utilities over GiNaC (0.6.4) It provides now gcd normalization, factor (square-free and rational roots), partfrac & integrate (assuming factor works on the denominator), for example: integrate '1/(x^4-1)^5' returns -1155/8192*log(1+x)-1155/4096*atan(x)+1155/8192*log(-1+x)+1/2048*(-893*x-1375*x^9+1755*x^5+385*x^13)*(-1+x^4)^(-4) # Time for (integrate)0.04 You can retrieve and test it at ftp://fourier.ujf-grenoble.fr/pub/hp48/giac.tgz Bye, Bernard Parisse

4 9

sqrfree
by Parisse Bernard 11 Mar '01

11 Mar '01

Hi, Doing some tests with square-free factorization, I believe I have found a bug. Let e by (x+y)^2*(x-y) expanded and call sqrfree(e,x). You do not get the original expression, it remains expanded. If y is replaced by 1 it seems it works. Bernard

3 3

Problem with CLN on Alpha
by maxa＠frodo.jia.czn.cz 13 Feb '01

13 Feb '01

Dear all! I am trying to compile CLN library on an Alpha machine, but I am getting a lot of undefined symbols (see below). Can anybody help? Thanks. Best regards, Ladislav Zejda zejdaz(a)mbox.vol.cz Problem: after running "make" I am getting: ... /bin/sh ../libtool --mode=link g++ -O6 -g -fno-exceptions exam.o exam_I.o exam_RA.o exam_SF.o exam_FF.o exam_DF.o exam_LF.o exam_I_gcd.o exam_I_sqrtp.o ../src/libcln.la -lm -o exam mkdir .libs g++ -O6 -g -fno-exceptions exam.o exam_I.o exam_RA.o exam_SF.o exam_FF.o exam_DF.o exam_LF.o exam_I_gcd.o exam_I_sqrtp.o ../src/.libs/libcln.so -lgmp -lm -o .libs/exam -Wl,--rpath -Wl,/usr/local/lib exam.o: In function `global constructors keyed to main': /usr/src/cln/cln-1.1/tests/../include/cln/io.h:91: undefined reference to `_GLOBAL_$I$cl_module__cl_prin_globals__firstglobalfun' /usr/src/cln/cln-1.1/tests/../include/cln/io.h:91: undefined reference to `_GLOBAL_$I$cl_module__cl_prin_globals__firstglobalfun' exam.o: In function `global constructors keyed to main': /usr/src/cln/cln-1.1/tests/../include/cln/random.h:36: undefined reference to `_GLOBAL_$I$cl_module__cl_random_def__firstglobalfun' /usr/src/cln/cln-1.1/tests/../include/cln/random.h:36: undefined reference to `_GLOBAL_$I$cl_module__cl_random_def__firstglobalfun' exam.o: In function `global constructors keyed to main': /usr/src/cln/cln-1.1/tests/../include/cln/dfloat_class.h:55: undefined reference to `_GLOBAL_$I$cl_module__cl_DF_globals__firstglobalfun' /usr/src/cln/cln-1.1/tests/../include/cln/dfloat_class.h:55: undefined reference to `_GLOBAL_$I$cl_module__cl_DF_globals__firstglobalfun' exam.o: In function `global constructors keyed to main': ... and so on My setup (Isn't it too old?): cat /proc/version Linux version 2.2.11 (root(a)zelva.jia.czn.cz) (gcc version egcs-2.91.60 19981201 (egcs-1.1.1 release)) #46 Mon Aug 2 13:19:49 PDT 1999 cat /proc/cpuinfo cpu : Alpha cpu model : EV56 cpu variation : 0 cpu revision : 0 cpu serial number : Linux_is_Great! system type : EB164 system variation : SX164 system revision : 0 system serial number : MILO-0000 cycle frequency [Hz] : 533171648 timer frequency [Hz] : 1024.00 page size [bytes] : 8192 phys. address bits : 40 max. addr. space # : 127 BogoMIPS : 528.48 kernel unaligned acc : 7 (pc=fffffc0000449014,va=fffffc000029ef32) user unaligned acc : 2007464 (pc=120185630,va=120ca2dac) platform string : N/A g++ -v Reading specs from /usr/local/lib/gcc-lib/alphapca56-unknown-linux-gnu/egcs-2.91.60/specs gcc version egcs-2.91.60 19981201 (egcs-1.1.1 release)

2 7

Re: GiNaC/ginac power.cpp
by Richard B. Kreckel 30 Jan '01

30 Jan '01

On Mon, 29 Jan 2001, Christian Bauer wrote: > Update of /home/cvs/GiNaC/ginac > > Modified Files: > power.cpp > Log Message: > expand() always expands the exponent and transforms x^(a+b) -> x^a*x^b Question, Sir: Now, that we have expand(x^((a+b)^2-a^2-b^2-2*a*b)) -> 1, wouldn't it be natural to also have normal(x^((x^2+2*x*y+y^2)/(x+y)^2)) -> x? Or does this contradict the "first make a polynomial out of it" principle? -richy. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

2 1

Install Problems
by Peter S. Christopher 28 Jan '01

28 Jan '01

Hi there, I hope that you folks aren't too busy to help me out. I am trying to compile GiNaC on a dual pentium III machine. The result of g++ -v is : gcc version 2.96 20000731 (Red Hat Linux 7.0) The part where the compiler seems to be troubled is at this stage: c++ -DHAVE_CONFIG_H -I. -I. -I.. -I/usr/local/include -g -O2 -c inifcns.cpp -fPIC -DPIC -o .libs/inifcns.lo The message that the machine keeps giving me is: In file included from inifcns.cpp:31: matrix.h:144: non-local function `GiNaC::ex GiNaC::determinant (const GiNaC::matrix &, unsigned int = GiNaC::determinant_algo::automatic)' uses anonymous type If there is more information I can provide please let me know. Any information would be much appreciated. Thank you in advance Pete Christopher Univ. of Chicago Chicago IL 60615

2 1

GiNaC with GUI
by Stefan Weinzierl 26 Jan '01

26 Jan '01

Dear all, we wrote a small program called gTybalt, which combines GiNaC, Root and TeXmacs. This yields a graphical user interface to GiNaC. Since the Root-package includes CINT, you basically have the same functionality as ginac-cint. In addition, gTybalt provides routines for plotting functions of one or two variables. The plotting routines are based on Root and are quite fancy. Furthermore by integrating TeXmacs, you get high-quality output, e.g. mathematical formulae are displayed using TeX fonts. TeXmacs is a what-you-see-is-what-you-get editor and the resulting formulae can easily be converted to LaTeX. At the present stage, the very first alpha-release of gTybalt is available. (This means: Don't expect too much !) You are welcome to visit the web-site http://www.fis.unipr.it/~stefanw/gtybalt.html for more information and the source code. Please note that at the moment gTybalt is based on OLD VERSIONS of GiNaC and CLN. You need GiNaC 0.6.4 and CLN 1.0.3. Best wishes, Stefan Weinzierl and Roberta Marani

2 2

GiNaC-0.7.0
by Eugene Gorshkov 12 Jan '01

12 Jan '01

Hi GiNaC's authors! I'm using GiNaC in project for simplifying and normalizing math expressions (and physics, hereafter). My program checks student's solution of math or physics problem without teacher help (teacher need to enter some variants of solution only once). This program will be used for Distant Learning. So, I need to compare two different expressions on identity. I've found out some bugs (or features?) in new version of GiNaC 0.7.0 There they are: When I'm entering command like this: > normal(a^b*a^c, 2); in ginsh or simple c++ code, GiNaC stops to work or exits by Segmentation fault (there is no output). If normalization level == 1, the result is: a^b*a^c. But I need a^(b+c)! Could you explain to me what must I do to correct this? And I have some questions about GiNaC. 1. If I need to normalize expression (a+b)^2 and result must be a^2+b^2+2*a*b. Have I need to call only normal() method or expand() also? 2. Do you planning to include some simple operations with vectors in GiNaC? Also, I'd like to help you in extending GiNaC. Therefore I'll examine GiNaC's sources and try to write some additional code... --------------------------------------- Machine configuration: Red Hat Linux release 7.0 (Guinness) Kernel 2.2.16-22 on an i586 Linux host.domain 2.2.16-22 #1 Tue Aug 22 16:16:55 EDT 2000 i586 unknown gcc version: 2.96 glibc-2.1.92-14 libc-5.3.12-31 GiNaC-0.7.0 cln-1.1 -- With best regards, Eugene Gorshkov aka Atany .. e-mail: atany(a)ic.kharkov.ua --- XEmacs 21.1 "Channel Islands" (linux)

2 1

bug
by Ralf Hoch 04 Jan '01

04 Jan '01

Hello ginac, I use gcc version 2.96 20000731 (Red Hat Linux 7.0) on Pent.III gcc version egcs-2.91.66 19990314/Linux (egcs-1.1.2 release) Whenever I want to compile ginac I got the following error: In file included from /usr/local/include/ginac/ginac.h:39, /usr/local/include/ginac/matrix.h:144: non-local function `GiNaC::ex GiNaC::determinant (const GiNaC::matrix &, unsigned int = GiNaC::determinant_algo::automatic)' uses anonymous type -- With best regards, Ralf Hoch University of Kassel Phone: (+49)-5542-98-1609 Dept. of Agriculture and Fax: (+49)-5542-98-1556 Environmental Protection (FB 11) email: hoch(a)wiz.uni-kassel.de D-37213 Witzenhausen, (FRG)

3 2