add createGeometry to discreteRegion, so that we can selectively choose which...

add createGeometry to discreteRegion, so that we can selectively choose which discrete volumes we want to remesh

Geo/discreteRegion.cpp

@@ -7,12 +7,27 @@
 #include "GModelIO_GEO.h"
 #include "Geo.h"
 
-discreteRegion::discreteRegion(GModel *model, int num) : GRegion(model, num)
+discreteRegion::discreteRegion(GModel *model, int num)
+  : GRegion(model, num), _geometry(false)
 {
   ::Volume *v = CreateVolume(num, MSH_VOLUME_DISCRETE);
   Tree_Add(model->getGEOInternals()->Volumes, &v);
 }
 
-discreteRegion::discreteRegion(GModel *model) : GRegion(model, 0)
+discreteRegion::discreteRegion(GModel *model)
+  : GRegion(model, 0), _geometry(false)
 {
 }
+
+int discreteRegion::createGeometry()
+{
+  if(l_faces.empty()) {
+    Msg::Error("Discrete volume %d has no surfaces on its boundary", tag());
+    return 1;
+  }
+
+  // we should check that the surfaces form one or more (closed) shell(s)
+  _geometry = true;
+
+  return 0;
+}

Geo/discreteRegion.h

@@ -11,11 +11,15 @@
 #include "MFace.h"
 
 class discreteRegion : public GRegion {
+private:
+  bool _geometry;
 public:
   discreteRegion(GModel *model, int num);
   discreteRegion(GModel *model);
   virtual ~discreteRegion() {}
   virtual GeomType geomType() const { return DiscreteVolume; }
+  virtual bool haveParametrization() { return _geometry; }
+  int createGeometry();
 };
 
 #endif