4 int pcg(
int n,
int (*atimes) (
int,
double *,
double *,
void *),
5 int (*msolve) (
int,
double *,
double *,
int (*atimes)(
int,
double *,
double *,
void *),
void *),
6 double *x,
double *b,
double *r,
double *z,
double *p,
int itmax,
7 double eps,
double *err,
void *,
double *valpha= 0,
double *vbeta= 0);