diff --git a/utils/api_demos/mainRemesh.cpp b/utils/api_demos/mainRemesh.cpp
index 69713b774e20ef7f1524d3af011d46a70bed4ccc..ca0169226db46969479d8cb4d4f5baa6bd465b9d 100644
--- a/utils/api_demos/mainRemesh.cpp
+++ b/utils/api_demos/mainRemesh.cpp
@@ -17,8 +17,20 @@ int main(int argc, char **argv)
   GmshSetOption("Mesh", "Optimize", 0.); // not yet: need boundary!
 
   GModel *m = new GModel();
+
   m->readMSH("cube.msh");
 
+  // discreteRegion *gr = new discreteRegion(m);
+  // MVertex *v0 = new MVertex(x, y, z, gr, tag);
+  // MVertex *v1 = new MVertex(x, y, z, gr, tag);
+  // ...
+  // gr->mesh_vertices.push_back(v0);
+  // ...
+  // MTetrahedron *t = new MTetrahedron(v0, v1, v2, v3, tag);
+  // gr->tetrahedra.push_back(t);
+
+  //MVertex *v2 = gr->mesh_vertices[2];
+
   for(GModel::riter it = m->firstRegion(); it != m->lastRegion(); ++it){
     discreteRegion *r = dynamic_cast<discreteRegion*>(*it);
     if(r){