diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-09-22 09:27:15 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-09-22 09:27:15 +0200 |
| commit | 79196db2ac3d31543056b3b7abe90ad609c51b4f (patch) | |
| tree | fb8a8bc480d8af824e00ae36cf20dd1ea882f499 /examples/one_particle_sedimentation_2d/step.hpp | |
| parent | 54dfad7ac102fe40d71d8818d77801160a7073c1 (diff) | |
| download | libs-lbm-79196db2ac3d31543056b3b7abe90ad609c51b4f.tar.gz | |
Fixing two way coupling
Diffstat (limited to 'examples/one_particle_sedimentation_2d/step.hpp')
| -rw-r--r-- | examples/one_particle_sedimentation_2d/step.hpp | 11 |
1 files changed, 6 insertions, 5 deletions
diff --git a/examples/one_particle_sedimentation_2d/step.hpp b/examples/one_particle_sedimentation_2d/step.hpp index c3a3b8a..4bc9dde 100644 --- a/examples/one_particle_sedimentation_2d/step.hpp +++ b/examples/one_particle_sedimentation_2d/step.hpp @@ -8,6 +8,7 @@ namespace lbm { template<typename T, typename Desc, typename Coll> saw::error_or<void> step( const converter<T>& conv, + const saw::data<sch::Pure<T>>& tau, saw::data<sch::Ptr<sch::ChunkStruct<T,Desc>>,encode::Sycl<saw::encode::Native>>& fields, saw::data<sch::Ptr<sch::MacroStruct<T,Desc>>,encode::Sycl<saw::encode::Native>>& macros, saw::data<sch::Ptr<sch::ParticleSpheroidGroup<T,Desc>>,encode::Sycl<saw::encode::Native>>& particles, @@ -47,7 +48,7 @@ saw::error_or<void> step( }).wait(); q.submit([&](acpp::sycl::handler& h){ - component<T,Desc,cmpt::Hlbm,encode::Sycl<saw::encode::Native>> collision{1.0}; + component<T,Desc,cmpt::Hlbm,encode::Sycl<saw::encode::Native>> collision{tau}; component<T,Desc,cmpt::BounceBack,encode::Sycl<saw::encode::Native>> bb; h.parallel_for(acpp::sycl::range<Desc::D>{dim_x,dim_y}, [=](acpp::sycl::id<Desc::D> idx){ @@ -99,14 +100,14 @@ saw::error_or<void> step( index.at({{i}}).set(idx[i]); } - fplbm_one_part.apply(fields,macros,particles,index,t_i,{16u}); + fplbm_one_part.apply(fields,macros,particles,index,t_i,{32u}); }); }).wait(); // auto coll_ev = q.submit([&](acpp::sycl::handler& h){ - component<T,Desc,cmpt::BGK,encode::Sycl<saw::encode::Native>> bgk{1.0}; - component<T,Desc,cmpt::FpLbm,encode::Sycl<saw::encode::Native>> collision{1.0}; + component<T,Desc,cmpt::BGK,encode::Sycl<saw::encode::Native>> bgk{tau}; + component<T,Desc,cmpt::FpLbm,encode::Sycl<saw::encode::Native>> collision{tau}; component<T,Desc,cmpt::BounceBack,encode::Sycl<saw::encode::Native>> bb; component<T,Desc,cmpt::ZouHeHorizontal<false>,encode::Sycl<saw::encode::Native>> flow_out{1.0}; @@ -182,7 +183,7 @@ saw::error_or<void> step( // auto coll_ev = q.submit([&](acpp::sycl::handler& h){ component<T,Desc,cmpt::BGK,encode::Sycl<saw::encode::Native>> bgk{1.0}; - component<T,Desc,cmpt::Psm,encode::Sycl<saw::encode::Native>> collision{1.0}; + component<T,Desc,cmpt::Psm,encode::Sycl<saw::encode::Native>> collision{tau}; component<T,Desc,cmpt::BounceBack,encode::Sycl<saw::encode::Native>> bb; component<T,Desc,cmpt::ZouHeHorizontal<false>,encode::Sycl<saw::encode::Native>> flow_out{1.0}; |