void dilupm(long,double[],double*,long[],double[],double[],long[],long[],long[],
      double[]);