Function MeshUtils_remesh

Code/Source/PythonAPI/MeshUtils_PyModule.cxx:119–185

------------ MMG_remesh ------------

//------------
//
static PyObject *
MeshUtils_remesh(PyObject* self, PyObject* args, PyObject* kwargs)
{
  auto api = PyUtilApiFunction("O|ddddds", PyRunTimeErr, __func__);
  static char *keywords[] = {"surface", "hmin", "hmax", "angle", "hgrad", "hausd", "log_file", NULL};
  PyObject* surfaceArg;
  double hmin = 0.1;
  double hmax = 0.1;
  double angle = 45.0;
  double hgrad = 1.1;
  double hausd = 0.01;
  char *logFileName = nullptr;

  if (!PyArg_ParseTupleAndKeywords(args, kwargs, api.format, keywords, &surfaceArg, &hmin, &hmax, &angle, 
          &hgrad, &hausd, &logFileName)) {
      return api.argsError();
  }

  // Redirect stdout to the 'mesh.log' file.
  //
  // [TODO:DaveP] this does not work for mmg or doesn't work on ubunutu?
  //
  // [TODO:DaveP] put this in py utils as a class.
  //
/*
  if (logFileName == nullptr) {
      logFileName = "/dev/null";
  }
  int stdout_dupe = dup(fileno(stdout));
  freopen(logFileName, "w", stdout);
*/

  // Get the vtkPolyData objectsfrom the Python object.
  //
  auto surfPolydata = GetVtkPolyData(api, surfaceArg);
  if (surfPolydata == nullptr) {
      return nullptr;
  }
  //cvPolyData cvSurfPolydata(surfPolydata);
  surfPolydata->BuildLinks();

  // Remesh the entire surface polydata.
  //
  int useSizingFunction = 0;
  int numAddedRefines = 0;
  vtkDoubleArray *meshSizingFunction = NULL;

  if (MMGUtils_SurfaceRemeshing(surfPolydata, hmin, hmax, hausd, angle, hgrad, useSizingFunction, meshSizingFunction,
          numAddedRefines) != SV_OK) {
      api.error("Remeshing failed.");
      return nullptr;
  }

  auto result = new cvPolyData(surfPolydata);
  if (result == NULL) {
      api.error("Error creating polydata from the remeshed object.");
      return nullptr;
  }