void dchol(double*,long,long,double[],double*,double,long*);
double ddot(long,double[],long,double[],long);