Changeset 4337

Timestamp:

Apr 14, 2023, 3:22:32 AM (5 months ago)

Author:

Peter

Message:

avoid a couple of multiplications

File:

• trunk/yat/statistics/ln_pdf.cc (modified) (2 diffs)

trunk/yat/statistics/ln_pdf.cc

@@ -69 +69 @@
     // pdf = 1/(s sqrt(2pi)) * exp(-1/2 * (x/s)^2)
 
-    return -std::log(s) - 0.5*M_LN2 - 0.5*M_LNPI - 0.5 * std::pow(x/s, 2);
+    return -std::log(s) - 0.5 * (M_LN2 + M_LNPI + std::pow(x/s, 2));
   }
 
@@ -126 +126 @@
     // pdf = 1/(sqrt(2pi)) * exp(-1/2 * x^2)
 
-    return - 0.5*M_LN2 - 0.5*M_LNPI - 0.5 * x*x;
+    return - 0.5 * (M_LN2 + M_LNPI + x*x);
   }