From 80deace4d4b2c93945feeec42a8269f3ba2443d2 Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Tue, 4 Aug 2026 19:52:08 +0200
Subject: Asking myself how I had a stable implementation before

---
 examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp | 11 ++++++-----
 1 file changed, 6 insertions(+), 5 deletions(-)

(limited to 'examples/moving_poiseulle_particles_2d_fplbm_gpu')

diff --git a/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp b/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp
index f6c09e9..09b2616 100644
--- a/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp
+++ b/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp
@@ -141,7 +141,7 @@ saw::error_or<void> init(
 	rad.at({}) = dim_y * 0.1f;
 	saw::data<sch::Scalar<T>> dense;
 	dense.at({}) = 1.0f;
-	particles = create_spheroid_particle_group<T,Desc::D,particle_amount>(rad,dense,{16u});
+	particles = create_spheroid_particle_group<T,Desc::D,particle_amount>(rad,dense,{1u});
 	{
 		auto& rb = particles.template get<"particles">().at({0u}).template get<"rigid_body">();
 		auto& rbp = rb.template get<"position">();
@@ -179,7 +179,7 @@ saw::error_or<void> step(
 			for(uint64_t i = 0u; i < Desc::D; ++i){
 				index.at({{i}}).set(idx[i]);
 			}
-
+			
 			fplbm_reset.apply(fields,macros,index,t_i);
 		});
 	}).wait();
@@ -193,12 +193,13 @@ saw::error_or<void> step(
 				index.at({{i}}).set(idx[i]);
 			}
 
-			fplbm_one_part.apply(fields,macros,particles,index,t_i,{1u});
+			fplbm_one_part.apply(fields,macros,particles,index,t_i,{16u});
 		});
 	}).wait();
 
 	// auto coll_ev = 
 	q.submit([&](acpp::sycl::handler& h){
+		component<T,Desc,cmpt::BGK,encode::Sycl<saw::encode::Native>> bgk{0.8};
 		component<T,Desc,cmpt::FpLbm,encode::Sycl<saw::encode::Native>> collision{0.8};
 		component<T,Desc,cmpt::BounceBack,encode::Sycl<saw::encode::Native>> bb;
 
@@ -232,11 +233,11 @@ saw::error_or<void> step(
 					break;
 				case 3u:
 					flow_in.apply(fields,index,t_i);
-					collision.apply(fields,macros,index,t_i);
+					bgk.apply(fields,macros,index,t_i);
 					break;
 				case 4u:
 					flow_out.apply(fields,index,t_i);
-					collision.apply(fields,macros,index,t_i);
+					bgk.apply(fields,macros,index,t_i);
 				break;
 				default:
 				break;
-- 
cgit v1.2.3


From a0767a86817055744230239284441fdbb024fab5 Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Wed, 5 Aug 2026 17:13:19 +0200
Subject: WOrking on psm 2way

---
 examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'examples/moving_poiseulle_particles_2d_fplbm_gpu')

diff --git a/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp b/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp
index 09b2616..c82f632 100644
--- a/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp
+++ b/examples/moving_poiseulle_particles_2d_fplbm_gpu/sim.cpp
@@ -185,7 +185,7 @@ saw::error_or<void> step(
 	}).wait();
 
 	q.submit([&](acpp::sycl::handler& h){
-		component<T,Desc,cmpt::FpLbmOneParticle,encode::Sycl<saw::encode::Native>> fplbm_one_part;
+		component<T,Desc,cmpt::FpLbmOneParticleImplicit,encode::Sycl<saw::encode::Native>> fplbm_one_part;
 
 		h.parallel_for(acpp::sycl::range<1u>{particle_amount}, [=](acpp::sycl::id<1u> idx){
 			saw::data<sch::FixedArray<sch::UInt64,1u>> index;
@@ -262,7 +262,7 @@ saw::error_or<void> lbm_main(int argc, char** argv){
 	}
 	auto& lbm_dir = eo_lbm_dir.get_value();
 
-	auto out_dir = lbm_dir / "moving_poiseulle_particles_2d_fplbm_gpu";
+	auto out_dir = lbm_dir / "moving_poiseulle_particles_2d_fplbm_imp_gpu";
 
 	{
 		std::error_code ec;
-- 
cgit v1.2.3