From 3af0e20773b21b519ea474e34e7c1df73d04a307 Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Tue, 24 Mar 2026 17:41:28 +0100
Subject: Dangling changes

---
 examples/settling_cubes_2d_ibm_gpu/sim.cpp |  7 +++++++
 lib/core/c++/collision.hpp                 | 15 +++++++++++++--
 lib/core/c++/particle/particle.hpp         |  5 ++++-
 3 files changed, 24 insertions(+), 3 deletions(-)

diff --git a/examples/settling_cubes_2d_ibm_gpu/sim.cpp b/examples/settling_cubes_2d_ibm_gpu/sim.cpp
index 4501a54..05a4a09 100644
--- a/examples/settling_cubes_2d_ibm_gpu/sim.cpp
+++ b/examples/settling_cubes_2d_ibm_gpu/sim.cpp
@@ -102,6 +102,11 @@ saw::error_or<void> setup_initial_conditions(
 	return saw::make_void();
 }
 
+template<typename T, typename Desc>
+void add_particles(saw::data<sch::Array<sch::Particle<T,Desc::D>>>& particles){
+	////// TODO
+}
+
 template<typename T, typename Desc>
 saw::error_or<void> step(
 		saw::data<sch::Ptr<sch::ChunkStruct<T,Desc>>,encode::Sycl<saw::encode::Native>>& fields,
@@ -148,6 +153,8 @@ saw::error_or<void> step(
 				break;
 			}
 		});
+
+
 	}).wait();
 
 	
diff --git a/lib/core/c++/collision.hpp b/lib/core/c++/collision.hpp
index f8312c0..df45bbe 100644
--- a/lib/core/c++/collision.hpp
+++ b/lib/core/c++/collision.hpp
@@ -95,10 +95,21 @@ class component<T, Descriptor, cmpt::BGKGuo, Encode> {
 private:
 	typename saw::native_data_type<T>::type relaxation_;
 	saw::data<T> frequency_;
+	saw::data<sch::Vector<T,Descriptor::D>> global_force_;
 public:
 	component(typename saw::native_data_type<T>::type relaxation__):
 		relaxation_{relaxation__},
-		frequency_{typename saw::native_data_type<T>::type(1) / relaxation_}
+		frequency_{typename saw::native_data_type<T>::type(1) / relaxation_},
+		global_force_{}
+	{}
+
+	component(
+		typename saw::native_data_type<T>::type relaxation__,
+		saw::data<sch::Vector<T,Descriptor::D>>& global_force__
+	):
+		relaxation_{relaxation__},
+		frequency_{typename saw::native_data_type<T>::type(1) / relaxation_},
+		global_force_{global_force__}
 	{}
 
 	using Component = cmpt::BGKGuo;
@@ -153,7 +164,7 @@ public:
 			saw::data<sch::Scalar<T>> w;
 			w.at({}).set(dfi::weights[i]);
 
-			auto F_i_d = saw::math::dot(force * w, (ci - vel * dfi_inv_cs2 + ci * ci_dot_u * dfi_inv_cs2 * dfi_inv_cs2 ));
+			auto F_i_d = saw::math::dot((force+global_force_) * w, (ci - vel * dfi_inv_cs2 + ci * ci_dot_u * dfi_inv_cs2 * dfi_inv_cs2 ));
 			/*
 			saw::data<sch::Scalar<T>> F_i_sum;
 			for(uint64_t d = 0u; d < Descriptor::D; ++d){
diff --git a/lib/core/c++/particle/particle.hpp b/lib/core/c++/particle/particle.hpp
index 8ef590d..f8104fd 100644
--- a/lib/core/c++/particle/particle.hpp
+++ b/lib/core/c++/particle/particle.hpp
@@ -43,12 +43,15 @@ using ParticleCollisionSpheroid = Struct<
 	Member<Scalar<T>, "radius">
 >;
 
+tem
+
 template<typename T, uint64_t D, typename CollisionType = ParticleCollisionSpheroid<T>>
 using Particle = Struct<
 	Member<ParticleRigidBody<T,D>, "rigid_body">,
 	Member<CollisionType, "collision">,
 	// Problem is that dynamic data would two layered
-	// Member<FixedArray<Float64,D,D,D>, "mask">, 
+	// Member<Array<Float64,D>, "mask">, 
+	Member<UInt64, "mask_id">,
 	Member<Scalar<T>, "mass">
 >;
 }
-- 
cgit v1.2.3