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
78
79
80
81
82
83
|
function CreateDiscMaterial(filename, data, mat_db, mesh)
% function CreateDiscMaterial(filename, data, mat_db, mesh)
%
% Create the discrete material hdf5 file (version 2) usable by AddDiscMaterial
%
% Note: This function currently requires Matlab. Octave is missing the
% necessary hdf5 write functions.
%
% arguments:
% filename: hdf5 file to create (must not exist)
% data: voxel based index data, index as used in mat_db-1
% mat_db: material database
% mesh: used voxel mesh. Note size is size(data)+[1 1 1]
%
% example:
% mat_db.epsR = [1 3 4]; %relative permittivity
% mat_db.kappa = [0 0.2 0.4]; %electric conductivity (S/m)
% mat_db.density = [0 1000 1010]; %material density (kg/m³)
% mat_db.Name = {'Background','mat2','mat3'};
%
% data = [0 1 0; 2 1 2; 0 1 0]; % 3x3x3 data
% mesh.x = [0 0.1 0.2 0.3]; % 4 mesh lines in x-dir (3 cells)
% mesh.y = [-0.1 0 0.1 0.2]; % 4 mesh lines in y-dir (3 cells)
% mesh.z = [0 0.4 0.8 1.2]; % 4 mesh lines in z-dir (3 cells)
%
% CreateDiscMaterial('test_mat.h5', data, mat_db, mesh);
%
% See also AddDiscMaterial
%
% CSXCAD matlab interface
% -----------------------
% author: Thorsten Liebig (2013)
if isOctave
error('CreateDiscMaterial currently does not support Octave, due to missing hdf5 functions!');
end
if (exist('h5create')==0)
error('CSXCAD:CreateDiscMaterial','Your matlab seems to be too old, h5create cannot be found!');
end
data_size = size(data);
mesh_size = [numel(mesh.x) numel(mesh.y) numel(mesh.z)];
if ((mesh_size-1)~=data_size)
error('data size and mesh size mismatch');
end
if (exist(filename,'file'))
error(['file "' filename '" already exist. Delete/rename first!']);
end
h5create(filename, '/DiscData',data_size, 'Datatype', 'uint8', 'ChunkSize',data_size, 'Deflate',9);
h5write(filename, '/DiscData', data);
clear data;
h5writeatt(filename, '/DiscData', 'DB_Size', int32(numel(mat_db.epsR)));
h5writeatt(filename, '/DiscData', 'epsR', mat_db.epsR);
h5writeatt(filename, '/DiscData', 'kappa', mat_db.kappa);
h5writeatt(filename, '/DiscData', 'density', mat_db.density);
h5writeatt(filename, '/DiscData', 'Name', strjoin(mat_db.Name,','));
h5create(filename, '/mesh/x', mesh_size(1));
h5write(filename, '/mesh/x', mesh.x);
h5create(filename, '/mesh/y', mesh_size(2));
h5write(filename, '/mesh/y', mesh.y);
h5create(filename, '/mesh/z', mesh_size(3));
h5write(filename, '/mesh/z', mesh.z);
h5writeatt(filename, '/', 'Version', 2);
function out = strjoin(names, delimitier)
out = names{1};
for n=2:numel(names)
out = [out delimitier names{n}];
end
|
