songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Upgrade example Customized_subdivision

This commit is contained in:
Andreas Fabri 2017-01-03 12:21:23 +01:00 committed by Mael Rouxel-Labbé
parent 20caf9a7c3
commit 9ecbd56009
1 changed files with 40 additions and 38 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

78
Subdivision_method_3/examples/Subdivision_method_3/Customized_subdivision.cpp
View File

@ -1,15 +1,15 @@
#include <CGAL/Simple_cartesian.h>
#include <CGAL/Subdivision_method_3.h>
#include <CGAL/Surface_mesh.h>
#include <CGAL/Polyhedron_3.h>
#include <boost/foreach.hpp>
#include <cstdio>
#include <iostream>
#include <CGAL/Polyhedron_3.h>
#include <CGAL/IO/Polyhedron_iostream.h>
typedef CGAL::Simple_cartesian<double> Kernel;
typedef CGAL::Polyhedron_3<Kernel> Polyhedron;
//typedef CGAL::Polyhedron_3<Kernel> Polyhedron;
typedef CGAL::Surface_mesh<Kernel::Point_3> Polyhedron;
using namespace std;
using namespace CGAL;
@ -18,43 +18,43 @@ template <class Poly>
class WLoop_mask_3 {
typedef Poly Polyhedron;
typedef typename Polyhedron::Vertex_iterator Vertex_iterator;
typedef typename Polyhedron::Halfedge_iterator Halfedge_iterator;
typedef typename Polyhedron::Facet_iterator Facet_iterator;
typedef typename boost::graph_traits<Polyhedron>::vertex_descriptor vertex_descriptor;
typedef typename boost::graph_traits<Polyhedron>::halfedge_descriptor halfedge_descriptor;
typedef typename Polyhedron::Halfedge_around_facet_circulator
Halfedge_around_facet_circulator;
typedef typename Polyhedron::Halfedge_around_vertex_circulator
Halfedge_around_vertex_circulator;
typedef Halfedge_around_target_circulator<Poly> Halfedge_around_target_circulator;
typedef typename Polyhedron::Traits Traits;
typedef typename Traits::Kernel Kernel;
typedef typename Kernel::FT FT;
typedef typename Kernel::Point_3 Point;
typedef typename Kernel::Vector_3 Vector;
typedef typename boost::property_map<Polyhedron, vertex_point_t>::type Vertex_pmap;
typedef typename boost::property_traits<Vertex_pmap>::value_type Point;
Polyhedron& polyhedron;
Vertex_pmap vpm;
public:
void edge_node(Halfedge_iterator eitr, Point& pt) {
Point& p1 = eitr->vertex()->point();
Point& p2 = eitr->opposite()->vertex()->point();
Point& f1 = eitr->next()->vertex()->point();
Point& f2 = eitr->opposite()->next()->vertex()->point();
WLoop_mask_3(Polyhedron& polyhedron)
: polyhedron(polyhedron), vpm(get(CGAL::vertex_point, polyhedron))
{}
void edge_node(halfedge_descriptor hd, Point& pt) {
Point& p1 = get(vpm, target(hd,polyhedron));
Point& p2 = get(vpm, target(opposite(hd,polyhedron),polyhedron));
Point& f1 = get(vpm, target(next(hd,polyhedron),polyhedron));
Point& f2 = get(vpm, target(next(opposite(hd,polyhedron),polyhedron),polyhedron));
pt = Point((3*(p1[0]+p2[0])+f1[0]+f2[0])/8,
(3*(p1[1]+p2[1])+f1[1]+f2[1])/8,
(3*(p1[2]+p2[2])+f1[2]+f2[2])/8 );
}
void vertex_node(Vertex_iterator vitr, Point& pt) {
void vertex_node(vertex_descriptor vd, Point& pt) {
double R[] = {0.0, 0.0, 0.0};
Point& S = vitr->point();
Point& S = get(vpm,vd);
Halfedge_around_vertex_circulator vcir = vitr->vertex_begin();
std::size_t n = circulator_size(vcir);
for (std::size_t i = 0; i < n; i++, ++vcir) {
Point& p = vcir->opposite()->vertex()->point();
std::size_t n = 0;
BOOST_FOREACH(halfedge_descriptor hd, halfedges_around_target(vd, polyhedron)){
++n;
Point& p = get(vpm, target(opposite(hd,polyhedron),polyhedron));
R[0] += p[0]; R[1] += p[1]; R[2] += p[2];
}
if (n == 6) {
pt = Point((10*S[0]+R[0])/16, (10*S[1]+R[1])/16, (10*S[2]+R[2])/16);
} else if (n == 3) {
@ -67,15 +67,17 @@ public:
pt = Point((A*S[0]+B*R[0]), (A*S[1]+B*R[1]), (A*S[2]+B*R[2]));
}
}
void border_node(Halfedge_iterator eitr, Point& ept, Point& vpt) {
Point& ep1 = eitr->vertex()->point();
Point& ep2 = eitr->opposite()->vertex()->point();
void border_node(halfedge_descriptor hd, Point& ept, Point& vpt) {
Point& ep1 = get(vpm, target(hd,polyhedron));
Point& ep2 = get(vpm, target(opposite(hd,polyhedron),polyhedron));
ept = Point((ep1[0]+ep2[0])/2, (ep1[1]+ep2[1])/2, (ep1[2]+ep2[2])/2);
Halfedge_around_vertex_circulator vcir = eitr->vertex_begin();
Point& vp1 = vcir->opposite()->vertex()->point();
Point& vp0 = vcir->vertex()->point();
Point& vp_1 = (--vcir)->opposite()->vertex()->point();
Halfedge_around_target_circulator vcir(hd,polyhedron);
Point& vp1 = get(vpm, target(opposite(*vcir,polyhedron),polyhedron));
Point& vp0 = get(vpm, target(*vcir,polyhedron));
--vcir;
Point& vp_1 = get(vpm,target(opposite(*vcir,polyhedron),polyhedron));
vpt = Point((vp_1[0] + 6*vp0[0] + vp1[0])/8,
(vp_1[1] + 6*vp0[1] + vp1[1])/8,
(vp_1[2] + 6*vp0[2] + vp1[2])/8 );
@ -95,7 +97,7 @@ int main(int argc, char **argv) {
Polyhedron P;
cin >> P; // read the .off
Subdivision_method_3::PTQ(P, WLoop_mask_3<Polyhedron>(), d);
Subdivision_method_3::PTQ(P, WLoop_mask_3<Polyhedron>(P), d);
cout << P; // write the .off