From 9a1ee9e0b348f75f5f4f8fb41ea901a0dea02f13 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Tue, 25 Aug 2026 16:49:40 +0200 Subject: Weird things are happening? --- .../schaefer_turek_durst_krause_rannbacher/sim.cpp | 4 ++-- .../step.hpp | 22 ++++++---------------- 2 files changed, 8 insertions(+), 18 deletions(-) (limited to 'examples/schaefer_turek_durst_krause_rannbacher') diff --git a/examples/schaefer_turek_durst_krause_rannbacher/sim.cpp b/examples/schaefer_turek_durst_krause_rannbacher/sim.cpp index 59f61d9..5a92681 100644 --- a/examples/schaefer_turek_durst_krause_rannbacher/sim.cpp +++ b/examples/schaefer_turek_durst_krause_rannbacher/sim.cpp @@ -48,7 +48,7 @@ saw::error_or lbm_main(const saw::data& args){ // delta_x {{0.41 / dim_y}}, // delta_t - {{1.067e-4 * 0.5}} + {{1.067e-4}} }; print_lbm_meta(conv,{1e-3},{1e-4},{0.5 * dim_y}); @@ -122,7 +122,7 @@ saw::error_or lbm_main(const saw::data& args){ auto lsdm_view = make_view(lbm_sycl_macro_data); auto lsdp_view = make_view(lbm_sycl_particle_data); - saw::data time_steps{128u*1024ul}; + saw::data time_steps{64u*1024ul}; auto& info_f = lsd_view.template get<"info">(); diff --git a/examples/schaefer_turek_durst_krause_rannbacher/step.hpp b/examples/schaefer_turek_durst_krause_rannbacher/step.hpp index beb123b..9a75e6f 100644 --- a/examples/schaefer_turek_durst_krause_rannbacher/step.hpp +++ b/examples/schaefer_turek_durst_krause_rannbacher/step.hpp @@ -46,7 +46,7 @@ saw::error_or step( }).wait(); q.submit([&](acpp::sycl::handler& h){ - component> collision{0.75}; + component> collision{1.0}; component> bb; component,encode::Sycl> flow_out{1.0}; @@ -72,19 +72,9 @@ saw::error_or step( component,encode::Sycl> flow_in{ [&]() -> saw::data> { saw::data> vel; - saw::data scale; - { - scale.set([&](){ - if(t_i.get() > 1024u){ - return 1.0f; - } - - return t_i.get() / 1024.0f; - }()); - } { auto y_si = conv.meter_lbm_to_si(index.at({{1u}}).template cast_to()).handle(); - vel.at({{0u}}) = conv.velocity_si_to_lbm({{static_cast::type>((1.2 / 0.41) * y_si.get() - (1.2 / (0.41*0.41)) * y_si.get() * y_si.get())}}).handle() * scale; + vel.at({{0u}}) = conv.velocity_si_to_lbm({{static_cast::type>((1.2 / 0.41) * y_si.get() - (1.2 / (0.41*0.41)) * y_si.get() * y_si.get())}}).handle(); vel.at({{1u}}) = 0.0f; } return vel; @@ -131,8 +121,8 @@ saw::error_or step( // auto coll_ev = q.submit([&](acpp::sycl::handler& h){ - component> bgk{0.75}; - component> collision{0.75}; + component> bgk{1.0}; + component> collision{1.0}; component> bb; component,encode::Sycl> flow_out{1.0}; @@ -207,8 +197,8 @@ saw::error_or step( // auto coll_ev = q.submit([&](acpp::sycl::handler& h){ - component> bgk{0.75}; - component> collision{0.75}; + component> bgk{1.0}; + component> collision{1.0}; component> bb; component,encode::Sycl> flow_out{1.0}; -- cgit v1.2.3