gaussian.cpp

Go to the documentation of this file.

/*
 *  SPDX-FileCopyrightText: Copyright 2021, Siavash Ameli <sameli@berkeley.edu>
 *  SPDX-License-Identifier: BSD-3-Clause
 *  SPDX-FileType: SOURCE
 *
 *  This program is free software: you can redistribute it and/or modify it
 *  under the terms of the license found in the LICENSE.txt file in the root
 *  directory of this source tree.
 */
 
 
// =======
// Headers
// =======
 
// Before including cmath, define _USE_MATH_DEFINES. This is only required to
// define the math constants like M_PI, etc, in win32 operating system.
#if defined(WIN32) || defined(_WIN32) || defined(__WIN32) && \
    !defined(__CYGWIN__)
    #define _USE_MATH_DEFINES
#endif
 
#include <cmath>  // exp, M_SQRT1_2, M_2_SQRTPI
#include "./gaussian.h"
 
 
// ========
// Gaussian
// ========
 
 
Gaussian::Gaussian(double mu_, double sigma_)
{
    this->mu = mu_;
    this->sigma = sigma_;
}
 
// =================
// Gaussian function (float)
// =================
 
 
float Gaussian::function(const float lambda_) const
{
    float mu_ = static_cast<float>(this->mu);
    float sigma_ = static_cast<float>(this->sigma);
    float x = (lambda_ - mu_) / sigma_;
    return (0.5 * M_SQRT1_2 * M_2_SQRTPI / sigma_) * exp(-0.5 * x * x);
}
 
 
// =================
// Gaussian function (double)
// =================
 
 
double Gaussian::function(const double lambda_) const
{
    double x = (lambda_ - this->mu) / this->sigma;
    return (0.5 * M_SQRT1_2 * M_2_SQRTPI / this->sigma) * exp(-0.5 * x * x);
}
 
 
// =================
// Gaussian function (long double)
// =================
 
 
long double Gaussian::function(const long double lambda_) const
{
    long double mu_ = static_cast<long double>(this->mu);
    long double sigma_ = static_cast<long double>(this->sigma);
    long double x = (lambda_ - mu_) / sigma_;
    return (0.5 * M_SQRT1_2 * M_2_SQRTPI / sigma_) * exp(-0.5 * x * x);
}