refine.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_mesh_simplicial_refine_h__)
#define __geom_mesh_simplicial_refine_h__

#include "coarsen.h"
#include "inc_opt.h"
#include "insert.h"
#include "set.h"

#include "../iss/ISS_SignedDistance.h"
#include "../iss/PointsOnManifold.h"

#include "../../../ads/functor/constant.h"

#include <set>

// CONTINUE: I use the following for debugging.
#if 0
#include "file_io.h"
#include <sstream>
#include <iomanip>
#endif

BEGIN_NAMESPACE_GEOM

// CONTINUE: Refining is really slow when the max edge length is relatively
// small.  This is due to deep recursion.  Fix this problem.


template<int N, int M, typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         int MM, int SD,
         class MaxEdgeLength>
int
refine(SimpMeshRed<N,M,T,Node,Cell,Cont>* mesh, 
       PointsOnManifold<N,MM,SD,T>* manifold,
       const MaxEdgeLength& f);


  

template<int N, int M, typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         class MaxEdgeLength>
inline
int
refine(SimpMeshRed<N,M,T,Node,Cell,Cont>* mesh, const MaxEdgeLength& f) {
  PointsOnManifold<N,M-1,1,T>* null = 0;
  // Call the above function with a null boundary manifold pointer.
  return refine(mesh, null, f);
}



  

template<int N, int M, typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         int MM, int SD,
         typename IntInputIterator>
int
refine(SimpMeshRed<N,M,T,Node,Cell,Cont>* mesh, 
       PointsOnManifold<N,MM,SD,T>* manifold,
       IntInputIterator begin, IntInputIterator end);




template<int N, int M, typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         typename IntInputIterator>
inline
int
refine(SimpMeshRed<N,M,T,Node,Cell,Cont>* mesh, 
       IntInputIterator begin, IntInputIterator end) {
  PointsOnManifold<N,M-1,1,T>* null = 0;
  // Call the above function with a null boundary manifold pointer.
  return refine(mesh, null, begin, end);
}




template<int N, int M, typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         int MM, int SD,
         typename IntInputIterator,
         class MaxEdgeLength>
int
refine(SimpMeshRed<N,M,T,Node,Cell,Cont>* mesh, 
       PointsOnManifold<N,MM,SD,T>* manifold,
       IntInputIterator begin, IntInputIterator end,
       const MaxEdgeLength& f);





template<int N, int M, typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         typename IntInputIterator,
         class MaxEdgeLength>
inline
int
refine(SimpMeshRed<N,M,T,Node,Cell,Cont>* mesh, 
       IntInputIterator begin, IntInputIterator end,
       const MaxEdgeLength& f) {
  PointsOnManifold<N,M-1,1,T>* null = 0;
  // Call the above function with a null boundary manifold pointer.
  return refine(mesh, null, begin, end, f);
}




  

template<typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         class ISS>
int
refineBoundary(SimpMeshRed<2,2,T,Node,Cell,Cont>* x, 
               const ISS& boundary,
               T maxAngle, 
               T minEdgeLength,
               int maxSweeps = 10);
  


template<typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         class Distance,
         class ClosestPoint,
         class MaxAngle,
         class MinEdgeLength>
int
refineBoundary(SimpMeshRed<2,2,T,Node,Cell,Cont>* x, 
               const Distance& distance,
               const ClosestPoint& closestPoint,
               const MaxAngle& maxAngle,
               const MinEdgeLength& minEdgeLength,
               int maxSweeps = 10);
  


template<typename T,
         template<class> class Node,
         template<class> class Cell,
         template<class,class> class Cont,
         class MaxEdgeLength>
int
refineAdjust(SimpMeshRed<2,2,T,Node,Cell,Cont>* mesh, const MaxEdgeLength& f);
  

END_NAMESPACE_GEOM

#define __geom_mesh_simplicial_refine2_ipp__
#include "refine2.ipp"
#undef __geom_mesh_simplicial_refine2_ipp__

#define __geom_mesh_simplicial_refine3_ipp__
#include "refine3.ipp"
#undef __geom_mesh_simplicial_refine3_ipp__

#define __geom_mesh_simplicial_refine_ipp__
#include "refine.ipp"
#undef __geom_mesh_simplicial_refine_ipp__

#endif

---