diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-09-07 15:08:20 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-09-07 15:08:20 +0200 |
| commit | 96619fc23312e91e9fbe453cb115f116784c032a (patch) | |
| tree | 9d931191117f54f2134b3ad3615e2092870a99ca /modules/core/c++/hlbm.hpp | |
| parent | d2a406aceeb935d4eac1b5ae02e268be4463b6b0 (diff) | |
| parent | 54dfad7ac102fe40d71d8818d77801160a7073c1 (diff) | |
| download | libs-lbm-96619fc23312e91e9fbe453cb115f116784c032a.tar.gz | |
Diffstat (limited to 'modules/core/c++/hlbm.hpp')
| -rw-r--r-- | modules/core/c++/hlbm.hpp | 15 |
1 files changed, 4 insertions, 11 deletions
diff --git a/modules/core/c++/hlbm.hpp b/modules/core/c++/hlbm.hpp index e8ea66a..0df5896 100644 --- a/modules/core/c++/hlbm.hpp +++ b/modules/core/c++/hlbm.hpp @@ -102,12 +102,6 @@ public: template<typename T, typename Desc, typename Encode> class component<T, Desc, cmpt::HlbmOneParticle, Encode> final { private: -/* - template<typename CellFieldSchema, typename MacroFieldSchema, typename ParticleSchema, uint64_t i> - void apply_i(const saw::data<CellFieldSchema, Encode>& field, const saw::data<MacroFieldSchema,Encode>& macros, const saw::data<ParticleSchema,Encode>& part_groups, saw::data<sch::FixedArray<sch::UInt64,1u>> index, saw::data<sch::UInt64> time_step) const { - // if constexpr ( i < ) - } -*/ public: template<typename CellFieldSchema, typename MacroFieldSchema, typename ParticleSchema> void apply(const saw::data<CellFieldSchema, Encode>& field, const saw::data<MacroFieldSchema,Encode>& macros, const saw::data<ParticleSchema,Encode>& part_group, saw::data<sch::FixedArray<sch::UInt64,1u>> index, saw::data<sch::UInt64> time_step, saw::data<sch::UInt64> sub_steps) const { @@ -167,10 +161,6 @@ public: auto vel_p_old = (pirb_pos-pirb_pos_old) / ts; iterator<Desc::D>::apply([&](const auto& index_f) -> void{ - // ask for the d_k value here. - // For every value im iterating over I need sth - // std::cout<<"Pos: "<<index.at({0u}).get()<<" "<<index.at({1u}).get()<<std::endl; - auto& dfs = dfs_old_f.at(index_f); saw::data<sch::Vector<T,Desc::D>> rel_dist = saw::math::vectorize_data(index_f).template cast_to<T>() - pirb_pos; saw::data<sch::Scalar<T>> eps; @@ -207,9 +197,12 @@ public: momentum = momentum + ei_dfs; } - // technically needs to adjust for rotation as well + + // TODO technically needs to adjust for rotation as well auto& rho = mrho_f.at(index_f); + + macros.template get<"force">().at(index_f) = momentum; force_p = force_p + momentum; },start,stop); |