linepath_d1.h

Go to the documentation of this file.

#ifndef LINEPATHD1_H
#define LINEPATHD1_H

//#include <cassert>
//#include <iostream>

//using namespace std;

#include <typeop.h>

template< typename FN, typename XI, typename X, typename T >
class linepath_d1
{
public:

  uintc dim;

  typename typeop<FN>::Tref fn;
  XI xi;
  X x0;
  X di;

  linepath_d1
  (
    uintc dim_, 
    typename typeop<FN>::Tref fn_, 
    XI xi_, 
    X x0_, 
    X di_ 
  );

  void eval(T const t)
  {
    for (uint i=0; i<dim; ++i)
      xi[i] = x0[i]+di[i]*t;
//cout << "***";
//cout << printvecfunc(xi,dim) << endl;
  }
  void eval(T& fval, T const t)
    { eval(t); fn(fval); }
  //void fneval(T& fval)
  //  { fn(fval); }

  ostream & print(ostream & os) const;

};


//-----------------------------------------------
//  Implementation


template< typename FN, typename XI, typename X, typename T >
ostream & linepath_d1<FN,XI,X,T>::print(ostream & os) const
{
  os << "dim=" << dim << endl;
  os << "xi[]=" << printvecfunc(xi,dim) << endl;
  os << "x0[]=" << printvecfunc(x0,dim) << endl;
  os << "di[]=" << printvecfunc(di,dim) << endl;

  return os;
}

template< typename FN, typename XI, typename X, typename T >
linepath_d1<FN,XI,X,T>::linepath_d1
(
  uintc dim_, 
  typename typeop<FN>::Tref fn_, 
  XI xi_, 
  X x0_, 
  X di_ 
)
  : dim(dim_), fn(fn_), xi(xi_), x0(x0_), di(di_)
{
}






#endif