extern double *c;
extern double c_det;
extern double *c_factor;
extern int c_factored;
extern double *zp_mean;