Codebase list octave-iso2mesh / debian/1.9.6+ds-2 readtetgen.m
debian/1.9.6+ds-2

Tree @debian/1.9.6+ds-2 (Download .tar.gz)

readtetgen.m @debian/1.9.6+ds-2 — raw · history · blame 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
function [node,elem,face]=readtetgen(fstub)
%
% [node,elem,face]=readtetgen(fstub)
%
% read tetgen output files
%
% author: Qianqian Fang, <q.fang at neu.edu>
% date: 2007/11/21
%
% input:
%    fstub: file name stub
%
% output:
%    node: node coordinates of the tetgen mesh
%    elem: tetrahedra element list of the tetgen mesh
%    face: surface triangles of the tetgen mesh
%
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%

% read node file
fp=fopen([fstub,'.node'],'rt');
if(fp==0) 
	error('node file is missing!'); 
end
[dim,count] = fscanf(fp,'%d',4);
if(count<4) error('wrong node file'); end
node=fscanf(fp,'%f',[4,dim(1)]);
idx=node(1,:);
node=node(2:4,:)';
fclose(fp);

% read element file
fp=fopen([fstub,'.ele'],'rt');
if(fp==0) 
        error('elem file is missing!'); 
end
[dim,count] = fscanf(fp,'%d',3);
if(count<3) error('wrong elem file'); end
elem=fscanf(fp,'%d',[dim(2)+dim(3)+1,dim(1)]);
elem=elem';
elem(:,1)=[];
elem(:,1:dim(2))=elem(:,1:dim(2))+(1-idx(1));
fclose(fp);

% read surface mesh file
fp=fopen([fstub,'.face'],'rt');
if(fp==0)
        error('surface data file is missing!');
end
[dim,count] = fscanf(fp,'%d',2);
if(count<2) error('wrong surface file'); end
face=fscanf(fp,'%d',[5,dim(1)]);
face=[face(2:end-1,:)+1;face(end,:)]';
fclose(fp);