src/INTERP_KERNEL/DirectedBoundingBox.hxx

Go to the documentation of this file.

// Copyright (C) 2009-2011  OPEN CASCADE
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License.
//
// This library 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
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
//
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//

#ifndef __DIRECTEDBOUNDINGBOX_HXX__
#define __DIRECTEDBOUNDINGBOX_HXX__

#include "INTERPKERNELDefines.hxx"

#include <vector>

namespace INTERP_KERNEL
{

  class INTERPKERNEL_EXPORT DirectedBoundingBox
  {
  public:

    DirectedBoundingBox();

    DirectedBoundingBox(const double* pts, const unsigned numPts, const unsigned dim);

    DirectedBoundingBox(const double** pts, const unsigned numPts, const unsigned dim);

    //~DirectedBoundingBox();

    void enlarge(const double tol);
    
    bool isDisjointWith(const DirectedBoundingBox& box) const;

    bool isDisjointWith(const double* box) const;

    bool isOut(const double* point) const;


    // return internal data
    std::vector<double> getData() const;

    // initialize with data returned by getData()
    void setData(const double* data);

    // return size of internal data
    static int dataSize(int dim);

  private:

    //void computeAxes3D(const std::vector<double>& tensor);

    //void computeAxes2D(const std::vector<double>& tensor);

    inline void addPointToBox(const double* coord);

    void toLocalCS(const double* p, double* pLoc) const;

    void fromLocalCS(const double* p, double* pGlob) const;

    inline bool isLocalOut(const double* pLoc) const;

    void getCorners(std::vector<double>& corners, const double* minmax) const;

    unsigned _dim;

    std::vector<double> _axes; 
    std::vector<double> _minmax; 

  };

  //================================================================================
  //================================================================================

  inline bool DirectedBoundingBox::isLocalOut(const double* pLoc) const
    {
      for ( int i = 0; i < (int)_dim; ++i )
        if ( pLoc[i] < _minmax[i*2] || pLoc[i] > _minmax[i*2+1] )
          return true;
      return false;
    }

  //================================================================================
  //================================================================================

  inline void DirectedBoundingBox::addPointToBox(const double* coord)
  {
    for ( int i = 0; i < (int)_dim; ++i )
      {
        double c = 0;
        for ( int j = 0; j < (int)_dim; ++j ) c += coord[j]*_axes[i*_dim+j];
        if ( c < _minmax[i*2] )   _minmax[i*2] = c;
        if ( c > _minmax[i*2+1] ) _minmax[i*2+1] = c;
      }
  }
}
#endif