gaussian_int.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>  // erf, M_SQRT1_2
#include "./gaussian_int.h"
 
 
// ============
// Gaussian Int
// ============
 
 
GaussianInt::GaussianInt(double mu_, double sigma_)
{
    this->mu = mu_;
    this->sigma = sigma_;
}
 
// =====================
// Gaussian Int function (float)
// =====================
 
 
float GaussianInt::function(const float lambda_) const
{
    float mu_ = static_cast<float>(this->mu);
    float sigma_ = static_cast<float>(this->sigma);
    float x = (lambda_ - mu_) / sigma_;
 
    // C++11 and later
    #if __cplusplus >= 201103L
        return 0.5f * (1.0f + std::erf(x * M_SQRT1_2));
    #else
        return 0.5f * (1.0f + static_cast<float>(erf(x * M_SQRT1_2)));
    #endif
}
 
 
// =====================
// Gaussian Int function (double)
// =====================
 
 
double GaussianInt::function(const double lambda_) const
{
    double x = (lambda_ - this->mu) / this->sigma;
 
    // C++11 and later
    #if __cplusplus >= 201103L
        return 0.5 * (1.0 + std::erf(x * M_SQRT1_2));
    #else
        return 0.5 * (1.0 + erf(x * M_SQRT1_2));
    #endif
}
 
 
// =====================
// Gaussian Int function (long double)
// =====================
 
 
long double GaussianInt::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_;
 
    // C++11 and later
    #if __cplusplus >= 201103L
        return 0.5l * (1.0l + std::erf(x * M_SQRT1_2));
    #else
        return 0.5l * (1.0l + static_cast<long double>(erf(x * M_SQRT1_2)));
    #endif
}