From 92f5645809449f56c39c0e4c6c29045b8a4acea6 Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Mon, 16 Feb 2026 14:39:38 +0100
Subject: Dangling changes

---
 lib/core/c++/geometry/poiseulle_channel.hpp | 7 +++++++
 lib/core/c++/hlbm.hpp                       | 4 +++-
 2 files changed, 10 insertions(+), 1 deletion(-)
 create mode 100644 lib/core/c++/geometry/poiseulle_channel.hpp

(limited to 'lib')

diff --git a/lib/core/c++/geometry/poiseulle_channel.hpp b/lib/core/c++/geometry/poiseulle_channel.hpp
new file mode 100644
index 0000000..f675a99
--- /dev/null
+++ b/lib/core/c++/geometry/poiseulle_channel.hpp
@@ -0,0 +1,7 @@
+#pragma once
+
+namespace kel {
+namespace lbm {
+
+}
+}
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);
 
-- 
cgit v1.2.3