Neighborhood_Criterion.h

#ifndef QNANO_NEW_NEIGHBORHOOD_CRITERION_DEFINED_H
#define QNANO_NEW_NEIGHBORHOOD_CRITERION_DEFINED_H

#include "structure/Atom.h"

class Neighborhood_Criterion{
public:

  virtual bool is_neighbor(const Atom &a1, const Atom &a2) const =0;
  bool operator() (const Atom &a1, const Atom &a2) const {return is_neighbor(a1,a2);}

  virtual double upper_bound_radius() const=0;
};

class Neighborhood_Criterion_Distance: public Neighborhood_Criterion{
public: 
  double distance;
  
  Neighborhood_Criterion_Distance(double distance_) :distance(distance_){
  }

  virtual bool is_neighbor(const Atom &a1, const Atom &a2) const{
    if(dist(a1,a2)<distance)return true;
    else return false;
  } 
  virtual double upper_bound_radius() const{ return distance;}

};


#endif