mapping.h

#ifndef MAPPING__H
#define MAPPING__H


#include <vector>
#include "eigenstuff.h"
#include "countedref.h"

class FatPoint;
class Point;

class Mapping  : public RefCount
{

    public :
    
    virtual unsigned Npar() const = 0;


    virtual void GetMappingIndices(std::vector<unsigned> &Indices) const = 0;


    virtual void ComputeTransformAndDerivatives(const FatPoint &Where,
                        FatPoint &OutPos,
                        Eigen::MatrixX2d &H) const = 0;
    virtual void TransformPosAndErrors(const FatPoint &Where,
                       FatPoint &OutPos) const = 0;


    virtual void  PosDerivative(Point &Where, Eigen::Matrix2d &Der, const double & Eps) const = 0;

    virtual ~Mapping() {};
    
};

// typedef std::list<CountedRef<Mapping> > MappingList;
// typedef MappingList::iterator MappingIterator;
// typedef MappingList::const_iterator MappingCIterator;


#endif /* MAPPING__H */