void daccum(double*,long,long,double*,long,long,long*,
	long,long*);
void dhtcc(long,long,long,long,double[],double*,double[],
	long,long);
double dnrm2(long,double[],long);
void drotg(double*,double*,double*,double*);
void ierm1(STRING,long,long,STRING,STRING,long,byte);
void ierv1(STRING,long,byte);