void silupm(long,float[],float*,long[],float[],float[],long[],long[],long[],
      float[]);