PrCellStructure.h

//===========================================================================
// GoTools - SINTEF Geometry Tools version 1.1
//
// GoTools module: parametrization
//
// Copyright (C) 2000-2005 SINTEF ICT, Applied Mathematics, Norway.
//
// This program is free software; you can redistribute it and/or          
// modify it under the terms of the GNU General Public License            
// as published by the Free Software Foundation version 2 of the License. 
//
// This program is distributed in the hope that it will be useful,        
// but WITHOUT ANY WARRANTY; without even the implied warranty of         
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the          
// GNU General Public License for more details.                           
//
// You should have received a copy of the GNU General Public License      
// along with this program; if not, write to the Free Software            
// Foundation, Inc.,                                                      
// 59 Temple Place - Suite 330,                                           
// Boston, MA  02111-1307, USA.                                           
//
// Contact information: E-mail: tor.dokken@sintef.no                      
// SINTEF ICT, Department of Applied Mathematics,                         
// P.O. Box 124 Blindern,                                                 
// 0314 Oslo, Norway.                                                     
//
// Other licenses are also available for this software, notably licenses
// for:
// - Building commercial software.                                        
// - Building software whose source code you wish to keep private.        
//===========================================================================
#ifndef PRCELLSTRUCTURE_H
#define PRCELLSTRUCTURE_H

#include "Array.h"
using Go::Vector3D;
#include <vector>
using std::vector;
using std::istream;
using std::ostream;

/*<PrCellStructure-syntax: */

class PrCellStructure
{
private:
    vector<Vector3D> xyz_;
    //VecSimplest(CgPoint3d) xyz_;
    int n_; // no points,
    double min_[3]; // min value in each coordinate
    int ncells_[3]; // number of cells in each coordinate
    double cell_size_; // dimension of cell, each cell is a cube
    vector<vector<int> > ind_;
    //VecSimplest(PrIntList) ind_;
    int max_no_cells_;

    void makeCellStructure();

public:
    PrCellStructure() {}
    PrCellStructure(int n, double* xyz_points, int max_no_cells = 10);

    virtual ~PrCellStructure() {};

    void      attach(int n, const double* xyz_points);

    void      setNumCells(int max_no_cells) {max_no_cells_ = max_no_cells; }

    int       getNumCells() const {return max_no_cells_; }

    int       getNumNodes() const {return xyz_.size(); }

    Vector3D get3dNode(int i) const {return xyz_[i]; }

    void       set3dNode(int i, const Vector3D& p) 
    {xyz_[i].x() = p.x(); xyz_[i].y() = p.y(); xyz_[i].z() = p.z();}

    void getBall(const Vector3D& p, double radius,
                 vector<int>& neighbours, int notP = 0) const;

    void getKNearest(const Vector3D& p, int k,
                     vector<int>& neighbours, int notP = 0) const;

    int getI(int i, int j, int k) const
    {return i + ncells_[0] * (j + ncells_[1] * k); }

    void getIJK(int ii, int& i, int& j, int& k) const;

    void whichCell(const Vector3D& xyz, int& i, int& j, int& k) const;

    //print and scan routines
    void print(ostream& os);

    void scan(istream& is);
};

/*>PrCellStructure-syntax: */

/*Class:PrCellStructure

Name:              PrCellStructure
Syntax:            @PrCellStructure-syntax
Keywords:
Description:       This class represents a set of points in three dimensions
                   and a cell structure: a 3D grid of cubical cells
                   each with a local list of those points which lie inside it.
                   This allows one to perform various queries, such as
                   finding all points in a ball, very efficiently.
                   Creating the cell structure requires only O(N) operations,
                   where $N$ is the number of points.
Member functions:

Constructors:
Files:
Example:


See also:
Developed by:      SINTEF Applied Mathematics, Oslo, Norway
Author:            Michael Floater, SINTEF
Date:              Sep. 99
*/

#endif // PRCELLSTRUCTURE_H

---