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