bessel_y

auto eve::bessel_y = eve::functor<bessel_y_t>

inlineconstexpr

#include <eve/functions.hpp>

Callable Signatures

namespace eve
{
   // Regular overload
  constexpr auto bessel_y(value auto n, floating_value auto z)                  noexcept; // 1
 
   // Semantic modifier
   constexpr auto bessel_y[cylindrical](value auto n, floating_value auto z)    noexcept; // 1
   constexpr auto bessel_y[spherical](value auto n, floating_value auto z)      noexcept; // 2
 
   // Lanes masking
   constexpr auto bessel_y[conditional_expr auto c](/*any previous overload*/)  noexcept; // 3
   constexpr auto bessel_y[logical_value auto m](/*any previous overload*/)     noexcept; // 3
}

Parameters

• n: order (integral or floating)
• z: Value to process.
• m: Logical value masking the operation.

Return value

1. returns \(Y_n\)(z) (cylindrical).
2. returns \(y_n\)(z) (spherical).
3. The operation is performed conditionally.

Note

if n is not integral nor flint negative z entries produce NaN outputs as the correct result is not real. You may want to use KYOSU to avoid this behavior and produce complex outputs (that are not available in EVE).

External references

• Wolfram MathWorld: Bessel Function of the Second Kind
• Wolfram MathWorld: Spherical Bessel Function of the Second Kind
• Wikipedia: Bessel function
• DLMF: Bessel functions

Example

// revision 1
#include <eve/module/bessel.hpp>
#include <iostream>
 
eve::wide<double> wdf([](auto i, auto c)->double{ return 2*(i-c/2)+2/3.0;});
eve::wide<std::uint64_t> wu([](auto i, auto )->std::uint64_t{ return 2*i+1;});
eve::wide wd{0.5, 1.5, 0.1, 1.0, 19.0, 0.0, 0.0, 10000.0};
eve::wide m{0.5, -1.0, 1.5, -2.0, 2.5, -2.6, 3.2, -12.0};
 
int main()
{
   std::cout << "<- wdf = " << wdf << "\n";
   std::cout << "<- wu = " << wu << "\n";
   std::cout << "<- m  = " << m << "\n";
   std::cout << "<- wd  = " << wd << "\n";
 
   std::cout << "-> bessel_y(wu, wdf)                = " << eve::bessel_y(wu, wdf) << "\n";
   std::cout << "-> bessel_y[ignore_last(2)](wu, wdf)= " << eve::bessel_y[eve::ignore_last(2)](wu, wdf) << "\n";
   std::cout << "-> bessel_y[wu != 3u](wu, wdf)      = " << eve::bessel_y[wu != 3u](wu, wdf) << "\n";
   std::cout << "-> bessel_y(m, wd)                  = " << eve::bessel_y(m, wd) << "\n";
   std::cout << "-> bessel_y[spherical](m, wd)       = " << eve::bessel_y[eve::spherical](m, wd) << "\n";
}