jacobi_oo_STL
ascii_write_mesh.m
Go to the documentation of this file.
1 function ascii_write_mesh( xc, ia, e, basename)
2 %
3 % Saves the 2D triangular mesh in the minimal way (only coordinates, vertex connectivity, minimal boundary edge info)
4 % in an ASCII file.
5 % Matlab indexing is stored (starts with 1).
6 %
7 % The output file format is compatible with Mesh_2d_3_matlab:Mesh_2d_3_matlab(std::string const &fname) in jacobi_oo_stl/geom.h
8 %
9 % IN:
10 % coordinates xc: [2][nnode]
11 % connectivity ia: [4][nelem] with t(4,:) are the subdomain numbers
12 % edges e: [7][nedges] boundary edges
13 % e([1,2],:) - start/end vertex of edge
14 % e([3,4],:) - start/end values
15 % e(5,:) - segment number
16 % e([6,7],:) - left/right subdomain
17 % basename: file name without extension
18 %
19 % Data have been generated via <https://de.mathworks.com/help/pde/ug/initmesh.html initmesh>.
20 %
21 fname = [basename, '.txt'];
22 
24 ndim = int32(size(xc,1));
27 
28 
29 dlmwrite(fname,nnode,'delimiter','\t','precision',16) % number of nodes
30 dlmwrite(fname,ndim,'-append','delimiter','\t','precision',16) % space dimension
31 dlmwrite(fname,nelem,'-append','delimiter','\t','precision',16) % number of elements
32 dlmwrite(fname,nvert_e,'-append','delimiter','\t','precision',16) % number of vertices per element
33 
34 % dlmwrite(fname,xc(:),'-append','delimiter','\t','precision',16) % coordinates
35 dlmwrite(fname,xc([1,2],:).','-append','delimiter','\t','precision',16) % coordinates
36 
37 % no subdomain info transferred
38 tmp=int32(ia(1:3,:));
39 % dlmwrite(fname,tmp(:),'-append','delimiter','\t','precision',16) % connectivity in Matlab indexing
40 dlmwrite(fname,tmp(:,:).','-append','delimiter','\t','precision',16) % connectivity in Matlab indexing
41 
42 % store only start and end point of boundary edges,
43 nbedges = size(e,2);
44 dlmwrite(fname,nbedges,'-append','delimiter','\t','precision',16) % number boundary edges
45 tmp=int32(e(1:2,:));
46 % dlmwrite(fname,tmp(:),'-append','delimiter','\t','precision',16) % boundary edges in Matlab indexing
47 dlmwrite(fname,tmp(:,:).','-append','delimiter','\t','precision',16) % boundary edges in Matlab indexing
48 
49 end
xc
xc
Definition: ascii_read_meshvector.m:30
tmp
tmp(:,:).'
int32
ia(1:3,: int32()
dlmwrite
dlmwrite(fname, nnode, 'delimiter','\t', 'precision', 16) % number of nodes dlmwrite(fname
Mesh_2d_3_matlab
Definition: geom.h:347
connectivity
function vertex connectivity
Definition: ascii_write_mesh.m:3
ascii_write_mesh
function ascii_write_mesh(xc, ia, e, basename) % % Saves the 2D triangular mesh in the minimal way(only coordinates
txt
id txt()
ia
ia
Definition: ascii_read_meshvector.m:35
size
e, 2 size()
ndim
ndim
Definition: ascii_write_mesh.m:24
nnode
function vertex minimal boundary edge info in an ASCII file Matlab indexing is stored(starts with 1). % % The output file format is compatible with Mesh_2d_3_matlab nnode
Definition: ascii_write_mesh.m:23
nelem
nelem
Definition: ascii_write_mesh.m:25
nvert_e
nvert_e
Definition: ascii_write_mesh.m:26