void cholesky ( double a[], int n, int nn, double u[], int *nullty, 
  int *ifault );
double r8_max ( double x, double y );
void subchl ( double a[], int b[], int n, double u[], int *nullty, 
  int *ifault, int ndim, double *det );
void timestamp ( );