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
56
57
58
59
60
61
62
63
|
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'], 'rb');
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'], 'rb');
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'], 'rb');
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);
elem(:, 1:4) = meshreorient(node(:, 1:3), elem(:, 1:4));
|
