- GiNaC-devel - ginac.de

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-406-g9ada7a7
by git＠ginac.de 15 Dec '13

15 Dec '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 9ada7a7f5c47e512b7bf6057d4c013612be9a33b (commit) from e5c76f659e2e882da3d5dba60502d6851f782bf3 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 9ada7a7f5c47e512b7bf6057d4c013612be9a33b Author: Matthias Dellweg <dellweg(a)tp1.uni-duesseldorf.de> Date: Tue Dec 10 10:55:28 2013 +0200 expairseq::match(): no side effects if match failed. Fixes spurious match failures. > match(sin(y)*exp(b)+sin(x)*exp(a), sin($0)*exp(a)+exp(b)*sin($1)) FAIL The reason is that expairseq::match() might assign a wildcard even if the match fails. The first attempted submatch is sin(y)*exp(b) with sin($0)*exp(a). It fails (as it should) but $0 == y gets assigned as a side effect (which is wrong). Next submatch is sin(x)*exp(a) with sin($0)*exp(a) (the same pattern as in the first submatch). This one fails because of spurious $0 == y assignment. Due to the unpredicatable term ordering the sequence of submatches might be different and the match might succeed (as it should). This makes debugging a bit more funny. Signed-off-by: Matthias Dellweg <dellweg(a)tp1.uni-duesseldorf.de> [Alexei Sheplyakov: figure out the cause of the problem, make a test case] ----------------------------------------------------------------------- Summary of changes: check/match_bug.cpp | 44 ++++++++++++++++++++++++++++++++++++++++++++ ginac/expairseq.cpp | 26 +++++++++++++++++++++----- 2 files changed, 65 insertions(+), 5 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

Derivative of conjugated is conjugated of derivative.
by Vladimir V. Kisil 09 Nov '13

09 Nov '13

Dear All, Further to my previous patch on functions derivatives, this small patch adds the rule conjugate(f)'=conjugate(f') to the GiNaC. The change can be seen here: cout << conjugate(log(z)).diff(z) << endl; // (before) -> z^(-1)*D[0](conjugate)(log(z)) // (now) -> conjugate(z)^(-1) Best wishes, Vladimir -- Vladimir V. Kisil email: kisilv(a)maths.leeds.ac.uk www: http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Transformations http://www.worldscientific.com/worldscibooks/10.1142/p835

3 8

Roots of polynomial
by Rekha Penmetsa 08 Nov '13

08 Nov '13

Hi, I am using GiNaC in my project. I want to find roots of a univariate polynomial. I am deriving a polynomial from a matrix determinant, which is coming to, 6-1/15*r-19/3138418483200*r^4-23/20713561989120000*r^5+1/18900*r^2-71/7264857600*r^3 I am using RealAlgebraicNumberFactory::realRoots() method as in the below code to find roots. I have installed GiNaCRA as well and example programs are compiling and running fine. RationalUnivariatePolynomial p(det,r); list<RealAlgebraicNumberPtr> roots = RealAlgebraicNumberFactory::realRoots( p ); This call is returning a list of zero size. Am I missing anything ? how to solve for roots for polynomials of this type? Thanks Rekha

1 0

Even powers of modulus are rationalized
by Vladimir V. Kisil 05 Nov '13

05 Nov '13

Dear All, I believe that |z|^2= z*conjugate(z) is valid for any z. Also, the automatic substitution abs(z)^(2n) -> z^n * conjugate(z)^n shall make an expression more rational and, thus, more accessible to GiNaC. Therefore, I propose the attached alteration of abs_power(). Best wishes, Vladimir -- Vladimir V. Kisil email: kisilv(a)maths.leeds.ac.uk www: http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Transformations http://www.worldscientific.com/worldscibooks/10.1142/p835

1 0

Releasing?
by Richard B. Kreckel 04 Nov '13

04 Nov '13

Hi! Does anybody have plans for releasing? Doing releases may have become unpopular among developers since the current version can always easily be pulled off from git. On the other hand, we've had quite a number of improvements and build fixes since 1.6.2, so it might be a good opportunity to officially get all this released. What do others think? And, in particular, @Alex and @Vladimir: do you have things in the pipeline? Cheers -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

2 1

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-405-ge5c76f6
by git＠ginac.de 03 Nov '13

03 Nov '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via e5c76f659e2e882da3d5dba60502d6851f782bf3 (commit) via 5a3641098c88d3c6ea0765df65b801312ecfb91b (commit) from 5bf87cea66bb2071222c2910ed68c2649a98906c (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit e5c76f659e2e882da3d5dba60502d6851f782bf3 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Sun Nov 3 22:21:00 2013 +0200 Expansion rules for exp(), log(), and abs(). exp(a + b) -> exp(a)*exp(b) log(a*b) -> log(a) + log(b) abs(z*w) -> abs(z)*abs(w) log and exp are transformed only if expand_options::expand_transcendental is given. Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit 5a3641098c88d3c6ea0765df65b801312ecfb91b Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Sat Nov 2 21:06:53 2013 +0200 Introduce expand_options::expand_transcendental. In general expand() treats the argument as a rational expression (and leave transcendental functions as they are). However sometimes it's convenient to expand transcendental functions too, like log(a*b) = log(a) + log(b) exp(a + b) = exp(a)*exp(b) Applying such transformation by default doesn't seem to be a smart idea (think of log(p^2/mu^2) transformed to 2*log(p) - 2*log(mu)). Therefore introduce expand_options::expand_transcendental. As the name implies expand() tries to transform transcendental functions only if this flag is set. Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> ----------------------------------------------------------------------- Summary of changes: check/exam_inifcns.cpp | 65 ++++++++++++++++++++++++++++++++++++++++++ doc/tutorial/ginac.texi | 53 ++++++++++++++++++++++++++++++++++ ginac/flags.h | 3 +- ginac/inifcns.cpp | 22 ++++++++++++++ ginac/inifcns_trans.cpp | 72 +++++++++++++++++++++++++++++++++++++++++++++++ 5 files changed, 214 insertions(+), 1 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-403-g5bf87ce
by git＠ginac.de 03 Nov '13

03 Nov '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 5bf87cea66bb2071222c2910ed68c2649a98906c (commit) from 0eb0cb55c76ceab744847d6e66533e292774e1cb (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 5bf87cea66bb2071222c2910ed68c2649a98906c Author: Richard Kreckel <kreckel(a)ginac.de> Date: Sun Nov 3 12:49:07 2013 +0000 Keep CLang/libc++ happy. ----------------------------------------------------------------------- Summary of changes: ginac/inifcns_nstdsums.cpp | 28 ++++++++++++++-------------- 1 files changed, 14 insertions(+), 14 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

Re: [GiNaC-devel] Simplifying powers
by Vladimir V. Kisil 23 Oct '13

23 Oct '13

Dear Alexei, Many thanks for careful revision of proposed patches. >>>>> On Wed, 16 Oct 2013 08:26:54 +0300, Alexei Sheplyakov <alexei.sheplyakov(a)gmail.com> said: ASh> I'm working on making the main part (power::expand) reasonably ASh> efficient (in particular when the basis has no positive terms ASh> at all). Just let me know if I can be useful in any way. ASh> possible a custom expand() for functions. However I don't quite ASh> like changes to log and exp functions. I think expand() should ASh> treat the argument as a rational expression (and leave ASh> transcendental functions as they are). Transforming ASh> log(p^2/mu^2) into 2*log(p) - 2*log(mu) as a side effect of ASh> expand() looks like a bad idea. On the other hand, if someone want automatic simplification for log(p/mu)+log(mu), this expansion is the simplest way. ASh> default behavior. Therefore I propose to add an extra expand ASh> option, say, expand_options::expand_transcendental, and apply ASh> log_expand (and exp_expand) only if that option specified. I have implemented this suggestion in the attached patch. I am not sure if expansion of the absolute value shall depend from the same flag. It seems, that abs() did not break rationality you mentioned above. Best wishes, Vladimir -- Vladimir V. Kisil email: kisilv(a)maths.leeds.ac.uk www: http://www.maths.leeds.ac.uk/~kisilv/ Book: Geometry of Mobius Transformations http://www.worldscientific.com/worldscibooks/10.1142/p835

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-402-g0eb0cb5
by git＠ginac.de 16 Oct '13

16 Oct '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 0eb0cb55c76ceab744847d6e66533e292774e1cb (commit) via 798d53ebb4da4e8e3865ed7bd7f31412fe2be3a7 (commit) via 0358859aa6bc5806805fa415db79c039af3960d6 (commit) via 2c886a86f1d020e927d64ab2c79b20437a6bd561 (commit) via 4e1e7c6547dafc0b0f8695f9eced2591ae816a31 (commit) via d916416984a857e80962239a0ee93e7216f803bb (commit) via 10365850aa3803337bfea1fc201b81b6752096d4 (commit) via 31fef4a6bc14afa46a297b3b87f94204d6c726d5 (commit) via 50ed15a65e858858999be45301b3f49821cb1e71 (commit) from 52c44ffaee28a70b6df3280c3354ee4241473b14 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 0eb0cb55c76ceab744847d6e66533e292774e1cb Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Wed Sep 11 09:26:30 2013 +0100 abs() and conjugate() provide some info() now. commit 798d53ebb4da4e8e3865ed7bd7f31412fe2be3a7 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Wed Sep 11 09:26:30 2013 +0100 Make it possible to override info() for functions. Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit 0358859aa6bc5806805fa415db79c039af3960d6 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Thu Aug 1 14:37:45 2013 +0100 Power with a real base and even exponent reports that it's non-negative Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit 2c886a86f1d020e927d64ab2c79b20437a6bd561 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Thu Aug 1 14:37:45 2013 +0100 Improve abs_eval(): abs(conjugate(z)) => abs(z), ... ... abs(step(z)) => step(z) Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit 4e1e7c6547dafc0b0f8695f9eced2591ae816a31 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Mon Oct 7 08:04:30 2013 +0300 Functions can define a custom expand() now. Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit d916416984a857e80962239a0ee93e7216f803bb Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Wed Jul 31 09:48:24 2013 +0100 More evaluation rules: abs(x^n) => abs(x)^n (x > 0, n is real). A simple test is added as well. Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit 10365850aa3803337bfea1fc201b81b6752096d4 Author: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> Date: Tue Jul 30 16:40:38 2013 +0100 Transform abs(x)^n => x^n if x is real and n is even. While at it add some .hold() to returns which cannot be processed further. Signed-off-by: Vladimir V. Kisil <kisilv(a)maths.leeds.ac.uk> commit 31fef4a6bc14afa46a297b3b87f94204d6c726d5 Author: Alexei Sheplyakov <Alexei.Sheplyakov(a)gmail.com> Date: Sun Oct 13 13:16:16 2013 +0300 mul::info(): cache the result of the positiveness/negativeness checks, ... ... and use the cached result if possible. commit 50ed15a65e858858999be45301b3f49821cb1e71 Author: Alexei Sheplyakov <Alexei.Sheplyakov(a)gmail.com> Date: Sat Aug 3 17:01:42 2013 +0300 Introduce status_flags::{is_positive,is_negative,purely_indefinite} Quite a number of transformations can be applied if the expression is known to be positive (negative) (for instance, sqrt(x^2*y^2) = x*y if x >= 0 and y >= 0). However the check itself might be quite expensive. Hence a few status_flags are introduced to cache the result. ----------------------------------------------------------------------- Summary of changes: check/CMakeLists.txt | 1 + check/check_mul_info.cpp | 14 +++++++ check/exam_inifcns.cpp | 44 +++++++++++++++++++++- ginac/flags.h | 12 ++++- ginac/function.cppy | 69 ++++++++++++++++++++++++++++----- ginac/function.hppy | 11 +++++ ginac/function.py | 2 +- ginac/inifcns.cpp | 95 +++++++++++++++++++++++++++++++++++++++++++--- ginac/mul.cpp | 27 ++++++++++++- ginac/power.cpp | 2 + ginac/version.h | 4 +- 11 files changed, 256 insertions(+), 25 deletions(-) create mode 100644 check/check_mul_info.cpp hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0

[SCM] GiNaC -- a C++ library for symbolic computations branch, master, updated. release_1-4-0-393-g52c44ff
by git＠ginac.de 11 Oct '13

11 Oct '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "GiNaC -- a C++ library for symbolic computations". The branch, master has been updated via 52c44ffaee28a70b6df3280c3354ee4241473b14 (commit) from d266da8d47161d02ddd068f7f5b9acfabe92a192 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 52c44ffaee28a70b6df3280c3354ee4241473b14 Author: Richard Kreckel <kreckel(a)ginac.de> Date: Fri Oct 11 23:11:35 2013 +0200 Add subdir-objects to Automake options. This option is needed to keep automake 1.14 happy. ----------------------------------------------------------------------- Summary of changes: configure.ac | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) hooks/post-receive -- GiNaC -- a C++ library for symbolic computations

1 0