void dtggrd(double[],double[],long,long[],double[],long[],long,long[][4],long,
      long*,long[]);
void dtgpd(double[],double[],double[],double[][2],long,long[],long,long[]);
void dtgprg(double[],double[],long,long[],long[][4],long,long);
void dtgrec(double[],double[],double[],double[][2],long,long[],long,long[][4],
      long,double[],long,long,double,double*,long,long,long,LOGICAL32,double*);