---

BRMeshRefine.H

Go to the documentation of this file.

/* _______              __
  / ___/ /  ___  __ _  / /  ___
 / /__/ _ \/ _ \/  ' \/ _ \/ _ \
 \___/_//_/\___/_/_/_/_.__/\___/ 
*/
//
// This software is copyright (C) by the Lawrence Berkeley
// National Laboratory.  Permission is granted to reproduce
// this software for non-commercial purposes provided that
// this notice is left intact.
// 
// It is acknowledged that the U.S. Government has rights to
// this software under Contract DE-AC03-765F00098 between
// the U.S.  Department of Energy and the University of
// California.
//
// This software is provided as a professional and academic
// contribution for joint exchange. Thus it is experimental,
// is provided ``as is'', with no warranties of any kind
// whatsoever, no support, no promise of updates, or printed
// documentation. By using this software, you acknowledge
// that the Lawrence Berkeley National Laboratory and
// Regents of the University of California shall have no
// liability with respect to the infringement of other
// copyrights by any part of this software.

#ifndef _BR_MESH_REFINE_H_
#define _BR_MESH_REFINE_H_
// Include files:

#include <climits>

#include "Vector.H"
#include "Box.H"
#include "IntVectSet.H"
#include "REAL.H"
#include "MeshRefine.H"
#include <list>

// Constants:

// Minimum acceptable ratio of tagged cells to total cells for 
// the Berger-Rigoutsos algorithm in \func{makeBoxes}.
// Used as default for \var{FillRatio} optional argument.
#ifndef _BR_MIN_BOX_FILL_RATIO_
#define _BR_MIN_BOX_FILL_RATIO_ ( 0.75 )
#endif


class BRMeshRefine : public MeshRefine
{
public:
  BRMeshRefine();
  
  BRMeshRefine(
               const Box& a_baseDomain,
               const Vector<int>& a_refRatios,
               const Real a_fillRatio,
               const int a_blockFactor,
               const int a_bufferSize,
               const int a_maxSize);

  BRMeshRefine(
               const ProblemDomain& a_baseDomain,
               const Vector<int>& a_refRatios,
               const Real a_fillRatio,
               const int a_blockFactor,
               const int a_bufferSize,
               const int a_maxSize);
  

  virtual ~BRMeshRefine();
  
  void define(
              const Box& a_baseDomain,
              const Vector<int>& a_refRatios,
              const Real a_fillRatio,
              const int a_blockFactor,
              const int a_bufferSize,
              const int a_maxSize);


  void define(
              const ProblemDomain& a_baseDomain,
              const Vector<int>& a_refRatios,
              const Real a_fillRatio,
              const int a_blockFactor,
              const int a_bufferSize,
              const int a_maxSize);
  


  void 
  makeBoxes(
            Vector<Box>&      a_mesh,        
            const IntVectSet& a_tags,        
            const IntVectSet& a_pnd,         
            const ProblemDomain& a_domain,
            const int         a_maxSize) const;
  

protected:
  void 
  makeBoxes(
            std::list<Box>&      a_mesh,        
            const IntVectSet& a_tags,        
            const IntVectSet& a_pnd,         
            const ProblemDomain& a_domain,
            const int         a_maxSize,
                        const int         a_depth
                        ) const;
  void
  splitBox(std::list<Box> & a_boxes ,
                   const std::list<Box>::iterator& a_boxindex,
           const int a_dimension ,const int a_maxboxsize ) const;
  
   Vector<int>
  makeTrace( const IntVectSet& a_Ivs ,int a_dir ) const;
  
  void
  makeTraces( const IntVectSet& a_Ivs ,Vector<int>* a_traces ) const;
   int
  findSplit( const Vector<int>& a_trace ) const;
  
   int
  findMaxInflectionPoint( const Vector<int>& a_trace ,int& a_maxVal ) const;
  
   void
  splitTags( const IntVectSet& a_tags,
             const int a_split_dir ,const int a_split_indx,
             IntVectSet& a_tags_lo ,IntVectSet& a_tags_hi ) const;
  
  void 
  breakBoxes(Vector<Box>& a_vboxin,  const int& a_maxSize, 
             const int& a_idir) const;
  
   int
  maxloc( const int* a_V ,const int a_Size ) const;


};


extern void domainSplit(const ProblemDomain& a_domain, Vector<Box>& a_vbox, 
                        int a_maxSize, int a_blockfactor=1);


extern void domainSplit(const Box& a_domain, Vector<Box>& a_vbox, 
                        int a_maxSize, int a_blockfactor=1);

#endif

---