sqrt.hpp

Go to the documentation of this file.

#pragma once
 
#include <cmath>
#include <complex>
 
#include "NumCpp/Core/Internal/StlAlgorithms.hpp"
#include "NumCpp/NdArray.hpp"
 
namespace nc
{
    //============================================================================
    // Method Description:
    template<typename dtype>
    auto sqrt(dtype inValue) noexcept
    {
        STATIC_ASSERT_ARITHMETIC_OR_COMPLEX(dtype);
 
        return std::sqrt(inValue);
    }
 
    //============================================================================
    // Method Description:
    template<typename dtype>
    auto sqrt(const NdArray<dtype>& inArray)
    {
        NdArray<decltype(sqrt(dtype{ 0 }))> returnArray(inArray.shape());
        stl_algorithms::transform(inArray.cbegin(),
                                  inArray.cend(),
                                  returnArray.begin(),
                                  [](dtype inValue) noexcept -> auto { return sqrt(inValue); });
 
        return returnArray;
    }
} // namespace nc