Mercurial > octave
diff libinterp/corefcn/__qp__.cc @ 21269:3c8a3d35661a
better use of templates for Cholesky factorization
* liboctave/numeric/chol.h, liboctave/numeric/chol.cc: New files
generated from CmplxCHOL.cc, fCmplxCHOL.cc, floatCHOL.cc, CmplxCHOL.h,
dbleCHOL.cc, dbleCHOL.h, fCmplxCHOL.h, and floatCHOL.h and converted
to templates.
* liboctave/numeric/module.mk: Update.
* __qp__.cc, chol.cc, CMatrix.cc, CMatrix.h, dMatrix.cc, dMatrix.h,
fCMatrix.cc, fCMatrix.h, fMatrix.cc, fMatrix.h, eigs-base.cc,
mx-defs.h, mx-ext.h: Use new classes.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 16 Feb 2016 02:47:29 -0500 |
parents | fcac5dbbf9ed |
children | cbced1c09916 |
line wrap: on
line diff
--- a/libinterp/corefcn/__qp__.cc Tue Feb 16 00:32:29 2016 -0500 +++ b/libinterp/corefcn/__qp__.cc Tue Feb 16 02:47:29 2016 -0500 @@ -26,7 +26,7 @@ #include <cfloat> -#include "dbleCHOL.h" +#include "chol.h" #include "dbleSVD.h" #include "mx-m-dm.h" #include "EIG.h" @@ -191,7 +191,7 @@ // factorization since the Hessian is positive // definite. - CHOL cholH (H); + chol<Matrix> cholH (H); R = cholH.chol_matrix (); @@ -250,7 +250,7 @@ // Computing the Cholesky factorization (pR = 0 means // that the reduced Hessian was positive definite). - CHOL cholrH (rH, pR); + chol<Matrix> cholrH (rH, pR); Matrix tR = cholrH.chol_matrix (); if (pR == 0) R = tR;