void daxpy(long,double,double[],long,double[],long);
void dcopy(long,double[],long,double[],long);
void drot(long,double[],long,double[],long,double,double);
void drotg(double*,double*,double*,double*);