Codebase list octave-iso2mesh / debian/1.9.6+ds-3 meshasphere.m
debian/1.9.6+ds-3

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

meshasphere.m @debian/1.9.6+ds-3 — 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
function [node,face,elem]=meshasphere(c0,r,tsize,maxvol)
%
% [node,face,elem]=meshasphere(c0,r,tsize,maxvol)
%
% create the surface and tetrahedral mesh of a sphere
%
% author: Qianqian Fang, <q.fang at neu.edu>
%
% input: 
%   c0:  center coordinates (x0,y0,z0) of the sphere
%   r:   radius of the sphere
%   tsize: maximum surface triangle size on the sphere
%   maxvol: maximu volume of the tetrahedral elements
%
% output:
%   node: node coordinates, 3 columns for x, y and z respectively
%   face: integer array with dimensions of NB x 3, each row represents
%         a surface mesh face element 
%   elem: integer array with dimensions of NE x 4, each row represents
%         a tetrahedron 
%
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
% 

if(nargin<3)
    error('you must at least provide c0, r and tsize');
end
if(nargin==3)
    maxvol=tsize*tsize*tsize;
end
if(nargout==3)
    [node,face,elem]=meshunitsphere(tsize/r,maxvol/(r*r*r));
else
    [node,face]=meshunitsphere(tsize/r);
end

node=node*r+repmat(c0(:)',size(node,1),1);