Changeset ac6d04 for src/samples/discretization_poisson_fv.cpp

Timestamp:

Apr 10, 2012, 1:55:49 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

a40eea

Parents:

d24c2f

Message:

Merge recent changes of the vmg library into ScaFaCos.

Includes a fix for the communication problems on Jugene.

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

File:

• src/samples/discretization_poisson_fv.cpp (modified) (8 diffs)

src/samples/discretization_poisson_fv.cpp

@@ -23 +23 @@
   Index i_c, i_f;
 
-  const vmg_float h2_inv = 0.5 / sol_f.MeshWidth();
+  const Vector h2_inv = 0.5 / sol_f.Extent().MeshWidth();
 
-  const Index b1_c = sol_c.Local().AlignmentBegin();
-  const Index b2_c = sol_c.Local().AlignmentEnd() - 1;
+  const Index b1_c = sol_c.Local().FinerBeginFoo();
+  const Index b2_c = sol_c.Local().FinerEndFoo() - 1;
   const Index b1_f = 0;
   const Index b2_f = rhs_f.Local().SizeTotal() - 1;
@@ -32 +32 @@
   const Index begin_f = sol_f.Local().Begin();
   const Index end_f = sol_f.Local().End();
-  const Index begin_c = sol_c.Local().AlignmentBegin();
+  const Index begin_c = sol_c.Local().FinerBeginFoo();
 
   const vmg_float c_1_3 = 1.0 / 3.0;
-  const vmg_float c_2_3_sp = 2.0 / 3.0 * sol_f.MeshWidth();
+  const Vector c_2_3_sp = 2.0 / 3.0 * sol_f.Extent().MeshWidth();
   const vmg_float c_4_3 = 4.0 / 3.0;
 
@@ -43 +43 @@
   for (i_f.Y()=begin_f.Y(), i_c.Y()=begin_c.Y(); i_f.Y()<end_f.Y(); i_f.Y()+=2, ++i_c.Y())
     for (i_f.Z()=begin_f.Z(), i_c.Z()=begin_c.Z(); i_f.Z()<end_f.Z(); i_f.Z()+=2, ++i_c.Z()) {
-      assert(b2_c.X()+1 >= 0 && b2_c.X() < sol_c.Local().SizeTotal().X());
-      assert(i_c.Y() >= 0 && i_c.Y() < sol_c.Local().SizeTotal().Y());
-      assert(i_c.Z() >= 0 && i_c.Z() < sol_c.Local().SizeTotal().Z());
-      assert(b2_f.X()-1 >= 0 && b2_f.X()-1 < sol_f.Local().SizeTotal().X());
-      assert(i_f.Y() >= 0 && i_f.Y() < sol_f.Local().SizeTotal().Y());
-      assert(i_f.Z() >= 0 && i_f.Z() < sol_f.Local().SizeTotal().Z());
-      rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) = (sol_c(b1_c.X()-1,i_c.Y(),i_c.Z()) - sol_f(b1_f.X()+1,i_f.Y(),i_f.Z())) * h2_inv;
-      rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) = (sol_c(b2_c.X()+1,i_c.Y(),i_c.Z()) - sol_f(b2_f.X()-1,i_f.Y(),i_f.Z())) * h2_inv;
+      rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) = (sol_c(b1_c.X()-1,i_c.Y(),i_c.Z()) - sol_f(b1_f.X()+1,i_f.Y(),i_f.Z())) * h2_inv.X();
+      rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) = (sol_c(b2_c.X()+1,i_c.Y(),i_c.Z()) - sol_f(b2_f.X()-1,i_f.Y(),i_f.Z())) * h2_inv.X();
     }
 
@@ -82 +76 @@
     for (i_f.Z()=begin_f.Z(); i_f.Z()<end_f.Z(); ++i_f.Z()) {
 
-      rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) = sol_f(b1_f.X(),i_f.Y(),i_f.Z()) = c_2_3_sp * rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) +
+      rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) = sol_f(b1_f.X(),i_f.Y(),i_f.Z()) = c_2_3_sp.X() * rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) +
         c_4_3 * sol_f(b1_f.X()+1,i_f.Y(),i_f.Z()) -
         c_1_3 * sol_f(b1_f.X()+2,i_f.Y(),i_f.Z());
 
-      rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) = sol_f(b2_f.X(),i_f.Y(),i_f.Z()) = c_2_3_sp * rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) +
+      rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) = sol_f(b2_f.X(),i_f.Y(),i_f.Z()) = c_2_3_sp.X() * rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) +
         c_4_3 * sol_f(b2_f.X()-1,i_f.Y(),i_f.Z()) -
         c_1_3 * sol_f(b2_f.X()-2,i_f.Y(),i_f.Z());
@@ -97 +91 @@
   for (i_f.X()=begin_f.X(), i_c.X()=begin_c.X(); i_f.X()<end_f.X(); i_f.X()+=2, ++i_c.X())
     for (i_f.Z()=begin_f.Z(), i_c.Z()=begin_c.Z(); i_f.Z()<end_f.Z(); i_f.Z()+=2, ++i_c.Z()) {
-      rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) = (sol_c(i_c.X(),b1_c.Y()-1,i_c.Z()) - sol_f(i_f.X(),b1_f.Y()+1,i_f.Z())) * h2_inv;
-      rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) = (sol_c(i_c.X(),b2_c.Y()+1,i_c.Z()) - sol_f(i_f.X(),b2_f.Y()-1,i_f.Z())) * h2_inv;
+      rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) = (sol_c(i_c.X(),b1_c.Y()-1,i_c.Z()) - sol_f(i_f.X(),b1_f.Y()+1,i_f.Z())) * h2_inv.Y();
+      rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) = (sol_c(i_c.X(),b2_c.Y()+1,i_c.Z()) - sol_f(i_f.X(),b2_f.Y()-1,i_f.Z())) * h2_inv.Y();
     }
 
@@ -130 +124 @@
     for (i_f.Z()=begin_f.Z(); i_f.Z()<end_f.Z(); ++i_f.Z()) {
 
-      rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) = sol_f(i_f.X(),b1_f.Y(),i_f.Z()) = c_2_3_sp * rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) +
+      rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) = sol_f(i_f.X(),b1_f.Y(),i_f.Z()) = c_2_3_sp.Y() * rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) +
         c_4_3 * sol_f(i_f.X(),b1_f.Y()+1,i_f.Z()) -
         c_1_3 * sol_f(i_f.X(),b1_f.Y()+2,i_f.Z());
 
-      rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) = sol_f(i_f.X(),b2_f.Y(),i_f.Z()) = c_2_3_sp * rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) +
+      rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) = sol_f(i_f.X(),b2_f.Y(),i_f.Z()) = c_2_3_sp.Y() * rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) +
         c_4_3 * sol_f(i_f.X(),b2_f.Y()-1,i_f.Z()) -
         c_1_3 * sol_f(i_f.X(),b2_f.Y()-2,i_f.Z());
@@ -145 +139 @@
   for (i_f.X()=begin_f.X(), i_c.X()=begin_c.X(); i_f.X()<end_f.X(); i_f.X()+=2, ++i_c.X())
     for (i_f.Y()=begin_f.Y(), i_c.Y()=begin_c.Y(); i_f.Y()<end_f.Y(); i_f.Y()+=2, ++i_c.Y()) {
-      rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) = (sol_c(i_c.X(),i_c.Y(),b1_c.Z()-1) - sol_f(i_f.X(),i_f.Y(),b1_f.Z()+1)) * h2_inv;
-      rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) = (sol_c(i_c.X(),i_c.Y(),b2_c.Z()+1) - sol_f(i_f.X(),i_f.Y(),b2_f.Z()-1)) * h2_inv;
+      rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) = (sol_c(i_c.X(),i_c.Y(),b1_c.Z()-1) - sol_f(i_f.X(),i_f.Y(),b1_f.Z()+1)) * h2_inv.Z();
+      rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) = (sol_c(i_c.X(),i_c.Y(),b2_c.Z()+1) - sol_f(i_f.X(),i_f.Y(),b2_f.Z()-1)) * h2_inv.Z();
     }
 
@@ -177 +171 @@
     for (i_f.Y()=begin_f.Y(); i_f.Y()<end_f.Y(); ++i_f.Y()) {
 
-      rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) = sol_f(i_f.X(),i_f.Y(),b1_f.Z()) = c_2_3_sp * rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) +
+      rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) = sol_f(i_f.X(),i_f.Y(),b1_f.Z()) = c_2_3_sp.Z() * rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) +
         c_4_3 * sol_f(i_f.X(),i_f.Y(),b1_f.Z()+1) -
         c_1_3 * sol_f(i_f.X(),i_f.Y(),b1_f.Z()+2);
 
-      rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) = sol_f(i_f.X(),i_f.Y(),b2_f.Z()) = c_2_3_sp * rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) +
+      rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) = sol_f(i_f.X(),i_f.Y(),b2_f.Z()) = c_2_3_sp.Z() * rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) +
         c_4_3 * sol_f(i_f.X(),i_f.Y(),b2_f.Z()-1) -
         c_1_3 * sol_f(i_f.X(),i_f.Y(),b2_f.Z()-2);