ViennaTools · tobre1 · Jan 2, 2026 · Dec 9, 2025 · Dec 10, 2025 · Dec 10, 2025
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -121,6 +121,7 @@ CPMFindPackage(
 CPMFindPackage(
   NAME ViennaLS
   VERSION 5.3.0
+  GIT_TAG master
   GIT_REPOSITORY "https://github.com/ViennaTools/ViennaLS"
   EXCLUDE_FROM_ALL ${VIENNAPS_BUILD_PYTHON})
 

diff --git a/README.md b/README.md
@@ -53,23 +53,32 @@ ViennaPS is also available on the [Python Package Index (PyPI)](https://pypi.org
 
 * C++17 Compiler with OpenMP support
 
-### Dependencies (installed automatically)
+### ViennaTools Dependencies (installed automatically)
 
-* [ViennaCore](https://github.com/ViennaTools/viennacore) 
+ViennaPS is part of the ViennaTools ecosystem and depends on several lightweight, header-only ViennaTools libraries. During configuration, CMake will look for them and fetch them automatically as part of the ViennaPS build. No separate installation step is required:
 
+* [ViennaCore](https://github.com/ViennaTools/viennacore) 
 * [ViennaLS](https://github.com/ViennaTools/viennals) 
-  * [ViennaHRLE](https://github.com/ViennaTools/viennahrle) 
-  * [VTK](https://vtk.org/) (9.0.0+)
-
+* [ViennaHRLE](https://github.com/ViennaTools/viennahrle) 
 * [ViennaRay](https://github.com/ViennaTools/viennaray) 
-  * [Embree](https://www.embree.org/) (4.0.0+)
-
 * [ViennaCS](https://github.com/ViennaTools/viennacs)
 
-* [pybind11](https://github.com/pybind/pybind11) (3.0.0+, only for building Python libs)
+### External Dependencies
+
+The following external dependencies are required to build ViennaPS. On most systems, installing them via a package manager (e.g. `apt`, `brew`, or `vcpkg`) is the fastest option:
+
+* [VTK](https://vtk.org/) (9.0.0+)
+* [Embree](https://www.embree.org/) (4.0.0+)
+
+CMake automatically checks for these dependencies during configuration. If they are not found, they can be built from source as part of the build.
+
+To prefer a specific local installation, point CMake to it via `VIENNAPS_LOOKUP_DIRS` (a semicolon-separated list of prefixes):
+
+```bash
+cmake -B build -DVIENNAPS_LOOKUP_DIRS="/path/to/vtk;/path/to/embree"
+```
 
-The CMake configuration automatically checks if the dependencies are installed.
-If the dependencies are not found on the system, they will be built from source. To use local installations of the dependencies, the `VIENNAPS_LOOKUP_DIRS` variable can be set to the installation path of the dependencies.
+Alternatively (or additionally), you can use `CMAKE_PREFIX_PATH` if that better matches your local setup.
 
 ## Installing
 

diff --git a/examples/DRAMWiggling/DRAMWiggling.cpp b/examples/DRAMWiggling/DRAMWiggling.cpp
@@ -60,8 +60,8 @@ int main(int argc, char **argv) {
 
   // Advection parameters
   AdvectionParameters advectionParams;
-  advectionParams.integrationScheme = util::convertIntegrationScheme(
-      params.get<std::string>("integrationScheme"));
+  advectionParams.spatialScheme =
+      util::convertSpatialScheme(params.get<std::string>("spatialScheme"));
 
   RayTracingParameters rayParams;
   rayParams.raysPerPoint = params.get<int>("raysPerPoint");

diff --git a/examples/DRAMWiggling/DRAMWiggling.py b/examples/DRAMWiggling/DRAMWiggling.py
@@ -57,8 +57,8 @@
 rayTracingParams.raysPerPoint = int(params["raysPerPoint"])
 
 advectionParams = ps.AdvectionParameters()
-advectionParams.integrationScheme = ps.util.convertIntegrationScheme(
-    params["integrationScheme"]
+advectionParams.spatialScheme = ps.util.convertSpatialScheme(
+    params["spatialScheme"]
 )
 
 # process setup

diff --git a/examples/DRAMWiggling/config.txt b/examples/DRAMWiggling/config.txt
@@ -17,7 +17,7 @@ ionExponent=1000
 meanEnergy=200 # eV
 sigmaEnergy=10 # eV
 
-integrationScheme=LF_2
+spatialScheme=LF_2
 
 numSteps=100
 raysPerPoint=1000

diff --git a/examples/blazedGratingsEtching/blazedGratingsEtching.cpp b/examples/blazedGratingsEtching/blazedGratingsEtching.cpp
@@ -32,8 +32,8 @@ int main(int argc, char **argv) {
   geometry->saveSurfaceMesh("initial");
 
   AdvectionParameters advParams;
-  advParams.integrationScheme =
-      viennals::IntegrationSchemeEnum::LAX_FRIEDRICHS_2ND_ORDER;
+  advParams.spatialScheme =
+      viennals::SpatialSchemeEnum::LAX_FRIEDRICHS_2ND_ORDER;
   advParams.timeStepRatio = 0.25;
 
   RayTracingParameters rayTracingParams;
@@ -83,4 +83,4 @@ int main(int argc, char **argv) {
     process.apply();
     geometry->saveSurfaceMesh("BlazedGratingsEtch_P" + std::to_string(i));
   }
-}
+}
diff --git a/examples/blazedGratingsEtching/blazedGratingsEtching.py b/examples/blazedGratingsEtching/blazedGratingsEtching.py
@@ -54,7 +54,7 @@
 
 # ----- Model Setup ----- #
 advectionParams = ps.AdvectionParameters()
-advectionParams.integrationScheme = ps.IntegrationScheme.LAX_FRIEDRICHS_2ND_ORDER
+advectionParams.spatialScheme = ps.SpatialScheme.LAX_FRIEDRICHS_2ND_ORDER
 advectionParams.timeStepRatio = 0.25
 
 rayTracingParams = ps.RayTracingParameters()

diff --git a/examples/boschProcess/boschProcessEmulate.py b/examples/boschProcess/boschProcessEmulate.py
@@ -38,9 +38,7 @@
 direction[args.dim - 1] = -1.0
 
 # Geometric advection model for deposition
-depoModel = ps.SphereDistribution(
-    radius=params["depositionThickness"], gridDelta=params["gridDelta"]
-)
+depoModel = ps.SphereDistribution(radius=params["depositionThickness"])
 
 # Define purely directional rate for depo removal
 etchDir = ps.RateSet(

diff --git a/examples/cantileverWetEtching/cantileverWetEtching.cpp b/examples/cantileverWetEtching/cantileverWetEtching.cpp
@@ -71,8 +71,8 @@ int main(int argc, char **argv) {
           {ps::Material::Si, -1.}});
 
   ps::AdvectionParameters advectionParams;
-  advectionParams.integrationScheme =
-      viennals::IntegrationSchemeEnum::STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
+  advectionParams.spatialScheme =
+      viennals::SpatialSchemeEnum::STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
 
   ps::Process<NumericType, D> process;
   process.setDomain(geometry);

diff --git a/examples/cantileverWetEtching/cantileverWetEtching.py b/examples/cantileverWetEtching/cantileverWetEtching.py
@@ -1,6 +1,6 @@
 # This example only works in 3D mode
 import viennaps.d3 as psd
-from viennaps import BoundaryType, Material, AdvectionParameters, IntegrationScheme
+from viennaps import BoundaryType, Material, AdvectionParameters, SpatialScheme
 
 maskFileName = "cantilever_mask.gds"
 
@@ -56,8 +56,8 @@
 )
 
 advectionParams = AdvectionParameters()
-advectionParams.integrationScheme = (
-    IntegrationScheme.STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER
+advectionParams.spatialScheme = (
+    SpatialScheme.STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER
 )
 
 process = psd.Process()

diff --git a/examples/emulation/FinFET.cpp b/examples/emulation/FinFET.cpp
@@ -220,8 +220,8 @@ int main() {
     auto model =
         SmartPointer<IsotropicProcess<NumericType, D>>::New(-1., masks);
     AdvectionParameters advectionParams;
-    advectionParams.integrationScheme =
-        viennals::IntegrationSchemeEnum::LAX_FRIEDRICHS_2ND_ORDER;
+    advectionParams.spatialScheme =
+        viennals::SpatialSchemeEnum::LAX_FRIEDRICHS_2ND_ORDER;
     Process<NumericType, D> process(domain, model, 5.);
     process.setParameters(advectionParams);
     process.apply();
@@ -236,8 +236,8 @@ int main() {
     domain->duplicateTopLevelSet(Material::SiGe);
     Logger::setLogLevel(LogLevel::INFO);
     AdvectionParameters advectionParams;
-    advectionParams.integrationScheme =
-        viennals::IntegrationSchemeEnum::STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
+    advectionParams.spatialScheme =
+        viennals::SpatialSchemeEnum::STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
     lsInternal::StencilLocalLaxFriedrichsScalar<NumericType, D,
                                                 1>::setMaxDissipation(100);
 

diff --git a/examples/emulation/FinFET.py b/examples/emulation/FinFET.py
@@ -183,7 +183,7 @@ def writeSurface(domain):
 masks = [ps.Material.PolySi, ps.Material.SiO2, ps.Material.Si3N4]
 model = ps.IsotropicProcess(rate=-1.0, maskMaterial=masks)
 advParams = ps.AdvectionParameters()
-advParams.integrationScheme = ps.IntegrationScheme.LAX_FRIEDRICHS_2ND_ORDER
+advParams.spatialScheme = ps.SpatialScheme.LAX_FRIEDRICHS_2ND_ORDER
 process = ps.Process(domain, model, 5.0)
 process.setParameters(advParams)
 process.apply()
@@ -195,8 +195,8 @@ def writeSurface(domain):
 print("S/D Epitaxy ...", end="", flush=True)
 domain.duplicateTopLevelSet(ps.Material.SiGe)
 advectionParams = ps.AdvectionParameters()
-advectionParams.integrationScheme = (
-    ps.IntegrationScheme.STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER
+advectionParams.spatialScheme = (
+    ps.SpatialScheme.STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER
 )
 ps.StencilLocalLaxFriedrichsScalar.setMaxDissipation(1000)
 material = [

diff --git a/examples/faradayCageEtching/faradayCageEtching.cpp b/examples/faradayCageEtching/faradayCageEtching.cpp
@@ -46,8 +46,8 @@ int main(int argc, char *argv[]) {
       cageParams, maskMaterials);
 
   ps::AdvectionParameters advectionParams;
-  advectionParams.integrationScheme =
-      ps::IntegrationScheme::LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
+  advectionParams.spatialScheme =
+      ps::SpatialScheme::LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
 
   ps::RayTracingParameters rayParams;
   rayParams.raysPerPoint = params.get<int>("raysPerPoint");

diff --git a/examples/holeEtching/config.txt b/examples/holeEtching/config.txt
@@ -28,7 +28,8 @@ A_O=2          # passivation layer yield constant
 A_Si=7         # Si yield constant
 
 etchStopDepth=-10   # maximum etching depth
-integrationScheme=EO_1
+spatialScheme=EO_1
+temporalScheme=RK3
 
 raysPerPoint=1000
 

diff --git a/examples/holeEtching/holeEtching.cpp b/examples/holeEtching/holeEtching.cpp
@@ -60,8 +60,10 @@ int main(int argc, char *argv[]) {
   rayTracingParams.raysPerPoint = params.get<unsigned>("raysPerPoint");
 
   AdvectionParameters advectionParams;
-  advectionParams.integrationScheme = util::convertIntegrationScheme(
-      params.get<std::string>("integrationScheme"));
+  advectionParams.spatialScheme =
+      util::convertSpatialScheme(params.get<std::string>("spatialScheme"));
+  advectionParams.temporalScheme =
+      util::convertTemporalScheme(params.get<std::string>("temporalScheme"));
 
   // process setup
   Process<NumericType, D> process(geometry, model);

diff --git a/examples/holeEtching/holeEtching.py b/examples/holeEtching/holeEtching.py
@@ -55,8 +55,8 @@
 rayParams.smoothingNeighbors = 2
 
 advParams = ps.AdvectionParameters()
-advParams.integrationScheme = ps.util.convertIntegrationScheme(
-    params["integrationScheme"]
+advParams.spatialScheme = ps.util.convertSpatialScheme(
+    params["spatialScheme"]
 )
 
 # process setup

diff --git a/examples/holeEtching/testFluxes.py b/examples/holeEtching/testFluxes.py
@@ -57,7 +57,7 @@
 processDuration = 3  # min
 
 advParams = ps.AdvectionParameters()
-advParams.integrationScheme = ps.IntegrationScheme.ENGQUIST_OSHER_2ND_ORDER
+advParams.spatialScheme = ps.SpatialScheme.ENGQUIST_OSHER_2ND_ORDER
 
 rayParams = ps.RayTracingParameters()
 rayParams.raysPerPoint = int(1000)

diff --git a/examples/ionBeamEtching/ionBeamEtching.cpp b/examples/ionBeamEtching/ionBeamEtching.cpp
@@ -57,8 +57,8 @@ int main(int argc, char *argv[]) {
   model->setPrimaryDirection(direction);
 
   AdvectionParameters advectionParams;
-  advectionParams.integrationScheme =
-      viennals::IntegrationSchemeEnum::LAX_FRIEDRICHS_2ND_ORDER;
+  advectionParams.spatialScheme =
+      viennals::SpatialSchemeEnum::LAX_FRIEDRICHS_2ND_ORDER;
 
   Process<NumericType, D> process(geometry, model);
   process.setProcessDuration(params.get("processTime"));

diff --git a/examples/ionBeamEtching/ionBeamEtching.py b/examples/ionBeamEtching/ionBeamEtching.py
@@ -59,7 +59,7 @@
 model.setPrimaryDirection(direction)
 
 advParams = ps.AdvectionParameters()
-advParams.integrationScheme = ps.IntegrationScheme.LAX_FRIEDRICHS_2ND_ORDER
+advParams.spatialScheme = ps.SpatialScheme.LAX_FRIEDRICHS_2ND_ORDER
 
 process = ps.Process(geometry, model)
 process.setProcessDuration(params["processTime"])

diff --git a/examples/selectiveEpitaxy/selectiveEpitaxy.cpp b/examples/selectiveEpitaxy/selectiveEpitaxy.cpp
@@ -48,8 +48,9 @@ int main(int argc, char *argv[]) {
       params.get("R111"), params.get("R100"));
 
   ps::AdvectionParameters advectionParams;
-  advectionParams.integrationScheme =
-      viennals::IntegrationSchemeEnum::STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
+  advectionParams.integrationScheme = viennals::IntegrationSchemeEnum::
+      // advectionParams.spatialScheme = viennals::SpatialSchemeEnum::
+      STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
   // advectionParams.velocityOutput = true;
   lsInternal::StencilLocalLaxFriedrichsScalar<NumericType, D,
                                               1>::setMaxDissipation(1000);

diff --git a/examples/selectiveEpitaxy/selectiveEpitaxy.py b/examples/selectiveEpitaxy/selectiveEpitaxy.py
@@ -48,8 +48,8 @@
 )
 
 advectionParams = ps.AdvectionParameters()
-advectionParams.integrationScheme = (
-    ps.IntegrationScheme.STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER
+advectionParams.spatialScheme = (
+    ps.SpatialScheme.STENCIL_LOCAL_LAX_FRIEDRICHS_1ST_ORDER
 )
 
 process = ps.Process(geometry, model, params["processTime"])

diff --git a/examples/stackEtching/stackEtching.cpp b/examples/stackEtching/stackEtching.cpp
@@ -70,8 +70,7 @@ int main(int argc, char *argv[]) {
       SmartPointer<FluorocarbonEtching<NumericType, D>>::New(parameters);
 
   AdvectionParameters advectionParams;
-  advectionParams.integrationScheme =
-      IntegrationScheme::LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
+  advectionParams.spatialScheme = SpatialScheme::LOCAL_LAX_FRIEDRICHS_1ST_ORDER;
   advectionParams.timeStepRatio = 0.25;
 
   CoverageParameters coverageParams;

diff --git a/examples/stackEtching/stackEtching.py b/examples/stackEtching/stackEtching.py
@@ -83,7 +83,7 @@
 covParams.tolerance = 1e-4
 
 advParams = vps.AdvectionParameters()
-advParams.integrationScheme = vps.IntegrationScheme.LOCAL_LAX_FRIEDRICHS_1ST_ORDER
+advParams.spatialScheme = vps.SpatialScheme.LOCAL_LAX_FRIEDRICHS_1ST_ORDER
 
 process = vps.d2.Process()
 process.setDomain(geometry)

diff --git a/examples/trenchDepositionGeometric/trenchDepositionGeometric.py b/examples/trenchDepositionGeometric/trenchDepositionGeometric.py
@@ -35,9 +35,7 @@
 # copy top layer to capture deposition
 geometry.duplicateTopLevelSet(ps.Material.SiO2)
 
-model = ps.SphereDistribution(
-    radius=params["layerThickness"], gridDelta=params["gridDelta"]
-)
+model = ps.SphereDistribution(radius=params["layerThickness"])
 
 geometry.saveHullMesh("initial")
 

diff --git a/include/viennaps/process/psALPStrategy.hpp b/include/viennaps/process/psALPStrategy.hpp
@@ -93,7 +93,9 @@ class ALPStrategy final : public ProcessStrategy<NumericType, D> {
 
     // Initialize translation field. Converts points ids from level set points
     // to surface points
-    const int translationMethod = context.needsExtendedVelocities() ? 2 : 1;
+    /// TODO: find out why translationMethod 1 does not work properly
+    const int translationMethod =
+        2; //  context.needsExtendedVelocities() ? 2 : 1;
     context.translationField =
         SmartPointer<TranslationField<NumericType, D>>::New(
             context.model->getVelocityField(), context.domain->getMaterialMap(),
@@ -148,7 +150,7 @@ class ALPStrategy final : public ProcessStrategy<NumericType, D> {
       VIENNACORE_LOG_INFO("Cycle: " + std::to_string(cycle + 1) + "/" +
                           std::to_string(numCycles));
 
-      // Prepare advection (expand level set based on integration scheme)
+      // Prepare advection (expand level set based on discretization scheme)
       advectionHandler_.prepareAdvection(context);
 
       updateState(context);
@@ -322,4 +324,4 @@ class ALPStrategy final : public ProcessStrategy<NumericType, D> {
   }
 };
 
-} // namespace viennaps
+} // namespace viennaps