bubble-dynamics

commit d5cd12233e9f3fe0497e3122a592816312dc9e9a
parent 568c00aa72a9cef6a69769409285fb88fcf922d4
Author: Ursula Anna Rasthofer <ursular@euler02.hpc-lca.ethz.ch>
Date:   Fri,  3 Feb 2017 19:08:11 +0100

changed linear solver for clouds from fullPivHouseholderQr() to colPivHouseholderQr()

Diffstat:
M
src/kernels/KMCluster_TY.h
 | 
4
+++-

1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/kernels/KMCluster_TY.h b/src/kernels/KMCluster_TY.h
@@ -129,7 +129,9 @@ public:
         }
 
         /* std::cout << A << std::endl << std::endl; */
-        Eigen::Matrix<Real, Eigen::Dynamic, 1> Rddot = A.fullPivHouseholderQr().solve(b);
+        //Eigen::Matrix<Real, Eigen::Dynamic, 1> Rddot = A.fullPivHouseholderQr().solve(b);
+        // this is notably faster and sufficiently accurate (rasthofer Feb 3, 2017)
+        Eigen::Matrix<Real, Eigen::Dynamic, 1> Rddot = A.colPivHouseholderQr().solve(b);
 
         for (size_t i = 0; i < _N; ++i)
         {