From c0549d71b2109f10c1238db8b22362e7826ba61b Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Sun, 5 Jul 2026 15:59:23 +0200 Subject: Just rename from lib to modules --- lib/core/c++/momentum_gather.hpp | 54 ---------------------------------------- 1 file changed, 54 deletions(-) delete mode 100644 lib/core/c++/momentum_gather.hpp (limited to 'lib/core/c++/momentum_gather.hpp') diff --git a/lib/core/c++/momentum_gather.hpp b/lib/core/c++/momentum_gather.hpp deleted file mode 100644 index 7b1fc96..0000000 --- a/lib/core/c++/momentum_gather.hpp +++ /dev/null @@ -1,54 +0,0 @@ -#pragma once - -#include "common.hpp" -#include "macroscopic.hpp" -#include "component.hpp" -#include "equilibrium.hpp" - -namespace kel { -namespace lbm { -namespace cmpt { -// Gather the force, Luke! -struct ForceGather {}; -} - -template -class component { -private: -public: - template - void apply(const saw::data& field, const saw::data& macros, - saw::data> index, saw::data time_step) const { - - using dfi = df_info; - bool is_even = ((time_step.get() % 2) == 0); - - auto& dfs_old_f = (is_even) ? field.template get<"dfs_old">() : field.template get<"dfs">(); - - auto& porous_f = macros.template get<"porosity">(); - auto& porous = porous_f.at(index); - - auto& force_f = macros.template get<"force">(); - auto& force = force_f.at(index); - - for(uint64_t k{0u}; k < Descriptor::D; ++k){ - force.at({{k}}).set(0); - } - - auto& dfs = dfs_old_f.at(index); - for(uint64_t i{0u}; i < Descriptor::Q; ++i){ - uint64_t i_opp = dfi::opposite_index[i]; - auto dfs_diff = dfs.at({i}) - dfs.at({i_opp}); - for(uint64_t k{0u}; k < Descriptor::D; ++k){ - force.at({{k}}) = force.at({{k}}) + dfs_diff * saw::data{dfi::directions[i][k]}; - } - } - - saw::data> one; - one.at({}) = 1.0; - auto flip_porous = one - porous; - force = force * flip_porous; - } -}; -} -} -- cgit v1.2.3