Mercurial > octave

comparison liboctave/array/CMatrix.cc @ 21269:3c8a3d35661a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 e69eaee28737
children cbced1c09916
comparison
equal deleted inserted replaced
21268:f08ae27289e4 21269:3c8a3d35661a
36 #include <sys/types.h> 36 #include <sys/types.h>
37 37
38 #include "Array-util.h" 38 #include "Array-util.h"
39 #include "boolMatrix.h" 39 #include "boolMatrix.h"
40 #include "chMatrix.h" 40 #include "chMatrix.h"
41 #include "chol.h"
41 #include "dMatrix.h" 42 #include "dMatrix.h"
42 #include "CMatrix.h" 43 #include "CMatrix.h"
43 #include "CNDArray.h" 44 #include "CNDArray.h"
44 #include "CRowVector.h" 45 #include "CRowVector.h"
45 #include "dRowVector.h" 46 #include "dRowVector.h"
46 #include "CDiagMatrix.h" 47 #include "CDiagMatrix.h"
47 #include "dDiagMatrix.h" 48 #include "dDiagMatrix.h"
48 #include "CmplxCHOL.h"
49 #include "schur.h" 49 #include "schur.h"
50 #include "CmplxSVD.h" 50 #include "CmplxSVD.h"
51 #include "DET.h" 51 #include "DET.h"
52 #include "f77-fcn.h" 52 #include "f77-fcn.h"
53 #include "functor.h" 53 #include "functor.h"
1107 ret = tinverse (mattype, info, rcon, force, calc_cond); 1107 ret = tinverse (mattype, info, rcon, force, calc_cond);
1108 else 1108 else
1109 { 1109 {
1110 if (mattype.is_hermitian ()) 1110 if (mattype.is_hermitian ())
1111 { 1111 {
1112 ComplexCHOL chol (*this, info, true, calc_cond); 1112 chol<ComplexMatrix> chol (*this, info, true, calc_cond);
1113 if (info == 0) 1113 if (info == 0)
1114 { 1114 {
1115 if (calc_cond) 1115 if (calc_cond)
1116 rcon = chol.rcond (); 1116 rcon = chol.rcond ();
1117 else 1117 else