void daxpy(long,double,double[],long,double[],long);
void dcopy(long,double[],long,double[],long);
void drotm(long,double[],long,double[],long,double[]);
void drotmg(double*,double*,double*,double,double[]);
void dvecp(double[],long,STRING);