ltpda: testing/utp_1.1/utps/ao/utp_ao

author	Daniele Nicolodi <nicolodi@science.unitn.it>
date	Tue, 06 Dec 2011 18:42:11 +0100
parents
children

rev	line source
44 409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	1 % UTP_AO_DSMEAN a set of UTPs for the ao/dsmean method
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	2 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	3 % M Hewitson 06-08-08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	4 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	5 % $Id: utp_ao_dsmean.m,v 1.11 2009/09/20 16:51:33 hewitson Exp $
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	6 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	7
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	8 % <MethodDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	9 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	10 % The dsmean method of the ao class performs downsampling by taking the
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	11 % mean of every N samples of a time-series AOs.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	12 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	13 % </MethodDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	14
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	15 function results = utp_ao_dsmean(varargin)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	16
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	17 % Check the inputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	18 if nargin == 0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	19
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	20 % Some keywords
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	21 class = 'ao';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	22 mthd = 'dsmean';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	23
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	24 results = [];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	25 disp('******************************************************');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	26 disp(['**** Running UTPs for ' class '/' mthd]);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	27 disp('******************************************************');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	28
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	29 % Test AOs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	30 [at1,at2,at3,at4,at5,at6,atvec,atmat] = eval(['get_test_objects_' class]);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	31
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	32 % Exception list for the UTPs:
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	33 [ple1,ple2,ple3,ple4,ple5,ple6] = get_test_ples();
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	34
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	35 % Run the tests
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	36 results = [results utp_01]; % getInfo call
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	37 results = [results utp_02]; % Vector input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	38 results = [results utp_03]; % Matrix input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	39 results = [results utp_04]; % List input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	40 results = [results utp_05]; % Test with mixed input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	41 results = [results utp_06]; % Test history is working
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	42 results = [results utp_07]; % Test the modify call works
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	43 results = [results utp_08]; % Test the data shape
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	44 results = [results utp_09]; % Test output of the data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	45 results = [results utp_11(mthd, at1, ple1)]; % Test plotinfo doesn't disappear
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	46 results = [results utp_12(mthd, at1, ple1)]; % Test errors are cleared
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	47
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	48 disp('Done.');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	49 disp('******************************************************');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	50
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	51 elseif nargin == 1 % Check for UTP functions
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	52 if strcmp(varargin{1}, 'isutp')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	53 results = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	54 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	55 results = 0;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	56 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	57 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	58 error('### Incorrect inputs')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	59 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	60
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	61 %% UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	62
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	63 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	64 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	65 % Tests that the getInfo call works for this method.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	66 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	67 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	68 function result = utp_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	69
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	70
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	71 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	72 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	73 % Test that the getInfo call works for no sets, all sets, and each set
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	74 % individually.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	75 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	76 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	77
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	78 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	79 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	80 % Call for no sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	81 io(1) = eval([class '.getInfo(''' mthd ''', ''None'')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	82 % Call for all sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	83 io(2) = eval([class '.getInfo(''' mthd ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	84 % Call for each set
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	85 for kk=1:numel(io(2).sets)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	86 io(kk+2) = eval([class '.getInfo(''' mthd ''', ''' io(2).sets{kk} ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	87 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	88 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	89 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	90 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	91 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	92 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	93 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	94
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	95 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	96 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	97 % 1) Check that getInfo call returned an minfo object in all cases.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	98 % 2) Check that all plists have the correct parameters.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	99 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	100 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	101
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	102 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	103 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	104 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	105 % check we have minfo objects
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	106 if isa(io, 'minfo')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	107 %%% SET 'None'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	108 if ~isempty(io(1).sets), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	109 if ~isempty(io(1).plists), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	110 %%% Check all Sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	111 if ~any(strcmpi(io(2).sets, 'Default')), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	112 if numel(io(2).plists) ~= numel(io(2).sets), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	113 %%%%%%%%%% SET 'Default'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	114 if io(3).plists.nparams ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	115 % Check key
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	116 if ~io(3).plists.isparam('fsout'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	117 % Check default value
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	118 if ~isequal(io(3).plists.find('fsout'), 10), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	119 % Check options
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	120 if ~isequal(io(3).plists.getOptionsForParam('fsout'), {10}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	121 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	122 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	123 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	124 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	125 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	126
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	127 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	128 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	129 end % END UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	130
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	131 %% UTP_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	132
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	133 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	134 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	135 % Tests that the dsmean method works with a vector of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	136 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	137 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	138 function result = utp_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	139
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	140 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	141 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	142 % Test that the dsmean method works for a vector of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	143 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	144 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	145
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	146 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	147 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	148 fsout = 3;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	149 avec = [at1 at5 at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	150 out = dsmean(avec, plist('fsout', fsout));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	151 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	152 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	153 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	154 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	155 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	156 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	157
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	158 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	159 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	160 % 1) Check that the number of elements in 'out' is the square of the
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	161 % number in the input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	162 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	163 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	164 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	165
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	166 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	167 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	168 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	169 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	170 if numel(out) ~= numel(avec), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	171 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	172 for kk = 1:numel(avec)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	173 dsf = round(avec(kk).data.fs/fsout);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	174 n = floor(length(avec(kk).data.y) / dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	175 y = avec(kk).data.getY(1:n*dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	176 y = mean(reshape(y, dsf, length(y)/dsf)).';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	177 % Check y-axis
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	178 if ~isequal(out(kk).y, y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	179 % Check t0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	180 t0 = avec(kk).data.t0 + dsf/(2*avec(kk).data.fs);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	181 if ~eq(out(kk).t0, t0), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	182 % Check fs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	183 if ~eq(out(kk).fs, fsout), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	184 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	185 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	186 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	187 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	188 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	189
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	190 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	191 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	192 end % END UTP_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	193
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	194 %% UTP_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	195
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	196 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	197 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	198 % Tests that the dsmean method works with a matrix of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	199 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	200 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	201 function result = utp_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	202
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	203 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	204 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	205 % Test that the dsmean method works for a matrix of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	206 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	207 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	208
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	209 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	210 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	211 fsout = 3;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	212 amat = [at1 at5 at6; at5 at6 at1];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	213 out = dsmean(amat, plist('fsout', fsout));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	214 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	215 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	216 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	217 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	218 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	219 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	220
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	221 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	222 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	223 % 1) Check that the number of elements in 'out' is the square of the
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	224 % number in the input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	225 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	226 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	227 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	228
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	229 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	230 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	231 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	232 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	233 if numel(out) ~= numel(amat), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	234 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	235 for kk = 1:numel(amat)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	236 dsf = round(amat(kk).data.fs/fsout);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	237 n = floor(length(amat(kk).data.y) / dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	238 y = amat(kk).data.getY(1:n*dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	239 y = mean(reshape(y, dsf, length(y)/dsf)).';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	240 % Check y-axis
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	241 if ~isequal(out(kk).y, y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	242 % Check t0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	243 t0 = amat(kk).data.t0 + dsf/(2*amat(kk).data.fs);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	244 if ~eq(out(kk).t0, t0), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	245 % Check fs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	246 if ~eq(out(kk).fs, fsout), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	247 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	248 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	249 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	250 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	251 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	252
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	253 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	254 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	255 end % END UTP_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	256
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	257 %% UTP_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	258
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	259 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	260 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	261 % Tests that the dsmean method works with a list of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	262 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	263 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	264 function result = utp_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	265
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	266 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	267 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	268 % Test that the dsmean method works for a list of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	269 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	270 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	271
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	272 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	273 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	274 fsout = 3;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	275 out = dsmean(at1,at5,at6, plist('fsout', fsout));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	276 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	277 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	278 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	279 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	280 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	281 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	282
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	283 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	284 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	285 % 1) Check that the number of elements in 'out' is the square of the
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	286 % number in the input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	287 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	288 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	289 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	290
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	291 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	292 aoin = [at1, at5, at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	293 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	294 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	295 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	296 if numel(out) ~= 3, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	297 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	298 for kk = 1:numel(aoin)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	299 dsf = round(aoin(kk).data.fs/fsout);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	300 n = floor(length(aoin(kk).data.y) / dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	301 y = aoin(kk).data.getY(1:n*dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	302 y = mean(reshape(y, dsf, length(y)/dsf)).';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	303 % Check y-axis
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	304 if ~isequal(out(kk).y, y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	305 % Check t0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	306 t0 = aoin(kk).data.t0 + dsf/(2*aoin(kk).data.fs);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	307 if ~eq(out(kk).t0, t0), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	308 % Check fs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	309 if ~eq(out(kk).fs, fsout), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	310 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	311 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	312 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	313 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	314 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	315
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	316 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	317 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	318 end % END UTP_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	319
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	320 %% UTP_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	321
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	322 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	323 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	324 % Tests that the dsmean method works with a mix of different shaped AOs as
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	325 % input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	326 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	327 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	328 function result = utp_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	329
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	330 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	331 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	332 % Test that the dsmean method works with an input of matrices and vectors
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	333 % and single AOs.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	334 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	335 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	336
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	337 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	338 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	339 fsout = 2;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	340 out = dsmean(at1,[at5 at6],at5,[at5 at1; at6 at1],at6, plist('fsout', fsout));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	341 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	342 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	343 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	344 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	345 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	346 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	347
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	348 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	349 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	350 % 1) Check that the number of elements in 'out' is the same as in
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	351 % input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	352 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	353 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	354 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	355
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	356 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	357 aoin = [at1, reshape([at5 at6], 1, []), at5, reshape([at5 at1; at6 at1], 1, []), at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	358 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	359 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	360 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	361 if numel(out) ~= 9, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	362 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	363 for kk = 1:numel(aoin)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	364 dsf = round(aoin(kk).data.fs/fsout);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	365 n = floor(length(aoin(kk).data.y) / dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	366 y = aoin(kk).data.getY(1:n*dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	367 y = mean(reshape(y, dsf, length(y)/dsf)).';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	368 % Check y-axis
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	369 if ~isequal(out(kk).y, y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	370 % Check t0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	371 t0 = aoin(kk).data.t0 + dsf/(2*aoin(kk).data.fs);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	372 if ~eq(out(kk).t0, t0), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	373 % Check fs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	374 if ~eq(out(kk).fs, fsout), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	375 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	376 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	377 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	378 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	379 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	380
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	381 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	382 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	383 end % END UTP_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	384
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	385 %% UTP_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	386
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	387 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	388 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	389 % Tests that the dsmean method properly applies history.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	390 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	391 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	392 function result = utp_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	393
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	394 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	395 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	396 % Test that the result of applying the dsmean method can be processed back
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	397 % to an m-file.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	398 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	399 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	400
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	401 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	402 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	403 out = dsmean(at5,plist('fsout', 2));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	404 mout = rebuild(out);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	405 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	406 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	407 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	408 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	409 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	410 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	411
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	412 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	413 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	414 % 1) Check that the last entry in the history of 'out' corresponds to
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	415 % 'dsmean'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	416 % 2) Check that the re-built object is the same object as 'out'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	417 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	418 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	419
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	420 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	421 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	422 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	423 % Check the last step in the history of 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	424 if ~strcmp(out.hist.methodInfo.mname, 'dsmean'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	425 % Check the re-built object
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	426 if ~eq(mout, out, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	427 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	428 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	429 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	430 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	431
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	432 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	433 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	434 end % END UTP_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	435
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	436 %% UTP_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	437
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	438 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	439 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	440 % Tests that the dsmean method can modify the input AO.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	441 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	442 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	443 function result = utp_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	444
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	445 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	446 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	447 % Test that the dsmean method can modify the input AO by calling
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	448 % with no output.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	449 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	450 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	451
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	452 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	453 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	454 % copy at1 to work with
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	455 fsout = 3;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	456 pl = plist('fsout', fsout);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	457 ain = ao(at1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	458 % modify ain
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	459 aout = ain.dsmean(pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	460 ain.dsmean(pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	461 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	462 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	463 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	464 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	465 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	466 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	467
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	468 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	469 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	470 % 1) Check that 'at1' and 'ain' are now different.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	471 % 2) Check that 'ain' is dsmean(at1).
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	472 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	473 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	474
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	475 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	476 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	477 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	478 % Check that dsmean modified the input by comparing to the copy
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	479 if eq(ao(at1), ain, ple1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	480 % Check that dsmean doesn't modified the input for the function notation
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	481 if ~eq(aout, ain, ple1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	482 % Check that the modified input is the dsmean of the copy
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	483 dsf = round(at1.data.fs/fsout);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	484 n = floor(length(at1.data.y) / dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	485 y = at1.data.getY(1:n*dsf);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	486 y = mean(reshape(y, dsf, length(y)/dsf)).';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	487 % Check y-axis
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	488 if ~isequal(ain.y, y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	489 % Check t0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	490 t0 = at1.data.t0 + dsf/(2*at1.data.fs);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	491 if ~eq(ain.t0, t0), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	492 % Check fs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	493 if ~eq(ain.fs, fsout), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	494 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	495 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	496 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	497 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	498
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	499 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	500 result = utp_prepare_result(atest, stest, dbstack, mfilename);

44

409a22968d5e Add unit tests