public void TryInvert()
{
// Regular, square
MatrixD m = new MatrixD(new double[,] {{1, 2, 3, 4},
{2, 5, 8, 3},
{7, 6, -1, 1},
{4, 9, 7, 7}});
MatrixD inverse = m.Clone();
Assert.AreEqual(true, m.TryInvert());
Assert.IsTrue(MatrixD.AreNumericallyEqual(MatrixD.CreateIdentity(4, 4), m * inverse));
// Full column rank, rectangular
m = new MatrixD(new double[,] {{1, 2, 3},
{2, 5, 8},
{7, 6, -1},
{4, 9, 7}});
inverse = m.Clone();
Assert.AreEqual(true, m.TryInvert());
Assert.IsTrue(MatrixD.AreNumericallyEqual(m, m * inverse * m));
// singular
m = new MatrixD(new double[,] {{1, 2, 3},
{2, 5, 8},
{3, 7, 11}});
inverse = m.Clone();
Assert.AreEqual(false, m.TryInvert());
}
