double dasum(long,double[],long);
void daxpy(long,double,double[],long,double[],long);
double dnrm2(long,double[],long);
void dscal(long,double,double[],long);
void dswap(long,double[],long,double[],long);
long idamax(long,double[],long);