void shtcc(long,long,long,long,float[],float*,float[],
	long,long);
float snrm2(long,float[],long);