- GiNaC-list - ginac.de

GiNaC 1.1.1 released
by Christian Bauer 18 Jun '03

18 Jun '03

Hi! With one last-minute bug-fix sneaked in by the ever-helpful Chris Dams, we'd like to present GiNaC 1.1.1 which is crunchy on the outside but soft and delicious on the inside: - lst (and exprseq) provide iterators for read-only element access. For sequential access this is one order faster than using op(). - Implemented relational::subs() (this was done in 1.0.9 but inadvertently omitted from the 1.1 branch). - pole_error and do_taylor are available to library users. - Added on-line help and Tab-completion for print(), iprint(), print_latex() and print_csrc() in ginsh. Get it from your favorite FTP server: ftp://ftpthep.physik.uni-mainz.de/pub/GiNaC There's nobody here to drink champagne with, so I can save to money for buying more C++ books... Happy computing, Christian -- / Physics is an algorithm \/ http://www.uni-mainz.de/~bauec002/

1 0

Atoms of an expression
by Dr. Vassilis S. Vassiliadis 10 Jun '03

10 Jun '03

How is it possible to get GiNaC to yield all the atoms occuring in an expression? (i.e. return a list of all the leaves of the expression tree). This is of interest in the generation of sparse jacobian / hessian matrices as one example of applications. Vassilis ----------------------------------------------------------------- Dr. Vassilis S. Vassiliadis, Lecturer, UNIVERSITY OF CAMBRIDGE, DEPARTMENT OF CHEMICAL ENGINEERING Pembroke Street, Cambridge CB2 3RA, UK. e-mail: vsv20(a)cheng.cam.ac.uk -----------------------------------------------------------------

2 1

simplify_indexed() didn't fully simplify the expressions
by Yong Xiao 09 Jun '03

09 Jun '03

Hi! I am using GiNaC1.1.0 It seems that for some expression(suppose it's defined as e1) simplify_indexed() doesn't do simplification fully. However, if I define another ex objec tmp = e1, and e1.simplify_indexed() will give a more simple expression which, although, is still not of the simplest form. Anyone knows a more efficient way? Another problem I met is that when I simplify an expression which is composed of 240 terms, where every term is a simple product of indexed objects, the complilation will stop, saying "out of memory". I don't think that GiNaC can't handle such expressions. Anything wrong?

2 3

cln and compilation!
by Jean-Sebastien Gagne 28 Apr '03

28 Apr '03

Hello, I installed cln in private a dir and am trying to compile some basic c++ code using cl_I class. The problems is I am not able to get an executable out my code. The installation seems to have worked ok. I did not have any error and I can execute the examples. What woul be the correct command to get exe from my code? I known the -I and -l flags are involved but cannot get the right combo. I am a little new at this! Thanks! Jean-Sebastien Gagne

2 1

GiNaC 1.1.0 released
by Christian Bauer 04 Apr '03

04 Apr '03

Lo! Receive now, mylords, version 1.1 of GiNaC, as a gift from one who has nothing, and as a memorial of macros and print_contexts. For now our fate is fulfilled, and-- Oops, I think I got the wrong script here. Anyway... Today, we'd like to present GiNaC 1.1.0. This is the first release since 1.0.0 to introduce incompatible API changes (hence the "1" after the first dot), but unless you are extending GiNaC with your own algebraic classes you probably won't notice (be sure, though, to check the change log for a list of things to watch out for). So rejoice, noble ladies and lords, and behold its mighty features: - Added stream manipulators dflt, latex, python, python_repr, tree, csrc, csrc_float, csrc_double, csrc_cl_N, index_dimensions and no_index_dimensions to control the output format. [That's, like, the coolest thing, if you ask me. Lets you, like, say "cout << latex", and then anything you print will be, like, in LaTeX code. Totally smooth, man...] - subs() accepts an optional "options" argument. The option subs_option::subs_algebraic enables "smart" substitutions in products and powers. [Brought to you by Chris Dams (the real one, not his substitute).] - The scalar_products mechanism allows the specification of an index dimension. [What could be beyond height and width, depth and time?] - Added ex::to_polynomial(). [ex::to_rational()'s evil twin brother.] - Removed deprecated macros is_ex_a, is_ex_exactly_a and friends for good. [Begone, macros! Begone!] - Removed dirac_gamma6/7(). [That's a feature?] - Made the hashing more simple and efficient. [Or so Richy claims...] 'Tis time, methinks, to join the feast: ftp://ftpthep.physik.uni-mainz.de/pub/GiNaC Good night, Christian -- / Coding on PowerPC and proud of it \/ http://www.uni-mainz.de/~bauec002/

1 0

GiNaC 1.1 API freeze
by Christian Bauer 19 Mar '03

19 Mar '03

Hi! We plan to release GiNaC 1.1.0 by the end of next week. It will become the new "stable" branch, and any new and exciting experimental stuff will then go into 1.2. Further development of 1.0.x will be limited to bugfixes. The idea is to keep the programming interfaces compatible across all versions of 1.1.x, as it has been the case for 1.0.x, so should there be anyone who has cause why the current 1.1 CVS should not become the 1.1 release version (suggestions for parameter/return type changes, useful additional methods etc.), they must speak now or forever hold their peace. Bye, Christian -- / Coding on PowerPC and proud of it \/ http://www.uni-mainz.de/~bauec002/

1 0

Problems with linking GiNaC
by paul81＠webb.de 17 Mar '03

17 Mar '03

Hello, I'm new to GiNaC. In fact I've just installed the Debian package. But when I tried to compile the examples given in the section "A Tour of GiNaC" in the info pages it won't work. The object file is correctly made; but the linking process just returns a lot of errors. I used the command g++ $(ginac-config --libs) example.c++ ld returns many error messages (too much to send them). Basicly there two types. First many messages like example.o(.text+0x47): In function `main': : undefined reference to `GiNaC::symbol::symbol[in-charge](std::basic_string<char, std::char_traits<char>, std::allocator< char> > const&)' And for (I think) every reference in the object file. And second there also messages like: /usr/lib/libginac.so: undefined reference to `cerr' /usr/lib/libginac.so: undefined reference to `cln::compare(cln::cl_R const &, cln::cl_R const &)' (Very much of them, too). I have checked ginac-config. It seems to report the right libraries. When I tried static linking some of the messages disappeared, but only a few. Can anyone tell me, how I can get those examples run? Thank in advance, Paul PS: My systems configuration: g++ (GCC) 3.2.3 20030309 (Debian prerelease) GiNaC 1.0.11-1 (the debian packages: libginac-dev 1.0.11-1 & libginac0 1.0.11-1) CLN 1.1.5-1 (Debian: libcln2 1.1.5-1) and Debian GNU/Linux 2.4.17

1 0

here we go...
by Richard B. Kreckel 01 Mar '03

01 Mar '03

This is to let folks know that Christian Bauer-san has just rolled the release of GiNaC-1.0.14 of which the only thing we know is that it works for him (tm). Can't ask him any more details. He just left the room. Hope he's fetching the champagne. Wait, I can do better and cite the NEWS file: * Improved the C-source output: complex numbers are printed correctly (using the STL complex<> template or cln::complex()), rational numbers use cl_RA() in the CLN output, and small integers are printed in a more compact format (e.g. "2.0" instead of "2.0000000e+00"). * function_options::set_return_type() and function_options::do_not_evalf_params() now actually work. Enjoy -richy. -- Richard B. Kreckel <Richard.Kreckel(a)GiNaC.DE> <http://www.ginac.de/~kreckel/>

1 0

matrices and expressions
by Johan Rutgeerts 27 Feb '03

27 Feb '03

Hi, I only just started using ginac and there's some things about matrices and their relation to expressions that I do not understand. As for instance: why does GiNaC::unity_matrix return a GiNaC::ex ? or why does A + B with A and B both a GiNaC::matrix return a GiNaC::ex? And how do I convert such a GiNaC::ex to a matrix? Can this be done? Greetings, Johan Rutgeerts

2 1

Differentiating matrix/vector expressions
by junk2＠microserf.org.uk 24 Feb '03

24 Feb '03

Hi I've had a trawl around the ginac site, and it's not obvious that ginac supports operations on matrix/vector expressions (particularly differentiation). Although it's clear that matrices and vectors are supported, can I perform a differentiation on a matrix or vector expression? Cheers

2 1