public void SymmetricRandomMatrixInverse() { for (int d = 1; d <= 100; d = d + 11) { SymmetricMatrix M = TestUtilities.CreateSymmetricRandomMatrix(d, 1); SymmetricMatrix MI = M.Inverse(); Assert.IsTrue(TestUtilities.IsNearlyEqual(MI * M, UnitMatrix.OfDimension(d))); } }
public void SymmetricRandomMatrixInverse() { for (int d = 1; d <= 100; d = d + 11) { Console.WriteLine("d={0}", d); SquareMatrix I = TestUtilities.CreateSquareUnitMatrix(d); SymmetricMatrix M = TestUtilities.CreateSymmetricRandomMatrix(d, 1); SymmetricMatrix MI = M.Inverse(); Assert.IsTrue(TestUtilities.IsNearlyEqual(M * MI, I)); } }