Changeset 952298

Timestamp:

May 11, 2012, 7:45:28 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

36d56c

Parents:

e271f0

Message:

Improve serial performance by 10 - 20 percent by explicit loop unrolling in GSRB smoother.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@1814 5161e1c8-67bf-11de-9fd5-51895aff932f

Location:

src

Files:

• Makefile.am (modified) (1 diff)
• smoother/gsrb_poisson_2.cpp (added)
• smoother/gsrb_poisson_2.hpp (added)
• smoother/gsrb_poisson_4.cpp (added)
• smoother/gsrb_poisson_4.hpp (added)
• units/particle/interface_fcs.cpp (modified) (2 diffs)

src/Makefile.am

@@ -117 +117 @@
         smoother/gsrb.cpp \
         smoother/gsrb.hpp \
+        smoother/gsrb_poisson_2.cpp \
+        smoother/gsrb_poisson_2.hpp \
+        smoother/gsrb_poisson_4.cpp \
+        smoother/gsrb_poisson_4.hpp \
         smoother/jacobi.cpp \
         smoother/jacobi.hpp \

src/units/particle/interface_fcs.cpp

@@ -52 +52 @@
 #include "samples/discretization_poisson_fv.hpp"
 #include "samples/techniques.hpp"
-#include "smoother/gsrb.hpp"
+#include "smoother/gsrb_poisson_2.hpp"
+#include "smoother/gsrb_poisson_4.hpp"
 #ifdef HAVE_LAPACK
 #include "solver/dgesv.hpp"
@@ -158 +159 @@
    * Use Gauss-Seidel Red-Black ordering
    */
-  Smoother* smoother = new GaussSeidelRB();
-  smoother->Register("SMOOTHER");
+  if (discretization_order == 2) {
+    Smoother* smoother = new GaussSeidelRBPoisson2();
+    smoother->Register("SMOOTHER");
+  }else {
+    Smoother* smoother = new GaussSeidelRBPoisson4();
+    smoother->Register("SMOOTHER");
+  }
 
   /*