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
treefac74455fa47d6575686cf5480d1d2a74858758f
parentafb0ccaa49a0cca001d854594e09125a58434123
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
ginac/matrix.cpp