void daxpy(long,double,double[],long,double[],long);
double ddot(long,double[],long,double[],long);
void dgbsl(double*,long,long,long,long,long[],double[],
	long);