public void LUDecompositionInvertTest() { for (Int32 matrixIndex = 0; matrixIndex < this.matrices.Length; matrixIndex++) { Matrix inverse = LUDecomposition.Invert(this.matrices[matrixIndex]); inverse.ShouldBe(this.expectedInverse[matrixIndex], 0.01); } Should.Throw <ArgumentNullException>(() => LUDecomposition.ComputeDeterminant(null)); }