prime.hpp

Go to the documentation of this file.

#pragma once
 
#ifndef NUMCPP_NO_USE_BOOST
 
#include <string>
 
#include "boost/math/special_functions/prime.hpp"
 
#include "NumCpp/Core/Internal/Error.hpp"
#include "NumCpp/Core/Internal/StlAlgorithms.hpp"
#include "NumCpp/Core/Types.hpp"
#include "NumCpp/NdArray.hpp"
 
namespace nc::special
{
    //============================================================================
    // Method Description:
    inline uint32 prime(uint32 n)
    {
        if (n > boost::math::max_prime)
        {
            THROW_INVALID_ARGUMENT_ERROR("input n must be less than or equal to " +
                                         std::to_string(boost::math::max_prime));
        }
 
        return boost::math::prime(n);
    }
 
    //============================================================================
    // Method Description:
    inline NdArray<uint32> prime(const NdArray<uint32>& inArray)
    {
        NdArray<uint32> returnArray(inArray.shape());
 
        stl_algorithms::transform(inArray.cbegin(),
                                  inArray.cend(),
                                  returnArray.begin(),
                                  [](uint32 inValue) -> uint32 { return prime(inValue); });
 
        return returnArray;
    }
} // namespace nc::special
 
#endif // #ifndef NUMCPP_NO_USE_BOOST