From d1b0e854b603c06b44c736df06ba6ee62711d83c Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Tue, 25 Aug 2026 10:21:35 +0200
Subject: two to two thirds

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

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

diff --git a/modules/core/c++/fplbm.hpp b/modules/core/c++/fplbm.hpp
index 1d2b2d4..d130c8d 100644
--- a/modules/core/c++/fplbm.hpp
+++ b/modules/core/c++/fplbm.hpp
@@ -92,10 +92,9 @@ public:
 		auto& vel_f = macros.template get<"momentum">();
 		saw::data<sch::Vector<T,Descriptor::D>> vel = vel_f.at(index);
 
-		auto eq = equilibrium<T,Descriptor>(rho,vel);
-		vel = vel + force * half / rho;
-
 		// compute_rho_u<T,Descriptor>(dfs,rho,vel);
+		vel = vel + force * half / rho;
+		auto eq = equilibrium<T,Descriptor>(rho,vel);
 
 		using dfi = df_info<T,Descriptor>;
 
@@ -176,6 +175,9 @@ public:
 		}
 		auto& aabb = eo_aabb.get_value();
 
+		saw::data<sch::Scalar<T>> three;
+		three.at({}).set(3.0f);
+
 		saw::data<sch::Scalar<T>> two;
 		two.at({}).set(2.0f);
 
@@ -204,7 +206,7 @@ public:
 
 			// vel_s is technically time the density of the particle?
 
-			force = ( vel_s * rho - vel * rho ) * two * flip_por;
+			force = ( vel_s * rho - vel * rho ) * (two / three) * flip_por;
 
 			force_p = force_p - force;
 		}, aabb.template get<"a">(), aabb.template get<"b">());
-- 
cgit v1.2.3