G4KDNode.icc

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
// Author: Mathieu Karamitros (kara (AT) cenbg . in2p3 . fr) 
//
// History:
// -----------
// 10 Oct 2011 M.Karamitros created
//
// -------------------------------------------------------------------
 
//______________________________________________________________________
template<typename PointT>
  G4KDNode<PointT>::G4KDNode(G4KDTree* tree,
                             PointT* point,
                             G4KDNode_Base* parent) :
      G4KDNode_Base(tree, parent)
  {
    fPoint = point;
    fValid = true;
  }
 
// Copy constructor should not be used
template<typename PointT>
  G4KDNode<PointT>::G4KDNode(const G4KDNode<PointT>& right) :
      G4KDNode_Base(right), fPoint(0)
  {
    fValid = false;
  }
 
// Assignement should not be used
template<typename PointT>
  G4KDNode<PointT>& G4KDNode<PointT>::operator=(const G4KDNode<PointT>& right)
  {
    if(this == &right) return *this;
    fPoint = right.fPoint;
    fTree = right.fTree;
    fLeft = right.fLeft;
    fRight = right.fRight;
    fParent = right.fParent;
    fSide = right.fSide;
    fAxis = right.fAxis;
    return *this;
  }
 
template<typename PointT>
  G4KDNode<PointT>::~G4KDNode()
  {
  }
 
template<typename Position>
  G4KDNode_Base* G4KDNode_Base::FindParent(const Position& x0)
  {
    G4KDNode_Base* aParent = 0;
    G4KDNode_Base* next = this;
    int split = -1;
    while(next)
    {
      split = next->fAxis;
      aParent = next;
 
      // works if node called "next" is valid
      if(x0[split] > (*next)[split]) next = next->fRight;
      else next = next->fLeft;
    }
    return aParent;
  }
 
template<typename PointT>
  G4KDNode_Base* G4KDNode_Base::Insert(PointT* point)
  {
    G4KDNode_Base* aParent = FindParent(*point);
    // TODO check p == aParent->pos
    // Exception
 
    G4KDNode_Base* newNode = new G4KDNode<PointT>(fTree, point, aParent);
 
    if((*point)[aParent->fAxis] > (*aParent)[aParent->fAxis])
    {
      aParent->fRight = newNode;
      newNode->fSide = 1;
    }
    else
    {
      aParent->fLeft = newNode;
      newNode->fSide = -1;
    }
 
    return newNode;
  }
 
template<typename PointT>
  G4KDNode_Base* G4KDNode_Base::Insert(const PointT& point)
  {
    G4KDNode_Base* aParent = FindParent(point);
    // TODO check p == aParent->pos
    // Exception
 
    G4KDNode_Base* newNode = new G4KDNodeCopy<PointT>(fTree, point, aParent);
 
    if(point[aParent->fAxis] > (*aParent)[aParent->fAxis])
    {
      aParent->fRight = newNode;
      newNode->fSide = 1;
    }
    else
    {
      aParent->fLeft = newNode;
      newNode->fSide = -1;
    }
 
    return newNode;
  }