Pane_connectivity.h

Go to the documentation of this file.

/* *******************************************************************
 * Rocstar Simulation Suite                                          *
 * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
 *                                                                   *
 * Illinois Rocstar LLC                                              *
 * Champaign, IL                                                     *
 * www.illinoisrocstar.com                                           *
 * sales@illinoisrocstar.com                                         *
 *                                                                   *
 * License: See LICENSE file in top level of distribution package or *
 * http://opensource.org/licenses/NCSA                               *
 *********************************************************************/
/* *******************************************************************
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
 * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
 * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
 * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
 * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
 *********************************************************************/
// $Id: Pane_connectivity.h,v 1.15 2008/12/06 08:43:20 mtcampbe Exp $

/* Author: Xiangmin Jiao
 * Date:   Aug. 26, 2002
 */

#ifndef _PANE_CONNECTIVITY_H_
#define _PANE_CONNECTIVITY_H_

#include <list>

#include "mapbasic.h"
#include "roccom_devel.h"
#include "Simple_manifold_2.h"
#include "commpi.h"

MAP_BEGIN_NAMESPACE

class KD_tree_3;

class Pane_connectivity {
  typedef std::vector<const COM::Pane*>  Pane_set;
  typedef MAP::Point_3<double>           Point_3;
public:
  explicit Pane_connectivity( const COM::Attribute *mesh,
                              MPI_Comm c=MPI_COMM_WORLD);

  explicit Pane_connectivity( const Simple_manifold_2 **mani2,
                              MPI_Comm c=MPI_COMM_WORLD);
  
  void compute_pconn( COM::Attribute *pconn_n,
                      COM::Attribute *pconn_f=NULL) throw(int);

  static void size_of_cpanes( const COM::Attribute *pconn, const int *pane_id,
                              int *npanes_total, int *npanes_ghost=NULL);

  static int pconn_nblocks( const COM::Attribute *pconn);

  static int pconn_offset() {
    return _pconn_offset;
  }

protected:

  void create_b2map( COM::Attribute *pconn,
                     bool for_facet) throw(int);

  void create_b2map( std::vector< std::vector<int> > &b2,
                     bool for_facet) throw(int);

  double get_local_boundary_nodes( std::vector<int> &nodes, 
                                   std::vector<Point_3> &pnts,
                                   bool for_facet) throw(int);
  
  double collect_boundary_nodes( std::vector<int>     &nodes, 
                                 std::vector<Point_3> &pnts,
                                 bool for_facet) throw(int);

  // Determine the nodes that coincide with given isolated points.
  // If the last argument is present, then it is mapped to KD_tree.
  void determine_coisolated_nodes( const COM::Pane &pn, 
                                   std::vector<Point_3> &pnts,
                                   const double tol, 
                                   std::vector<bool> &is_co,
                                   KD_tree_3 *tree=NULL) throw(int);
private:
  void collect_nodes( const std::vector<std::vector<int> > &ns,
                      std::vector<int> &nodes, 
                      std::vector<Point_3> &pnts, 
                      bool for_facet);

  void collect_points( const std::vector<std::vector<int> > &ns,
                       std::vector<Point_3> &pnts);

private:
  const COM::Window * const         _win;
  const Pane_set                    _panes;
  const Simple_manifold_2         **_mani2;
  MPI_Comm                          _comm;
  static const int                  _pconn_offset;


};

MAP_END_NAMESPACE

#endif /* _PANE_CONNECTIVITY_H_ */