ginac.git
3 years ago[build] Install compiler.h header file. c++11
Richard Kreckel [Mon, 4 Apr 2016 18:42:12 +0000 (20:42 +0200)]
[build] Install compiler.h header file.

Some header files need the deprecated macro. (Looking forward to C++14.)

3 years agoInstall ginac-excompiler in $LIBEXECDIR, not in $BINDIR.
Richard Kreckel [Mon, 4 Apr 2016 07:05:50 +0000 (09:05 +0200)]
Install ginac-excompiler in $LIBEXECDIR, not in $BINDIR.

...and make the GiNaC library aware of where it is installed.

The ginac-excompiler script is only invoked by GiNaC::compile_ex(...)
and serves no purpose on its own. In compliance with the FHS, it should
be installed in $LIBEXECDIR, not in $BINDIR. This also disburdens
distribution packagers from having to provide a manpage (which may be
required for all programs in $BINDIR).

The location for $LIBEXECDIR defaults to ${prefix}/libexec/. It may be
overwritten at configuration time. (Many distributions want to set it to
${prefix}/lib/ginac/.)

3 years ago[build] Ship FindLibDL.cmake, it's needed by CMakeLists.txt.
Richard Kreckel [Sun, 3 Apr 2016 15:46:47 +0000 (17:46 +0200)]
[build] Ship FindLibDL.cmake, it's needed by CMakeLists.txt.

3 years agoRemove ex::to_rational(lst) and ex::to_polynomial(lst).
Richard Kreckel [Sat, 12 Mar 2016 23:02:57 +0000 (00:02 +0100)]
Remove ex::to_rational(lst) and ex::to_polynomial(lst).

Those were GiNaC 1.1 compatibility functions. Use the corresponding
functions with exmap arguments instead.

3 years agoUse C++11 'auto' and range-based for loops where possible.
Richard Kreckel [Thu, 25 Feb 2016 21:12:28 +0000 (22:12 +0100)]
Use C++11 'auto' and range-based for loops where possible.

3 years agoUse C++11 'nullptr' where applicable.
Richard Kreckel [Mon, 22 Feb 2016 21:51:40 +0000 (22:51 +0100)]
Use C++11 'nullptr' where applicable.

3 years ago[build] Force C++11 compilation when using CMake.
Richard Kreckel [Thu, 11 Feb 2016 22:10:55 +0000 (23:10 +0100)]
[build] Force C++11 compilation when using CMake.

3 years ago[build] Fix CMake build.
Richard Kreckel [Thu, 11 Feb 2016 22:08:12 +0000 (23:08 +0100)]
[build] Fix CMake build.

3 years agominor optimization in power::expand_add
Richard Kreckel [Sat, 6 Feb 2016 22:19:35 +0000 (23:19 +0100)]
minor optimization in power::expand_add

3 years agoFix a comment.
Richard Kreckel [Sat, 6 Feb 2016 11:53:14 +0000 (12:53 +0100)]
Fix a comment.

3 years agospeed up multiply_lcm(e, lcm) a bit
Richard Kreckel [Sat, 6 Feb 2016 00:35:07 +0000 (01:35 +0100)]
speed up multiply_lcm(e, lcm) a bit

and comment it a little

3 years ago[bugfix] fix elusive bug in quo, rem,...
Richard Kreckel [Fri, 5 Feb 2016 23:47:08 +0000 (00:47 +0100)]
[bugfix] fix elusive bug in quo, rem,...

The power of two rational numeric objects needs not be rational. As a
result, some care is required when transforming (b^e)*l -> (b*l^(1/e))^e
for some rational e and l. This is a common transformation in order to
convert a polynomial over Q into a polynomial over Z when computing
their quotient, remainder, etc. Failure to be careful can potentially
introduce spurious non-rational numbers into rational polynomials and
make those operations fail. This patch avoids this transformation when
l^(1/e) is not a rational number.

3 years ago[bugfix] fix sqrfree(poly) for zero polynomials in disguise.
Richard Kreckel [Wed, 3 Feb 2016 07:13:58 +0000 (08:13 +0100)]
[bugfix] fix sqrfree(poly) for zero polynomials in disguise.

Yun's algorithm now handles polynomials which would become zero after
expanding.

3 years agoRemove extra cases for missing <stdint.h>.
Richard Kreckel [Sun, 31 Jan 2016 09:15:50 +0000 (10:15 +0100)]
Remove extra cases for missing <stdint.h>.

Use <cstdint> unconditionally -- it's been standardized in C++11.

3 years ago[C++17] Remove inheritance from std::binary_function, std::unary_function.
Richard Kreckel [Fri, 29 Jan 2016 22:45:21 +0000 (23:45 +0100)]
[C++17] Remove inheritance from std::binary_function, std::unary_function.

This inheritance is unneeded since C++11 and the two classes may be
removed in C++-17.

3 years ago[C++17] Replace std::bind2nd with a lambda.
Richard Kreckel [Fri, 29 Jan 2016 21:57:52 +0000 (22:57 +0100)]
[C++17] Replace std::bind2nd with a lambda.

3 years agoRemove global variable max_recursion_level.
Richard Kreckel [Thu, 28 Jan 2016 23:11:14 +0000 (00:11 +0100)]
Remove global variable max_recursion_level.

With the removal of the last 'level' parameter, this variable has
become obsolete.

3 years agoRemove 'level' argument of normal().
Richard Kreckel [Thu, 28 Jan 2016 21:45:56 +0000 (22:45 +0100)]
Remove 'level' argument of normal().

The 'level' argument was modeled after that of the eval() methods
(removed in 6c946d4c). It has never been very useful except for
confusing developers and it hasn't been documented in the tutorial.
Moreover, I have no indication that it has ever been used.

3 years agoRemove 'level' argument of evalf().
Richard Kreckel [Thu, 28 Jan 2016 21:11:46 +0000 (22:11 +0100)]
Remove 'level' argument of evalf().

The 'level' argument was modeled after that of the eval() methods
(removed in 6c946d4c). It has never been very useful except for
confusing developers. Moreover, I have no indication that it has
ever been used.

3 years agoAvoid x^0 and Order(x^0) terms together in series expansion.
Richard Kreckel [Wed, 6 Jan 2016 21:34:23 +0000 (22:34 +0100)]
Avoid x^0 and Order(x^0) terms together in series expansion.

At the branch cut, the series expansions of log(), atan(), and atanh()
assembled pseries objects which contained both x^0 and Order(x^0) terms.
This patch removes the extra x^0 term if the order is 0. It also adds
a GINAC_ASSERT for these kinds of invariants to the pseries ctor and
simplifies the loops in pseries' degree(), ldegree(), eval() and evalf()
member functions.

3 years agoHappy New Year!
Richard Kreckel [Fri, 1 Jan 2016 01:05:06 +0000 (02:05 +0100)]
Happy New Year!

3 years agoDe-document eval() in ginsh help and manpage.
Richard Kreckel [Fri, 1 Jan 2016 01:02:49 +0000 (02:02 +0100)]
De-document eval() in ginsh help and manpage.

This was forgotten in 6c946d4c.

3 years agoMake relational::rhs() and lhs() nonvirtual.
Richard Kreckel [Thu, 31 Dec 2015 18:34:54 +0000 (19:34 +0100)]
Make relational::rhs() and lhs() nonvirtual.

I am unable to see a point in ever overriding these two member functions.

3 years agoClarification on symmetries of metric of clifford object.
Vladimir V. Kisil [Thu, 31 Dec 2015 15:44:16 +0000 (16:44 +0100)]
Clarification on symmetries of metric of clifford object.

The metric of a clifford object may be non-symmetric. Even if
a metric is defined by a symmetric tensor, clifford object may
not be aware of the symmetry and it needs to be explicitly declared.

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>
3 years agoAdd specific dbgprinttree() method to clifford.
Vladimir V. Kisil [Thu, 31 Dec 2015 15:41:48 +0000 (16:41 +0100)]
Add specific dbgprinttree() method to clifford.

Metric is an important element of a clifford object, thus it is useful
to have the metric to be printed in the detailed debug printing.

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>
3 years agocache pseries coeff accesses in pseries::mul_series.
Ralf Stephan [Thu, 31 Dec 2015 15:36:46 +0000 (16:36 +0100)]
cache pseries coeff accesses in pseries::mul_series.

Coeff is too general as long as we have only ints as exponents.
This patch fixes a potentially cubic slowdown.

3 years agoRemove reference to ftpthep.physik.uni-mainz.de.
Richard Kreckel [Sun, 20 Dec 2015 16:56:09 +0000 (17:56 +0100)]
Remove reference to ftpthep.physik.uni-mainz.de.

Higgs is down, alas...

3 years agoHoist typeid(*this) from loops.
Richard Kreckel [Sun, 20 Dec 2015 16:50:07 +0000 (17:50 +0100)]
Hoist typeid(*this) from loops.

The compilers don't seem to reliably recognize the invariance of this.

3 years agoRemove iter_swap specializations.
Richard Kreckel [Sun, 20 Dec 2015 16:30:34 +0000 (17:30 +0100)]
Remove iter_swap specializations.

It's a myth that iter_swap helps the compiler produce better code.
Let's just use specialized swap instead.

3 years agoAdd support for power::info(info_flags::real).
Richard Kreckel [Sun, 20 Dec 2015 03:44:01 +0000 (04:44 +0100)]
Add support for power::info(info_flags::real).

(And, while at it, clean up the surrounding code a bit.)

3 years agoRemove info_flags::algebraic.
Richard Kreckel [Sun, 20 Dec 2015 03:36:42 +0000 (04:36 +0100)]
Remove info_flags::algebraic.

It doesn't work, never has, and it is quire unclear how to fix it.

3 years agoRemove expairseq::construct_from_2_ex_via_exvector() member function.
Richard Kreckel [Thu, 17 Dec 2015 23:06:31 +0000 (00:06 +0100)]
Remove expairseq::construct_from_2_ex_via_exvector() member function.

This function is much slower than the optimized construct_from_2_ex()
member function, which explains why it wasn't used.

3 years agoinclude <algorithm> for use of std::min().
Richard Kreckel [Thu, 17 Dec 2015 21:57:04 +0000 (22:57 +0100)]
include <algorithm> for use of std::min().

Thanks to Jan Rheinländer for reporting this.

3 years agoRemove useless code in add::eval().
Richard Kreckel [Thu, 17 Dec 2015 21:45:50 +0000 (22:45 +0100)]
Remove useless code in add::eval().

With commit ae6c004b, we have actually a more rigorous solution for the
bug fixed first (and wrongly) on September 23 2010 with commit 89d5356b.
This patch removes the now useless code and adds the new regression check
from master, just in case.

3 years agoAdd rvalue reference ctors from epvectors for add, mul.
Richard Kreckel [Thu, 17 Dec 2015 08:32:15 +0000 (09:32 +0100)]
Add rvalue reference ctors from epvectors for add, mul.

3 years agoMerge some cosmetic patches.
Richard Kreckel [Wed, 16 Dec 2015 21:02:25 +0000 (22:02 +0100)]
Merge some cosmetic patches.

This is a collection of some minor optimizations and indentation fixes.

3 years agoMake .eval() evaluate top-level only.
Richard Kreckel [Wed, 16 Dec 2015 20:22:36 +0000 (21:22 +0100)]
Make .eval() evaluate top-level only.

With the previous patch, some old workarounds have become unnecessary:
If all expressions are evaluated, any object which is an aggregate of
expressions will only ever have to evaluate the top level. As such, it
has become pointless to evaluate child objects of an expression prior
to doing its own term-rewriting. This patch removes the evaluation of
children from all GiNaC objects. It also removes the now superfluous
parameter 'level' of the eval methods.

3 years agoMake add::eval(), mul::eval() work without compromise.
Richard Kreckel [Wed, 16 Dec 2015 12:00:30 +0000 (13:00 +0100)]
Make add::eval(), mul::eval() work without compromise.

If a user asks to evaluate an object, the expectation is that the
library returns an evaluated, canonical expression. Everything else is
a bug. (It doesn't matter whether the user asks explicitly or by
assigning to an ex.) It turns out that it was possible to construct
objects which didn't fully evaluate. This patch fixes this by making
eval() a little bit more careful.

This obsoletes the need to go through combine_ex_with_coeff_to_pair()
when constructing an epvector that is then used to construct an add
or mul. (Alas, this was omitted frequently enough...)

3 years agoImprove method of setting status_flags::dynallocated.
Richard Kreckel [Mon, 30 Nov 2015 22:06:53 +0000 (23:06 +0100)]
Improve method of setting status_flags::dynallocated.

There seems to be no way to obsolete the need to mark an object derived
from basic and handled by ex as being 'on the heap', at least none that
doesn't have significant performance impact. Having said that, this
patch aims at making this process simpler and more intuitive.

Where, before, one would return from a function returning an ex with
        return (new mul(a, b))->setflag(status_flags::dynallocated);
this patch lets us return with
        return dynallocate<mul>(a, b);
which should be much clearer. In any case, it involves less typing.

The two points where the status_flags::dynallocated are set are now
  * the dynallocate<B>(args...) template function and
  * the virtual duplicate() member functions.

This patch rolls out the new functionality throughout the library.

3 years agoRemove accidentally pushed binary.
Richard Kreckel [Sun, 29 Nov 2015 11:42:59 +0000 (12:42 +0100)]
Remove accidentally pushed binary.

In commit 79cad108, a compiled binary from the test suite was accidentally
added and pushed. Sorry.

3 years agoRename array of precomputed data in test suite.
Richard Kreckel [Sun, 29 Nov 2015 11:33:59 +0000 (12:33 +0100)]
Rename array of precomputed data in test suite.

Reason: C++17 may introduce a std::data<> template. Right now, the GCC 6.0
prerelease bails out at this code, when compiler with -std=c++17.
<http://en.cppreference.com/w/cpp/iterator/data>

3 years agoUse neseted initializer lists to construct matrix objects.
Richard Kreckel [Sat, 28 Nov 2015 14:43:46 +0000 (15:43 +0100)]
Use neseted initializer lists to construct matrix objects.

Add constructor of initializer_list<initializer_list<ex>> to matrix.
Use this syntax where, previously, ctor from comma-separated list of
elements was used. Deprecate the ctor from comma-separated list.

Note: The output format '[[a,b],[c,d]]' and ginsh syntax are
unchanged because lists are printed '{a,b,c}' and a matrix is not a
list of lists.

3 years ago[build] Update scripts for correct rpath for linking with CLN.
Richard Kreckel [Fri, 27 Nov 2015 10:49:16 +0000 (11:49 +0100)]
[build] Update scripts for correct rpath for linking with CLN.

Use up-to-date files from GNU gettext 0.19.6.

3 years agoUse initializer lists to construct container<>, lst.
Richard Kreckel [Thu, 26 Nov 2015 20:36:46 +0000 (21:36 +0100)]
Use initializer lists to construct container<>, lst.

Add constructor of initializer_list<ex> to container<C<ex>>. In
particular, this means that we can finally declare lst objects as
lst{a,2*b,2}. Convert GiNaC to this syntax throughout. Deprecate
the old constructors taking 1..16 ex parameters and the ones from
comma seaparated lists (without braces).

3 years agoInfrastructure: Update baseline language check M4 macros.
Richard Kreckel [Thu, 26 Nov 2015 14:04:52 +0000 (15:04 +0100)]
Infrastructure: Update baseline language check M4 macros.

Use the more general ax_cxx_compile_stdcxx.m4 macro from gnu.org
instead of the old ax_cxx_compile_stdcxx_11.m4.

3 years agoMake ample use of the contextual keyword 'override'.
Richard Kreckel [Thu, 26 Nov 2015 11:20:20 +0000 (12:20 +0100)]
Make ample use of the contextual keyword 'override'.

This patch adds the C++11 contextual keyword 'override' to every overridden
virtual function declaration except where that would incur macro duplications.
Along the way, it fixes some comments about member functions 'virtuality'.

3 years agoStandard header cleanup.
Richard Kreckel [Wed, 25 Nov 2015 14:17:20 +0000 (15:17 +0100)]
Standard header cleanup.

Sanitize some standard headers, remove include directives of unused
headers, complete configure check for used standard headers.

3 years agoReplace idx_is_not functor by a C++11 lambda.
Richard Kreckel [Wed, 25 Nov 2015 11:49:31 +0000 (12:49 +0100)]
Replace idx_is_not functor by a C++11 lambda.

3 years agoIn power::expand_add(), don't reserve excess monomial sizes.
Richard Kreckel [Wed, 25 Nov 2015 10:28:10 +0000 (11:28 +0100)]
In power::expand_add(), don't reserve excess monomial sizes.

There is no need to reserve n terms in each of the monomials of the
result of power(+(x,y,z...;0),n): We can compute it exactly as the
number of nonzero exponents in the multinomial expansion. The good
thing is that this counting is the same for each composition of a
partition, so it can be hoisted out of the loop over compositions.

3 years agoMake specialized power::expand() helpers static member fuctions.
Richard Kreckel [Wed, 25 Nov 2015 10:22:34 +0000 (11:22 +0100)]
Make specialized power::expand() helpers static member fuctions.

The class power member functions expand_add(), expand_add_2(), and
expand_mul() do not access any member variable of class power. (In
fact, the only reason not to hide them entirely is that they invoke
protected members of classs expairseq, add, and mul which requires
them to be friends of these classes.)

3 years agoTutorial: how to create noncommutative symbols?
Richard Kreckel [Tue, 24 Nov 2015 12:02:46 +0000 (13:02 +0100)]
Tutorial: how to create noncommutative symbols?

This seems to be frequently asked:
<http://www.ginac.de/pipermail/ginac-devel/2015-November/002213.html>

3 years agoAvoid some explicit copying of STL containers.
Richard Kreckel [Tue, 24 Nov 2015 10:57:35 +0000 (11:57 +0100)]
Avoid some explicit copying of STL containers.

3 years agomul::can_make_flat(): Finish ancient unfinished code.
Richard Kreckel [Tue, 24 Nov 2015 10:14:51 +0000 (11:14 +0100)]
mul::can_make_flat(): Finish ancient unfinished code.

This makes mul(expairseq({a*b,2},{b})) construct the canonical
mul(expairseq({a,3},{b,2})). Normally, such term rewriting is done
by the .eval() methods, but this is not done when a plain object of
a class derived from basic is constructed.

3 years agounit_matrix(r,c) can return an evaluated matrix object.
Richard Kreckel [Mon, 23 Nov 2015 20:38:30 +0000 (21:38 +0100)]
unit_matrix(r,c) can return an evaluated matrix object.

3 years agoMake (a+b+c+...)^n work for huuuge terms.
Richard Kreckel [Mon, 23 Nov 2015 20:21:18 +0000 (21:21 +0100)]
Make (a+b+c+...)^n work for huuuge terms.

This fix is neccessary for results exceeding 2^31 terms. While at it,
also restructured power::expand_add_() a bit to look more like the
more general power::expand_add().

3 years agoConvert some more ctors from copying to moving STL containers.
Richard Kreckel [Mon, 23 Nov 2015 18:36:53 +0000 (19:36 +0100)]
Convert some more ctors from copying to moving STL containers.

3 years agoReplace static const variables with constexpr, where applicable.
Richard Kreckel [Sun, 8 Nov 2015 22:33:33 +0000 (23:33 +0100)]
Replace static const variables with constexpr, where applicable.

3 years agoFix some pedantic compiler warnings.
Richard Kreckel [Sun, 8 Nov 2015 22:18:19 +0000 (23:18 +0100)]
Fix some pedantic compiler warnings.

Found with CLang++ -Wall. (There remain some more, but these are
more ugly to fix.)

3 years agoReplace GiNaC::ToString() with std::to_string().
Richard Kreckel [Sun, 8 Nov 2015 21:19:08 +0000 (22:19 +0100)]
Replace GiNaC::ToString() with std::to_string().

And, while at it, also make CMake exercise the new real_imag exam
introduced last week and document successful build with GCC-4.8.

3 years agoUse C++11 'nullptr' where applicable.
Richard Kreckel [Sat, 7 Nov 2015 21:24:29 +0000 (22:24 +0100)]
Use C++11 'nullptr' where applicable.

3 years agoRemove C++89 workaround for closing angle brackets in nested templates.
Richard Kreckel [Sat, 7 Nov 2015 12:23:26 +0000 (13:23 +0100)]
Remove C++89 workaround for closing angle brackets in nested templates.

3 years agoRemove 'discardable' option from ctors of container and derived classes.
Richard Kreckel [Sat, 7 Nov 2015 12:03:05 +0000 (13:03 +0100)]
Remove 'discardable' option from ctors of container and derived classes.

The whole idea of this was to allow the ctor to pilfer the data from the
constructed-from object, which is precisely the move semantics which
C++11 supports with rvalue references.

3 years agoAdd trivial shortcuts in expair plumbing of class add.
Richard Kreckel [Sat, 7 Nov 2015 00:50:00 +0000 (01:50 +0100)]
Add trivial shortcuts in expair plumbing of class add.

These little shortcuts in add::split_ex_to_pair(ex) and in
add::combine_ex_with_coeff_to_pair(ex, ex) avoid the creation of
new ex objects. These pieces of code are executed so often that this
patch speeds up GiNaC by 5-10%, depending on test.

3 years agoFix little missed optimization.
Richard Kreckel [Fri, 6 Nov 2015 14:04:49 +0000 (15:04 +0100)]
Fix little missed optimization.

3 years agoReplace dynamic throw() with static noexcept operator.
Richard Kreckel [Fri, 6 Nov 2015 00:00:55 +0000 (01:00 +0100)]
Replace dynamic throw() with static noexcept operator.

3 years agoAdd some minor optimizations.
Richard Kreckel [Thu, 5 Nov 2015 12:41:18 +0000 (13:41 +0100)]
Add some minor optimizations.

3 years agoSpeed up power::real_part() and power::imag_part().
Richard Kreckel [Thu, 5 Nov 2015 12:30:48 +0000 (13:30 +0100)]
Speed up power::real_part() and power::imag_part().

Add special case for real base and exponent (not just integer exponent).
Improve special case for integer exponent by explicitly constructing
the result using the Binomial expansion.
Add a test case for real/imaginary part expansion.

3 years agoUse AX_CXX_COMPILE_STDCXX_11 macro to force C++11 compilation.
Richard Kreckel [Tue, 3 Nov 2015 20:23:16 +0000 (21:23 +0100)]
Use AX_CXX_COMPILE_STDCXX_11 macro to force C++11 compilation.

3 years agoFix most remaining GCC compiler warnings.
Richard Kreckel [Tue, 3 Nov 2015 19:53:07 +0000 (20:53 +0100)]
Fix most remaining GCC compiler warnings.

3 years agoMove ctors for classes matrix and pseries.
Richard Kreckel [Tue, 3 Nov 2015 15:46:06 +0000 (16:46 +0100)]
Move ctors for classes matrix and pseries.

3 years agoUse C++11 range-based foor loops and auto, where possible.
Richard Kreckel [Tue, 3 Nov 2015 13:15:57 +0000 (14:15 +0100)]
Use C++11 range-based foor loops and auto, where possible.

3 years agoUpdate documentation of prerequisites.
Richard Kreckel [Sun, 1 Nov 2015 23:00:48 +0000 (00:00 +0100)]
Update documentation of prerequisites.

3 years agoReplace use of NULL by C++11 nullptr.
Richard Kreckel [Sun, 1 Nov 2015 22:48:05 +0000 (23:48 +0100)]
Replace use of NULL by C++11 nullptr.

3 years agoRemove dependence on depreacted std::auto_ptr<T>.
Richard Kreckel [Sun, 1 Nov 2015 19:16:12 +0000 (20:16 +0100)]
Remove dependence on depreacted std::auto_ptr<T>.

Most of the auto_ptr were introduced to avoid copying. C++11 supports
move semantics for this. Some other can be replaced by std::unique_ptr<T>.

3 years agoRemove support for Cint.
Richard Kreckel [Sat, 31 Oct 2015 16:18:55 +0000 (17:18 +0100)]
Remove support for Cint.

3 years agoRemove unfinished code for hash table-based expairseq.
Richard Kreckel [Sat, 31 Oct 2015 16:02:14 +0000 (17:02 +0100)]
Remove unfinished code for hash table-based expairseq.

This is unfinished code by A.F. Later generations don't ever
seem to have tinkered much with it.  ;-)

3 years agoFix make check with g++ -std=c++11.
Richard Kreckel [Thu, 15 Oct 2015 06:59:20 +0000 (08:59 +0200)]
Fix make check with g++ -std=c++11.

Our template function log() returning an object of type GiNaC::function
turns out to compete with a template from a libstdc++ header. We have
to help the compiler a bit.

3 years agoFix some apparent typos and misindentations 2.
Richard Kreckel [Sat, 10 Oct 2015 20:24:32 +0000 (22:24 +0200)]
Fix some apparent typos and misindentations 2.

(Formatting should comply with http://www.cebix.net/indentation.html.)

3 years agoFix algebraic power::has() for larger integer exponents.
Richard Kreckel [Sat, 10 Oct 2015 09:02:17 +0000 (11:02 +0200)]
Fix algebraic power::has() for larger integer exponents.

There is no need to overflow integers in (x^2147483648).has(x^2).

4 years agoInclude crc32.h only if needed.
Richard Kreckel [Mon, 21 Sep 2015 21:04:20 +0000 (23:04 +0200)]
Include crc32.h only if needed.

This shuts up a compiler warning about unused static function crc32.

4 years agoRemove unused variable.
Richard Kreckel [Mon, 21 Sep 2015 20:49:44 +0000 (22:49 +0200)]
Remove unused variable.

This shuts up a compiler warning.

4 years agoFix parser example in tutorial.
Richard Kreckel [Fri, 31 Jul 2015 22:52:11 +0000 (00:52 +0200)]
Fix parser example in tutorial.

Thanks to Jean-Michaël Celerier for pointing this out.

4 years agoFix some apparent typos and misindentations.
Richard Kreckel [Fri, 31 Jul 2015 22:41:30 +0000 (00:41 +0200)]
Fix some apparent typos and misindentations.

4 years agoFinalizing release 1.6.5. ginac_1-6-5
Richard Kreckel [Mon, 27 Jul 2015 19:28:44 +0000 (21:28 +0200)]
Finalizing release 1.6.5.

4 years agoFix bug in PolyLogs of rational numbers.
Stefan Weinzierl [Fri, 24 Jul 2015 18:50:57 +0000 (20:50 +0200)]
Fix bug in PolyLogs of rational numbers.

Thanks to Alexander Hasselhuhn for reporting this bug.

4 years agoSpeed up pow(+(...),n).expand() where n>2.
Richard Kreckel [Sat, 18 Jul 2015 22:04:49 +0000 (00:04 +0200)]
Speed up pow(+(...),n).expand() where n>2.

This new algorithm is faster than the old one because it re-uses
multinomial coefficients as much as possible.

4 years agoFix pow(+(...),2).expand().
Richard Kreckel [Sat, 18 Jul 2015 21:56:55 +0000 (23:56 +0200)]
Fix pow(+(...),2).expand().

Due to a failure to expand result terms, expand((sqrt(1+x)+y*sqrt(1+x))^2)
returned 1+y^2+x+x*y^2+2*y*(1+x). Note that 2*y*(1+x) was not expanded
to 2*y+2*y*x.

4 years agoFinalized 1.6.4 release. ginac_1-6-4
Richard Kreckel [Fri, 8 May 2015 07:35:24 +0000 (09:35 +0200)]
Finalized 1.6.4 release.

4 years agoFix pow(+(...),2).expand().
Richard Kreckel [Thu, 7 May 2015 20:33:13 +0000 (22:33 +0200)]
Fix pow(+(...),2).expand().

Due to a failure to recombine coeffs and rests to expairs,
expand((x+sqrt(2)*x)^2) returned x^2+2*x^2+2*sqrt(2)*x^2. The
2*x^2 term was not combined with the x^2 term to 3*x^2 because it
was not the canonical expair [[x^2,2]] but rather [[2*x^2,1]].

Thanks to Isuru Fernando for the bugreport.

4 years agoAdd test cases for improved abs(arg).
Richard Kreckel [Thu, 23 Apr 2015 20:33:52 +0000 (22:33 +0200)]
Add test cases for improved abs(arg).

Exam the evaluations abs(±sqrt(p)) -> sqrt(p) for p > 0.

4 years agoImprove abs(arg).
Vladimir Kisil [Wed, 22 Apr 2015 21:21:53 +0000 (23:21 +0200)]
Improve abs(arg).

Make abs(arg) return -arg if arg is negative or if -arg is known to be
nonnegative.

4 years agoFix pow::info(info_flags::nonnegative).
Vladimir Kisil [Wed, 22 Apr 2015 21:06:56 +0000 (23:06 +0200)]
Fix pow::info(info_flags::nonnegative).

This function call was missing the case where a positive base is raised
to a real power (where pow::info(info_flags::positive) correctly returned
true).

4 years agoFix mul::info(info_flags::nonnegative).
Vladimir Kisil [Wed, 22 Apr 2015 19:52:15 +0000 (21:52 +0200)]
Fix mul::info(info_flags::nonnegative).

Due to an obvious Boolean inversion, this function call always
returned the wrong answer.

4 years agoMake power::expand() (x*p)^c -> x^c * p^c, if p>0.
Vladimir Kisil [Tue, 14 Apr 2015 21:14:09 +0000 (23:14 +0200)]
Make power::expand() (x*p)^c -> x^c * p^c, if p>0.

This expansion seems to be helpful in many cases.

4 years agoMake compile_ex() work with CMake.
Martin Vala [Tue, 14 Apr 2015 20:50:31 +0000 (22:50 +0200)]
Make compile_ex() work with CMake.

Infrastructure fixes to make CMake find libdl and make ginac-excompiler
listen to the $CXXFLAGS environment variable.

4 years agoImprove normalization with nested functions.
Ladislav Zejda [Sat, 14 Mar 2015 10:12:11 +0000 (11:12 +0100)]
Improve normalization with nested functions.

normal() fails to fully normalize expressions where nested functions occur
multiple times because replace_with_symbol() searches repl for an
original expression but inserts into repl the substitued expression. Then
the repeated search does not work and a new symbol is introduced.

This patch fixes that by searching for the expression with the repl
substituted instead of for the original expression.

4 years agoFixing index in tutorial.
Vladimir V. Kisil [Sun, 8 Feb 2015 19:53:55 +0000 (20:53 +0100)]
Fixing index in tutorial.

The index entry for clifford::get_metric() did not work, change it to get_metric().

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>
4 years agoExplicit derivation of functions.
Vladimir V. Kisil [Sun, 8 Feb 2015 19:50:51 +0000 (20:50 +0100)]
Explicit derivation of functions.

Some function cannot be cleanly differentiated through the chain rule.
For example, it is natural to define derivative of the absolute value as

(abs(f))'=(f'*f.conjugate()+f*f'.conjugate())/2/abs(f)

This patch adds a possibility to define derivatives of functions in this way.
In particular the derivative of abs(), Order(), real_part(), imag_part() and
conjugate() are defined.

For example, conjugate of a derivative with respect of a real symbol
If x is real then U.diff(x)-I*V.diff(x) represents both
conjugate(U+I*V).diff(x) and conjugate((U+I*V).diff(x))
Thus in this patch we use the rule

conjugate(f)'=conjugate(f')

for a derivative with respect to the real symbol.

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>
4 years agoPut a reference on eta function into tutorial.
Vladimir V. Kisil [Sat, 31 Jan 2015 12:04:13 +0000 (13:04 +0100)]
Put a reference on eta function into tutorial.

4 years agoHappy new year!
Richard Kreckel [Sun, 25 Jan 2015 10:50:43 +0000 (11:50 +0100)]
Happy new year!