]> www.ginac.de Git - ginac.git/commit
git://www.ginac.de / ginac.git / commit
? search:
summary | shortlog | log | commit | commitdiff | tree
(parent: afb0cca) | patch
Make matrix::solve() work with non-normal zeros.
authorRichard Kreckel <kreckel@ginac.de>
Sun, 28 Jan 2018 23:07:17 +0000 (00:07 +0100)
committerRichard Kreckel <kreckel@ginac.de>
Sun, 28 Jan 2018 23:07:17 +0000 (00:07 +0100)
commitd0ff428fb5b7bb565a0aea69e05e5705d840c16d
treefac74455fa47d6575686cf5480d1d2a74858758ftree | snapshot
parentafb0ccaa49a0cca001d854594e09125a58434123commit | diff
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/exam_matrices.cpp diff | blob | history
ginac/matrix.cpp diff | blob | history
GiNaC -- a C++ library for symbolic computations