forked from fangq/iso2mesh
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plottetra.m
77 lines (71 loc) · 2.16 KB
/
plottetra.m
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
function hm=plottetra(node,elem,varargin)
%
% hm=plottetra(node,elem,opt)
%
% plot 3D surface meshes
%
% author: Qianqian Fang <fangq at nmr.mgh.harvard.edu>
%
% input:
% node: a node coordinate list, 3 columns for x/y/z; if node has a
% 4th column, it will be used to set the color at each node.
% elem: a tetrahedral element list; if elem has a 5th column,
% and the 5th column are all integers, it will be treated
% as labels of sub-domains and display them in different colors.
% if the 5th column contains non-integer values, it will be
% used to map to the color of triangles.
% opt: additional options for a patch object, see plotmesh
%
% output:
% hm: handle or handles (vector) to the plotted surfaces
%
% example:
%
% h=plottetra(node,elem);
% h=plottetra(node,elem,'facealpha',0.5);
%
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%
rngstate = rand ('state');
randseed=hex2dec('623F9A9E'); % "U+623F U+9A9E"
if(isoctavemesh) randseed=randseed+3; end
if(~isempty(getvarfrom({'caller','base'},'ISO2MESH_RANDSEED')))
randseed=getvarfrom({'caller','base'},'ISO2MESH_RANDSEED');
end
rand('state',randseed);
if(~iscell(elem))
if(size(elem,2)>4)
tag=elem(:,5);
types=unique(tag);
hold on;
h=[];
if(all(isinteger(tag)))
for i=1:length(types)
idx=find(tag==types(i));
face=volface(elem(idx,1:4));
if(size(node,2)==3)
h=[h plotsurf(node,face,'facecolor',rand(3,1),varargin{:})];
else
h=[h plotsurf(node,face,varargin{:})];
end
end
else
[face,eid]=volface(elem(:,1:4));
if(size(elem,2)>=5)
h=[h plotsurf(node,face,'FaceVertexCData',tag(eid), varargin{:})];
else
h=[h plotsurf(node,face,varargin{:})];
end
end
else
face=volface(elem(:,1:4));
h=plotsurf(node,face,varargin{:});
end
end
if(~isempty(h))
axis equal;
end
if(~isempty(h) & nargout>=1)
hm=h;
end
rand ('state',rngstate);