if (o) FREE(o);
}
-static real conjugate_gradient(Operator A, Operator precon, int n, real *x, real *rhs, real tol, int maxit, int *flag){
+static real conjugate_gradient(Operator A, Operator precon, int n, real *x, real *rhs, real tol, int maxit){
real *z, *r, *p, *q, res = 10*tol, alpha;
real rho = 1.0e20, rho_old = 1, res0, beta;
real* (*Ax)(Operator o, real *in, real *out) = A->Operator_apply;
b[i] = rhs[i*dim+k];
}
- res += conjugate_gradient(Ax, precond, n, x, b, tol, maxit, flag);
+ res += conjugate_gradient(Ax, precond, n, x, b, tol, maxit);
for (i = 0; i < n; i++) {
rhs[i*dim+k] = x[i];
}