Mercurial > octave

diff liboctave/numeric/eigs-base.cc @ 21271:7e67c7f82fc1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
better use of templates for lu factorization classes * liboctave/numeric/lu.h, liboctave/numeric/lu.cc: New files generated from base-lu.h, base-lu.cc, CmplxLU.cc, CmplxLU.h, dbleLU.cc, dbleLU.h, fCmplxLU.cc, fCmplxLU.h, floatLU.cc, and floatLU.h and converted to templates. * liboctave/numeric/module.mk: Update. * lu.cc, mx-defs.h, mx-ext.h, eigs-base.cc: Use new classes.
author John W. Eaton <jwe@octave.org>
date Tue, 16 Feb 2016 12:58:32 -0500
parents 3c8a3d35661a
children 40de9f8f23a6
line wrap: on
line diff
--- a/liboctave/numeric/eigs-base.cc	Tue Feb 16 11:13:55 2016 -0500
+++ b/liboctave/numeric/eigs-base.cc	Tue Feb 16 12:58:32 2016 -0500
@@ -30,11 +30,10 @@
 #include <iostream>
 
 #include "CSparse.h"
-#include "CmplxLU.h"
+#include "lu.h"
 #include "MatrixType.h"
 #include "chol.h"
 #include "dSparse.h"
-#include "dbleLU.h"
 #include "eigs-base.h"
 #include "f77-fcn.h"
 #include "mx-ops.h"
@@ -510,7 +509,7 @@
         p[i*(n+1)] -= sigma;
     }
 
-  LU fact (AminusSigmaB);
+  lu<Matrix> fact (AminusSigmaB);
 
   L = fact.P ().transpose () * fact.L ();
   U = fact.U ();
@@ -674,7 +673,7 @@
         p[i*(n+1)] -= sigma;
     }
 
-  ComplexLU fact (AminusSigmaB);
+  lu<ComplexMatrix> fact (AminusSigmaB);
 
   L = fact.P ().transpose () * fact.L ();
   U = fact.U ();