Changeset 894a5f for src/level/level_operator_cs.cpp

Timestamp:

Feb 2, 2012, 1:58:12 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

32ff22

Parents:

01be70

Message:

Parallel performance update.

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

File:

• src/level/level_operator_cs.cpp (modified) (2 diffs)

src/level/level_operator_cs.cpp

@@ -70 +70 @@
 
     for (stencil_iter=op.begin(); stencil_iter!=op.end(); ++stencil_iter)
-      res += stencil_iter->Val() * temp->GetVal(*iter_f + stencil_iter->Disp());
+      res += stencil_iter->Val() * temp->GetVal(iter_f->X() + stencil_iter->Disp().X(),
+                                                iter_f->Y() + stencil_iter->Disp().Y(),
+                                                iter_f->Z() + stencil_iter->Disp().Z());
 
     rhs_c(*iter_c) = res;
@@ -129 +131 @@
 
   for (iter_f=bounds_f.Begin(), iter_c=bounds_c.Begin(); iter_c!=bounds_c.End(); iter_f+=2, ++iter_c) {
-
     val = sol_c.GetVal(*iter_c);
-
     sol_f(*iter_f) += op.GetDiag() * val;
-
-    for (stencil_iter = op.begin(); stencil_iter != op.end(); ++stencil_iter) {
-
-      sol_f(*iter_f + stencil_iter->Disp()) += stencil_iter->Val() * val;
-
-    }
-
+    for (stencil_iter = op.begin(); stencil_iter != op.end(); ++stencil_iter)
+      sol_f(iter_f->X() + stencil_iter->Disp().X(),
+            iter_f->Y() + stencil_iter->Disp().Y(),
+            iter_f->Z() + stencil_iter->Disp().Z()) += stencil_iter->Val() * val;
   }