Feature.hpp

/* $Id$ */

#ifndef _SPAF_FEATURE_HPP_
#define _SPAF_FEATURE_HPP_

#include "jafarConfig.h"
#include <values.h>

#include <iostream>
#include <geom/Declarations.hpp>
#include <geom/Atom.hpp>

namespace jafar {

  namespace spaf {
    class ComparisonResult;
    class Location;
    class LocationsBag;
    class DisplayContext;
    class MapsTransformation;
    class Feature {
      public:
        struct Range {
          Range() : a(FLT_MIN), aCov(0), b(FLT_MAX), bCov(0)
          {
          }
          Range(double _a, double _aCov, double _b, double _bCov) : a(_a), aCov(_aCov), b(_b), bCov(_bCov)
          {
          }
          double a, aCov, b, bCov;
        };
      public:
        Feature();
        Feature(const Feature&);
        virtual ~Feature();
        virtual int id() const = 0;
        const LocationsBag& locations() const;
        virtual std::list<int> maps() const = 0;
        virtual const std::map<int, geom::Atom3D*>& atoms() const = 0;
        virtual geom::Atom3D* atom(int _map) =0;
        virtual const geom::Atom3D* atom(int _map) const =0;
        virtual void setAtom(const geom::Atom3D* _atom, int _map) = 0;
        virtual void fillDisplayContext( DisplayContext*, MapsTransformation* mapsTransfo ) const = 0;
        virtual void computeCharacteristicVector( const Feature* feature, geom::Atom3D::HomogenousVecD& vec, geom::Atom3D::HomogenousSymMatrixD& vecCov ) const = 0;
        virtual Feature::Range computeOccupancyRange( const Feature* _feature ) const = 0;
        virtual void computeCosAngle( const Feature* feature, double& angle, double& angleCov ) const = 0;
        virtual void computeDistance( const Feature* feature, double& distance, double& distanceCov ) const = 0;
        virtual void mergeWith( Feature* feature) = 0;
        virtual geom::BoundingBox3D globalBoundingBox(MapsTransformation* mapsTransfo) const = 0;
      private:
        struct Private;
        Private* const d;
    };
    std::ostream& operator<<(std::ostream& s, const Feature&);
  }
}

#endif