fail("*** cholesky: Error in method elmult() or elsolve()");
}
double[][] ainv = buildMatrix(a);
ach.inverse(ainv);
sbeps = 5.e-15;
localflag = maxel(matsub(matmul(ainv,aposdef),ident(aposdef.length,1.))) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {
fail("*** cholesky: Error in method inverse()");