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_MPOWER a set of UTPs for the ao/mpower 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_mpower.m,v 1.6 2011/02/28 15:55:13 mauro 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 mpower method of the ao class computes the matrix power of the y data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	11 % of the two inputs.
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_mpower(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 = 'mpower';
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 input data shape == output data shape
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	44
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	45 disp('Done.');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	46 disp('******************************************************');
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 elseif nargin == 1 % Check for UTP functions
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	49 if strcmp(varargin{1}, 'isutp')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	50 results = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	51 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	52 results = 0;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	53 end
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 error('### Incorrect inputs')
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
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	58 %% UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	59
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	60 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	61 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	62 % Tests that the getInfo call works for this method.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	63 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	64 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	65 function result = utp_01
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
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	68 % <SyntaxDescription>
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 % 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	71 % individually.
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 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	74
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	75 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	76 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	77 % Call for no sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	78 io(1) = eval([class '.getInfo(''' mthd ''', ''None'')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	79 % Call for all sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	80 io(2) = eval([class '.getInfo(''' mthd ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	81 % Call for each set
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	82 for kk=1:numel(io(2).sets)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	83 io(kk+2) = eval([class '.getInfo(''' mthd ''', ''' io(2).sets{kk} ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	84 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	85 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	86 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	87 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	88 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	89 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	90 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	91
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	92 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	93 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	94 % 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	95 % 2) Check that all plists have the correct parameters.
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 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	98
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	99 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	100 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	101 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	102 % check we have minfo objects
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	103 if isa(io, 'minfo')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	104 %%% SET 'None'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	105 if ~isempty(io(1).sets), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	106 if ~isempty(io(1).plists), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	107 %%% Check all Sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	108 if ~any(strcmpi(io(2).sets, 'Default')), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	109 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	110 %%%%%%%%%% SET 'Default'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	111 if io(3).plists.nparams ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	112 % Check key
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	113 if ~io(3).plists.isparam('exponent'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	114 % Check default value
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	115 if ~isEmptyDouble(io(3).plists.find('exponent')), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	116 % Check options
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	117 if ~isequal(io(3).plists.getOptionsForParam('exponent'), {[]}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	118 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	119 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	120 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	121 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	122 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	123
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	124 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	125 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	126 end % END UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	127
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	128 %% UTP_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	129
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	130 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	131 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	132 % Tests that the mpower method works with a vector of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	133 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	134 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	135 function result = utp_02
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 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	138 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	139 % Test that the mpower method works for a vector of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	140 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	141 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	142
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	143 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	144 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	145 % We need a vector of AOs all the same size.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	146 % The shape of the data shouldn't have an effect to the algorithm.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	147 avec = [at4 ao(3) ao(1)];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	148 out = mpower(avec);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	149 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	150 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	151 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	152 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	153 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	154 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	155
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	156 % <AlgoDescription>
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 % 1) Check that the number of elements in 'out' is 1
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	159 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	160 % 3) Check the y-units
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	161 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	162 % </AlgoDescription>
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 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	165 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	166 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	167 % Check we have the one output
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	168 if numel(out) ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	169 % Check the output against the matrix power of the inputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	170 s = avec(1).data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	171 for jj=2:numel(avec)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	172 s = s ^ avec(jj).data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	173 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	174 if ~isequal(s, out.data.getY), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	175 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	176 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	177 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	178 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	179
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	180 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	181 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	182 end % END UTP_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	183
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	184 %% UTP_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	185
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	186 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	187 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	188 % Tests that the mpower method works with a matrix of AOs as input.
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 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	191 function result = utp_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	192
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	193 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	194 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	195 % Test that the mpower method works for a matrix of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	196 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	197 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	198
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	199 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	200 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	201 % We need a matrix of AOs all the same size
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	202 amat = [at4 ao(1) ao(2); ao(3) ao(4) ao(5)];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	203 out = mpower(amat);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	204 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	205 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	206 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	207 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	208 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	209 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	210
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	211 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	212 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	213 % 1) Check that the number of elements in 'out' is 1
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	214 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	215 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	216 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	217
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	218 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	219 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	220 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	221 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	222 if numel(out) ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	223 % Check the output against the matrix power of the inputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	224 s = amat(1).data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	225 for kk=2:numel(amat)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	226 s = s ^ amat(kk).data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	227 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	228 if ~isequal(s, out.data.getY), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	229 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	230 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	231 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	232 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	233
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	234 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	235 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	236 end % END UTP_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	237
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	238 %% UTP_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	239
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	240 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	241 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	242 % Tests that the mpower method works with a list of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	243 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	244 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	245 function result = utp_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	246
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	247 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	248 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	249 % Test that the mpower method works for a list of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	250 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	251 % </SyntaxDescription>
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 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	254 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	255 out = mpower(at4,ao(7),ao(4));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	256 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	257 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	258 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	259 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	260 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	261 end
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 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	264 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	265 % 1) Check that the number of elements in 'out' is 1.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	266 % 2) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	267 % 3) Check the y-units
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	268 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	269 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	270
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	271 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	272 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	273 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	274 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	275 if numel(out) ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	276 % Check the output against the matrix power of the inputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	277 s = at4.data.getY ^ 7 ^ 4;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	278 if ~isequal(s, out.data.getY), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	279 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	280 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	281 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	282 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	283
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	284 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	285 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	286 end % END UTP_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	287
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	288 %% UTP_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	289
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	290 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	291 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	292 % Tests that the mpower method works with a mix of different shaped AOs as
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	293 % input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	294 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	295 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	296 function result = utp_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	297
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	298 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	299 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	300 % Test that the mpower method works with an input of matrices and vectors
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	301 % and single AOs.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	302 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	303 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	304
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	305 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	306 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	307 % This UTP doen't make sense.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	308 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	309 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	310 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	311 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	312 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	313 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	314
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	315 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	316 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	317 % 1) Nothing to check
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	318 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	319 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	320
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	321 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	322 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	323 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	324 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	325 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	326 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	327 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	328
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	329 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	330 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	331 end % END UTP_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	332
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	333 %% UTP_06
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 % <TestDescription>
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 % Tests that the mpower method properly applies history.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	338 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	339 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	340 function result = utp_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	341
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	342 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	343 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	344 % Test that the result of applying the mpower method can be processed back
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	345 % to an m-file.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	346 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	347 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	348
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	349 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	350 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	351 out1 = mpower(at4,5);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	352 out2 = mpower(at4,ao(5));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	353 mout1 = rebuild(out1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	354 mout2 = rebuild(out2);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	355 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	356 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	357 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	358 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	359 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	360 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	361
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	362 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	363 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	364 % 1) Check that the last entry in the history of 'out1' and 'out2' corresponds to
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	365 % 'mpower'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	366 % 2) Check that the re-built objects are the same objectd as 'out1' and 'out2'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	367 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	368 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	369
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	370 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	371 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	372 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	373 % Check the last step in the history of 'out1' and 'out2'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	374 if ~strcmp(out1.hist.methodInfo.mname, 'mpower'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	375 if ~strcmp(out2.hist.methodInfo.mname, 'mpower'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	376 % Check the re-built objects
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	377 if ~eq(mout1, out1, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	378 if ~eq(mout2, out2, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	379 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	380 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	381 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	382 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	383
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	384 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	385 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	386 end % END UTP_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	387
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	388 %% UTP_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	389
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	390 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	391 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	392 % Tests that the mpower method can modify the input AO.
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 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	395 function result = utp_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	396
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	397 % <SyntaxDescription>
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 % Test that the mpower method can modify the input AO by calling with no
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	400 % output and that the method doesn't change the input of the function
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	401 % notation (with a equal sign).
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	402 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	403 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	404
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	405 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	406 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	407 % copy at4 to work with
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	408 ain = ao(at4);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	409 % modify ain
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	410 aout = ain.mpower(5);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	411 ain ^ 5;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	412 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	413 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	414 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	415 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	416 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	417 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	418
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	419 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	420 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	421 % 1) Check that 'at4' and 'ain' are now different.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	422 % 2) Check that 'ain' is mpower(at4).
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	423 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	424 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	425
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	426 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	427 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	428 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	429 % Check that mpower modified the input by comparing to the copy
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	430 if eq(ao(at4), ain, ple1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	431 % Check that mpower doesn't modified the input for the function notation
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	432 if ~eq(aout, ain, ple1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	433 % Check that the modified input is the mpower value of the copy
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	434 if ~isequal(mpower(at4.y, 5), ain.y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	435 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	436 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	437 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	438 end
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 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	441 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	442 end % END UTP_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	443
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	444 %% UTP_08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	445
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	446 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	447 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	448 % Test the shape of the 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 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	451 function result = utp_08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	452
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	453 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	454 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	455 % Test that the mpower method keeps the data shape of the input object. The
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	456 % input AO must be an AO with row data and an AO with column data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	457 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	458 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	459
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	460 try
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 % This UTP doesn't make sense because the AO must contain a square matrix
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	463 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	464 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	465 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	466 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	467 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	468 end
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 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	471 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	472 % 1) Nothinf to check
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	473 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	474 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	475
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	476 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	477 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	478 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	479 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	480 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	481 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	482 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	483
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	484 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	485 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	486 end % END UTP_08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	487
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	488 end

44

409a22968d5e Add unit tests