[DOC] Explaining in tutorial sign ambiguity for numer() and denom().

Based on a patch provided by Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>.

Consider solve_algo::markowitz in automatic elimination algorithm selection.

Cf. <https://www.ginac.de/pipermail/ginac-list/2018-June/002211.html>.

Add optional matrix::rank() algorighm selection.

Before, matrix::rank() was hardcoded to use Bareiss elimination,
which for some matrices is way too slow (and I had to work around
this limitation in my own code). I propose adding an argument that
will allow users to select which elimination scheme to use. The
argument is the same as matrix::solve() takes (so, solve_algo::*).

I've tried to mimic how a similar argument was added to
matrix::inverse() by keeping the current set of rank() functions
and defining separate ones with the 'algo' argument, with the
old ones calling the new ones with solve_algo::automatic. This
is instead of using default argument values.

I've also put the code that makes the automatic selection of
the elimination method into a separate function so that both
matrix::solve() and matrix::rank() could share it. There's a bit
of similar code in matrix::determinant(), but since the self
of determinant algorithms is different from the elimination
algorithms, it's not clear if these two could be unified.

Correct wording in tutorial regarding degree(), ldegree().

Allow expreseq in the arguments of lsolve(eqns, vars).

See <http://www.cebix.net/pipermail/ginac-list/2018-May/002201.html>.

Add Markowitz-ordered Gaussian elimination algorithm.

This algorithm avoids the 'fill-in' problem of Gaussian elimination
and significantly improves the times for solving large sparse systems.
See: <https://www.ginac.de/pipermail/ginac-list/2018-May/002202.html>.

Remove MSVC-specific definition of __func__ and __alignof__.

They are not required any more in C++11.

Avoid multiple filling of partitions in ::get() member functions.

Finalize 1.7.4 release.

Fix ABI compatibility with so-version 6.

Commit 1be7026a unintentionally broke the ABI since g++ doesn't appear
to generate symbols for inline methods. Ugh.

Cf. https://bugzilla.redhat.com/show_bug.cgi?id=1541187

Finalize 1.7.3 release.

Prepare for release 1.7.3 (yeah, again).

Fix two comments in ginac/normal.cpp.

And, while at it, make some disabled debugging output code compile when
it is enabled.

Fix out-of-bounds error in multivar_diophant().

This could trigger an assertion in the factorization exams introduced
in afb0ccaa49 if the library was compiled with -DDO_GINAC_ASSERT.

Shortcut conversion ex->numeric->ex in expairseq::make_flat().

Of course, passing an ex by reference is always faster than converting
it to a numeric first (even if it is a numeric).

Improve gcd(a, b) where one argument is a power of a symbol.

The already implemented recursion
  gcd(x^n, x*p(x)) -> x*gcd(x^(n-1), p(x))
is not ambitious enough: If p(x) has a factor of x, it just goes through
the same step again, and if p(x) has a factor which is a power of x, this
is reapeted many times.

Example:
  gcd(x^n, expand(x^n*(1+x)))
used to go recurse through the gcd routine n times, which could
easily lead to a stack overflow for n about 10^5.

To improve the situation, introduce a special case for gcd where one
argument is a power of a symbol and just cancel the common factor.

This turned out to be the root cause of segfaults in matrix elimination
algorithms, as reported by Patrick Schulz and Vitaly Magerya:
Cf. <https://www.ginac.de/pipermail/ginac-list/2018-January/thread.html>

Make matrix::solve() work with non-normal zeros.

Normalize elements of augmented matrix before checking if they are zero.
This ensures that we don't divide by a non-normal zero. Also added a test.

This bug was reported by Vitaly Magerya <vmagerya@gmail.com>.

[CHECK] Add some more factorization exams.

Add the 15 multivariate polynomials from P. S. Wang's paper "An
Improved Multivariate Polynomial Factoring Algorithm" as exams.

Don't add them as benchmarks since timings vary considerably
depending on internal choice of variables. In any case, this should
never take hours. (But before 630db9a0b0 it did, occasionally.)

Add optional algorithm selection to matrix::inverse().

Following a proposal by Vitaly Magerya, a signature with a solve_algo
is added to class matrix. (Not using a default argument in order to
not break the ABI.)

Fix previous commit 'Speed up special cases of square-free...'.

Huh, I swear 4ffb3cbb3 worked for me! Sorry for that.

Speed up special cases of square-free factorization.

Square-free factorization of polynomials containing a factor which is
a high power P of a symbol x did scale like O(P) in space and time.
This patch introduces a shortcut in Yun's algorithm, such that the
computation is only O(1) in space and time.

This makes it possible to compute, say sqrfree(x^P + x^(P+1)) =>
(1+x)*x^P with P=123456789.

Found this to be a bottleneck while debugging one of Vitaly Magerya's
examples.

Prepare for release 1.7.3.

Happy New Year!

[BUGFIX] Speed up some cases of polynomial factorization.

Speed up a naive O(N) modulus operation in multivar_diophant().
This loop caused some polynomials to seemingly hang forever.

Do not expand pre-factored polynomials.

See <https://www.ginac.de/pipermail/ginac-list/2017-December/002155.html>.

Minor sreamlining of some internal code using class lst...

...without any external effect at all.

Add LaTeX pretty-print for function derivatives.

Update m4/ax_cxx_compile_stdcxx.m4 to current upstream version.

[PATCH] Fix bug in evaluation of real and imaginary parts of powers.

Be more careful when simplifying Re(a^c) -> a^c and Im(a^c) -> 0.

Thanks to Jan Rheinländer for reporting this.

Introducing method symbol::get_TeX_name() to match existing symbol::get_name().

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>

Fix excessive backslashes in the LaTeX function names.

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>

checks: prefer initializer lists over push_back on vectors.

Clean up combinatorial helpers.

* Convert all signatures from int to unsiged.
* Renamed partition_generator to partition_with_zero_parts_generator...
* ...and add class partition_generator which does not include zero parts.
*  Rename ::current() to ::get() and make them compute vectors on first use.

Move combinatorial helpers from power.cpp to utils.h.

This way, these can be used by other modules.

Finalize 1.7.2 release.

Happy New Year!

Update copyright statements.

Fix memory leaks in operators returning relationals.

These memory leaks were introduced in 82df71852.

Thanks to Thomas Luthe for reporting this.

GINAC_DECLARE_REGISTERED_CLASS(): add namespace of status_flags.

Fix MSC compiler warning.

Thanks to Jan Rheinländer for reporting this.

Clean up some header files and fix compilation with MSC.

On MSC there is no <unistd.h>. It turns out that we need to #include <io.h>
in order to use close(3). This fixes the build failure introduced in 8305ec38.

Thanks to Jan Rheinländer for reporting this.

Remove configure check for type 'long long'.

This type with at least 64 bits is mandatory since C++11.

Replace mktemp with mkstemp.

As the manpage says, mktemp is deprecated ("Never use mktemp") and gcc
warns about the usage of mktemp. Replace mktemp with its safe variant
mkstemp.

Finalize 1.7.1 release.

Suppress compiler warnings about inconsistent override.

The introduction of the 'override' keyword in GiNaC in d5b86dd10 was
deliberately incomplete: Duplication of code in macros was avoided.
It turns out that some compilers (e.g. g++ -Wsuggest-override or clang++
-Wall) complain if the 'override' keyword is used inconsistently. Let's
make this consistent now, even if it leads to the introduction of two
more macros.

Fix compilation of programs using __attribute__((deprecated)).

Ouu, we mustn't redefine 'deprecated' in order for this to work...

Thanks to Alexander Proskurin for reporting this.

Update ax_cxx_compile_stdcxx.m4.

Serial 3 fixes make CXXFLAGS=..., c.f. the macro's upstream revision history.

Improve performance of clifford_to_lst().

Rewrite get_clifford_comp() so it is less obscure and more efficient.

Fix convergence of Li_projection at ±I.

For arguments near ±I, direct summation does not practically converge.
Use the Bernoulli transformation instead.

Thanks to Peter Banks for reporting this:
<http://www.cebix.net/pipermail/ginac-list/2016-August/002085.html>

Add proper functions to make clifford_bar() and clifford_star().

Previously operations clifford_bar() and clifford_star() called
conjugate() method. This results in reversion of all non-commutative
entries (not only Clifford units) and produced complex conjugation
of all non-real items. The new routine operates on Clifford units only.

[bugfix] Fix crash in basic::subs().

A regression in 1.7.0 was introduced with 1b8bcb06 in function
basic::subs_one_level(): implicitly casting *this to an ex first for
finding *this in m and later in the function's return statement caused
a crash in the second cast because *this was deleted in the first one.
After all, *this was dynamically allocated in basic::subs().

This bug was reported and hunted down by Mario Prausa.

Finalize 1.7.0 release.

More shortcuts for trivial cases in mul::combine_*_to_pair().

Update reference to Robert C. Martin's "Acyclic Visitor" paper.

The source at objectmentor.com has been unavailble for a long while.

Add fderivative::derivatives() method.

Fix typos in comments.

Update compiler list and compile times.

Prefer C library wrappers over C headers, where possible.

Merge branch 'c++11' in preparation for version 1.7.0.

[build] Install compiler.h header file.

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

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/.)

[build] Ship FindLibDL.cmake, it's needed by CMakeLists.txt.

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.

Finalize 1.6.7 release.

Use C++11 'auto' and range-based for loops where possible.

Use C++11 'nullptr' where applicable.

[build] Force C++11 compilation when using CMake.

[build] Fix CMake build.

minor optimization in power::expand_add

Fix a comment.

speed up multiply_lcm(e, lcm) a bit

and comment it a little

[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.

[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.

[bugfix] fix sqrfree(poly) for zero polynomials in disguise.

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

[bugfix] fix sqrfree(poly) for zero polynomials in disguise.

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

Remove extra cases for missing <stdint.h>.

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

[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.

[C++17] Replace std::bind2nd with a lambda.

Remove global variable max_recursion_level.

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

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.

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.

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.

Happy New Year!

Happy New Year!

De-document eval() in ginsh help and manpage.

This was forgotten in 6c946d4c.

Make relational::rhs() and lhs() nonvirtual.

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

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>

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>

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.

Finalize 1.6.6 release.

Remove reference to ftpthep.physik.uni-mainz.de.

Higgs is down, alas...

Remove reference to ftpthep.physik.uni-mainz.de.

Higgs is down, alas...

Hoist typeid(*this) from loops.

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

Remove iter_swap specializations.

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

Add support for power::info(info_flags::real).

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

Remove info_flags::algebraic.

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

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.

include <algorithm> for use of std::min().

Thanks to Jan Rheinländer for reporting this.

Cherry-picked from a56196f3.

include <algorithm> for use of std::min().

Thanks to Jan Rheinländer for reporting this.

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.

[bugfix] Make sure add::eval() collects all numeric terms *correctly*.

This fixes a bug introduced with in 2010 with 89d5356b.