diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2025-08-04 16:16:17 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2025-08-04 16:16:17 +0200 |
| commit | 51a044e3160df05ad56102d3b8b1e0087c60d111 (patch) | |
| tree | 334f0ba1f7bb02872661207dfcc9f9d0c7cfa587 /c++/particle/geometry | |
| parent | 9ff280692f5cc78c0b006f8b6f1b66f7f36ba2a6 (diff) | |
emits weird behaviour on particles
Diffstat (limited to 'c++/particle/geometry')
| -rw-r--r-- | c++/particle/geometry/circle.hpp | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/c++/particle/geometry/circle.hpp b/c++/particle/geometry/circle.hpp index 77fa9d8..069a1a5 100644 --- a/c++/particle/geometry/circle.hpp +++ b/c++/particle/geometry/circle.hpp @@ -23,15 +23,15 @@ public: uint64_t size = rad_i + 2*boundary_nodes; grid = {{{size,size}}}; - saw::data<T> delta_x{2.0 / rad_i}; + saw::data<T> delta_x{static_cast<typename saw::native_data_type<T>::type>(2.0 / static_cast<double>(rad_i))}; saw::data<T> center = (saw::data<T>{1.0} + saw::data<T>{2.0} * saw::data<T>{static_cast<saw::native_data_type<T>::type>(boundary_nodes)/rad_i}); for(uint64_t i = 0; i < size; ++i){ for(uint64_t j = 0; j < size; ++j){ grid.at({{i,j}}).set(0); if(i >= boundary_nodes and j >= boundary_nodes and i < (rad_i + boundary_nodes) and j < (rad_i + boundary_nodes) ){ - saw::data<T> fi = saw::data<T>{0.5+static_cast<saw::native_data_type<T>::type>(i)} * delta_x - center; - saw::data<T> fj = saw::data<T>{0.5+static_cast<saw::native_data_type<T>::type>(j)} * delta_x - center; + saw::data<T> fi = saw::data<T>{static_cast<saw::native_data_type<T>::type>(0.5+i)} * delta_x - center; + saw::data<T> fj = saw::data<T>{static_cast<saw::native_data_type<T>::type>(0.5+j)} * delta_x - center; auto norm_f_ij = fi*fi + fj*fj; if(norm_f_ij.get() <= 1){ |