powerseries.cpp

Go to the documentation of this file.

#include <print.h>

#include <powerseries.h>


void powerseries::operator()(double & val, doublec t) const
{
  // ai[0] + t(ai[1] + t(ai[2] + t(...
  val = 0.0;
  for (uint i=0; i<N; ++i)
  {
    val *= t;
    val += ai[N-1-i];
  }
}

doublec powerseries::operator()(doublec t) const
{
  // ai[0] + t(ai[1] + t(ai[2] + t(...
  double val = 0.0;
  for (uint i=0; i<N; ++i)
  {
    val *= t;
    val += ai[N-1-i];
  }

  return val;
}

void powerseries::differentiate()
{
  for (uint i=1; i<N; ++i)
  {
    ai[i]*=i;
    ai[i-1] = ai[i];
  }
  ai[N-1] = 0.0;
}