1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143 | #include "MeshRefiner.hpp"
#include "EdgeSizeEvaluator.hpp"
#include "EntityRefiner.hpp"
#include "moab/Interface.hpp"
#include "RefinerTagManager.hpp"
#include "MeshOutputFunctor.hpp"
#ifdef MOAB_HAVE_MPI
#include "moab/ParallelComm.hpp"
#include "moab_mpi.h"
#else // MOAB_HAVE_MPI
typedef int MPI_Comm;
#endif // MOAB_HAVE_MPI
namespace moab
{
/**\brief Construct a mesh refiner.
* The input and output mesh pointers may be identical.
* Existing elements will <b>not</b> be removed from the input mesh
* as they are refined, so the adjacencies for entitities may appear
* strange after refinement.
*/
MeshRefiner::MeshRefiner( Interface* imesh, Interface* omesh )
{
this->mesh_in = imesh;
this->mesh_out = omesh;
this->tag_manager = new RefinerTagManager( this->mesh_in, this->mesh_out );<--- Class 'MeshRefiner' does not have a copy constructor which is recommended since it has dynamic memory/resource allocation(s).<--- Class 'MeshRefiner' does not have a operator= which is recommended since it has dynamic memory/resource allocation(s).
this->output_functor = new MeshOutputFunctor( this->tag_manager );
this->entity_refiner = 0;
this->comm = ParallelComm::get_pcomm( this->mesh_out, 0 );
}
/**\brief Destroy a mesh refiner.
*
* Note that any EntityRefiner object associated with the mesh refiner will be deleted inside this
* destructor. Destruction is virtual so subclasses may clean up after refinement.
*/
MeshRefiner::~MeshRefiner()
{
delete this->tag_manager;
delete this->output_functor;
if( this->entity_refiner ) delete this->entity_refiner;
}
/**\brief Specify which techniqe will be used to refine entities in the mesh.
* The entity refiner object you specify is ``owned'' by the mesh refiner after this call;
* the entity refiner will be deleted when this mesh refiner is destroyed.
* You should not delete the entity refiner yourself.
*/
bool MeshRefiner::set_entity_refiner( EntityRefiner* er )
{
if( !er || er == this->entity_refiner ) return false;
this->entity_refiner = er;
return true;
}
/**\brief A convenience method to reset the list of tags to be copied to output vertices.
* This simply calls the method of the same name on the tag manager.
*/
void MeshRefiner::reset_vertex_tags()
{
this->tag_manager->reset_vertex_tags();
}
/**\brief A convenience method to add a tag to be copied/interpolated from input vertices to output
* vertices. This simply calls the method of the same name on the tag manager.
*/
int MeshRefiner::add_vertex_tag( Tag tag_handle )
{
return this->tag_manager->add_vertex_tag( tag_handle );
}
struct MeshRefinerIterator
{
Range subset;
EntityHandle destination_set;
};
/**\brief Refine entities in a mesh set.
* This will recursively descend any mesh sets contained in the \a range.
* It returns false when not able to refine (because no entity refiner is
* set or no edge size evaluator has been set on the entity refiner) and
* true otherwise.
*/
bool MeshRefiner::refine( Range& range )
{
this->tag_manager->create_output_tags();
if( !this->entity_refiner->prepare( this->tag_manager, this->output_functor ) )
{ // Oops, perhaps the edge_size_evaluator was not set?
return false;
}
MeshRefinerIterator entry;
std::vector< MeshRefinerIterator > work;
entry.subset = range;
entry.destination_set = 0;
work.push_back( entry );
while( !work.empty() )
{
entry = work.back();
work.pop_back();
this->output_functor->destination_set = entry.destination_set;
for( Range::const_iterator it = entry.subset.begin(); it != entry.subset.end(); ++it )
{
EntityType etyp = this->mesh_in->type_from_handle( *it );
if( etyp == MBENTITYSET )
{
Range set_ents;
if( this->mesh_in->get_entities_by_handle( *it, set_ents, false ) == MB_SUCCESS )
{
// Create a matching set on the output mesh.
MeshRefinerIterator set_work;
unsigned int set_work_opts;
this->mesh_in->get_meshset_options( *it, set_work_opts );
this->mesh_out->create_meshset( set_work_opts, set_work.destination_set );
set_work.subset = set_ents;
work.push_back( set_work );
// Copy any per-element tag values the user has requested to the output set.
this->tag_manager->set_element_tags_from_ent( *it );
this->tag_manager->assign_element_tags( set_work.destination_set );
// Copy the global ID to the new set (assuming it exists).
this->tag_manager->copy_gid( *it, set_work.destination_set );
}
}
else
{
this->tag_manager->set_element_tags_from_ent( *it );
this->tag_manager->set_element_procs_from_ent( *it );
this->entity_refiner->refine_entity( etyp, *it );
}
}
}
this->output_functor->assign_global_ids( this->comm );
return true;
}
} // namespace moab
|