void scalcr(long,long,float[],long*,float[]);
void scalcj(long,long,float[],long*,float[]);
void sivset(long,long[],long,long,float[]);
void snlafu(long,long,float[],void(*)(long,long,float[],long*,float[]),long[],
      long,long,float[]);