Function Geom_point_inside

Code/Source/PythonAPI/Geometry_PyModule.cxx:276–312

");

static PyObject *
Geom_point_inside(PyObject* self, PyObject* args, PyObject* kwargs)
{
  //std::cout << "========== Geom_point_inside ==========" << std::endl;
  auto api = PyUtilApiFunction("OO!", PyRunTimeErr, __func__);
  static char *keywords[] = {"polydata", "point", NULL};
  PyObject* pdObj;
  PyObject* pointArg;

  if (!PyArg_ParseTupleAndKeywords(args, kwargs, api.format, keywords, &pdObj, &PyList_Type, &pointArg)) {
      return api.argsError();
  }

  // Get the vtkPolyData object from the Python object.
  auto polydata = GetVtkPolyData(api, pdObj);
  if (polydata == nullptr) {
      return nullptr;
  }

  // Get the point data.
  double point[3];
  std::string emsg;
  if (!PyUtilGetPointData<double>(pointArg, emsg, point)) {
      api.error("The point argument " + emsg);
      return nullptr;
  }

  // Result is -1 if the point is outside, 1 if it is inside.
  cvPolyData cvPolydata(polydata);
  int result;
  if (sys_geom_Classify(&cvPolydata, point, &result) != SV_OK) {
      api.error("Error classifying a point for the geometry.");
      return nullptr;
  }

  return Py_BuildValue("N", PyBool_FromLong(result+1));
}

//-------------------------------
// Geom_interpolate_closed_curve