public void CheckCountXVectorForJacobi() { double[,] matrix1 = { { 1, 2 }, { 2, 1 } }; double[,] matrix2 = { { 4, -1, -0.2, 2 }, { -1, 5, 0, -2 }, { 0.2, 1, 10, -1 }, { 0, -2, -1, 4 } }; double[] xVector1 = IterationMethodsOperations.SetDefaultVector(2); double[] xVector2 = IterationMethodsOperations.SetDefaultVector(4); double[] bVector1 = { 1, 1 }; double[] bVector2 = { 30, 0, -10, 5 }; JacobiOperations.CountXVector(xVector1, matrix1, bVector1); JacobiOperations.CountXVector(xVector2, matrix2, bVector2); double[] expect1 = { 1, 1 }; double[] expect2 = { 6.825, 2, -1.025, 1 }; //Assert.Equal(expect1, xVector1); Assert.Equal(expect2, xVector2); }
public double[] Jacobi(double[] bVector, int numberOfIterations) { double[] xVector = IterationMethodsOperations.SetDefaultVector(numberOfColumns); JacobiOperations.CountXVector(xVector, matrix, bVector); return(xVector); }