diff options
| author | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-06-30 23:52:02 +0200 |
|---|---|---|
| committer | Claudius "keldu" Holeksa <mail@keldu.de> | 2026-06-30 23:52:02 +0200 |
| commit | 2bae8b93faf614b1acb5df7269087b4f3786da5a (patch) | |
| tree | a6db1328ed035c6130b6b29110e0891f6c78c7e4 /lib/core/c++/particle | |
| parent | b88d59477a7973bdee102aaf0e26c13c9059048b (diff) | |
| download | libs-lbm-2bae8b93faf614b1acb5df7269087b4f3786da5a.tar.gz | |
Final run
Diffstat (limited to 'lib/core/c++/particle')
| -rw-r--r-- | lib/core/c++/particle/porosity.hpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/lib/core/c++/particle/porosity.hpp b/lib/core/c++/particle/porosity.hpp index 11185c5..2f39d2e 100644 --- a/lib/core/c++/particle/porosity.hpp +++ b/lib/core/c++/particle/porosity.hpp @@ -48,17 +48,17 @@ public: auto rad_high = rad + eps_h; if(s_dist.at({}) <= rad_low.at({})){ - por.at({}).set(1); + por.at({}).set(0); return por; } if(s_dist.at({}) >= rad_high.at({})){ - por.at({}).set(0); + por.at({}).set(1); return por; } { - typename saw::native_data_type<T>::type inner = (std::numbers::pi / ( eps.at({}).get() )) * (s_dist - rad_low).at({}).get(); - auto cos_inner = std::cos(inner); + typename saw::native_data_type<T>::type inner = (std::numbers::pi / ( 2 * eps.at({}).get() )) * (s_dist - rad_low).at({}).get(); + auto cos_inner = std::sin(inner); auto cos_i_2 = cos_inner * cos_inner; por.at({}).set(cos_i_2); } |
