summaryrefslogtreecommitdiff
path: root/examples
diff options
context:
space:
mode:
authorClaudius "keldu" Holeksa <mail@keldu.de>2026-07-08 18:55:16 +0200
committerClaudius "keldu" Holeksa <mail@keldu.de>2026-07-08 18:55:16 +0200
commit14dd342782540607cdc08f718700cd7917c017fc (patch)
treea06599af8e0d2f6cda57d267e3c7c858c5d07891 /examples
parente67393d7fcc6e9aa1dfe10e3a6a0d087b9621678 (diff)
downloadlibs-lbm-14dd342782540607cdc08f718700cd7917c017fc.tar.gz
Need that, well. no clue how to fix it though :/
Diffstat (limited to 'examples')
-rw-r--r--examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp1
-rw-r--r--examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp1
-rw-r--r--examples/stokes_drag_particle_2d_psm_gpu/sim.cpp1
3 files changed, 3 insertions, 0 deletions
diff --git a/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp
index 0d6ebec..9e96219 100644
--- a/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp
+++ b/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp
@@ -135,6 +135,7 @@ saw::error_or<void> setup_initial_conditions(
por_f.at(index).at({}) = {0};
rho.at({}) = {1};
auto& vel = vel_f.at(index);
+ vel.at({{0u}}) = 0.0002f;
auto eq = equilibrium<T,Desc>(rho,vel);
df = eq;
diff --git a/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp
index 6166aa1..9c28d88 100644
--- a/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp
+++ b/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp
@@ -132,6 +132,7 @@ saw::error_or<void> setup_initial_conditions(
por_f.at(index).at({}) = {1};
rho.at({}) = {1};
auto& vel = vel_f.at(index);
+ vel.at({{0u}}) = 0.0002f;
auto eq = equilibrium<T,Desc>(rho,vel);
df = eq;
diff --git a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp
index 3220a63..75cbe54 100644
--- a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp
+++ b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp
@@ -134,6 +134,7 @@ saw::error_or<void> setup_initial_conditions(
por_f.at(index).at({}) = {1};
rho.at({}) = {1};
auto& vel = vel_f.at(index);
+ vel.at({{0u}}) = 0.0002f;
auto eq = equilibrium<T,Desc>(rho,vel);