iter 0 residual = 2.666847220715e-15
Converged

Solution vector:

i = 0  a_sol[i] = -1.4414414414e-01
i = 1  a_sol[i] = 6.6666666667e-01
i = 2  a_sol[i] = 8.6486486486e-01

Matrix times solution vector vs. rhs vector:

i = 0  matrix * a_sol[i] = 2.000000e+00  rhs[i] = 2.000000e+00
i = 1  matrix * a_sol[i] = 6.000000e+00  rhs[i] = 6.000000e+00
i = 2  matrix * a_sol[i] = -4.000000e+00  rhs[i] = -4.000000e+00