VectorMath.hpp

Go to the documentation of this file.

#ifndef SOLIDUTILS_INCLUDE_VECTORMATH_HPP
#define SOLIDUTILS_INCLUDE_VECTORMATH_HPP


#include "Array.hpp"

#include <type_traits>


namespace sl
{


class VectorMath
{
  public:
    template <typename T>
    static T sum(
        T * const data,
        size_t const size) noexcept
    {
      T sum = 0;
      for (size_t i = 0; i < size; ++i) {
        sum += data[i];
      }
      return sum;
    }


    template <typename T>
    static void increment(
        T * const data,
        size_t const size,
        T const start = 0,
        T const inc = 1) noexcept
    {
      T val = start;
      for (size_t i = 0; i < size; ++i) {
        data[i] = val;
        val += inc;
      }
    }


    template <typename T>
    static void prefixSumExclusive(
        T * const data,
        size_t const size) noexcept
    {
      prefixSumExclusive(data, data+size);
    }


    template <typename T>
    static void prefixSumExclusive(
        T const begin,
        T const end) noexcept
    {
      using V = typename std::remove_reference<decltype(*begin)>::type;
      V sum = 0;
      for (T iter = begin; iter != end; ++iter) {
        V const val = *iter;
        *iter = sum;
        sum += val;
      }
    }
};




}




#endif