void dmpdrv(double[],long,double[],long*);
void ierm1(STRING,long,long,STRING,STRING,long,byte);