From 62a266f1040a01f42790c18533d2addc2417d86b Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Sat, 18 Mar 2017 08:55:00 +0000
Subject: [PATCH] save embedded entities in unrolled volumes
---
Geo/GRegion.cpp | 12 ++++++++++++
1 file changed, 12 insertions(+)
diff --git a/Geo/GRegion.cpp b/Geo/GRegion.cpp
index 24b2bc0bbc..bca79da37c 100644
--- a/Geo/GRegion.cpp
+++ b/Geo/GRegion.cpp
@@ -269,6 +269,18 @@ void GRegion::writeGEO(FILE *fp)
fprintf(fp, "Volume(%d) = {%d};\n", tag(), tag());
}
+ for(std::list<GFace*>::iterator it = embedded_faces.begin();
+ it != embedded_faces.end(); it++)
+ fprintf(fp, "Surface {%d} In Volume {%d};\n", (*it)->tag(), tag());
+
+ for(std::list<GEdge*>::iterator it = embedded_edges.begin();
+ it != embedded_edges.end(); it++)
+ fprintf(fp, "Line {%d} In Volume {%d};\n", (*it)->tag(), tag());
+
+ for(std::list<GVertex*>::iterator it = embedded_vertices.begin();
+ it != embedded_vertices.end(); it++)
+ fprintf(fp, "Point {%d} In Volume {%d};\n", (*it)->tag(), tag());
+
if(meshAttributes.method == MESH_TRANSFINITE){
fprintf(fp, "Transfinite Volume {%d}", tag());
if(meshAttributes.corners.size()){
--
GitLab