From d8321ed7a23f7554cd0144c4f5e8dc625bbee4fe Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Thu, 26 Feb 2026 14:35:57 +0100 Subject: Clean repo for run and make benchmarkable hlbm case --- examples/meta/meta.cpp | 20 ++++++++++++++++++-- examples/poiseulle_particles_2d_gpu/sim.cpp | 19 ++++++++++--------- 2 files changed, 28 insertions(+), 11 deletions(-) diff --git a/examples/meta/meta.cpp b/examples/meta/meta.cpp index c638021..a3ed7e7 100644 --- a/examples/meta/meta.cpp +++ b/examples/meta/meta.cpp @@ -50,6 +50,9 @@ saw::error_or meta_main(int argc, char** argv){ saw::data args_dat; { + { + + } saw::data args_enc{argc,argv}; saw::codec args_codec; @@ -66,8 +69,13 @@ saw::error_or meta_main(int argc, char** argv){ switch(args_pos.get<0u>().get()){ case 1u: { - assert(false); - return saw::make_error(); + switch(args_pos.get<1u>().get()){ + case 3u: + meta_lbm_info>(args_str); + break; + default: + return saw::make_error("Second positional needs to be a valid velocity set"); + } } break; case 2u: @@ -87,6 +95,14 @@ saw::error_or meta_main(int argc, char** argv){ case 3u: { switch(args_pos.get<1u>().get()){ + case 7u: + return saw::make_error("D3Q7 Velocity Set is not implemented"); + // meta_lbm_info>(args_str); + break; + case 19u: + return saw::make_error("D3Q17 Velocity Set is not implemented"); + // meta_lbm_info>(args_str); + break; case 27u: meta_lbm_info>(args_str); break; diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index 01fec89..fffcb31 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -10,9 +10,10 @@ namespace kel { namespace lbm { -constexpr uint64_t dim_x = 1024ul; constexpr uint64_t dim_y = 512ul; -constexpr uint64_t particle_size = 1ul; +constexpr uint64_t dim_x = dim_y * 20ul; + +constexpr uint64_t particle_amount = 1ul; namespace sch { using namespace saw::schema; @@ -60,7 +61,7 @@ template saw::error_or setup_initial_conditions( saw::data>& fields, saw::data>& macros, - saw::data, particle_size>>& particles + saw::data, particle_amount>>& particles ){ auto& info_f = fields.template get<"info">(); auto& porous_f = macros.template get<"porosity">(); @@ -158,14 +159,14 @@ saw::error_or setup_initial_conditions( df_f.get_dims() ); - for(saw::data i{0u}; i < saw::data{particle_size}; ++i){ + for(saw::data i{0u}; i < saw::data{particle_amount}; ++i){ auto& part = particles.at(i); saw::data> pos; - pos.at({{0u}}) = dim_x * 0.5; + pos.at({{0u}}) = dim_x * 0.25; pos.at({{1u}}) = dim_y * 0.5; saw::data> rad, dense, dt; - rad.at({}) = dim_y * 0.2; + rad.at({}) = dim_y * 0.1; dense.at({}) = 1.0; dt.at({}) = 1.0; part = create_spheroid_particle( @@ -182,7 +183,7 @@ template saw::error_or step( saw::data>,encode::Sycl>& fields, saw::data>,encode::Sycl>& macros, - saw::data, particle_size>>& particles, + saw::data, particle_amount>>& particles, saw::data t_i, device& dev ){ @@ -321,7 +322,7 @@ saw::error_or lbm_main(int argc, char** argv){ // saw::data> meta{{dim_x,dim_y}}; auto lbm_data_ptr = saw::heap>>(); auto lbm_macro_data_ptr = saw::heap>>(); - auto lbm_particle_data_ptr = saw::heap, particle_size>>>(); + auto lbm_particle_data_ptr = saw::heap, particle_amount>>>(); std::cout<<"Estimated Bytes: "<,sch::MacroStruct>().get()< lbm_main(int argc, char** argv){ saw::data, encode::Sycl> lbm_sycl_data{sycl_q}; saw::data, encode::Sycl> lbm_sycl_macro_data{sycl_q}; - saw::data, particle_size>, encode::Sycl> lbm_sycl_particle_data{sycl_q}; + saw::data, particle_amount>, encode::Sycl> lbm_sycl_particle_data{sycl_q}; sycl_q.wait(); auto lsd_view = make_chunk_struct_view(lbm_sycl_data); -- cgit v1.2.3