void dcopy(long,double[],long,double[],long);
void ddaslx(void(*)(),long,double*,double[],double[],double,long[],double[],
      double[],long*,double[],long,long[],long);
void djacg(long*,long,long,double[],double[],double*,long,double[],double[],
      long[],double[],long,long[],long);