def pointInPolygon3(p = "const Dim<3>::Vector&",
                    vertices = "const std::vector<Dim<3>::Vector>&",
                    ipoints = "const std::vector<unsigned>&",
                    normal = "const Dim<3>::Vector&",
                    countBoundary = ("const bool", "false"),
                    tol = ("const double", "1.0e-10")):
    """3D point in polygon test -- the point and all vertices of the polygon must be coplanar.
  p             : point we're testing
  vertices      : contains the vertices the polygon as a subset
  ipoints       : the subset of vertices that define the polygon
  normal        : normal to the plane of the polygon
  countBoundary : should points on the boundary count as internal?
  tol           : tolerance for on boundary"""
    return "bool"