poly.h

Go to the documentation of this file.

/*
  Copyright (c) Mircea Neacsu (2014-2025) Licensed under MIT License.
  This file is part of MLIB project. See LICENSE file for full license terms.
*/

#pragma once
 
#if __has_include("defs.h")
#include "defs.h"
#endif
 
#include <array>
#include <vector>
 
namespace mlib {
 
template <typename T>
T poly (T x, const T* coeff, int n)
{
  T val = coeff[n - 1];
  for (int i = n - 2; i >= 0; i--)
  {
    val *= x;
    val += coeff[i];
  }
  return val;
}

template <typename T, size_t N>
T poly (T x, std::array<T, N> coeff)
{
  return poly (x, coeff.data (), (int)N);
}

template <typename T>
T poly (T x, std::vector<T> coeff)
{
  return poly (x, coeff.data (), (int)coeff.size ());
}
 
} // namespace mlib