gael/bulkMatlab/averagesFields.m

function []=averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,varargin);
global averagesFields_count;

%algorithm:
%       first call (doInitFiles=1 doWriteMean=0) -> init files, arrays and set counter to 0     
%       loop calls (doInitFiles=0 doWriteMean=0) -> increment counter and arrays
%       dump calls (doInitFiles=0 doWriteMean=1) -> write one record, re-initialize arrays and counter
%       final call (doInitFiles=-1 doWriteMean=0) -> close files

%input: varargin is a series of structures

if doInitFiles==1&doWriteMean==0;
averagesFields_count=1;
for icur=1:length(varargin);
eval(['global struct' num2str(icur) '; struct' num2str(icur) '=varargin{icur}; ']);
vars_struct=fieldnames(varargin{icur});
eval(['global fid_struct' num2str(icur) '; fid_struct' num2str(icur) '=zeros(1,length(vars_struct));']); 
for vcur=1:length(vars_struct);
tmp1=cell2mat(vars_struct(vcur)); tmp2=[rep_out tmp1 suff_out];
eval(['fid_struct' num2str(icur) '(vcur)=fopen(''' tmp2 ''',''w'',''b'');']);
end
end

elseif doInitFiles==0&doWriteMean==0;
averagesFields_count=averagesFields_count+1;
for icur=1:length(varargin);
eval(['global struct' num2str(icur) ';']);
vars_struct=fieldnames(varargin{icur});
for vcur=1:length(vars_struct);
tmp1=cell2mat(vars_struct(vcur)); 
eval(['struct' num2str(icur) '.' tmp1 '=struct' num2str(icur) '.' tmp1 '+varargin{icur}.' tmp1 ';']);
end
end

elseif doInitFiles==0&doWriteMean==1;
for icur=1:length(varargin);
eval(['global fid_struct' num2str(icur) ' struct' num2str(icur) ';']);
vars_struct=fieldnames(varargin{icur});
for vcur=1:length(vars_struct);
tmp1=cell2mat(vars_struct(vcur));
eval(['fwrite(fid_struct' num2str(icur) '(vcur),struct' num2str(icur) '.' tmp1 '/averagesFields_count,''float32'');']);
eval(['struct' num2str(icur) '.' tmp1 '(:)=0;']);
end
end
averagesFields_count=0;

elseif doInitFiles==-1&doWriteMean==0;
for icur=1:length(varargin);
eval(['global fid_struct' num2str(icur) ' struct' num2str(icur) ';']);
vars_struct=fieldnames(varargin{icur});
for vcur=1:length(vars_struct); eval(['fclose(fid_struct' num2str(icur) '(vcur));']); end;
eval(['clear struct' num2str(icur) ';']);
end
averagesFields_count=[];

else
fprintf('error in averagesFields \n');
averagesFields_count=[];
end

1	function []=averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,varargin);
2	global averagesFields_count;
3
4	%algorithm:
5	% first call (doInitFiles=1 doWriteMean=0) -> init files, arrays and set counter to 0
6	% loop calls (doInitFiles=0 doWriteMean=0) -> increment counter and arrays
7	% dump calls (doInitFiles=0 doWriteMean=1) -> write one record, re-initialize arrays and counter
8	% final call (doInitFiles=-1 doWriteMean=0) -> close files
9
10	%input: varargin is a series of structures
11
12	if doInitFiles==1&doWriteMean==0;
13	averagesFields_count=1;
14	for icur=1:length(varargin);
15	eval(['global struct' num2str(icur) '; struct' num2str(icur) '=varargin{icur}; ']);
16	vars_struct=fieldnames(varargin{icur});
17	eval(['global fid_struct' num2str(icur) '; fid_struct' num2str(icur) '=zeros(1,length(vars_struct));']);
18	for vcur=1:length(vars_struct);
19	tmp1=cell2mat(vars_struct(vcur)); tmp2=[rep_out tmp1 suff_out];
20	eval(['fid_struct' num2str(icur) '(vcur)=fopen(''' tmp2 ''',''w'',''b'');']);
21	end
22	end
23
24	elseif doInitFiles==0&doWriteMean==0;
25	averagesFields_count=averagesFields_count+1;
26	for icur=1:length(varargin);
27	eval(['global struct' num2str(icur) ';']);
28	vars_struct=fieldnames(varargin{icur});
29	for vcur=1:length(vars_struct);
30	tmp1=cell2mat(vars_struct(vcur));
31	eval(['struct' num2str(icur) '.' tmp1 '=struct' num2str(icur) '.' tmp1 '+varargin{icur}.' tmp1 ';']);
32	end
33	end
34
35	elseif doInitFiles==0&doWriteMean==1;
36	for icur=1:length(varargin);
37	eval(['global fid_struct' num2str(icur) ' struct' num2str(icur) ';']);
38	vars_struct=fieldnames(varargin{icur});
39	for vcur=1:length(vars_struct);
40	tmp1=cell2mat(vars_struct(vcur));
41	eval(['fwrite(fid_struct' num2str(icur) '(vcur),struct' num2str(icur) '.' tmp1 '/averagesFields_count,''float32'');']);
42	eval(['struct' num2str(icur) '.' tmp1 '(:)=0;']);
43	end
44	end
45	averagesFields_count=0;
46
47	elseif doInitFiles==-1&doWriteMean==0;
48	for icur=1:length(varargin);
49	eval(['global fid_struct' num2str(icur) ' struct' num2str(icur) ';']);
50	vars_struct=fieldnames(varargin{icur});
51	for vcur=1:length(vars_struct); eval(['fclose(fid_struct' num2str(icur) '(vcur));']); end;
52	eval(['clear struct' num2str(icur) ';']);
53	end
54	averagesFields_count=[];
55
56	else
57	fprintf('error in averagesFields \n');
58	averagesFields_count=[];
59	end
60