From 8551a96ce0a7753cac5c7f3037f16b40ebad9b34 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Wed, 8 Jul 2026 17:02:58 +0200 Subject: PSM wrong force. oops :) --- examples/stokes_drag_particle_2d_psm_gpu/sim.cpp | 2 ++ modules/core/c++/psm.hpp | 16 ---------------- 2 files changed, 2 insertions(+), 16 deletions(-) diff --git a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp index 572fab0..2c75ef1 100644 --- a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp +++ b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp @@ -178,6 +178,7 @@ saw::error_or step( component> collision{0.8}; component> bb; component,encode::Sycl> abb; + component> one_part_moment; saw::data> eps; eps.at({}).set(1.5f); @@ -222,6 +223,7 @@ saw::error_or step( break; case 2u: opa.apply(macros,index,t_i); + one_part_moment.apply(fields,macros,index,t_i); collision.apply(fields,macros,index,t_i); fg.apply(fields,macros,index,t_i); break; diff --git a/modules/core/c++/psm.hpp b/modules/core/c++/psm.hpp index 6dc146f..10e8f6e 100644 --- a/modules/core/c++/psm.hpp +++ b/modules/core/c++/psm.hpp @@ -60,22 +60,6 @@ public: uint64_t i_opp = dfi::opposite_index[i]; dfs.at({i}) = dfs_cpy.at({i}) + frequency_ * (eq.at(i) - dfs_cpy.at({i})) * porous.at({}) + (dfs_cpy.at({i_opp}) - dfs_cpy.at({i}) ) * flip_porous; } - - auto& force_f = macros.template get<"force">(); - auto& force = force_f.at(index); - for(uint64_t k{0u}; k < Descriptor::D; ++k){ - force.at({{k}}).set(0); - } - - for(uint64_t i{0u}; i < Descriptor::Q; ++i){ - uint64_t i_opp = dfi::opposite_index[i]; - auto dfs_diff = dfs.at({i}) - dfs.at({i_opp}); - for(uint64_t k{0u}; k < Descriptor::D; ++k){ - force.at({{k}}) = force.at({{k}}) + dfs_diff; - } - } - - force = force * porous; } }; -- cgit v1.2.3 From 51783ae88f88b2d70cd4e21961b1b3e9bda0d3e6 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Wed, 8 Jul 2026 17:03:41 +0200 Subject: time step increase --- examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp | 2 +- examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp | 2 +- examples/stokes_drag_particle_2d_psm_gpu/sim.cpp | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp b/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp index c24f808..9e96219 100644 --- a/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp +++ b/examples/stokes_drag_particle_2d_fplbm_gpu/sim.cpp @@ -366,7 +366,7 @@ saw::error_or 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 time_steps{8u*4096ul}; + saw::data time_steps{16u*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 4d6029e..9c28d88 100644 --- a/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp +++ b/examples/stokes_drag_particle_2d_hlbm_gpu/sim.cpp @@ -368,7 +368,7 @@ saw::error_or 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 time_steps{8u*4096ul}; + saw::data time_steps{16u*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 2c75ef1..135bde4 100644 --- a/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp +++ b/examples/stokes_drag_particle_2d_psm_gpu/sim.cpp @@ -351,7 +351,7 @@ saw::error_or 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 time_steps{8u*4096ul}; + saw::data time_steps{16u*4096ul}; auto& info_f = lsd_view.template get<"info">(); -- cgit v1.2.3