summaryrefslogtreecommitdiff
path: root/lib/core/c++/particle
diff options
context:
space:
mode:
authorClaudius "keldu" Holeksa <mail@keldu.de>2026-06-30 23:52:02 +0200
committerClaudius "keldu" Holeksa <mail@keldu.de>2026-06-30 23:52:02 +0200
commit2bae8b93faf614b1acb5df7269087b4f3786da5a (patch)
treea6db1328ed035c6130b6b29110e0891f6c78c7e4 /lib/core/c++/particle
parentb88d59477a7973bdee102aaf0e26c13c9059048b (diff)
downloadlibs-lbm-2bae8b93faf614b1acb5df7269087b4f3786da5a.tar.gz
Final run
Diffstat (limited to 'lib/core/c++/particle')
-rw-r--r--lib/core/c++/particle/porosity.hpp8
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);
}