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