divmod.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <cmath>
#include <type_traits>
 
#include "NumCpp/Core/Internal/StaticAsserts.hpp"
#include "NumCpp/Core/Internal/StlAlgorithms.hpp"
#include "NumCpp/Core/Types.hpp"
#include "NumCpp/NdArray.hpp"
 
namespace nc
{
    //===========================================================================
    // Method Description:
    template<typename dtype>
    std::pair<NdArray<dtype>, NdArray<dtype>> divmod(const NdArray<dtype>& inArray1, const NdArray<dtype>& inArray2)
    {
        STATIC_ASSERT_ARITHMETIC(dtype);
 
        if (inArray1.size() != inArray2.size())
        {
            THROW_INVALID_ARGUMENT_ERROR("Arrays must have the same size.");
        }
 
        auto div = NdArray<dtype>(inArray1.shape());
        auto mod = NdArray<dtype>(inArray1.shape());
 
        for (auto i = 0u; i < inArray1.size(); ++i)
        {
            if constexpr (std::is_floating_point_v<dtype>)
            {
                div[i] = std::floor(inArray1[i] / inArray2[i]);
                mod[i] = std::fmod(inArray1[i], inArray2[i]);
            }
            else
            {
                div[i] = inArray1[i] / inArray2[i];
                mod[i] = inArray1[i] % inArray2[i];
            }
        }
 
        return std::make_pair(div, mod);
    }
} // namespace nc