powerseriesWD.h

Go to the documentation of this file.

#ifndef POWERSERIESWD_H
#define POWERSERIESWD_H

#include <cassert>
using namespace std;

#include <powerseries1D.h>

template< typename T >
class powerseriesWD
{
public:

  T * const mat;

  uintc W;
  uintc N;

  powerseries1D<T> currentP;
  
  void currentSet(uintc i);

  powerseriesWD( T * const mat_, uintc W_, uintc N_ );

  void operator()(T & val, T const t) const
    { currentP(val,t); }

  T const operator()(T const t) const
    { return currentP(t); }

  
};

template< typename T >
void powerseriesWD<T>::currentSet(uintc i)
{
  assert( (i==0) || (i<N) );
  currentP.construct(mat+i*N,N);
}



template< typename T >
powerseriesWD<T>::powerseriesWD
( 
  T * const mat_, 
  uintc W_, 
  uintc N_
)
  : mat(mat_), W(W_), N(N_)
{ 
  currentP.construct(mat_,N);
}


#endif