40
40
Float tol = std::numeric_limits<Float>::epsilon();
42
42
vec<Float> x (a.col_ownership(), 0);
43
int status = pcg (a, x, b, eye<Float>(), max_iter, tol, &derr);
43
int status = pcg (a, x, b, eye(), max_iter, tol, &derr);
46
46
Float err = sqrt(dot (x-x_ex,x-x_ex));