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);
}
