public void LUDecompositionComputeTest() { for (Int32 matrixIndex = 0; matrixIndex < this.matrices.Length; matrixIndex++) { LUDecomposition decomposition = new LUDecomposition(this.matrices[matrixIndex]); decomposition.Compute(); decomposition.L.ShouldBe(this.expectedL[matrixIndex], 0.01); decomposition.U.ShouldBe(this.expectedU[matrixIndex], 0.01); decomposition.P.ShouldBe(this.expectedP[matrixIndex], 0.01); } }
public void LUDecompositionComputeTest() { for (Int32 i = 0; i < _matrix.Length; i++) { LUDecomposition decomposition = new LUDecomposition(_matrix[i]); decomposition.Compute(); for (Int32 j = 0; j < _matrix[i].NumberOfRows; j++) { for (Int32 k = 0; k < _matrix[i].NumberOfColumns; k++) { Assert.AreEqual(_lExpected[i][j, k], decomposition.L[j, k], 0.01); Assert.AreEqual(_uExpected[i][j, k], decomposition.U[j, k], 0.01); Assert.AreEqual(_pExpected[i][j, k], decomposition.P[j, k], 0.01); } } } }