Increasing stabilization time and hoping for the best :>

3 files changed, 7 insertions, 40 deletions

diff --git a/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp
index 98c72a2..b860fb3 100644
--- a/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp
+++ b/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp
@@ -224,18 +224,7 @@ saw::error_or<void> step(
 				return vel;
 			}()
 		};
-		component<T,Desc,cmpt::ZouHeVelocityX<false>,encode::Sycl<saw::encode::Native>> flow_out{
-			[&]() -> saw::data<sch::Vector<T,Desc::D>>{
-				saw::data<sch::Vector<T,Desc::D>> vel;
-				uint64_t target_t_i = 8u;
-				if(t_i.get() < target_t_i){
-					vel.at({{0u}}) = (0.0002 / target_t_i) * t_i.get();
-				}else{
-					vel.at({{0u}}) = 0.0002;
-				}
-				return vel;
-			}()
-		};
+		component<T,Desc,cmpt::ZouHeHorizontal<false>,encode::Sycl<saw::encode::Native>> flow_out{1.0};
 
 		h.parallel_for(acpp::sycl::range<Desc::D>{dim_x,dim_y}, [=](acpp::sycl::id<Desc::D> idx){
 			saw::data<sch::FixedArray<sch::UInt64,Desc::D>> index;
@@ -377,7 +366,7 @@ saw::error_or<void> lbm_main(int argc, char** argv){
 	auto lsd_view = make_view(lbm_sycl_data);
 	auto lsdm_view = make_view(lbm_sycl_macro_data);
 	
-	saw::data<sch::UInt64> time_steps{32u*4096ul};
+	saw::data<sch::UInt64> time_steps{64u*4096ul};
 	
 	auto& info_f = lsd_view.template get<"info">();
 
diff --git a/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp
index e0cbda2..f8a9ded 100644
--- a/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp
+++ b/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp
@@ -224,18 +224,7 @@ saw::error_or<void> step(
 				return vel;
 			}()
 		};
-		component<T,Desc,cmpt::ZouHeVelocityX<false>,encode::Sycl<saw::encode::Native>> flow_out{
-			[&]() -> saw::data<sch::Vector<T,Desc::D>>{
-				saw::data<sch::Vector<T,Desc::D>> vel;
-				uint64_t target_t_i = 8u;
-				if(t_i.get() < target_t_i){
-					vel.at({{0u}}) = (0.0002 / target_t_i) * t_i.get();
-				}else{
-					vel.at({{0u}}) = 0.0002;
-				}
-				return vel;
-			}()
-		};
+		component<T,Desc,cmpt::ZouHeHorizontal<false>,encode::Sycl<saw::encode::Native>> flow_out{1.0};
 
 		h.parallel_for(acpp::sycl::range<Desc::D>{dim_x,dim_y}, [=](acpp::sycl::id<Desc::D> idx){
 			saw::data<sch::FixedArray<sch::UInt64,Desc::D>> index;
@@ -379,7 +368,7 @@ saw::error_or<void> lbm_main(int argc, char** argv){
 	auto lsd_view = make_view(lbm_sycl_data);
 	auto lsdm_view = make_view(lbm_sycl_macro_data);
 	
-	saw::data<sch::UInt64> time_steps{32u*4096ul};
+	saw::data<sch::UInt64> time_steps{64u*4096ul};
 	
 	auto& info_f = lsd_view.template get<"info">();
 
diff --git a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp
index 9d678a8..e4f468c 100644
--- a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp
+++ b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp
@@ -202,19 +202,8 @@ saw::error_or<void> step(
 				return vel;
 			}()
 		};
-		component<T,Desc,cmpt::ZouHeVelocityX<false>,encode::Sycl<saw::encode::Native>> flow_out{
-			[&]() -> saw::data<sch::Vector<T,Desc::D>>{
-				saw::data<sch::Vector<T,Desc::D>> vel;
-				uint64_t target_t_i = 8u;
-				if(t_i.get() < target_t_i){
-					vel.at({{0u}}) = (0.0002 / target_t_i) * t_i.get();
-				}else{
-					vel.at({{0u}}) = 0.0002;
-				}
-				return vel;
-			}()
-		};
-
+		component<T,Desc,cmpt::ZouHeHorizontal<false>,encode::Sycl<saw::encode::Native>> flow_out{1.0};
+	
 		component<T,Desc,cmpt::ForceGather, encode::Sycl<saw::encode::Native>> fg;
 
 		h.parallel_for(acpp::sycl::range<Desc::D>{dim_x,dim_y}, [=](acpp::sycl::id<Desc::D> idx){
@@ -360,7 +349,7 @@ saw::error_or<void> lbm_main(int argc, char** argv){
 	auto lsd_view = make_view(lbm_sycl_data);
 	auto lsdm_view = make_view(lbm_sycl_macro_data);
 
-	saw::data<sch::UInt64> time_steps{32u*4096ul};
+	saw::data<sch::UInt64> time_steps{64u*4096ul};
 
 	auto& info_f = lsd_view.template get<"info">();