---

func.cpp

//$Header: /home/ben/Mapper/c++/RCS/func.cpp,v 6.3 2002/06/14 22:28:47 ben Exp $
// Copyright Benedict Adamson 2002.
// This file is part of Mapper.

// Mapper is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.

// Mapper is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.

// You should have received a copy of the GNU General Public License
// along with Mapper; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

#include "func.h"
#include "vector_math.h"



func::projection::projection(
   const vector<double> &ix0,
   const vector<double> &idxdt,
   const multi_dimensional *f
   )
   :
   x0( ix0 ),
   dxdt( idxdt ),
   multi_dimensional_function( f )
{
   {//postconditions:
      assert( this->x0 == ix0 );
      assert( this->dxdt == idxdt );
      assert( this->multi_dimensional_function == f );
   };
}



func::projection::~projection(void)
{
   //Do nothing
}



vector< double > func::projection::line(
   const double t
   ) const
{
   return x0 + t*dxdt;
}




double func::projection::operator() (
   const double t
   ) const
{
   {//preconditions:
      assert( this->multi_dimensional_function );
   };
   //(*this)( t ) == (*(this->multi_dimensional_function))( x0 + t*dxdt );
   const multi_dimensional &f(
      *(this->multi_dimensional_function)
      );
   vector<double> xt( this->line( t ) );
   return f( xt );
}

---