File: ReadHDF5Mesh.m

package info (click to toggle)
openems 0.0.35+dfsg.1-3
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 6,288 kB
  • sloc: cpp: 40,259; yacc: 580; lex: 350; makefile: 258; sh: 169; ruby: 19
file content (80 lines) | stat: -rw-r--r-- 2,280 bytes parent folder | download | duplicates (2)
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
function hdf_mesh = ReadHDF5Mesh(file)
% function hdf_mesh = ReadHDF5Mesh(file)
%
%   Get the raw mesh data stored in the hdf5 dump file created by openEMS
%
% returns:
% hdf_mesh.type     (0-> cartesian, 1-> cylindrical mesh type)
% hdf_mesh.names    (e.g. 'Mesh/y')
% hdf_mesh.lines    (e.g. [0,1,2,3,4])
%
% openEMS matlab interface
% -----------------------
% author: Thorsten Liebig
%
% See also ReadHDF5FieldData

isOctave = exist('OCTAVE_VERSION','builtin') ~= 0;
if isOctave
    hdf_mesh = ReadHDF5Mesh_octave(file);
    return
end

info = hdf5info(file);

for n=1:numel(info.GroupHierarchy.Groups)
    if strcmp(info.GroupHierarchy.Groups(n).Name,'/Mesh')
        for m=1:numel(info.GroupHierarchy.Groups(n).Datasets)
            names{m} = info.GroupHierarchy.Groups(n).Datasets(m).Name;
        end
    end
end

hdf_mesh.names = names;
for n=1:numel(names)
    hdf_mesh.lines{n} = double(hdf5read(file,names{n}));
end

if (strcmp(names{1},'/Mesh/alpha'))
    % alpha and rho are in the wrong order, flip to have rho, alpha, z
    hdf_mesh.names(1:2) = fliplr(hdf_mesh.names(1:2));
    hdf_mesh.lines(1:2) = fliplr(hdf_mesh.lines(1:2));
    hdf_mesh.type=1;
    return
end
if (strcmp(names{1},'/Mesh/phi'))
    % reorder coordinates
    hdf_mesh.names = hdf_mesh.names([2 3 1]);
    hdf_mesh.lines = hdf_mesh.lines([2 3 1]);
    hdf_mesh.type=2;
    return
end

hdf_mesh.type=0;


function hdf_mesh = ReadHDF5Mesh_octave(file)
hdf = load( '-hdf5', file );
hdf_mesh.names = fieldnames(hdf.Mesh);
hdf_mesh.type = 0; % cartesian mesh
for n=1:numel(hdf_mesh.names)
    hdf_mesh.lines{n} = hdf.Mesh.(hdf_mesh.names{n});
    hdf_mesh.names{n} = ['/Mesh/' hdf_mesh.names{n}];
    if strcmp(hdf_mesh.names{n},'/Mesh/alpha')
        hdf_mesh.type = 1; % cylindrical mesh
    end
    if strcmp(hdf_mesh.names{n},'/Mesh/phi')
        hdf_mesh.type = 2; % cylindrical mesh
    end
end

if (hdf_mesh.type==1)
    % alpha and rho are in the wrong order, flip to have rho, alpha, z
    hdf_mesh.names(1:2) = fliplr(hdf_mesh.names(1:2));
    hdf_mesh.lines(1:2) = fliplr(hdf_mesh.lines(1:2));
end
if (hdf_mesh.type==2)
    % alpha and rho are in the wrong order, flip to have rho, alpha, z
    hdf_mesh.names = hdf_mesh.names([2 3 1]);
    hdf_mesh.lines = hdf_mesh.lines([2 3 1]);
end