OctTreeOps.hpp

/*------------------------------------------------------------------------*/
/*      phdMesh : Parallel Heterogneous Dynamic unstructured Mesh         */
/*                Copyright (2007) Sandia Corporation                     */
/*                                                                        */
/*  Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive   */
/*  license for use of this work by or on behalf of the U.S. Government.  */
/*                                                                        */
/*  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, or (at your option) any later version.                       */
/*                                                                        */
/*  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                                                                   */
/*------------------------------------------------------------------------*/
#ifndef pddgeom_OctTreeOps_hpp
#define pddgeom_OctTreeOps_hpp

#include <utility>
#include <vector>
#include <util/Parallel.hpp>
#include <util/OctTree.hpp>

namespace phdmesh {

//----------------------------------------------------------------------
struct IdentProcBox : public IdentProc {
  unsigned part ;
  float box[6] ;

  IdentProcBox() {}

  IdentProcBox( const IdentProcBox & rhs )
    : IdentProc( rhs ), part( rhs.part ) { Copy<6>( box , rhs.box ); }

  IdentProcBox & operator = ( const IdentProcBox & rhs )
    {
      IdentProc::operator =( rhs );
      part = rhs.part ;
      Copy<6>( box , rhs.box );
      return *this ;
    }
};

//----------------------------------------------------------------------
void box_global_bounds(
  ParallelMachine            arg_comm ,
  const unsigned             arg_domain_boxes_number ,
  const IdentProcBox * const arg_domain_boxes ,
  const unsigned             arg_range_boxes_number ,
  const IdentProcBox * const arg_range_boxes ,
        float        * const arg_global_box );

//----------------------------------------------------------------------
bool oct_tree_proximity_search(
  ParallelMachine            arg_comm ,
  const float        * const arg_global_box ,
  const unsigned             arg_domain_boxes_number ,
  const IdentProcBox * const arg_domain_boxes ,
  const unsigned             arg_range_boxes_number ,
  const IdentProcBox * const arg_range_boxes ,
  const OctTreeKey   * const arg_cuts ,
  std::vector< std::pair<IdentProc,IdentProc> > & arg_relation ,
  unsigned * const arg_search_tree_stats = NULL );

//----------------------------------------------------------------------
bool hsfc_box_covering( const float * const global_box ,
                        const float * const small_box ,
                        OctTreeKey  * const covering ,
                        unsigned &          number );

//----------------------------------------------------------------------
void oct_tree_partition_course(
  ParallelMachine          comm ,
  const unsigned           length ,
  const OctTreeKey * const keys ,
  const float      * const weights ,
        OctTreeKey * const cut_keys ,
  const unsigned           override_ncuts = 0 );

//----------------------------------------------------------------------
void oct_tree_partition_fine(
  ParallelMachine          comm ,
  const unsigned           length ,
  const OctTreeKey * const keys ,
  const float      * const weights ,
        OctTreeKey * const cut_keys );

//----------------------------------------------------------------------
void oct_tree_partition_private(
  const unsigned p_first ,
  const unsigned p_end ,
  const unsigned depth ,
  const double   tolerance ,
  float * const weights ,
  const unsigned cuts_length ,
  OctTreeKey * const cuts );

//----------------------------------------------------------------------

}

#endif