public void TestProfileFormatJacobi() { int predCoeff = 0; Vector expectedAnswer = new Vector(dim); for (int i = 0; i < dim; i++) { expectedAnswer[i] = i + 1; } var matrix = generateProfileMatrix(predCoeff); var rightPart = generateRightPart(matrix, expectedAnswer); Jacobi solver = new Jacobi(); EmptyPreconditioner precond = EmptyPreconditioner.Create(matrix); Vector answer = solver.Solve(precond, rightPart, new Vector(dim), Logger.Instance, Logger.Instance, new JacobiParametrs(eps, maxIter)); Console.Write("Generate complete...\n"); for (int i = 0; i < dim; i++) { Console.Write(expectedAnswer[i].ToString() + "\t" + answer[i].ToString() + "\n"); Assert.AreEqual(expectedAnswer[i], answer[i], 0.001, "Not equal!"); } }