lib_bessel.cpp

#include "../util/math_functions.hpp"
 
#include <boost/math/constants/constants.hpp>
 
namespace NiHu
{
namespace bessel
{
    template <>
    std::complex<double> K<0>(std::complex<double> const &z)
    {
        using namespace std::complex_literals;
        using namespace boost::math::double_constants;
        return half_pi * 1.i * (std::imag(z) < 0. ? H<0, 1>(1.i*z) : -H<0, 2>(-1.i*z));
    }
 
    template <>
    std::complex<double> K<1>(std::complex<double> const &z)
    {
        using namespace std::complex_literals;
        using namespace boost::math::double_constants;
        return -half_pi * (std::imag(z) < 0. ? H<1, 1>(1.i*z) : H<1, 2>(-1.i*z));
    }
}
}