vifa_gaussian.cxx

Go to the documentation of this file.

#include <vcl_cmath.h>
#include <vifa/vifa_gaussian.h>

float vifa_gaussian::pdf(float  x)
{
  if (sigma_ < StatEPSILON)
  {
    // Degenerate distribution:
    // The variance is zero, so the pdf is a unit impulse at the mean.
    // Return 1.0 if x == mean, and 0.0 otherwise.

    return vcl_fabs(x - mu_) < StatEPSILON ? 1.0f : 0.0f;
  }
  else
  {
    float  xp = (x - mu_) / sigma_;
    return norm_dens(xp) / sigma_;
  }
}


float vifa_gaussian::norm_dens(float  x)
{
  // Check to see if the magnitude of x is large enough to
  // warrant clipping the pdf to 0.0
  return x < -EXPLIMIT || x > EXPLIMIT ? 0.0f :
                                         float(I_SQRT_2PI * vcl_exp(-0.5 * x * x));
}

---