vectorfield.h

Go to the documentation of this file.

#ifndef __VECTORFIELD_H
#define __VECTORFIELD_H

#include <boost/numeric/ublas/vector.hpp>
#include <boost/multi_array.hpp>
#include <functional>

namespace mimas {

template< class Tensor >
class vectorfield: public std::unary_function
   < boost::numeric::ublas::vector< double >,
     boost::numeric::ublas::vector< double > >
{
 public:
  typedef boost::numeric::ublas::vector< double > Vector;
  vectorfield( int width, int height, const Tensor &aTensor );
  Vector operator()( const Vector &v ) const {
    int
      height = data.shape()[ 0 ],
      width = data.shape()[ 1 ],
      x = (int)rint( v( 0 ) + width * 0.5 ),
      y = (int)rint( v( 1 ) + height * 0.5 );
    assert( x >= 0 && x < width );
    assert( y >= 0 && y < height );
    return data[ y ][ x ];
  };
 protected:
  boost::multi_array< Vector, 2 > data;
};

};

#include "vectorfield.tcc"

#endif

---