void daxpy(long,double,double[],long,double[],long);
double ddot(long,double[],long,double[],long);
void dhtgen(long,long,long,long,double[],long,LOGICAL32,
	double*,double[],long,long,LOGICAL32);
double dnrm2(long,double[],long);