dupl.h

Go to the documentation of this file.

#ifndef MIDAS_DUPL_H
#define MIDAS_DUPL_H

#include "alias.h"

#include "gelib.h"
#include "mathlib.h"
#include "arrays.h"


namespace midaspace {

class DuplicatePoints
{
  static const bool PRINT = false;
  
  long np;                              // number of points
  const PoinT **points;                 // array of points
  double averdist;                      // average distance between points
  double Zero;                          // pocitacova nula--zaokrouhlovani == priblizne 1e-10
  double Zero_dist;                     // maximal distance between duplicite points
  
  double mincoo[3], maxcoo[3], dcoo[3]; // boundary box
  long cc[3];                           // count of cells in xyz axis directions
  int  ****ccellpoints;                 // count of points in one cell
  long ****cellpoints;                  //    id of points in one cell
  
 public:
  DuplicatePoints (long n, double ad, double zero) {
    np = n;
    points = new const PoinT* [np];
    
    averdist = ad;
    Zero = zero;
    Zero_dist = averdist/500.0;
  }
  ~DuplicatePoints () { delete [] points; }
  
  void set_point (long i, const PoinT* p) { points[i] = p; }
  void assign_cellpoints (void);
  void find_duplicitys (long &nd, long *Cduplicity, long **duplicity);
  
 private:
  void add_point_to_cell (long id, long pos[]);
  bool point_is_in_cell (long id, long pos[], double offset);
  
};

} // namespace midaspace

#endif // MIDAS_DUPL_H