- GiNaC-devel - ginac.de

NTL goes GPL
by Richard B. Kreckel 10 Apr '00

10 Apr '00

Hello all, For those interested in factorization over finite fields, Victor Shoup's library NTL (which is used as a module for some "world records" established by MuPAD (if you can call it "established" (MuPAD as a mere frontend))) has just been released by the author under the GPL license. Of course, it is still a long way from there to factorization of multivariate polynomials over ZZ but it may be start. Regards -rbk. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 1

eps bug
by Dr K.M. Briggs 03 Apr '00

03 Apr '00

Dear ginac-bugs(a)ginac.de, Installing GiNaC-0.5.4 on my redhat 6.0 system fails at the doc step: make[3]: Entering directory `/home/kbriggs/GiNaC-0.5.4/doc/tutorial' Running /usr/X11R6/bin/fig2dev -L eps -m 0.9 classhierarchy.fig classhierarchy.eps... Unknown graphics language eps Known languages are: acad box epic eepic eepicemu gif ibmgl jpeg latex mf pcx pic pictex png ppm ps pstex pstex_t textyl tiff tk tpic xbm xpm I upgraded my transfig package to the latest from redhat 6.1 (Release 3.2 Patchlevel 1), but it has the same problem. Actually, changing the command to fig2dev -L ps -m 0.9 classhierarchy.fig classhierarchy.eps gives a perfectly good tutorial.ps. Maybe the configure script should check whether the fig2dev version installed knows about eps? Best wishes, Keith -- +-------------------------------------------------------------------+ | Dr. Keith M. Briggs, Complexity Research Group, BT Research Labs. | | Snail mail: Adastral Park adm2 pp5, Martlesham Heath, IP5 3RE, UK | | Tel. 01473 641 911, Fax. 01473 647 410. Home: 10 Saddler's Place, | | Martlesham Heath, Suffolk IP5 3SS. Home telephone: 01473 625 972 | | NEW! personal homepage: http://www.labs.bt.com/people/briggsk2/ | +-------------------------------------------------------------------+

2 1

MT safety
by Richard B. Kreckel 28 Mar '00

28 Mar '00

Since people here have been pondering about this issue lately: the design document <http://www.sgi.com/Technology/STL/thread_safety.html> applies to GCC's templates as well. IIUC this looks rather sound, the worst thing thus being our refcounting. How do I do a portable spin-lock? Regards -rbk. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 0

The Gamma function.
by Richard B. Kreckel 23 Mar '00

23 Mar '00

Hi folks, I have just checked with Martin v. Löwis about that issue. Martin is quite C-standard-literal. The question arouse how to call the (already present) Gamma function in GiNaC once and for all. The standard ISO/IEC 9899:1999, frequently dubbed as C99, states: > 7.12.8.3 The lgamma functions > Synopsis > #include <math.h> > double lgamma(double x); > float lgammaf(float x); > long double lgammal(long double x); > > Description > > The lgamma functions compute the natural logarithm of the absolute > value of gamma of x: log-e | G(x) |. A range error occurs if x is too > large or if x is a negative. > > Returns > > The lgamma functions return the value of the natural logarithm of the > absolute value of gamma of x. > > 7.12.8.4 The tgamma functions > > Synopsis > #include <math.h> > double tgamma(double x); > float tgammaf(float x); > long double tgammal(long double x); > > Description > > The tgamma functions compute the gamma function of x: G(x). A domain > error occurs if x is a negative integer or zero. A range error may > occur if the magnitude of x is too large or too small. > > Returns > > The tgamma functions return the gamma function value. The next revision of ISO/IEC 14882:1998 (the C++ standard) will follow this definition from the C standard. Since it is kind of a policy not to work against the language C++ I vote for forgetting about esthetic qualms and adopting this behaviour and call it ex tgamma(ex) in GiNaC 0.x.y | x>5. A function ex lgamma(ex) may be introduced as well to make the distinction clear. Any serious objections? -richy. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 0

JaRoG
by Richard B. Kreckel 15 Mar '00

15 Mar '00

...means "Just another Release of GiNaC". Version 0.5.4 has just been uploaded to the usual places. See the file NEWS to learn about what has changed. Enjoy! -rbk. PS: No, I didn't forget to tag... PPS: ...and no, config.h.in is not empty. ;-) -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 0

evalf() and collect()...
by Richard B. Kreckel 02 Mar '00

02 Mar '00

...don't commute. But there is more: squark:~$ ginsh ginsh - GiNaC Interactive Shell (GiNaC V0.5.3) __, _______ Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, (__) * | Germany. This is free software with ABSOLUTELY NO WARRANTY. ._) i N a C | You are welcome to redistribute it under certain conditions; <-------------' see the file COPYING for details. Type ?? for a list of help topics. > series(gamma(x),x,0,2); -EulerGamma+Order(x^2)+(1/12*Pi^2+1/2*EulerGamma^2)*x+x^(-1) > collect(",x); -EulerGamma+Order(x^2)+(1/12*Pi^2+1/2*EulerGamma^2)*x+x^(-1) > evalf("); -0.5772156649015328606L0+(0.98905599532797255544L0)*x+Order(x^(2.0L0))+x^(-1.0L0) > quit; squark:~$ ginsh ginsh - GiNaC Interactive Shell (GiNaC V0.5.3) __, _______ Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, (__) * | Germany. This is free software with ABSOLUTELY NO WARRANTY. ._) i N a C | You are welcome to redistribute it under certain conditions; <-------------' see the file COPYING for details. Type ?? for a list of help topics. > series(gamma(x),x,0,2); -EulerGamma+Order(x^2)+(1/12*Pi^2+1/2*EulerGamma^2)*x+x^(-1) > evalf("); -0.5772156649015328606L0+(0.98905599532797255544L0)*x+Order(x^(2.0L0))+x^(-1.0L0) > collect(",x); Not a 32-bit integer: -9223372028264841217 squark:~$ In any case, exponents should remain integer and even rational. Cheers -rbk. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

2 1

license
by Daniel Webb 01 Mar '00

01 Mar '00

You should have a link from the main page stating the license of your software. I could not find it. ____________________________________________________________________ Get free email and a permanent address at http://www.netaddress.com/?N=1

1 0

Version 0.5.3 on the wire
by Richard B. Kreckel 24 Feb '00

24 Feb '00

i.e. on <ftp://ftpthep.physik.uni-mainz.de/pub/GiNaC/>. There has been a bug in the last version where .normal() was returning wrong results for certain rational functions like 1/(5-4*a)*(-5-2*a-(3+2*a)*(2-5*a)/(-2+a)) which is now fixed. Also, this version features an improved and more flexible scheme for registering functions which makes extensions easier and also allows writing remembering functions in a clean way. See <http://www.ginac.de/reference/class_remember_table.html>. GiNaC-cint is making slow progress---we have it running with the current cint 5.14.31 now. Sex and drugs and Unix -rbk. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

1 0

Build problems
by Patrick K Notz 22 Feb '00

22 Feb '00

Hi, I'me trying to build ginac and am having problems getting 'configure' to find the readline and cln libraries. I've installed them in directories (non-standard since I'm not root) under my account. I just can't figure out how to tell configure/make to look there. I've tried the suggestion in Section 3.2 of the tutorial but the syntax given doesn't work. I've tried adding "LDFLAGS = $(LDFLAGS) /my/lib/dir" to all of the Makefile.in files but no luck again. I'm not familiar with autoconf so I'm not sure which files I should try tweeking. Can you offer any suggestions? Thanks ~ pat <notz(a)purdue.edu>

2 1

Today's bugs
by Richard B. Kreckel 21 Feb '00

21 Feb '00

Consider this: higgsino:~/projects/GiNaC/ginsh$ ./ginsh ginsh - GiNaC Interactive Shell (GiNaC V0.5.2) __, _______ Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, (__) * | Germany. This is free software with ABSOLUTELY NO WARRANTY. ._) i N a C | You are welcome to redistribute it under certain conditions; <-------------' see the file COPYING for details. Type ?? for a list of help topics. > series(gamma(x),x,0,2); collect(series(gamma(x),x,0,2),x); -EulerGamma+x*(1/12*Pi^2+1/2*EulerGamma^2)+x^(-1)+Order(x^2) -EulerGamma+x*(1/12*Pi^2+1/2*EulerGamma^2)+x^2*Order(1)+x^(-1) The other bug is with the fraction free elimination in ginac/matrix.cpp: If I try to solve the system [[ [[0,2-a]], [[5-4*a,3+2*a]] ]] * [[ [[x]], [[y]] ]] == [[ [[-2+5*a]], [[-5-2*a]] ]] I run into trouble within assert-blocks. I will have a look at the second problem. -rbk. -- Richard Kreckel <Richard.Kreckel(a)Uni-Mainz.DE> <http://wwwthep.physik.uni-mainz.de/~kreckel/>

3 3