- GiNaC-list - ginac.de

build problems with ver.1.3.3
by Joerg Arndt 06 Dec '05

06 Dec '05

attached is the log of a failed build of ginac 1.3.3 gcc used is: gcc-4.1.0_20051123-3 gcc-c++-4.1.0_20051123-3 The log is virtually identical for i386, ia64, ppc, ppc64, s390, s390x, x86_64. As the problems look rather generic I'd recommend these errors and warnings be fixed in the original package (instead of producing another SuSE patch). [resent because uncompressed log is >40KB and rejected] attached is the bzip2 compressed log -- p=2^q-1 prime <== q>2, cosh(2^(q-2)*log(2+sqrt(3)))%p=0 Life is hard and then you die.

3 12

bug in reposition_dummy_indices
by varg＠theor.jinr.ru 22 Nov '05

22 Nov '05

Hello! This simple program fails: #include <iostream> #include <cassert> #include <ginac/ginac.h> using namespace std; using namespace GiNaC; int main(int argc, char** argv) { varidx mu1(symbol("mu1"), 4); varidx mu2(symbol("mu2"), 4); varidx mu3(symbol("mu3"), 4); varidx mu4(symbol("mu4"), 4); varidx mu5(symbol("mu5"), 4); varidx mu6(symbol("mu6"), 4); symbol A("A"); exvector ev2; ev2.push_back(mu3.toggle_variance()); ev2.push_back(mu6); ev2.push_back(mu5.toggle_variance()); ev2.push_back(mu6.toggle_variance()); ev2.push_back(mu5); ev2.push_back(mu3); ex test_cycl = indexed(A, sy_cycl(), ev2); test_cycl = test_cycl.simplify_indexed(); assert(test_cycl.get_free_indices().size() == 0); return 0; } The reason is that reposition_dummy_indices() (indexed.cpp:639) ignores symmetry properties of expression[s]. To be more specific, if the expression e has some symmetry properties, than the code like this (indexed.cpp:654): e = e.subs(lst( *it2 == ex_to<varidx>(*it2).toggle_variance(), ex_to<varidx>(*it2).toggle_variance() == *it2 ), subs_options::no_pattern); can yield some non-trivial permutation of indices (apart from swapping covariant and contravariant instances of some index), so subsequent manipulations become incorrect. However, I don't know how to fix this. Any ideas? -- All science is either physics or stamp collecting.

1 0

GiNaC and Geometric Algebra
by Alan Bromborsky 18 Nov '05

18 Nov '05

Considering some of the applications you are developing GiNaC for I wonder if you familiar with the book "Geometric Algebra for Physicists," by Chis Doran and Anthony Lasenby, especially with regard to the geometric algebra formulation of tensor analysis for gravitation and quantum mechanics? The geometric algebra/calculus formulation leads to simplifications and unification of several areas of mathematical physics. Question, in your representation of Clifford algebra do you use a multiplication table or a matrix representation?

2 1

Clifford Metric Array
by Alan Bromborsky 17 Nov '05

17 Nov '05

You probably handle this correctly internally, but the documentation example is not consistent. In section 4.15.1.2, "A Generic Clifford Algebra", you state that e[i]*e[j]+e[j]*e[i] = M(i,j) and then give an example where ex M = diag_matrix(lst(1,-1,0,s)); ex e = clifford_unit(nu,M); and say M(0,0) = 1 -> e[0]^2 = 1. I get 2*e[0]^2 = 1 so that e[0]^2 = 1/2. Is M internally multiplied by 2?

2 1

small portability fix
by varg＠theor.jinr.ru 08 Nov '05

08 Nov '05

Hello! The attached patch removes non-portable code from ginsh and checks, so it possible to compile ginsh and run checks on MinGW. -- All science is either physics or stamp collecting.

2 6

ginac 1.3.3 fails make check under Mac OS X 10.4.x
by Marius Schamschula 08 Nov '05

08 Nov '05

Hello all, I've been running into a problem building ginac under Mac OS X 10.4.x (Tiger)/Darwin 8.x. I can compile w/o any problems, but when I run make check I get the following error: /bin/sh ../libtool --tag=CXX --mode=link g++ -g -O2 -o exams exam_paranoia.o exam_numeric.o exam_powerlaws.o exam_inifcns.o exam_inifcns_nstdsums.o exam_differentiation.o exam_polygcd.o exam_normalization.o exam_pseries.o exam_matrices.o exam_lsolve.o exam_indexed.o exam_color.o exam_clifford.o exam_archive.o exam_structure.o exam_hashmap.o exam_misc.o exams.o ../ginac/ libginac.la -L/usr/local/lib -lcln -lgmp g++ -g -O2 -o .libs/exams exam_paranoia.o exam_numeric.o exam_powerlaws.o exam_inifcns.o exam_inifcns_nstdsums.o exam_differentiation.o exam_polygcd.o exam_normalization.o exam_pseries.o exam_matrices.o exam_lsolve.o exam_indexed.o exam_color.o exam_clifford.o exam_archive.o exam_structure.o exam_hashmap.o exam_misc.o exams.o -Wl,-bind_at_load ../ginac/.libs/ libginac.dylib -L/usr/local/lib /usr/local/lib/libcln.dylib /usr/ local/lib/libgmp.dylib /usr/bin/ld: Undefined symbols: GiNaC::structure<sprod_s, GiNaC::compare_std_less>::reg_info collect2: ld returned 1 exit status make[3]: *** [exams] Error 1 make[2]: *** [check-am] Error 2 make[1]: *** [check-recursive] Error 1 make: *** [check] Error 2 Note: I don't see this problem under Mac OS 10.3.x (Panther)/Darwin 7.x. I tried building using both gcc 3.3 and gcc 4.0. gcc 3.3 dies with some linking errors. Any ideas? TIA Marius -- Marius Schamschula, Alabama A & M University, Department of Physics The Center for Hydrology Soil Climatology and Remote Sensing http://optics.physics.aamu.edu/ - http://www.physics.aamu.edu/ http://wx.aamu.edu/ - http://www.aamu.edu/hscars/

2 4

bug(?) in clifford_moebius_map
by varg＠theor.jinr.ru 31 Oct '05

31 Oct '05

Hello! I've noticed that example code from cs.MS/0410044 ("An Example of Clifford Algebras Calculations with GiNaC") does not work. $ g++ -O0 -g -o parabolic parabolic1-bin.cpp `ginac-config --libs` $ ./parabolic Metric is: e. *** Got problem in vector fields: Expression is a non-scalar Clifford number! *** Got problem in vector fields: Expression is a non-scalar Clifford number! *** Got problem in vector fields: Expression is a non-scalar Clifford number! Vect field Direct In Cayley In Cayley1 Segmentation fault (core dumped) $ ginac-config --version 1.3.3 $ g++ --version g++ (GCC) 4.0.3 20051023 (prerelease) (Debian 4.0.2-3) Copyright (C) 2005 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. I think the reason is that clifford_moebius_map mixes matrices and indexed expressions. $ gdb --silent ./parabolic [some boring stuff snipped] (gdb) break clifford.cpp:1023 (gdb) bt #0 GiNaC::clifford_moebius_map (a=@0x8088548, b=@0x808854c, c=@0x8088550, d=@0x8088554, v=@0x5ca926d4, G=@0x5ca926d0, rl=0 '\0') at clifford.cpp:1203 #1 0x2e89292e in GiNaC::clifford_moebius_map (M=@0x5ca9205c, v=@0x5ca926d4, G=@0x5ca926d0, rl=0 '\0') at clifford.cpp:1210 #2 0x2e8929fb in GiNaC::clifford_moebius_map (M=@0x5ca9205c, v=@0x5ca926d4, G=@0x5ca926d0) at clifford.cpp:1217 #3 0x08052acf in main (argc=1, argv=0x5ca94034) at parabolic1-bin.cpp:604 (gdb) frame 0 #0 GiNaC::clifford_moebius_map (a=@0x8088548, b=@0x808854c, c=@0x8088550, d=@0x8088554, v=@0x5ca926d4, G=@0x5ca926d0, rl=0 '\0') at clifford.cpp:1203 1203 return clifford_to_lst(e, cu, false); (gdb) call e.dbgprint() e~symbol13*[[x],[y]].symbol13*exp(t)*exp(-t)^(-1) (gdb) Due to this meaningless combination of matrix and indexed expression clifford_to_lst fails, and since the example code is not really bullet-proof, the end result is a segfault. Any ideas? Best regards, Alexei -- All science is either physics or stamp collecting.

2 3

Unsubscribe notifications
by Christian Bauer 25 Oct '05

25 Oct '05

Hi! If you've been getting "You have been unsubscribed from ginac-list@..." messages, these are from the old, defunct lists at thep.physik.uni-mainz.de. I noticed they were still active and, lacking the means to properly delete them, unsubscribed everyone (and configured them to discard all incoming messages). I guess I should have disabled the "Send good-bye message" option first... If you are reading this, you are still subscribed to the real GiNaC mailing lists at ginac.de, and you do not have to take any action. I'm sorry if I caused any confusion. Bye, Christian -- / Physics is an algorithm \/ www.cebix.net

1 0

GiNaC 1.3.3 release
by Richard B. Kreckel 24 Oct '05

24 Oct '05

Hi, GiNaC 1.3.3 is out and available. The changes are: - Fixed bug occuring in algebraic substitutions with expressions involving indices. - Fixed bug that selected the wrong algorithm for determinants of purely numeric matrices with at least one floating point element. - Fixed bug in function H that caused an infinite recursion for arguments around +-I. - Fixed const-correctness in printing handler for GCC 4.0.2. - Made lookup in adaptivesimpson precision-aware. - Added series expansion for functions (classical) Li and S around x==0. - Added fsolve() numerical univariate real-valued function solver. - Added functions sub_matrix() and reduced_matrix(). - Small cleanups. Less warnings with latest GCC. As always, this release can be downloaded from http://www.ginac.de/Download.html Enjoy, -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

1 0

ViewCVS frontend
by Christian Bauer 15 Oct '05

15 Oct '05

Hi! The GiNaC CVS can also be browsed at http://www.ginac.de/viewcvs/GiNaC/ Bye, Christian -- / Physics is an algorithm \/ www.cebix.net

1 0