From a71826a08d3d357fa2aa7a6198756693c5101f57 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Thu, 27 Oct 2005 08:25:40 +0000
Subject: [PATCH] *** empty log message ***
---
utils/converters/matlab/load_gmsh.m | 21 ++++++++++++++++-----
1 file changed, 16 insertions(+), 5 deletions(-)
diff --git a/utils/converters/matlab/load_gmsh.m b/utils/converters/matlab/load_gmsh.m
index b6df2927a8..f46fd75d8b 100644
--- a/utils/converters/matlab/load_gmsh.m
+++ b/utils/converters/matlab/load_gmsh.m
@@ -47,10 +47,12 @@ function mesh = load_gmsh ( filename )
mesh.nbPoints = 0;
mesh.nbTriangles = 0;
mesh.nbQuads = 0;
- mesh.POINTS=zeros(mesh.nbElm,2);
- mesh.BEAMS=zeros(mesh.nbElm,3);
- mesh.TRIANGLES=zeros(mesh.nbElm,4);
- mesh.QUADS=zeros(mesh.nbElm,5);
+ mesh.nbTet = 0;
+ mesh.POINTS=zeros(mesh.nbElm,2);
+ mesh.BEAMS=zeros(mesh.nbElm,3);
+ mesh.TRIANGLES=zeros(mesh.nbElm,4);
+ mesh.QUADS=zeros(mesh.nbElm,5);
+ mesh.TETS=zeros(mesh.nbElm,5);
for(I=1:mesh.nbElm)
mesh.ELE_INFOS(I,:) = fscanf (fid,'%d',5);
NODES_ELEM = fscanf (fid,'%d',mesh.ELE_INFOS(I,5));
@@ -72,7 +74,7 @@ function mesh = load_gmsh ( filename )
mesh.TRIANGLES ( mesh.nbTriangles , 3) = IDS (NODES_ELEM ( 3) );
mesh.TRIANGLES ( mesh.nbTriangles , 4) = I;
end
- if (mesh.ELE_INFOS(I,2) == 4) %% quadrangle
+ if (mesh.ELE_INFOS(I,2) == 3) %% quadrangle
mesh.nbQuads = mesh.nbQuads + 1;
mesh.QUADS ( mesh.nbQuads , 1) = IDS (NODES_ELEM ( 1) );
mesh.QUADS ( mesh.nbQuads , 2) = IDS (NODES_ELEM ( 2) );
@@ -80,6 +82,15 @@ function mesh = load_gmsh ( filename )
mesh.QUADS ( mesh.nbQuads , 4) = IDS (NODES_ELEM ( 4) );
mesh.QUADS ( mesh.nbQuads , 5) = I;
end
+ if (mesh.ELE_INFOS(I,2) == 4) %% tet
+ mesh.nbTet = mesh.nbTet + 1;
+ mesh.QUADS ( mesh.nbTet , 1) = IDS (NODES_ELEM ( 1) );
+ mesh.QUADS ( mesh.nbTet , 2) = IDS (NODES_ELEM ( 2) );
+ mesh.QUADS ( mesh.nbTet , 3) = IDS (NODES_ELEM ( 3) );
+ mesh.QUADS ( mesh.nbTet , 4) = IDS (NODES_ELEM ( 4) );
+ mesh.QUADS ( mesh.nbTet , 5) = mesh.ELE_INFOS(I,3);
+ end
+
end
tline = fgetl(fid);
disp('elements have been read')
--
GitLab