--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_ao_multiply.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,217 @@
+mc
+
+%% Some AOs
+
+c1 = ao(1); c1.setName;
+c2 = ao(2); c2.setName;
+c3 = ao(3); c3.setName;
+c4 = ao(randn(3,3)); c4.setName;
+c5 = ao(randn(3,3)); c5.setName;
+c6 = ao(randn(3,1)); c6.setName;
+c7 = ao(randn(1,3)); c7.setName;
+
+
+t1 = ao(plist('tsfcn', 't', 'fs', 10, 'nsecs', 10)); t1.setName;
+t2 = ao(plist('tsfcn', 'sqrt(t)', 'fs', 10, 'nsecs', 10)); t2.setName;
+t3 = ao(plist('tsfcn', 't.^2', 'fs', 10, 'nsecs', 10)); t3.setName;
+
+f1 = ao(plist('fsfcn', 'f', 'f', 'logspace(-4,1,100)')); f1.setName;
+f2 = ao(plist('fsfcn', '1./f', 'f', 'logspace(-4,1,100)')); f2.setName;
+f3 = ao(plist('fsfcn', 'sqrt(f)', 'f', 'logspace(-4,1,100)')); f3.setName;
+
+x1 = ao(plist('xyfcn', 'x', 'x', '1:100')); x1.setName;
+x2 = ao(plist('xyfcn', '1./x', 'x', '1:100')); x2.setName;
+x3 = ao(plist('xyfcn', 'sqrt(x)', 'x', '1:100')); x3.setName;
+
+
+%% Rule 1
+
+%---------- CDATA
+% cdata  .* cdata
+r1 = times(c1,c3);
+% cdata  .* tsdata
+r2 = times(c1,t1);
+% cdata  .* fsdata
+r3 = times(c1,f1);
+% cdata  .* xydata
+r4 = times(c1,x1);
+
+%---------- TSDATA
+% tsdata .* tsdata
+r5 = times(t1, t2);
+% tsdata .* fsdata
+try
+  r6 = times(t1, f2);
+  error('Data types should be incompatible: booo!');
+catch
+  lasterr
+  warning('Data types not compatible: fsdata.*tsdata: hoorah!');
+end
+% tsdata .* xydata
+r7 = times(t1, x2);
+% tsdata .* cdata
+r8 = times(t1, c2);
+
+%---------- FSDATA
+% fsdata .* tsdata
+try
+  r9 = times(f1, t2);
+  error('Data types should be incompatible: booo!');
+catch
+  lasterr
+  warning('Data types not compatible: fsdata.*tsdata: hoorah!');
+end
+% fsdata .* fsdata
+r10 = times(f1, f2);
+% fsdata .* xydata
+r11 = times(f1, x2);
+% fsdata .* cdata
+r12 = times(f1, c2);
+
+%---------- XYDATA
+% xydata .* tsdata
+r13 = times(x1, t2);
+% xydata .* fsdata
+r14 = times(x1, f2);
+% xydata .* xydata
+r15 = times(x1, x2);
+% xydata .* cdata
+r16 = times(x1, c2);
+
+
+
+
+%% Rule 2: vector .* single
+
+v = [t1 t2];
+s = c2;
+
+% one output
+r1 = times(v,s)
+r2 = times(s,v)
+
+% multiple outputs
+[r1 r2] = times(v,s)
+
+
+%% Rule 3: V_N .* V_M
+
+v1 = [t1 t2];
+v2 = [c1 c2 c3];
+
+try
+  r = times(v1,v2);
+  error('Rule 3 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 3 fails: hoorah!');
+end
+
+%% Rule 4: V_N .* U_N
+
+v1 = [t1 t2 t3];
+v2 = [x1 x2 x3];
+
+r = times(v1,v2);
+[r1 r2 r3] = times(v1,v2);
+
+%% Rule 5: matrix .* single
+
+m = [t1 t2 t3; x1 x2 x3];
+s = c2;
+
+r1 = times(m,s);
+r2 = times(s,m);
+
+%% Rule 6: matrix .* column vector
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1; c2];
+
+r1 = times(m,v)
+r2 = times(v,m)
+
+%% Rule 7: matrix .* row vector
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1 c2 c3];
+
+r1 = times(m,v)
+r2 = times(v,m)
+
+%% Rule 8: M_NP .* V_1Q (or V_Q1)
+
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1 c2];
+try
+  r = times(m,v)
+  error('Rule 8 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 8 fails: hoorah!');
+end
+try
+  r = times(v,m)
+  error('Rule 8 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 8 fails: hoorah!');
+end
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1; c2; c3];
+try
+  r = times(m,v)
+  error('Rule 8 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 8 fails: hoorah!');
+end
+try
+  r = times(v,m)
+  error('Rule 8 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 8 fails: hoorah!');
+end
+
+%% Rule 9: M_NP .* H_PQ
+
+m = [t1 t2; x1 x2];
+h = [c1 c2 c3; c1 c2 c3];
+
+try
+  r = times(m,h)
+  error('Rule 9 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 9 fails: hoorah!');
+end
+
+%% Rule 10: matrix .* matrix
+
+m = [t1 t2 t3; x1 x2 x3];
+h = [c1 c2 c3; c1 c2 c3];
+
+r = times(m,h)
+
+%% Rule 11: more than two inputs
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1; c2];
+
+r = times(m,v, c3, v, m)
+
+%% Check units
+
+
+% empty .* m = m
+t1.setYunits('m');
+t2.setYunits('');
+r = times(t1,t2)
+
+% V .* m = mV
+t2.setYunits('V');
+r = times(t1,t2)
+