app-radu/tree03.hh

Go to the documentation of this file.

#ifndef appRaduTree03_hh
#define appRaduTree03_hh

#include "cluster/tree03.hh"

// ************************************************************ CebyCenter3d **

template<class Node>
class CebyCenter3d {
public:
  typedef cluster::TreeTraits03<Node> Traits;
  typedef typename Traits::F F;

  inline CebyCenter3d(uint nvtx);
  inline ~CebyCenter3d() {delete[] vtx_;}

  void operator()(const Node* lfset, concepts::Real3d& c, concepts::Real& r,
      const concepts::Real3d* (&e)[4], uint& ne);

private:
  static concepts::Real eps_;
  uint nvtx_;
  const concepts::Real3d** vtx_;
  __gnu_cxx::hash_map<uint, uint> vtxtag_;
  uint tag_;
};

template<class Node>
inline CebyCenter3d<Node>::CebyCenter3d(uint nvtx)
  : nvtx_(nvtx), vtxtag_(), tag_(0) {
  vtx_ = new const concepts::Real3d*[nvtx_];
}

// **************************************************************** Tree3d03 **

template<class CNode>
class Tree3d03 : public cluster::Tree<CNode> {
public:
  typedef cluster::TreeTraits03<CNode> Traits;
  typedef typename Traits::F  F;

  Tree3d03(const concepts::Space<F>& spc, const cluster::BBall<F>& ball,
     uint nlfmin = 1);
  inline ~Tree3d03() {destructor_(root_);  delete root_;}

  inline uint nclst() const {return nclst_;}
  inline uint nleaf() const {return nlf_;}
  inline const CNode* root() const {return root_;}

  void sketch(std::ostream& os, const CNode* lnk = 0, uint l = 0);

private:
  uint nlfmin_;
  uint nlf_;
  uint nclst_;
  CNode* root_;

  CNode* constructor_(CNode* lfset, CebyCenter3d<CNode>& cby);
  void destructor_(CNode* clst);
  CNode* concat_(CNode* a, CNode* b);
};

#endif