text: /// @file
  /// @author David Pilger <dpilger26@gmail.com>
  /// [GitHub Repository](https://github.com/dpilger26/NumCpp)
  /// @version 1.3
  ///
  /// @section License
  /// Copyright 2020 David Pilger
  ///
  /// Permission is hereby granted, free of charge, to any person obtaining a copy of this
  /// software and associated documentation files(the "Software"), to deal in the Software
  /// without restriction, including without limitation the rights to use, copy, modify,
  /// merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
  /// permit persons to whom the Software is furnished to do so, subject to the following
  /// conditions :
  ///
  /// The above copyright notice and this permission notice shall be included in all copies
  /// or substantial portions of the Software.
  ///
  /// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
  /// INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR
  /// PURPOSE AND NONINFRINGEMENT.IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE
  /// FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR
  /// OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
  /// DEALINGS IN THE SOFTWARE.
  ///
  /// @section Description
  /// Special Functions
  ///
  #pragma once
  
  #include "NumCpp/NdArray.hpp"
  #include "NumCpp/Core/StlAlgorithms.hpp"
  
  #include "boost/math/special_functions/zeta.hpp"
  
  namespace nc
  {
      namespace special
      {
          //============================================================================
          // Method Description:
          ///	The Riemann Zeta function
          /// https://en.wikipedia.org/wiki/Riemann_zeta_function
          ///
          /// @param
          ///				inValue
          /// @return
          ///				double
          ///
          template<typename dtype>
          double riemann_zeta(dtype inValue) noexcept
          {
              return boost::math::zeta(static_cast<double>(inValue));
          }
  
          //============================================================================
          // Method Description:
          ///	The Riemann Zeta function
          /// https://en.wikipedia.org/wiki/Riemann_zeta_function
          ///
          /// @param
          ///				inArray
          /// @return
          ///				NdArray<double>
          ///
          template<typename dtype>
          NdArray<double> riemann_zeta(const NdArray<dtype>& inArray) noexcept
          {
              NdArray<double> returnArray(inArray.shape());
  
              stl_algorithms::transform(inArray.cbegin(), inArray.cend(), returnArray.begin(),
                  [](dtype inValue) noexcept -> double
                  { 
                      return riemann_zeta(inValue); 
                  });
  
              return returnArray;
          }
      }
  }
  