void daxpy(long,double,double[],long,double[],long);
void dcopy(long,double[],long,double[],long);
void ddasco(double*,long,long,double[],double[]);
double ddot(long,double[],long,double[],long);
double dnrm2(long,double[],long);
void drot(long,double[],long,double[],long,double,double);
void drotg(double*,double*,double*,double*);
void dscal(long,double,double[],long);
void dswap(long,double[],long,double[],long);
long idamax(long,double[],long);