summaryrefslogtreecommitdiff
path: root/modules/core/c++/particle
diff options
context:
space:
mode:
authorClaudius "keldu" Holeksa <mail@keldu.de>2026-07-23 17:28:51 +0200
committerClaudius "keldu" Holeksa <mail@keldu.de>2026-07-23 17:28:51 +0200
commitc8bea2d6a8c6cd6b0950d75646eef0becf1fce1e (patch)
tree81e00012524b5b303770e0977769c24716d1f181 /modules/core/c++/particle
parentda8cb1c7dd40ef18b85685f99369a31ee36370a1 (diff)
parent3489ac1d725c0eadb7d75329f9e6074ed177d125 (diff)
downloadlibs-lbm-c8bea2d6a8c6cd6b0950d75646eef0becf1fce1e.tar.gz
Merge branch 'dev'
Diffstat (limited to 'modules/core/c++/particle')
-rw-r--r--modules/core/c++/particle/particle.hpp12
1 files changed, 11 insertions, 1 deletions
diff --git a/modules/core/c++/particle/particle.hpp b/modules/core/c++/particle/particle.hpp
index f759c8d..a3669b4 100644
--- a/modules/core/c++/particle/particle.hpp
+++ b/modules/core/c++/particle/particle.hpp
@@ -74,9 +74,19 @@ saw::data<sch::ParticleGroup<T, D, PartAmount, coll::Spheroid<T>>> create_sphero
}
total_mass.at({}) = total_mass.at({}) + dpi;
-
},{},mask_dims);
+
+ // TODO ? Do mass properly
+ if constexpr ( D == 2u ) {
+ total_mass.at({}) = rad_d * rad_d * density.at({}) * 3.141592;
+ }else if constexpr ( D == 3u ){
+
+ }else if constexpr ( D== 1u ){
+
+ }
+ std::cout<<"Total Mass: "<<total_mass.at({}).get()<<std::endl;
+
for(uint64_t i = 0u; i < D; ++i){
com.at({{i}}) = com.at({{i}}) / total_mass.at({});
}