From 2b335de452a4cbb2c8859fecd29e1f4069cf193a Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Wed, 8 Jul 2026 19:46:29 +0200 Subject: Just compare against BGK ._. --- examples/stokes_drag_particle_2d_bgk_gpu/sim.cpp | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) (limited to 'examples/stokes_drag_particle_2d_bgk_gpu') diff --git a/examples/stokes_drag_particle_2d_bgk_gpu/sim.cpp b/examples/stokes_drag_particle_2d_bgk_gpu/sim.cpp index 9426371..3520fc3 100644 --- a/examples/stokes_drag_particle_2d_bgk_gpu/sim.cpp +++ b/examples/stokes_drag_particle_2d_bgk_gpu/sim.cpp @@ -188,15 +188,16 @@ saw::error_or step( p_pos.at({{0u}}) = dim_x * 0.5; p_pos.at({{1u}}) = dim_y * 0.5; } - component> opa{p_pos,rad,eps}; - - component,encode::Sycl> flow_in{ - [&](){ - uint64_t target_t_i = 64u; + component,encode::Sycl> flow_in{ + [&]() -> saw::data>{ + saw::data> vel; + uint64_t target_t_i = 8u; if(t_i.get() < target_t_i){ - return 1.0 + (0.0002 / target_t_i) * t_i.get(); + vel.at({{0u}}) = (0.0002 / target_t_i) * t_i.get(); + }else{ + vel.at({{0u}}) = 0.0002; } - return 1.0002; + return vel; }() }; component,encode::Sycl> flow_out{1.0}; @@ -217,7 +218,6 @@ saw::error_or step( abb.apply(fields,index,t_i); break; case 2u: - opa.apply(macros,index,t_i); collision.apply(fields,macros,index,t_i); break; case 3u: -- cgit v1.2.3