April 2006 - GiNaC-list

[patch] fix template template declarations [GCC 4.2 compilation fix]
by varg＠theor.jinr.ru 11 May '06

11 May '06

Hello! The attached patch corrects template template declarations, so GiNaC can be compiled with G++ 4.2. Without this patch compilation fails, see http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=362220 for details. Note that `make check' still fails with G++ 4.2. Basically, ex e("-1 + 2*I", symbol()); triggers an exception. My patch *does not* fix this error. Best regards, Alexei. -- All science is either physics or stamp collecting.

2 4

Announcement: SyFi 0.3.0
by Kent-Andre Mardal 19 Apr '06

19 Apr '06

Dear GiNaC users, I would like to announce the 0.3.0 version of SyFi (syfi.sf.net) which is a finite element library based on GiNaC. At present, we have implemented several finite elements, like e.g., Lagrangian, Raviart-Thomas, and Nedelec elements of arbitrary order in both 2D and 3D. The latest release adds several examples with code generation which demonstrate the use of libraries like e.g. the Dolfin (www.fenics.org) PyCC (heim.ifi.uio.no/~skavhaug/heart_simulations.html) and Trilinos (software.sandia.gov/trilinos) The name SyFi stands for Symbolic Finite Elements. The package provides polygonal domains, polynomial spaces, and degrees of freedom as symbolic expressions that are easily manipulated. This makes it easy to define finite elements. SyFi is open source. Thanks for GiNaC!! Regards, Kent-Andre Mardal Simula Research Laboratory

1 0

GiNaC 1.3.4 released
by Jens Vollinga 13 Apr '06

13 Apr '06

Hi everybody, GiNaC 1.3.4 is out and available. The changes are: - More general index dimensions are now possible. - Improved algebraic substitutions. - Fixed wrong matching in .has(). - Fixed bug in differentiation of a power. - collect_common_factors also works with negative powers of common factors now. - Code clean-up and preparations for upcoming cln 1.2. As always, this release can be downloaded from http://www.ginac.de/Download.html Best wishes, Jens

1 0

Rot, div, grad operator will be implemented in feelfem FEM code generator
by Fujio Hidehiro 10 Apr '06

10 Apr '06

Dear GiNaC developers, This is Hidehiro Fujio of NEC Corp. Japan. At first, I would like to say thank you very much for your significant contribution to the formula manipulation area by GiNaC library in the world, and I also have to say that thank you very much for your link to my small document in the Related Software page. (My system is feelfem, an FEM code generator.) Yesterday, I met my supervisor, and got his advice. I wll expand my feelfem system to deal with 3D Maxwell equation, and introduce the differential operator 'rot' and 'div' and 'grad' (rotation, divergence, gradient, respectively), and my professors suggestion is it may be possible to expand these differential operators into just a set of basic differential operators (like dx,dy,dz) inside GiNaC. I would like to investigate your documentation agaian, and I would like to make a reort, how to make an FEM code generator by GiNaC. Again, thank you very much for GiNaC. And I'm afraid it is not an appropriate mail............ Best regards, Hidehiro Fujio # mail to both ( h-fujio(a)ab.jp.nec.com , oimo(a)world.email.ne.jp ,if needed)

2 1

[patch] small improvement of collect_common_factors
by varg＠theor.jinr.ru 04 Apr '06

04 Apr '06

Hello! This ginsh code snippet: collect_common_factors(x/y + z/y^2); returns the original expression unchanged: z*y^(-2)+x*y^(-1) The patch below makes it return (x+z*y^(-1))*y^(-1) Index: normal.cpp =================================================================== RCS file: /home/cvs/GiNaC/ginac/normal.cpp,v retrieving revision 1.94.2.8 diff -u -r1.94.2.8 normal.cpp --- normal.cpp 30 Jun 2005 14:03:16 -0000 1.94.2.8 +++ normal.cpp 10 Mar 2006 13:40:26 -0000 @@ -2388,6 +2388,8 @@ { if (exponent.info(info_flags::posint)) return power(basis.to_rational(repl), exponent); + else if (exponent.info(info_flags::negint)) + return power(replace_with_symbol(power(basis, _ex_1), repl), -exponent); else return replace_with_symbol(*this, repl); } Best regards, Alexei. -- All science is either physics or stamp collecting.

2 1