dfer/matlab_stuff/concatstruct.m

function NewField = concatstruct(OldField1,OldField2);

name1=fieldnames(OldField1);
name2=fieldnames(OldField2);

if size(name1,1) ~= size(name2,1)
   disp('The 2 structures must have the same number of fields')
   NewField=[];
   return
end

NewField=[];

for i=1:size(name1,1)
 if char(name1(i)) ==  char(name2(i))
   nameloc = char(name1(i));
   disp(['Field : ',nameloc])
     %size(getfield(OldField1,nameloc));

   if ndims(getfield(OldField2,nameloc)) == 4
      disp(['    concatenating, 4d'])
      NewField=setfield(NewField,nameloc, ...
        cat(4,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));

   elseif ndims(getfield(OldField2,nameloc)) == 3
     if length(getfield(OldField2,'iters_read_from_file')) == 1
      disp(['    concatenating, 4d'])
      NewField=setfield(NewField,nameloc, ...
        cat(4,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
    else
      disp(['    concatenating, 3d'])
      NewField=setfield(NewField,nameloc, ...
        cat(3,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
    end
   elseif ndims(getfield(OldField2,nameloc)) == 2
       disp(['    concatenating, 2d'])
      NewField=setfield(NewField,nameloc, ...
        cat(3,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
   elseif strcmp(nameloc,'T')
      disp(['    concatenating, 1d'])
      NewField=setfield(NewField,nameloc, ...
        cat(1,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));

   elseif strcmp(nameloc,'iter')
      disp(['    concatenating, 1d'])
      NewField=setfield(NewField,nameloc, ...
        cat(1,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));

   elseif strcmp(nameloc,'iters_read_from_file')
      disp(['    concatenating, 1d'])
      NewField=setfield(NewField,nameloc, ...
        cat(1,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
   else
      disp(['    replacing'])
      NewField=setfield(NewField,nameloc,getfield(OldField1,nameloc));
 
   end

 end
end


return
1	dfer	1.1	function NewField = concatstruct(OldField1,OldField2);
2
3			name1=fieldnames(OldField1);
4			name2=fieldnames(OldField2);
5
6			if size(name1,1) ~= size(name2,1)
7			disp('The 2 structures must have the same number of fields')
8			NewField=[];
9			return
10			end
11
12			NewField=[];
13
14			for i=1:size(name1,1)
15			if char(name1(i)) == char(name2(i))
16			nameloc = char(name1(i));
17			disp(['Field : ',nameloc])
18			%size(getfield(OldField1,nameloc));
19
20			if ndims(getfield(OldField2,nameloc)) == 4
21			disp([' concatenating, 4d'])
22			NewField=setfield(NewField,nameloc, ...
23			cat(4,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
24
25			elseif ndims(getfield(OldField2,nameloc)) == 3
26			if length(getfield(OldField2,'iters_read_from_file')) == 1
27			disp([' concatenating, 4d'])
28			NewField=setfield(NewField,nameloc, ...
29			cat(4,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
30			else
31			disp([' concatenating, 3d'])
32			NewField=setfield(NewField,nameloc, ...
33			cat(3,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
34			end
35			elseif ndims(getfield(OldField2,nameloc)) == 2
36			disp([' concatenating, 2d'])
37			NewField=setfield(NewField,nameloc, ...
38			cat(3,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
39			elseif strcmp(nameloc,'T')
40			disp([' concatenating, 1d'])
41			NewField=setfield(NewField,nameloc, ...
42			cat(1,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
43
44			elseif strcmp(nameloc,'iter')
45			disp([' concatenating, 1d'])
46			NewField=setfield(NewField,nameloc, ...
47			cat(1,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
48
49			elseif strcmp(nameloc,'iters_read_from_file')
50			disp([' concatenating, 1d'])
51			NewField=setfield(NewField,nameloc, ...
52			cat(1,getfield(OldField1,nameloc),getfield(OldField2,nameloc)));
53			else
54			disp([' replacing'])
55			NewField=setfield(NewField,nameloc,getfield(OldField1,nameloc));
56
57			end
58
59			end
60			end
61
62
63			return