void dhtcc(long,long,long,long,double[],double*,double[],
	long,long);
double dnrm2(long,double[],long);