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_DIFF a set of UTPs for the ao/diff 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_diff.m,v 1.14 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 diff method of the ao class computes the derivative of the input data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	11 % using different methods.
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_diff(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 = 'diff';
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 atvec = [at1 at5 at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	37 atmat = [at1 at5 at6; at5 at6 at1];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	38 pli = plist('method', '3point', 'neval', true);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	39
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	40 results = [results utp_01]; % getInfo call
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	41 results = [results utp_02(mthd, atvec, @algo_test_y, pli, ple3)]; % Vector input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	42 results = [results utp_03(mthd, atmat, @algo_test_y, pli, ple3)]; % Matrix input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	43 results = [results utp_04(mthd, at1, at5, at6, @algo_test_y, pli, ple3)]; % List input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	44 results = [results utp_05(mthd, at1, atvec, atmat, @algo_test_y, pli, ple3)]; % Test with mixed input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	45 results = [results utp_06(mthd, at1, pli, ple2)]; % Test history is working
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	46 results = [results utp_07(mthd, at1, pli, ple2)]; % Test the modify call works
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	47 results = [results utp_09(mthd, at5, at6)]; % Test input data shape == output data shape
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	48 results = [results utp_10(mthd, at1, at5, ple2)]; % Test output of the data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	49 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	50 results = [results utp_12(mthd, at1, ple1)]; % Test errors are cleared
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	51
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	52 results = [results utp_13]; % Test with plist: method = 'ORDER2'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	53 results = [results utp_14]; % Test with plist: method = 'ORDER2SMOOTH'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	54 results = [results utp_15]; % Test with plist: method = '5POINT'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	55
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	56
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	57 disp('Done.');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	58 disp('******************************************************');
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 elseif nargin == 1 % Check for UTP functions
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	61 if strcmp(varargin{1}, 'isutp')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	62 results = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	63 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	64 results = 0;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	65 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	66 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	67 error('### Incorrect inputs')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	68 end
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 %% Algorithm test for UTP 02,03,04,05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	71
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	72 function atest = algo_test_y(in, out, pli)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	73 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	74 % 3 point derivative
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	75 x = in.data.getX;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	76 dx = diff(x);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	77 y = in.data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	78 z = zeros(size(y));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	79 z(2:end-1) = (y(3:end)-y(1:end-2)) ./ (dx(2:end)+dx(1:end-1));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	80 z(1) = (y(2)-y(1)) ./ (dx(1));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	81 z(end) = 2*z(end-1)-z(end-2);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	82 if ~isequal(out.y, z), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	83 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	84
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	85 %% UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	86
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	87 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	88 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	89 % Tests that the getInfo call works for this method.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	90 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	91 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	92 function result = utp_01
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
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	95 % <SyntaxDescription>
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 % 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	98 % individually.
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 % </SyntaxDescription>
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 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	103 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	104 % Call for no sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	105 io(1) = eval([class '.getInfo(''' mthd ''', ''None'')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	106 % Call for all sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	107 io(2) = eval([class '.getInfo(''' mthd ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	108 % Call for each set
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	109 for kk=1:numel(io(2).sets)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	110 io(kk+2) = eval([class '.getInfo(''' mthd ''', ''' io(2).sets{kk} ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	111 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	112 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	113 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	114 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	115 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	116 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	117 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	118
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	119 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	120 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	121 % 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	122 % 2) Check that all plists have the correct parameters.
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 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	125
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	126 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	127 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	128 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	129 % check we have minfo objects
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	130 if isa(io, 'minfo')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	131 %%% SET 'None'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	132 if ~isempty(io(1).sets), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	133 if ~isempty(io(1).plists), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	134 %%% Check all Sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	135 if ~any(strcmpi(io(2).sets, 'Default')), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	136 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	137 %%%%%%%%%% SET 'Default'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	138 if io(3).plists.nparams ~= 4, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	139 % Check key
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	140 if ~io(3).plists.isparam('method'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	141 if ~io(3).plists.isparam('f0'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	142 if ~io(3).plists.isparam('order'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	143 if ~io(3).plists.isparam('coeff'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	144 % Check default value
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	145 if ~isequal(io(3).plists.find('method'), '2POINT'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	146 if ~isequal(io(3).plists.find('f0'), '1/Nsecs'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	147 if ~isequal(io(3).plists.find('order'), 'ZERO'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	148 if ~isEmptyDouble(io(3).plists.find('coeff')), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	149 % Check options
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	150 if ~isequal(io(3).plists.getOptionsForParam('method'), {'2POINT', '3POINT', '5POINT', 'ORDER2', 'ORDER2SMOOTH', 'FILTER', 'FPS'}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	151 if ~isequal(io(3).plists.getOptionsForParam('f0'), {'1/Nsecs'}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	152 if ~isequal(io(3).plists.getOptionsForParam('order'), {'ZERO', 'FIRST', 'SECOND'}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	153 if ~isequal(io(3).plists.getOptionsForParam('coeff'), {[]}), atest = false; end
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 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	156 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	157 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	158 end
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 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	161 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	162 end % END UTP_01
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 %% UTP_13
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 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	167 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	168 % Control the method with a plist.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	169 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	170 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	171 function result = utp_13
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	172
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	173 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	174 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	175 % Test the computation of derivative using a 2nd order
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	176 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	177 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	178
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	179 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	180 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	181 pl = plist('method', 'ORDER2');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	182 out = diff(at5, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	183 mout = rebuild(out);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	184 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	185 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	186 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	187 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	188 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	189 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	190
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	191 % <AlgoDescription>
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 % 1) Check that the diff method uses the 2nd order derivative.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	194 % 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	195 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	196 % </AlgoDescription>
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 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	199 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	200 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	201 % Compute derivative using a 2nd order
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	202 x = at5.data.getX;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	203 dx = diff(x);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	204 y = at5.data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	205 z = zeros(size(y));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	206 m = length(y);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	207 % y'(x1)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	208 z(1) = (1/dx(1)+1/dx(2))*(y(2)-y(1))+...
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	209 dx(1)/(dx(1)dx(2)+dx(2)^2)(y(1)-y(3));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	210 % y'(xm)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	211 z(m) = (1/dx(m-2)+1/dx(m-1))*(y(m)-y(m-1))+...
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	212 dx(m-1)/(dx(m-1)dx(m-2)+dx(m-2)^2)(y(m-2)-y(m));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	213 % y'(xi) (i>1 & i<m)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	214 dx1 = repmat(dx(1:m-2),1,1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	215 dx2 = repmat(dx(2:m-1),1,1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	216 y1 = y(1:m-2); y2 = y(2:m-1); y3 = y(3:m);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	217 z(2:m-1) = 1./(dx1.dx2.(dx1+dx2)).*...
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	218 (-dx2.^2.y1+(dx2.^2-dx1.^2).y2+dx1.^2.*y3);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	219 % Check the 2nd oder derivative
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	220 if ~isequal(out.y, z), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	221 % Check the re-built object
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	222 if ~eq(mout, out, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	223 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	224 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	225 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	226 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	227
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	228 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	229 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	230 end % END UTP_13
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	231
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	232 %% UTP_14
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 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	235 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	236 % Control the method with a plist.
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 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	239 function result = utp_14
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	240
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	241 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	242 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	243 % Test the computation of derivative using a 2nd order with a parabolic fit
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	244 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	245 % </SyntaxDescription>
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 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	248 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	249 pl = plist('method', 'ORDER2SMOOTH');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	250 out = diff(at5, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	251 mout = rebuild(out);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	252 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	253 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	254 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	255 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	256 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	257 end
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 % <AlgoDescription>
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 % 1) Check that the diff method uses the 2nd order derivative with a
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	262 % parabolic fit
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	263 % 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	264 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	265 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	266
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	267 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	268 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	269 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	270 % Compute derivative using a 2nd order with a parabolic fit
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	271 x = at5.data.getX;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	272 y = at5.data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	273 dx = diff(x);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	274 m = length(y);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	275 h = mean(dx);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	276 z = zeros(size(y));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	277 % y'(x1)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	278 z(1) = sum(y(1:5).*[-54; 13; 40; 27; -26])/70/h;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	279 % y'(x2)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	280 z(2) = sum(y(1:5).*[-34; 3; 20; 17; -6])/70/h;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	281 % y'(x{m-1})
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	282 z(m-1) = sum(y(end-4:end).*[6; -17; -20; -3; 34])/70/h;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	283 % y'(xm)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	284 z(m) = sum(y(end-4:end).*[26; -27; -40; -13; 54])/70/h;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	285 % y'(xi) (i>2 & i<(N-1))
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	286 Dc = [2 1 0 -1 -2];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	287 tmp = convn(Dc,y)/10/h;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	288 z(3:m-2) = tmp(5:m);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	289 % Check the 2nd oder derivative
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	290 if ~isequal(out.y, z), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	291 % Check the re-built object
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	292 if ~eq(mout, out, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	293 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	294 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	295 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	296 end
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 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	299 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	300 end % END UTP_14
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	301
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	302 %% UTP_15
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	303
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	304 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	305 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	306 % Control the method with a plist.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	307 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	308 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	309 function result = utp_15
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	310
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	311 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	312 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	313 % Test the 5 point derivative.
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 % </SyntaxDescription>
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 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	318 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	319 pl = plist('method', '5POINT');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	320 out = diff(at5, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	321 mout = rebuild(out);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	322 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	323 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	324 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	325 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	326 stest = 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 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	330 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	331 % 1) Check that the diff method uses the 5 point derivative.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	332 % 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	333 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	334 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	335
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	336 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	337 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	338 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	339 % Compute 5 point derivative
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	340 x = at5.data.getX;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	341 dx = diff(x);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	342 y = at5.data.getY;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	343 z = zeros(size(y));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	344 z(1) = (y(2)-y(1)) ./ (dx(1));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	345 z(2) = (y(3)-y(1))./(dx(2)+dx(1));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	346 z(3:end-2) = (-y(5:end) + 8.y(4:end-1) - 8.y(2:end-3) + y(1:end-4)) ./ (3.*(x(5:end)-x(1:end-4)));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	347 z(end-1) = 2*z(end-2)-z(end-3);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	348 z(end) = 2*z(end-1)-z(end-2);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	349 % Check the 5 point derivative
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	350 if ~isequal(out.y, z), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	351 % Check the re-built object
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	352 if ~eq(mout, out, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	353 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	354 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	355 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	356 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	357
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	358 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	359 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	360 end % END UTP_15
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 end

44

409a22968d5e Add unit tests