feat: use lognormal distribution

yoctoyotta1024 · Apr 29, 2024 · 5484126 · 5484126
1 parent a66607b
commit 5484126
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/libs/cartesiandomain/add_supers_at_domain_top.cpp b/libs/cartesiandomain/add_supers_at_domain_top.cpp
@@ -165,7 +165,8 @@ double CreateSuperdrop::droplet_numconc_distribution(const double radius,
   return dnumconc_dlnr * dlnr;  // number of droplets per unit volume
 }
 
-double lognormal_pdf(const double radius) {
+/* normalised lognormal distribution returns the probability density of a given radius */
+double CreateSuperdrop::lognormal_pdf(const double radius) const {
   const auto inverse_norm = radius * sigtilda * std::sqrt(2.0 * std::numbers::pi);
   const auto expo = (std::log(radius) - mutilda) / sigtilda;
   return std::exp(-0.5 * expo * expo) / inverse_norm;