我如何写一个简单的CGAL代码来对一个二维矩形进行网格划分？

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Delaunay_triangulation_2.h>
#include <CGAL/Triangulation_vertex_base_with_info_2.h>

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Triangulation_vertex_base_with_info_2<int, K> Vb;
typedef CGAL::Triangulation_data_structure_2<Vb> Tds;
typedef CGAL::Delaunay_triangulation_2<K, Tds> Delaunay;

int main()
{
  // 定义矩形的四个顶点
  std::vector<K::Point_2> points = {
    K::Point_2(0, 0),
    K::Point_2(1, 0),
    K::Point_2(1, 1),
    K::Point_2(0, 1)
  };

  // 创建Delaunay三角剖分对象
  Delaunay dt;
  dt.insert(points.begin(), points.end());

  // 遍历三角剖分的三角形
  for (auto it = dt.finite_faces_begin(); it != dt.finite_faces_end(); ++it)
  {
    // 获取三角形的三个顶点
    K::Point_2 p1 = it->vertex(0)->point();
    K::Point_2 p2 = it->vertex(1)->point();
    K::Point_2 p3 = it->vertex(2)->point();

    // 进行网格划分的操作，例如输出三角形的顶点坐标
    std::cout << "Triangle: (" << p1 << "), (" << p2 << "), (" << p3 << ")" << std::endl;
  }

  return 0;
}