Changeset 8180d8 for src/samples/discretization_poisson_fv.cpp

Timestamp:

Apr 5, 2013, 12:39:30 PM (13 years ago)

Author:

Julian Iseringhausen <julian.iseringhausen@…>

Children:

4a709e

Parents:

f57182

Message:

Merge stashed open boundary stuff.

File:

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

src/samples/discretization_poisson_fv.cpp

@@ -33 +33 @@
 #endif
 
+#include "grid/grid_index_translations.hpp"
 #include "samples/discretization_poisson_fv.hpp"
 
@@ -98 +99 @@
 void DiscretizationPoissonFV::SetInnerBoundaryCompute(Grid& sol_f, Grid& rhs_f, Grid& sol_c) const
 {
+/*
+  //FIXME
   Index i_c, i_f;
 
-  const Vector h2_inv = 0.5 / sol_f.Extent().MeshWidth();
-
+  Index b1_c, b2_c, b1_f, b2_f;
+  GridIndexTranslations::GetInnerBoundary(sol_c, b1_c, b2_c, sol_f, b1_f, b2_f);
+*/
+/*
   const Index b1_c = sol_c.Local().FinerBegin();
   const Index b2_c = sol_c.Local().FinerEnd() - 1;
   const Index b1_f = 0;
   const Index b2_f = rhs_f.Local().SizeTotal() - 1;
-
-  const Index& begin_f = sol_f.Local().Begin();
-  const Index& end_f = sol_f.Local().End();
-  const Index& begin_c = sol_c.Local().FinerBegin();
-
+*/
+
+/*
+  const Index begin_f = sol_f.Local().Begin();
+  const Index end_f = sol_f.Local().End();
+  const Index begin_c = sol_c.Local().FinerBegin();
+
+  const Vector h2_inv = 0.5 / sol_f.Extent().MeshWidth();
   const vmg_float c_1_3 = 1.0 / 3.0;
   const Vector c_2_3_sp = 2.0 / 3.0 * sol_f.Extent().MeshWidth();
@@ -145 +153 @@
 
       rhs_f(b1_f.X(),i_f.Y(),i_f.Z()) = 0.25 * (rhs_f(b1_f.X(),i_f.Y()-1,i_f.Z()-1) +
-                                                rhs_f(b1_f.X(),i_f.Y()+1,i_f.Z()-1) +
-                                                rhs_f(b1_f.X(),i_f.Y()-1,i_f.Z()+1) +
-                                                rhs_f(b1_f.X(),i_f.Y()+1,i_f.Z()+1));
+          rhs_f(b1_f.X(),i_f.Y()+1,i_f.Z()-1) +
+          rhs_f(b1_f.X(),i_f.Y()-1,i_f.Z()+1) +
+          rhs_f(b1_f.X(),i_f.Y()+1,i_f.Z()+1));
 
       rhs_f(b2_f.X(),i_f.Y(),i_f.Z()) = 0.25 * (rhs_f(b2_f.X(),i_f.Y()-1,i_f.Z()-1) +
-                                                rhs_f(b2_f.X(),i_f.Y()+1,i_f.Z()-1) +
-                                                rhs_f(b2_f.X(),i_f.Y()-1,i_f.Z()+1) +
-                                                rhs_f(b2_f.X(),i_f.Y()+1,i_f.Z()+1));
+          rhs_f(b2_f.X(),i_f.Y()+1,i_f.Z()-1) +
+          rhs_f(b2_f.X(),i_f.Y()-1,i_f.Z()+1) +
+          rhs_f(b2_f.X(),i_f.Y()+1,i_f.Z()+1));
     }
 
@@ -159 +167 @@
 
       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());
+          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.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());
+          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());
     }
 
@@ -194 +202 @@
 
       rhs_f(i_f.X(),b1_f.Y(),i_f.Z()) = 0.25 * (rhs_f(i_f.X()-1,b1_f.Y(),i_f.Z()-1) +
-                                                rhs_f(i_f.X()+1,b1_f.Y(),i_f.Z()-1) +
-                                                rhs_f(i_f.X()-1,b1_f.Y(),i_f.Z()+1) +
-                                                rhs_f(i_f.X()+1,b1_f.Y(),i_f.Z()+1));
+          rhs_f(i_f.X()+1,b1_f.Y(),i_f.Z()-1) +
+          rhs_f(i_f.X()-1,b1_f.Y(),i_f.Z()+1) +
+          rhs_f(i_f.X()+1,b1_f.Y(),i_f.Z()+1));
 
       rhs_f(i_f.X(),b2_f.Y(),i_f.Z()) = 0.25 * (rhs_f(i_f.X()-1,b2_f.Y(),i_f.Z()-1) +
-                                                rhs_f(i_f.X()+1,b2_f.Y(),i_f.Z()-1) +
-                                                rhs_f(i_f.X()-1,b2_f.Y(),i_f.Z()+1) +
-                                                rhs_f(i_f.X()+1,b2_f.Y(),i_f.Z()+1));
+          rhs_f(i_f.X()+1,b2_f.Y(),i_f.Z()-1) +
+          rhs_f(i_f.X()-1,b2_f.Y(),i_f.Z()+1) +
+          rhs_f(i_f.X()+1,b2_f.Y(),i_f.Z()+1));
     }
 
@@ -208 +216 @@
 
       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());
+          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.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());
+          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());
     }
 
@@ -242 +250 @@
     for (i_f.Y()=begin_f.Y()+1; i_f.Y()<end_f.Y(); i_f.Y()+=2) {
       rhs_f(i_f.X(),i_f.Y(),b1_f.Z()) = 0.25 * (rhs_f(i_f.X()-1,i_f.Y()-1,b1_f.Z()) +
-                                                rhs_f(i_f.X()+1,i_f.Y()-1,b1_f.Z()) +
-                                                rhs_f(i_f.X()-1,i_f.Y()+1,b1_f.Z()) +
-                                                rhs_f(i_f.X()+1,i_f.Y()+1,b1_f.Z()));
+          rhs_f(i_f.X()+1,i_f.Y()-1,b1_f.Z()) +
+          rhs_f(i_f.X()-1,i_f.Y()+1,b1_f.Z()) +
+          rhs_f(i_f.X()+1,i_f.Y()+1,b1_f.Z()));
 
       rhs_f(i_f.X(),i_f.Y(),b2_f.Z()) = 0.25 * (rhs_f(i_f.X()-1,i_f.Y()-1,b2_f.Z()) +
-                                                rhs_f(i_f.X()+1,i_f.Y()-1,b2_f.Z()) +
-                                                rhs_f(i_f.X()-1,i_f.Y()+1,b2_f.Z()) +
-                                                rhs_f(i_f.X()+1,i_f.Y()+1,b2_f.Z()));
+          rhs_f(i_f.X()+1,i_f.Y()-1,b2_f.Z()) +
+          rhs_f(i_f.X()-1,i_f.Y()+1,b2_f.Z()) +
+          rhs_f(i_f.X()+1,i_f.Y()+1,b2_f.Z()));
     }
 
@@ -256 +264 @@
 
       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);
+          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.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);
-    }
-
+          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);
+    }
+*/
 #ifdef DEBUG_MATRIX_CHECKS
   rhs_f.IsConsistent();