animaGaussLaguerreQuadrature.h

Go to the documentation of this file.

#pragma once
#include "AnimaSpecialFunctionsExport.h"
#include <vector>

#include <boost/math/quadrature/gauss.hpp>

namespace anima
{

class ANIMASPECIALFUNCTIONS_EXPORT GaussLaguerreQuadrature
{
public:
    GaussLaguerreQuadrature();
    virtual ~GaussLaguerreQuadrature() {}

    void SetInterestZone(double minVal, double maxVal);

    template <class FunctionType>
    double GetIntegralValue(FunctionType integrand)
    {
        double deltaPart = 0.0;
        if (m_DeltaValue > 0.0)
            deltaPart = boost::math::quadrature::gauss <double, 15>::integrate(integrand, 0.0, m_DeltaValue);

        double laguerrePart = 0.0;
        for (unsigned int i = 0;i < m_Abcissas.size();++i)
            laguerrePart += m_ValueWeights[i] * integrand(m_ScaleValue * m_Abcissas[i] + m_DeltaValue) * m_ScaleValue;

        return deltaPart + laguerrePart;

    }

private:
    double m_DeltaValue;

    double m_ScaleValue;

    const std::vector <double> m_ValueWeights =
    {
        2.25030314864247252e-1,
        5.25836052762342454e-1,
        8.31961391687087088e-1,
        1.14609924096375170,
        1.47175131696680859,
        1.81313408738134816,
        2.17551751969460745,
        2.56576275016502921,
        2.99321508637137516,
        3.47123448310209029,
        4.02004408644466886,
        4.67251660773285420,
        5.48742065798615247,
        6.58536123328921366,
        8.27635798436423448,
        1.18242775516584348e1
    };

    const std::vector <double> m_Abcissas =
    {
        8.76494104789278403e-2,
        4.62696328915080832e-1,
        1.14105777483122686,
        2.12928364509838062,
        3.43708663389320665,
        5.07801861454976791,
        7.07033853504823413,
        9.43831433639193878,
        1.22142233688661587e1,
        1.54415273687816171e1,
        1.91801568567531349e1,
        2.35159056939919085e1,
        2.85787297428821404e1,
        3.45833987022866258e1,
        4.19404526476883326e1,
        5.17011603395433184e1
    };
};

} // end of namespace anima