From 573bfdc0bbffdb1d68259bc348e36d57471d56be Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Tue, 24 Jun 2025 16:49:31 +0200 Subject: Broken cavity2d. why? --- c++/collision.hpp | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'c++/collision.hpp') diff --git a/c++/collision.hpp b/c++/collision.hpp index 7ca7c9f..128a499 100644 --- a/c++/collision.hpp +++ b/c++/collision.hpp @@ -40,18 +40,18 @@ public: template void apply(saw::data& field, saw::data> index, uint64_t time_step){ bool is_even = ((time_step % 2) == 0); - auto& cell = field.at(index); + auto& cell = field(index); auto& dfs_old = is_even ? cell.template get<"dfs_old">() : cell.template get<"dfs">(); auto& dfs = (!is_even) ? cell.template get<"dfs_old">() : cell.template get<"dfs">(); - typename saw::native_data_type::type rho; - std::array::type, Descriptor::D> vel; + saw::data rho; + saw::data> vel; compute_rho_u(dfs,rho,vel); - auto eq = equilibrium(rho,vel); + auto eq = equilibrium(rho,vel); for(uint64_t i = 0u; i < Descriptor::Q; ++i){ - dfs({i}).set(dfs({i}).get() + (1.0 / relaxation_) * (eq[i] - dfs({i}).get())); + dfs({i}).set(dfs({i}).get() + (1.0 / relaxation_) * (eq.at(i).get() - dfs({i}).get())); } } -- cgit v1.2.3