void cholesky ( double a[], int n, int nn, double u[], int *nullty, 
  int *ifault );
void syminv ( double a[], int n, double c[], double w[], int *nullty, 
  int *ifault );