void saxpy(long,float,float[],long,float[],long);
float sdot(long,float[],long,float[],long);
void sgesl(float*,long,long,long[],float[],long);