Newer
Older
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
84
85
86
87
88
89
90
91
92
93
94
%% export points to inp
samples = {'out_out_CD','out_out_MD','in_in_CD','in_in_MD'};
degree = {'0','45','90','180'};
flag_save = 1;
load_folder_layers = 'Results\layer\transfer\align_0dg_rest\';
save_folder = 'Results\export\';
if ~exist(save_folder,'dir')
mkdir(save_folder)
end
% if you want to run the measure separately you can use the for loop
%for i = 1
sample = samples{i};
disp(['Load Data ',sample])
load([load_folder_layers,sample,'_aligned'])
disp('Export mesh ...')
f = parfeval(@elements_mesh,1,aligned_inner_dg0);
ele = fetchOutputs(f);
sz0 = size(aligned_inner_dg0);
%0 dg
no = [reshape(permute(aligned_inner_dg0(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg0(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg0(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'inner_',sample,'_0']);
no = [reshape(permute(aligned_outer_dg0(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg0(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg0(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'outer_',sample,'_0']);
% 45
no = [reshape(permute(aligned_inner_dg45(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg45(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg45(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'inner_',sample,'_45']);
no = [reshape(permute(aligned_outer_dg45(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg45(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg45(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'outer_',sample,'_45']);
% 90
no = [reshape(permute(aligned_inner_dg90(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg90(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg90(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'inner_',sample,'_90']);
no = [reshape(permute(aligned_outer_dg90(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg90(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg90(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'outer_',sample,'_90']);
% 180
no = [reshape(permute(aligned_inner_dg180(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg180(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_inner_dg180(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'inner_',sample,'_180']);
no = [reshape(permute(aligned_outer_dg180(:,:,1,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg180(:,:,2,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1]), ...
reshape(permute(aligned_outer_dg180(:,:,3,:),[2,1,3,4]),[sz0(1)*sz0(2)*sz0(4),1])];
no = [(1:size(no,1))',no];
disp('Write mesh ...')
parfeval(@write_mesh,0,no,ele,[save_folder,'outer_',sample,'_180']);
%end