From 80371071aae12bfd19631f2352e65eda35ae0d37 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Tue, 11 Apr 2006 12:19:02 +0000
Subject: [PATCH] *** empty log message ***
---
Common/DefaultOptions.h | 7 +-
utils/converters/scilab/load_gmsh.sci | 121 ++++++++++++++++++++++++++
2 files changed, 122 insertions(+), 6 deletions(-)
create mode 100644 utils/converters/scilab/load_gmsh.sci
diff --git a/Common/DefaultOptions.h b/Common/DefaultOptions.h
index 5fb5ec2f4f..5f63875377 100644
--- a/Common/DefaultOptions.h
+++ b/Common/DefaultOptions.h
@@ -672,12 +672,7 @@ StringXNumber GeneralOptions_Number[] = {
"Display size of points (in pixels)" },
{ F|O, "PolygonOffsetAlwaysOn" , opt_general_polygon_offset_always , 0. ,
"Always apply polygon offset, instead of trying to detect when it is required" },
- { F|O, "PolygonOffsetFactor" , opt_general_polygon_offset_factor ,
-#if defined(__APPLE__)
- 0.5 ,
-#else
- 1.0 ,
-#endif
+ { F|O, "PolygonOffsetFactor" , opt_general_polygon_offset_factor , 1.0 ,
"Polygon offset factor (offset = factor * DZ + r * units)" },
{ F|O, "PolygonOffsetUnits" , opt_general_polygon_offset_units , 1. ,
"Polygon offset units (offset = factor * DZ + r * units)" },
diff --git a/utils/converters/scilab/load_gmsh.sci b/utils/converters/scilab/load_gmsh.sci
new file mode 100644
index 0000000000..3ff8d0cd54
--- /dev/null
+++ b/utils/converters/scilab/load_gmsh.sci
@@ -0,0 +1,121 @@
+function gmesh = load_gmsh ( filename )
+ gmesh = [];
+ gmesh.MIN = zeros(3,1);
+ gmesh.MAX = zeros(3,1);
+ fid = mopen ( filename, 'r');
+
+ while 1
+ endoffile = 0;
+
+// while (1)
+// tline = mgetl(fid,1);
+// if meof(fid), endoffile=1,break, end
+// if (tline == '$'), break, end
+// end
+
+ tline = mgetl(fid,1);
+ if meof(fid), endoffile=1,break, end
+ if (tline == '$'), break, end
+
+ if (endoffile == 1), break, end
+
+ if tline == '$NOD'
+ disp('reading nodes')
+ gmesh.nbNod = mfscanf (1, fid, '%d');
+ gmesh.POS = zeros (gmesh.nbNod, 3);
+ for(I=1:gmesh.nbNod)
+ iNod = mfscanf(1,fid,'%d');
+ X = mfscanf(3,fid,'%g');
+ IDS(iNod) = I;
+ if (I == 1)
+ gmesh.MIN = X;
+ gmesh.MAX = X;
+ else
+ if (gmesh.MAX(1) < X(1)), gmesh.MAX(1)=X(1);end
+ if (gmesh.MAX(2) < X(2)), gmesh.MAX(2)=X(2);end
+ if (gmesh.MAX(3) < X(3)), gmesh.MAX(3)=X(3);end
+ if (gmesh.MIN(1) > X(1)), gmesh.MIN(1)=X(1);end
+ if (gmesh.MIN(2) > X(2)), gmesh.MIN(2)=X(2);end
+ if (gmesh.MIN(3) > X(3)), gmesh.MIN(3)=X(3);end
+ end
+ gmesh.POS(I,1) = X(1);
+ gmesh.POS(I,2) = X(2);
+ gmesh.POS(I,3) = X(3);
+ end
+ // Here we have to read two lines, instead of one, in order to continue.
+ // The reason may be due to an unread end-of-line character from the previous
+ // scan.
+ tline = mgetl(fid,2); // read 2 dummy lines
+ disp('nodes have been read')
+
+ else
+ if (tline == '$ELM')
+ disp('reading elements')
+ gmesh.nbElm = mfscanf (1,fid, '%d');
+ gmesh.ELE_INFOS = zeros (gmesh.nbElm,5);
+ gmesh.nbLines = 0;
+ gmesh.nbPoints = 0;
+ gmesh.nbTriangles = 0;
+ gmesh.nbQuads = 0;
+ gmesh.nbPrism=0;
+ gmesh.nbPyr=0;
+ gmesh.nbTet = 0;
+ gmesh.nbHex = 0;
+ gmesh.nbQTriangles = 0;
+ gmesh.nbQQuads = 0;
+ gmesh.POINTS=[];
+ gmesh.LINES=[];
+ gmesh.TRIANGLES=[];
+ gmesh.QUADS=[];
+ gmesh.TETS=[];
+ gmesh.HEXA=[];
+ gmesh.PRISM=[];
+ gmesh.PYRAMID=[];
+ gmesh.QTRIANGLES=[];
+ gmesh.QQUADS=[];
+ for (I=1:gmesh.nbElm)
+ gmesh.ELE_INFOS(I,:) = mfscanf(5,fid,'%d')';
+ NODES_ELEM = mfscanf(gmesh.ELE_INFOS(I,5),fid,'%d')';
+ select gmesh.ELE_INFOS(I,2)
+ case 15 then // point
+ gmesh.nbPoints = gmesh.nbPoints + 1;
+ gmesh.POINTS =[gmesh.POINTS; IDS(NODES_ELEM (1))' I];
+ case 1 then // beam
+ gmesh.nbLines = gmesh.nbLines + 1;
+ gmesh.LINES =[gmesh.LINES; IDS(NODES_ELEM (1:2))' I];
+ case 2 then // triangle
+ gmesh.nbTriangles = gmesh.nbTriangles + 1;
+ gmesh.TRIANGLES = [gmesh.TRIANGLES; IDS(NODES_ELEM (1:3))' I];
+ case 3 then // quadrangle
+ gmesh.nbQuads = gmesh.nbQuads + 1;
+ gmesh.QUADS =[gmesh.QUADS; IDS(NODES_ELEM (1:4))' I];
+ case 4 then // tetrahedron (4 node)
+ gmesh.nbTet = gmesh.nbTet + 1;
+ gmesh.TETS = [gmesh.TETS; IDS(NODES_ELEM (1:4))' I];
+ case 5 then // hexahedron (8 nodes)
+ gmesh.nbHex = gmesh.nbHex + 1;
+ gmesh.HEXA =[gmesh.HEXA; IDS(NODES_ELEM (1:8))' I];
+ case 6 then // prism (6 nodes)
+ gmesh.nbPrism = gmesh.nbPrism + 1;
+ gmesh.PRISM =[gmesh.PRISM; IDS(NODES_ELEM (1:6))' I];
+ case 7 then // pyramid (5 nodes)
+ gmesh.nbPyr = gmesh.nbPyr + 1;
+ gmesh.PYRAMID =[gmesh.PYRAMID; IDS(NODES_ELEM (1:5))' I];
+ case 9 then // second order triangle (6 nodes)
+ gmesh.nbQTriangles = gmesh.nbQTriangles + 1;
+ gmesh.QTRIANGLES = [gmesh.QTRIANGLES; IDS(NODES_ELEM (1:6))' I];
+ case 10 then // second order quadrangle (9 nodes)
+ gmesh.nbQQuads = gmesh.nbQQuads + 1;
+ gmesh.QQUADS = [gmesh.QQUADS; IDS(NODES_ELEM (1:9))' I];
+ else
+ disp(' ')
+ warning('Unknown element type '+string(gmesh.ELE_INFOS(I,2)) +' !')
+ end
+ end
+ tline = mgetl(fid,1);
+ disp('elements have been read')
+ end
+ end
+ end
+mclose (fid);
+endfunction
--
GitLab