--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_ao_lincom.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,89 @@
+% Test script for ao/lincom
+%
+% $Id: test_ao_lincom.m,v 1.5 2010/03/19 16:56:24 ingo Exp $
+%
+
+%% Make some data
+ple = plist('Exceptions', {'created', 'proctime', 'UUID', 'methodInvars', 'plistUsed', 'ao/name'});
+
+fs    = 10;
+nsecs = 10;
+x1 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'T'));
+x2 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'm'));
+x3 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'C'));
+n  = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'm'));
+c = [ao(1,plist('yunits','m/T')) ao(2,plist('yunits','m/m')) ao(3,plist('yunits','m C^-1'))];
+y = c(1)*x1 + c(2)*x2 + c(3)*x3 + n;
+y.simplifyYunits;
+% Get a fit for c
+p_m = lscov(x1, x2, x3, y);
+
+% do linear combination: list + plist with pest
+yfit_01 = lincom(x1, x2, x3, plist('coeffs', p_m));
+r01 = yfit_01.rebuild();
+eq(yfit_01, r01, ple)
+
+
+% do linear combination: list + plist with vector of cdata aos
+yfit_02 = lincom(x1, x2, x3, plist('coeffs', p_m.find));
+r02 = yfit_02.rebuild();
+eq(yfit_02, r02, ple)
+
+% do linear combination: vector + plist with pest
+yfit_03 = lincom([x1, x2, x3], plist('coeffs', p_m));
+r03 = yfit_03.rebuild();
+eq(yfit_03, r03, ple)
+
+% do linear combination: vector + plist with vector of cdata aos
+yfit_04 = lincom([x1, x2, x3], plist('coeffs', p_m.find));
+r04 = yfit_04.rebuild();
+eq(yfit_04, r04, ple)
+
+% do linear combination: list + pest
+yfit_05 = lincom(x1, x2, x3, p_m);
+r05 = yfit_05.rebuild();
+eq(yfit_05, r05, ple)
+
+% do linear combination: list + vector of cdata aos
+yfit_06 = lincom(x1, x2, x3, p_m.find);
+r06 = yfit_06.rebuild();
+eq(yfit_06, r06, ple)
+
+% do linear combination: vector pest
+yfit_07 = lincom([x1, x2, x3], p_m);
+r07 = yfit_07.rebuild();
+eq(yfit_07, r07, ple)
+
+% % do linear combination: vector + vector of cdata aos
+yfit_08 = lincom([x1, x2, x3], p_m.find);
+r08 = yfit_08.rebuild();
+eq(yfit_08, r08, ple)
+
+%% Make some data
+fs    = 10;
+nsecs = 10;
+B1 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'T'));
+B2 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'T'));
+B3 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'T'));
+n  = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs, 'yunits', 'm'));
+c = ao([1 2 3], plist('yunits','m/T')); 
+
+% do linear combination: list + cdata ao with coefficients
+yfit_11 = lincom(B1, B2, B3, c);
+r11 = yfit_11.rebuild();
+eq(yfit_11, r11, ple)
+
+% do linear combination: vector + cdata ao with coefficients
+yfit_12 = lincom([B1, B2, B3], c);
+r12 = yfit_12.rebuild();
+eq(yfit_12, r12, ple)
+
+% do linear combination: list + plist with vector of double
+yfit_13 = lincom(B1, B2, B3, plist('coeffs', p_m.y));
+r13 = yfit_13.rebuild();
+eq(yfit_13, r13, ple)
+
+% do linear combination: vector + plist with vector of double
+yfit_14 = lincom([B1, B2, B3], plist('coeffs', p_m.y));
+r14 = yfit_14.rebuild();
+eq(yfit_14, r14, ple)