From 16700f013886b8d8cb02e322c2d66c00179d9c89 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Tue, 3 Mar 2026 19:46:31 +0100 Subject: Fixing ZouHe? --- examples/poiseulle_particles_2d_gpu/sim.cpp | 3 ++- lib/core/c++/boundary.hpp | 4 ++-- 2 files changed, 4 insertions(+), 3 deletions(-) diff --git a/examples/poiseulle_particles_2d_gpu/sim.cpp b/examples/poiseulle_particles_2d_gpu/sim.cpp index 6d012e0..bb81383 100644 --- a/examples/poiseulle_particles_2d_gpu/sim.cpp +++ b/examples/poiseulle_particles_2d_gpu/sim.cpp @@ -131,7 +131,7 @@ saw::error_or setup_initial_conditions( rho.at({}) = {1}; auto& vel = vel_f.at(index); if(info_f.at(index).get() == 2u){ - vel.at({{0u}}) = 0.01; + vel.at({{0u}}) = 0.0; } auto eq = equilibrium(rho,vel); @@ -433,6 +433,7 @@ saw::error_or lbm_main(int argc, char** argv){ std::cout<<"Status: "<().get() * 100 / time_steps.get())<<"%"< k{0u}; k < saw::data{Descriptor::Q}; ++k){ auto c_k = dfi::directions[k.get()]; - if(c_k[0u]*known_dir <= 0){ + if(c_k[0u]*known_dir >= 0){ sum_df += dfs_old.at({k}); } } @@ -147,7 +147,7 @@ public: for(saw::data k{0u}; k < saw::data{Descriptor::Q}; ++k){ auto c_k = dfi::directions[k.get()]; - if(c_k[0u]*known_dir < 0){ + if(c_k[0u]*known_dir > 0){ sum_unknown_dfs += dfs_old.at({k}) * c_k[0u]; } } -- cgit v1.2.3