ellint_3.hpp

Go to the documentation of this file.

#pragma once
 
#include <cmath>
 
#if defined(__cpp_lib_math_special_functions) || !defined(NUMCPP_NO_USE_BOOST)
 
#include "NumCpp/Core/Internal/Error.hpp"
#include "NumCpp/Core/Internal/StaticAsserts.hpp"
#include "NumCpp/Core/Internal/StlAlgorithms.hpp"
#include "NumCpp/NdArray.hpp"
 
#ifndef __cpp_lib_math_special_functions
#include "boost/math/special_functions/ellint_3.hpp"
#endif
 
#include <type_traits>
 
namespace nc::special
{
    //============================================================================
    // Method Description:
    template<typename dtype1, typename dtype2, typename dtype3>
    auto ellint_3(dtype1 inK, dtype2 inV, dtype3 inP)
    {
        STATIC_ASSERT_ARITHMETIC(dtype1);
        STATIC_ASSERT_ARITHMETIC(dtype2);
        STATIC_ASSERT_ARITHMETIC(dtype3);
 
#ifdef __cpp_lib_math_special_functions
        return std::ellint_3(inK, inV, inP);
#else
        return boost::math::ellint_3(inK, inV, inP);
#endif
    }
 
    //============================================================================
    // Method Description:
    template<typename dtype1, typename dtype2, typename dtype3>
    auto ellint_3(const NdArray<dtype1>& inArrayK, const NdArray<dtype2>& inArrayV, const NdArray<dtype3>& inArrayP)
    {
        if (inArrayK.size() != inArrayV.size() || inArrayK.size() != inArrayP.size())
        {
            THROW_INVALID_ARGUMENT_ERROR("Shapes of inArrayK, inArrayV, and inArrayP must match.");
        }
 
        NdArray<decltype(ellint_3(dtype1{ 0 }, dtype2{ 0 }, dtype3{ 0 }))> returnArray(inArrayK.shape());
 
        for (uint32 i = 0; i < inArrayK.size(); ++i)
        {
            returnArray[i] = ellint_3(inArrayK[i], inArrayV[i], inArrayP[i]);
        }
 
        return returnArray;
    }
} // namespace nc::special
 
#endif // #if defined(__cpp_lib_math_special_functions) || !defined(NUMCPP_NO_USE_BOOST)