Catmull-Clark細分割法(CGALライブラリ)

  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

 
 
 
 
-
!
 
 
 
 
 
-
-
|
|
|
|
!
|
|
|
|
|
|
|
|
-
!
|
|
!
 

#include <CGAL/Polyhedron_3.h>
#include <CGAL/Simple_cartesian.h>
#include <CGAL/Subdivision_method_3.h>
 
// IO
#include <CGAL/IO/Polyhedron_iostream.h>
 
typedef CGAL::Simple_cartesian<double> Kernel;
typedef CGAL::Polyhedron_3<Kernel>     Polyhedron;
 
int main(int argc, char* argv[])
{
  if (argc != 2) {
    std::cerr << "Usage: CatmullClark_subdivision d < filename" << std::endl;
    std::cerr << "       d: the depth of the subdivision (0 < d < 10)" << std::endl;
    std::cerr << "       filename: the input mesh (.off)" << std::endl;
    return EXIT_FAILURE;
  }
 
  int d = argv[1][0] - '0';
 
  Polyhedron P;
  std::cin >> P; // read the .off
 
  CGAL::Subdivision_method_3::CatmullClark_subdivision(P,d);
 
  // Output
  std::cout << P;
 
  return EXIT_SUCCESS;
}