Set oxide regrowth Dim 2

ViennaTools · Nov 6, 2023 · 8f106e5 · 8f106e5
1 parent 099b890
commit 8f106e5
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 4 deletions.
diff --git a/Examples/OxideRegrowth/OxideRegrowth.cpp b/Examples/OxideRegrowth/OxideRegrowth.cpp
@@ -10,7 +10,7 @@ int main(int argc, char **argv) {
   using NumericType = double;
 
   omp_set_num_threads(12);
-  constexpr int D = 3;
+  constexpr int D = 2;
 
   psLogger::setLogLevel(psLogLevel::INTERMEDIATE);
 
@@ -49,8 +49,11 @@ int main(int argc, char **argv) {
   auto &cellSet = domain->getCellSet();
   cellSet->addScalarData("byproductSum", 0.);
   cellSet->writeVTU("initial.vtu");
-  if constexpr (D == 3)
-    cellSet->setPeriodicBoundary({false, true, false});
+  if constexpr (D == 3) {
+    std::array<bool, D> boundaryConds = {false};
+    boundaryConds[1] = true;
+    cellSet->setPeriodicBoundary(boundaryConds);
+  }
   // we need neighborhood information for solving the
   // convection-diffusion equation on the cell set
   cellSet->buildNeighborhood();

diff --git a/Examples/OxideRegrowth/OxideRegrowth.py b/Examples/OxideRegrowth/OxideRegrowth.py
@@ -1,4 +1,4 @@
-DIM = 3
+DIM = 2
 
 if DIM is 3:
     stabFac = 0.145
@@ -43,10 +43,13 @@
 cellSet = domain.getCellSet()
 cellSet.addScalarData("byproductSum", 0.0)
 cellSet.writeVTU("initial.vtu")
+if DIM is 3:
+    cellSet.setPeriodicBoundary([False, True, False])
 # we need neighborhood information for solving the
 # convection-diffusion equation on the cell set
 cellSet.buildNeighborhood()
 
+
 # The redeposition model captures byproducts from the selective etching
 # process in the cell set. The byproducts are then distributed by solving a
 # convection-diffusion equation on the cell set.