def test_solve(self, backend, mat, b, x, f):
        mat.assemble()
        op2.solve(mat, x, b)
        eps = 1.e-8
        assert_allclose(x.data, f.data, eps)
Solve a linear system where the solution is equal to the right-hand
dataset/ETHPy150Open OP2/PyOP2/test/unit/test_matrices.py/TestMatrices.test_solve