void dcov3(double*,long,long,double[],double,double[],long*);
double ddot(long,double[],long,double[],long);
double dnrm2(long,double[],long);
void dsvdrs(double*,long,long,long,double*,long,long,double[],double[]);