diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-02-16 14:39:38 +0100 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-02-16 14:39:48 +0100 |
| commit | 92f5645809449f56c39c0e4c6c29045b8a4acea6 (patch) | |
| tree | 36e6a698dc76196cc74fb1450e2216ce7650f4ef /lib/core/c++/hlbm.hpp | |
| parent | 6e0412d790cb5990364531f65e287f9867696da2 (diff) | |
| download | libs-lbm-92f5645809449f56c39c0e4c6c29045b8a4acea6.tar.gz | |
Dangling changes
Diffstat (limited to 'lib/core/c++/hlbm.hpp')
| -rw-r--r-- | lib/core/c++/hlbm.hpp | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/lib/core/c++/hlbm.hpp b/lib/core/c++/hlbm.hpp index 0373125..ea57ef4 100644 --- a/lib/core/c++/hlbm.hpp +++ b/lib/core/c++/hlbm.hpp @@ -50,7 +50,9 @@ public: auto& N = particle_N_f.at(index); auto& D = particle_D_f.at(index); // Convex combination of velocities - vel = vel * porosity + N * flip_porosity / D; + vel = vel * porosity + [&]() -> saw::data<sch::Vector<T,Descriptor::D>> { + return (D.at({}).get() > 0.0 ? N * flip_porosity / D : N); + }(); // Equilibrium auto eq = equilibrium<T,Descriptor>(rho,vel); |