ellint_rf

auto eve::ellint_rf = functor<ellint_rf_t>

inlineconstexpr

Header file

#include <eve/module/elliptic.hpp>

Callable Signatures

namespace eve
{
   // Regular overload
   constexpr auto ellint_rf(floating_value auto x, floating_value auto y, floating_value auto z)   noexcept; // 1
 
   // Lanes masking
   constexpr auto ellint_rf[conditional_expr auto c](floating_value auto x, floating_value auto y,
                                                     floating_value auto z)                        noexcept; // 2
   constexpr auto ellint_rf[logical_value auto m](floating_value auto x, floating_value auto y,
                                                  floating_value auto z)                           noexcept; // 2
}

Parameters

• x, y, z: floating values. All those arguments must be non-negative and at most one zero or the the result is NaN.
• c: Conditional expression masking the operation.
• m: Logical value masking the operation.

Return value

1. the value of the \(\mathbf{R}_\mathbf{F}\) Carlson elliptic integral: \( \mathbf{R}_\mathbf{F}(x, y) = \frac32 \int_{0}^{\infty} \scriptstyle[(t+x)(t+y)]^{-1/2} (t+z)^{-3/2}\scriptstyle\;\mathrm{d}t\). is returned
2. The operation is performed conditionnaly

External references

• DLMF: Elliptic Integral
• Wolfram MathWorld: Elliptic Integral

Example

// revision 1
#include <eve/module/elliptic.hpp>
#include <iostream>
 
eve::wide pf{1.0f, 0.0f, 1.5f, 3.0f};
eve::wide qf{1.0f, 4.0f, 0.2f, 0.5f};
eve::wide rf{2.0f, 1.0f, 0.1f, 0.4f};
 
int main()
{
  std::cout << "<- pf = " << pf << "\n";
  std::cout << "<- qf = " << qf << "\n";
  std::cout << "<- rf = " << rf << "\n";
 
  std::cout << "-> ellint_rf(pf, qf, rf)                = " << eve::ellint_rf(pf, qf, rf) << "\n";
  std::cout << "-> ellint_rf[ignore_last(2)](pf, qf, rf)= " << eve::ellint_rf[eve::ignore_last(2)](pf, qf, rf) << "\n";
  std::cout << "-> ellint_rf[qf != 4.0f](pf, qf, rf)    = " << eve::ellint_rf[qf != 4.0f](pf, qf, rf) << "\n";
}