void daxpy(long,double,double[],long,double[],long);
void dcopy(long,double[],long,double[],long);
double ddot(long,double[],long,double[],long);