52 #ifndef OPENMESH_SUBDIVIDER_UNIFORM_LOOPT_HH
53 #define OPENMESH_SUBDIVIDER_UNIFORM_LOOPT_HH
58 #include <OpenMesh/Core/System/config.hh>
60 #include <OpenMesh/Core/Utils/vector_cast.hh>
63 #if defined(OM_CC_MIPS)
73 namespace Subdivider {
87 template <
typename MeshType,
typename RealType =
float>
92 typedef RealType real_t;
93 typedef MeshType mesh_t;
96 typedef std::pair< real_t, real_t > weight_t;
97 typedef std::vector< std::pair<real_t,real_t> > weights_t;
102 LoopT(
void) :
parent_t(), _1over8( 1.0/8.0 ), _3over8( 3.0/8.0 )
106 LoopT( mesh_t& _m ) :
parent_t(_m), _1over8( 1.0/8.0 ), _3over8( 3.0/8.0 )
116 const char *
name()
const {
return "Uniform Loop"; }
122 weights_.resize(_max_valence);
123 std::generate(weights_.begin(), weights_.end(), compute_weight());
132 _m.add_property( vp_pos_ );
133 _m.add_property( ep_pos_ );
140 _m.remove_property( vp_pos_ );
141 _m.remove_property( ep_pos_ );
146 bool subdivide( mesh_t& _m,
size_t _n,
const bool _update_points =
true)
151 typename mesh_t::FaceIter fit, f_end;
152 typename mesh_t::EdgeIter eit, e_end;
153 typename mesh_t::VertexIter vit;
156 for (
size_t i=0; i < _n; ++i)
161 for (vit = _m.vertices_begin(); vit != _m.vertices_end(); ++vit) {
162 smooth(_m, vit.handle());
167 for (eit=_m.edges_begin(); eit != _m.edges_end(); ++eit)
168 compute_midpoint( _m, eit.handle() );
174 e_end = _m.edges_end();
175 for (eit=_m.edges_begin(); eit != e_end; ++eit)
176 split_edge(_m, eit.handle() );
182 f_end = _m.faces_end();
183 for (fit = _m.faces_begin(); fit != f_end; ++fit)
184 split_face(_m, fit.handle() );
188 for ( vit = _m.vertices_begin();
189 vit != _m.vertices_end(); ++vit) {
190 _m.set_point(vit, _m.property( vp_pos_, vit ) );
195 #if defined(_DEBUG) || defined(DEBUG)
208 struct compute_weight
210 compute_weight() : valence(-1) { }
211 weight_t operator() (
void)
213 #if !defined(OM_CC_MIPS)
222 double inv_v = 1.0/double(valence);
223 double t = (3.0 + 2.0 * cos( 2.0 * M_PI * inv_v) );
224 double alpha = (40.0 - t * t)/64.0;
226 return weight_t( 1.0-alpha, inv_v*alpha);
228 return weight_t(0.0, 0.0);
235 void split_face(
mesh_t& _m,
const typename mesh_t::FaceHandle& _fh)
237 typename mesh_t::HalfedgeHandle
238 heh1(_m.halfedge_handle(_fh)),
239 heh2(_m.next_halfedge_handle(_m.next_halfedge_handle(heh1))),
240 heh3(_m.next_halfedge_handle(_m.next_halfedge_handle(heh2)));
243 corner_cutting( _m, heh1 );
244 corner_cutting( _m, heh2 );
245 corner_cutting( _m, heh3 );
249 void corner_cutting(
mesh_t& _m,
const typename mesh_t::HalfedgeHandle& _he)
252 typename mesh_t::HalfedgeHandle
255 heh6(_m.next_halfedge_handle(heh1));
258 for (; _m.next_halfedge_handle(_m.next_halfedge_handle(heh5)) != heh1;
259 heh5 = _m.next_halfedge_handle(heh5))
263 vh1 = _m.to_vertex_handle(heh1),
264 vh2 = _m.to_vertex_handle(heh5);
266 typename mesh_t::HalfedgeHandle
267 heh2(_m.next_halfedge_handle(heh5)),
268 heh3(_m.new_edge( vh1, vh2)),
269 heh4(_m.opposite_halfedge_handle(heh3));
285 typename mesh_t::FaceHandle fh_old(_m.face_handle(heh6));
286 typename mesh_t::FaceHandle fh_new(_m.new_face());
290 _m.set_next_halfedge_handle(heh4, heh6);
291 _m.set_next_halfedge_handle(heh5, heh4);
293 _m.set_face_handle(heh4, fh_old);
294 _m.set_face_handle(heh5, fh_old);
295 _m.set_face_handle(heh6, fh_old);
296 _m.set_halfedge_handle(fh_old, heh4);
299 _m.set_next_halfedge_handle(heh1, heh3);
300 _m.set_next_halfedge_handle(heh3, heh2);
302 _m.set_face_handle(heh1, fh_new);
303 _m.set_face_handle(heh2, fh_new);
304 _m.set_face_handle(heh3, fh_new);
306 _m.set_halfedge_handle(fh_new, heh1);
310 void split_edge(
mesh_t& _m,
const typename mesh_t::EdgeHandle& _eh)
312 typename mesh_t::HalfedgeHandle
313 heh = _m.halfedge_handle(_eh, 0),
314 opp_heh = _m.halfedge_handle(_eh, 1);
316 typename mesh_t::HalfedgeHandle new_heh, opp_new_heh, t_heh;
319 typename mesh_t::Point midP(_m.point(_m.to_vertex_handle(heh)));
320 midP += _m.point(_m.to_vertex_handle(opp_heh));
327 _m.property( vp_pos_, vh ) = _m.property( ep_pos_, _eh );
331 if (_m.is_boundary(_eh))
334 _m.next_halfedge_handle(t_heh) != opp_heh;
335 t_heh = _m.opposite_halfedge_handle(_m.next_halfedge_handle(t_heh)))
340 for (t_heh = _m.next_halfedge_handle(opp_heh);
341 _m.next_halfedge_handle(t_heh) != opp_heh;
342 t_heh = _m.next_halfedge_handle(t_heh) )
346 new_heh = _m.new_edge(vh, vh1);
347 opp_new_heh = _m.opposite_halfedge_handle(new_heh);
348 _m.set_vertex_handle( heh, vh );
350 _m.set_next_halfedge_handle(t_heh, opp_new_heh);
351 _m.set_next_halfedge_handle(new_heh, _m.next_halfedge_handle(heh));
352 _m.set_next_halfedge_handle(heh, new_heh);
353 _m.set_next_halfedge_handle(opp_new_heh, opp_heh);
355 if (_m.face_handle(opp_heh).is_valid())
357 _m.set_face_handle(opp_new_heh, _m.face_handle(opp_heh));
358 _m.set_halfedge_handle(_m.face_handle(opp_new_heh), opp_new_heh);
361 _m.set_face_handle( new_heh, _m.face_handle(heh) );
362 _m.set_halfedge_handle( vh, new_heh);
363 _m.set_halfedge_handle( _m.face_handle(heh), heh );
364 _m.set_halfedge_handle( vh1, opp_new_heh );
367 _m.adjust_outgoing_halfedge( vh );
368 _m.adjust_outgoing_halfedge( vh1 );
373 void compute_midpoint(
mesh_t& _m,
const typename mesh_t::EdgeHandle& _eh)
375 #define V( X ) vector_cast< typename mesh_t::Normal >( X )
376 typename mesh_t::HalfedgeHandle heh, opp_heh;
378 heh = _m.halfedge_handle( _eh, 0);
379 opp_heh = _m.halfedge_handle( _eh, 1);
382 pos(_m.point(_m.to_vertex_handle(heh)));
384 pos += V( _m.point(_m.to_vertex_handle(opp_heh)) );
387 if (_m.is_boundary(_eh) )
394 pos += V(_m.point(_m.to_vertex_handle(_m.next_halfedge_handle(heh))));
395 pos += V(_m.point(_m.to_vertex_handle(_m.next_halfedge_handle(opp_heh))));
398 _m.property( ep_pos_, _eh ) = pos;
406 if (_m.is_boundary(_vh) )
408 typename mesh_t::HalfedgeHandle heh, prev_heh;
409 heh = _m.halfedge_handle( _vh );
411 if ( heh.is_valid() )
413 assert( _m.is_boundary( _m.edge_handle( heh ) ) );
415 prev_heh = _m.prev_halfedge_handle( heh );
418 to_vh = _m.to_vertex_handle( heh ),
419 from_vh = _m.from_vertex_handle( prev_heh );
422 pos = _m.point( _vh );
439 for (vvit=_m.vv_iter(_vh); vvit; ++vvit) {
443 pos *= weights_[valence].second;
444 pos += weights_[valence].first
448 _m.property( vp_pos_, _vh ) = pos;
458 const real_t _1over8;
459 const real_t _3over8;
469 #endif // OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH defined