From 616b0011c3638eceaae8890c3a66037625587f1d Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Wed, 3 Sep 2025 18:04:36 +0200
Subject: Error in guo forcing seems to persist

Some debug output left there, but there's no force. Might be that the
higher precision is impacting this
---
 examples/cavity_2d_gpu.cpp                  | 1 -
 examples/poiseulle_particles_channel_2d.cpp | 6 ++++--
 2 files changed, 4 insertions(+), 3 deletions(-)

(limited to 'examples')

diff --git a/examples/cavity_2d_gpu.cpp b/examples/cavity_2d_gpu.cpp
index b24ca38..19563e2 100644
--- a/examples/cavity_2d_gpu.cpp
+++ b/examples/cavity_2d_gpu.cpp
@@ -199,7 +199,6 @@ void lbm_step(
 	*/
 	component<sch::T, sch::D2Q9, cmpt::BGK> coll{0.59};
 	component<sch::T, sch::D2Q9, cmpt::BounceBack> bb;
-
 	component<sch::T, sch::D2Q9, cmpt::MovingWall> bb_lid;
 	bb_lid.lid_vel = {0.1,0.0};
 
diff --git a/examples/poiseulle_particles_channel_2d.cpp b/examples/poiseulle_particles_channel_2d.cpp
index e495043..fed78b2 100644
--- a/examples/poiseulle_particles_channel_2d.cpp
+++ b/examples/poiseulle_particles_channel_2d.cpp
@@ -643,7 +643,9 @@ int main(int argc, char** argv){
 				compute_rho_u<sch::T,sch::D2Q9>(dfs,rho,vel);
 				rho = rho * saw::data<sch::T>{dfi::cs2};
 
-				force =
+				for(uint64_t d = 0u; d < sch::D2Q9::D; ++d){
+					force.at({d}) = cell.template get<"force">()({d});
+				}
 
 				part_mask.set(0u);
 			}, {{0u,0u}}, meta);
@@ -683,7 +685,7 @@ int main(int argc, char** argv){
 					// Cast down to get lower corner.
 					// Before casting shift by 0.5 for closest pick
 					saw::data<sch::FixedArray<sch::UInt64,2u>> p_cell_pos {{
-						static_cast<uint64_t>(p_pos_lie.at(saw::data<T>{{0u}}).get()+0.5),
+						static_cast<uint64_t>(p_pos_lie.at({{0u}}).get()+0.5),
 						static_cast<uint64_t>(p_pos_lie.at({{1u}}).get()+0.5)
 					}};
 
-- 
cgit v1.2.3