save elem=numnode to test point datasets

Geo/GModel.cpp

@@ -785,9 +785,10 @@ void GModel::_storeVerticesInEntities(std::map<int, MVertex*> &vertices)
     MVertex *v = it->second;
     GEntity *ge = v->onWhat();
     if(ge){
-      if(ge->dim() || ge->mesh_vertices.empty()) // special case for points
+      if(ge->dim() || ge->mesh_vertices.empty()){ // special case for points
 	ge->mesh_vertices.push_back(v);
       }
+    }
     else
       delete v; // we delete all unused vertices
   }
@@ -800,9 +801,10 @@ void GModel::_storeVerticesInEntities(std::vector<MVertex*> &vertices)
     if(v){ // the vector is allowed to have null entries
       GEntity *ge = v->onWhat();
       if(ge) {
-	if(ge->dim() || ge->mesh_vertices.empty()) // special case for points
+	if(ge->dim() || ge->mesh_vertices.empty()){ // special case for points
 	  ge->mesh_vertices.push_back(v);
 	}
+      }
       else
         delete v; // we delete all unused vertices
     }

Post/PViewDataListIO.cpp

@@ -442,7 +442,7 @@ static void writeElementMSH(FILE *fp, int num, int nbnod, pVertex nod[8],
 {
   switch(dim){
   case 0:
-    fprintf(fp, "%d 15 0 %d\n", num, nod[0].Num);
+    fprintf(fp, "%d 15 2 0 %d %d\n", num, num, nod[0].Num);
     break;
   case 1:
     fprintf(fp, "%d 1 0 %d %d\n", num, nod[0].Num, nod[1].Num);