From 62c7dcd0c508ca9559e923a023181600196bf5ca Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Thu, 12 Feb 2026 21:43:44 +0100
Subject: HLBM kinda broken again?

---
 lib/core/c++/hlbm.hpp | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

(limited to 'lib/core/c++')

diff --git a/lib/core/c++/hlbm.hpp b/lib/core/c++/hlbm.hpp
index 8f267d5..0373125 100644
--- a/lib/core/c++/hlbm.hpp
+++ b/lib/core/c++/hlbm.hpp
@@ -40,15 +40,17 @@ public:
 		saw::data<sch::Scalar<T>>& rho = rho_f.at(index);
 		saw::data<sch::Vector<T,Descriptor::D>> vel = vel_f.at(index);
 
-		compute_rho_u<T,Descriptor::D>(dfs_old_f.at(index), rho, vel);
+		compute_rho_u<T,Descriptor>(dfs_old_f.at(index), rho, vel);
 
 		auto& porosity = porosity_f.at(index);
-		auto flip_porosity = saw::data<sch::Scalar<T>>{1.0} - porosity;
+		saw::data<sch::Scalar<T>> one;
+		one.at({}) = 1.0;
+		auto flip_porosity = one - porosity;
 
 		auto& N = particle_N_f.at(index);
 		auto& D = particle_D_f.at(index);
 		// Convex combination of velocities
-		vel = vel * flip_porosity + N * porosity / D;
+		vel = vel * porosity + N * flip_porosity / D;
 		// Equilibrium
 		auto eq = equilibrium<T,Descriptor>(rho,vel);
 
@@ -57,7 +59,7 @@ public:
 		}
 
 		porosity.at({}) = 1.0;
-
+		D.at({}) = 0.0;
 	}
 };
 
-- 
cgit v1.2.3