diff options
Diffstat (limited to 'lib/core/c++/particle/particle.hpp')
| -rw-r--r-- | lib/core/c++/particle/particle.hpp | 13 |
1 files changed, 12 insertions, 1 deletions
diff --git a/lib/core/c++/particle/particle.hpp b/lib/core/c++/particle/particle.hpp index 5110893..938131b 100644 --- a/lib/core/c++/particle/particle.hpp +++ b/lib/core/c++/particle/particle.hpp @@ -79,7 +79,8 @@ saw::data<sch::ParticleGroup<T,D, sch::ParticleCollisionSpheroid<T,radius>>> cre mask_dims.at({i}) = mask_resolution; } saw::data<sch::Scalar<T>> mask_step; - saw::data<T> dia_d{radius*2}; + saw::data<T> rad_d{radius}; + saw::data<T> dia_d = rad_d * 2; mask_step.at({}) = dia_d / mask_resolution.template cast_to<T>(); mask = {mask_dims}; @@ -89,6 +90,16 @@ saw::data<sch::ParticleGroup<T,D, sch::ParticleCollisionSpheroid<T,radius>>> cre com.at({{i}}) = {}; } saw::data<sch::UInt64> ele_ctr{0u}; + + // Radius ^ 2 + saw::data<sch::Scalar<T>> rad_2_d; + rad_2_d.at({}) = rad_d * rad_d; + + saw::data<sch::Vector<T,D>> center; + for(uint64_t i = 0u; i < D; ++i){ + com.at({{i}}) = ; + } + iterator<D>::apply([&](const auto& index){ ++ele_ctr; |