void daxpy(long,double,double[],long,double[],long);
void dgei(double*,long,long,long[],double[]);
void dscal(long,double,double[],long);
void dswap(long,double[],long,double[],long);
void ermsg(STRING,long,long,STRING,byte);