//-------------------------------------------- void NRLib::CholeskyInvert(SymmetricMatrix& A) //-------------------------------------------- { //NBNB legge til feilmld her, se metode ovenfor int infof = flens::potrf(A.upLo(), A.dim(), A.data(), A.leadingDimension()); int infoi = flens::potri(A.upLo(), A.dim(), A.data(), A.leadingDimension()); if(infof > 0 || infoi > 0) throw NRLib::Exception("Error in Cholesky inversion"); }