perf(geo model): only adjust bottom layers if etching

Author: tobre1

include/viennaps/models/psGeometricDistributionModels.hpp

@@ -20,6 +20,7 @@ class SphereDistribution : public ProcessModelCPU<NumericType, D> {
 
     auto geomModel =
         SmartPointer<GeometricModel<NumericType, D>>::New(dist, mask);
+    geomModel->setDeposition(radius > 0);
 
     this->setGeometricModel(geomModel);
     this->setProcessName("SphereDistribution");
@@ -45,6 +46,13 @@ class BoxDistribution : public ProcessModelCPU<NumericType, D> {
 
     auto geomModel =
         SmartPointer<GeometricModel<NumericType, D>>::New(dist, mask);
+    geomModel->setDeposition(true);
+    for (const auto &halfAxis : halfAxes) {
+      if (halfAxis < 0) {
+        geomModel->setDeposition(false);
+        break;
+      }
+    }
 
     this->setGeometricModel(geomModel);
     this->setProcessName("BoxDistribution");
@@ -73,6 +81,13 @@ class CustomSphereDistribution : public ProcessModelCPU<NumericType, D> {
 
     auto geomModel =
         SmartPointer<GeometricModel<NumericType, D>>::New(dist, mask);
+    geomModel->setDeposition(true);
+    for (const auto &radius : radii) {
+      if (radius < 0) {
+        geomModel->setDeposition(false);
+        break;
+      }
+    }
 
     this->setGeometricModel(geomModel);
     this->setProcessName("CustomSphereDistribution");

include/viennaps/process/psGeometricModel.hpp

@@ -15,6 +15,7 @@ VIENNAPS_TEMPLATE_ND(NumericType, D) class GeometricModel {
       nullptr;
   SmartPointer<viennals::Domain<NumericType, D>> mask = nullptr;
   std::vector<Material> maskMaterials;
+  bool isDepo = false;
 
 public:
   GeometricModel() = default;
@@ -46,6 +47,8 @@ VIENNAPS_TEMPLATE_ND(NumericType, D) class GeometricModel {
   auto &getDistribution() const { return dist; }
   auto &getMask() const { return mask; }
   auto &getMaskMaterials() const { return maskMaterials; }
+  auto isDeposition() const { return isDepo; }
+  void setDeposition(bool deposition) { isDepo = deposition; }
 };
 
 } // namespace viennaps

include/viennaps/process/psGeometricProcessStrategy.hpp

@@ -6,6 +6,8 @@ namespace viennaps {
 
 VIENNAPS_TEMPLATE_ND(NumericType, D)
 class GeometricProcessStrategy : public ProcessStrategy<NumericType, D> {
+  using GeometricKernel = viennals::GeometricAdvect<NumericType, D>;
+
 public:
   DEFINE_CLASS_NAME(GeometricProcessStrategy)
 
@@ -94,17 +96,19 @@ class GeometricProcessStrategy : public ProcessStrategy<NumericType, D> {
       }
     }
 
-    viennals::GeometricAdvect<NumericType, D>(
-        context.domain->getLevelSets().back(), dist, mask)
-        .apply();
-
-    // Intersect all other level sets with the last one to keep them consistent
-    for (int i = context.domain->getNumberOfLevelSets() - 1; i >= 0; --i) {
-      viennals::BooleanOperation<NumericType, D>(
-          context.domain->getLevelSets()[i],
-          context.domain->getLevelSets().back(),
-          viennals::BooleanOperationEnum::INTERSECT)
-          .apply();
+    GeometricKernel(context.domain->getSurface(), dist, mask).apply();
+
+    if (!geometricModel->isDeposition()) {
+      // If the top level set was etched, we need to make sure that the
+      // other level sets are still consistent by intersecting them with the
+      // last level set.
+      for (int i = context.domain->getNumberOfLevelSets() - 1; i >= 0; --i) {
+        viennals::BooleanOperation<NumericType, D>(
+            context.domain->getLevelSets()[i],
+            context.domain->getLevelSets().back(),
+            viennals::BooleanOperationEnum::INTERSECT)
+            .apply();
+      }
     }
 
     return ProcessResult::SUCCESS;