Operations on Surface_mesh in CGAL

http://doc.cgal.org/latest/BGL/group__PkgBGLIterators.html#ga7746cb7b28520877375dfff6ec7bd32c

#include <vector>

#include <boost/foreach.hpp>

#include <CGAL/Simple_cartesian.h>
#include <CGAL/Surface_mesh.h>

typedef CGAL::Simple_cartesian<double> K;
typedef CGAL::Surface_mesh<K::Point_3> Mesh;
typedef Mesh::Vertex_index vertex_descriptor;
typedef Mesh::Face_index face_descriptor;
int main()
{
  Mesh m;

  // u            x
  // +------------+
  // |            |
  // |            |
  // |      f     |
  // |            |
  // |            |
  // +------------+
  // v            w

  // Add the points as vertices
  vertex_descriptor u = m.add_vertex(K::Point_3(0,1,0));
  vertex_descriptor v = m.add_vertex(K::Point_3(0,0,0));
  vertex_descriptor w = m.add_vertex(K::Point_3(1,0,0));
  vertex_descriptor x = m.add_vertex(K::Point_3(1,1,0));

  face_descriptor f = m.add_face(u,v,w,x);
 
  {
    std::cout << "vertices around vertex " << v << std::endl;
    CGAL::Vertex_around_target_circulator<Mesh> vbegin(m.halfedge(v),m), done(vbegin);

    do {
      std::cout << *vbegin++ << std::endl;
    } while(vbegin != done);
  }
  {
	std::cout << "edges around face " << f << std::endl;
    CGAL::Halfedge_around_face_iterator<Mesh> hbegin, hend;
    for(boost::tie(hbegin, hend) = halfedges_around_face(m.halfedge(f), m);
        hbegin != hend; 
        ++hbegin){
      std::cout << *hbegin << std::endl;
	  std::cout << m.face(*hbegin) << std::endl;
	  std::cout << m.opposite(*hbegin) << std::endl;
	  std::cout << m.face(m.opposite(*hbegin)) << std::endl;
    }
  }
  { 
    std::cout << "vertices around face " << f << std::endl;
    CGAL::Vertex_around_face_iterator<Mesh> vbegin, vend;
    for(boost::tie(vbegin, vend) = vertices_around_face(m.halfedge(f), m);
        vbegin != vend; 
        ++vbegin){
      std::cout << *vbegin << std::endl;
    }
  }

  // or the same again, but directly with a range based loop
  BOOST_FOREACH(vertex_descriptor vd,vertices_around_face(m.halfedge(f), m)){
    std::cout << vd << std::endl;
  } 
    

  return 0;
}

note the difference between Surface_mesh and Polyhedron when referring to data from descriptor

In Polyhedron

Halfedge_handle h ;

Halfedge_handle g = h->next()->opposite()->next();

In Surface_mesh

CGAL::Halfedge_around_face_iterator<Mesh> hbegin;

m.face(m.opposite(*hbegin))