vecinterp.h

Go to the documentation of this file.

#ifndef VECINTERP_H
#define VECINTERP_H

#include <print.h>

#include <powerseries1D.h>

template< typename T >
class vecInterp2D0
{
public:
  
  T p0;
  T p1;

  vecInterp2D0(T const & _p0, T const & _p1 ) 
    : p0(_p0), p1(_p1) {}

  void operator()( T & val, doublec t) const
    { val = p0*(1.0-t) + p1*t; }

  T const operator()(doublec t) const
    { T val = p0*(1.0-t) + p1*t; return val; }

};


template< typename T >
class vecInterp2D1
{
  static double A0[];
  static double B0[];
  static double A1[];
  static double B1[];
  static powerseries1D<double> fA0;
  static powerseries1D<double> fB0;
  static powerseries1D<double> fA1;
  static powerseries1D<double> fB1;
public:
  
  T p0;
  T dp0;
  T p1;
  T dp1;

  vecInterp2D1
  (
    T const & _p0,
    T const & _dp0,
    T const & _p1,
    T const & _dp1
  ) : p0(_p0), dp0(_dp0), p1(_p1), dp1(_dp1) 
    {} 
  
  void operator()(T & val, doublec t) const
    { val = p0*fA0(t) + dp0*fB0(t) + p1*fA1(t) + dp1*fB1(t); }

  T const operator()(doublec t) const
    { return p0*fA0(t) + dp0*fB0(t) + p1*fA1(t) + dp1*fB1(t); }
    
};

template< typename T >
double vecInterp2D1<T>::A0[] = { 1.0, 0.0, -3.0, 2.0 };

template< typename T >
double vecInterp2D1<T>::B0[] = { 0.0, 1.0, -2.0, 1.0 };

template< typename T >
double vecInterp2D1<T>::A1[] = { 0.0, 0.0, 3.0, -2.0 };

template< typename T >
double vecInterp2D1<T>::B1[] = { 0.0, 0.0, -1.0, 1.0 };

template< typename T >
powerseries1D<double> vecInterp2D1<T>::fA0(A0,4);

template< typename T >
powerseries1D<double> vecInterp2D1<T>::fB0(B0,4);

template< typename T >
powerseries1D<double> vecInterp2D1<T>::fA1(A1,4);

template< typename T >
powerseries1D<double> vecInterp2D1<T>::fB1(B1,4);



#endif