--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_ao_minus.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,223 @@
+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 = minus(c1,c3);
+% cdata  + tsdata
+r2 = minus(c1,t1);
+% cdata  + fsdata
+r3 = minus(c1,f1);
+% cdata  + xydata
+r4 = minus(c1,x1);
+
+%---------- TSDATA
+% tsdata + tsdata
+r5 = minus(t1, t2);
+% tsdata + fsdata
+try
+  r6 = minus(t1, f2);
+  error('Data types should be incompatible: booo!');
+catch
+  lasterr
+  warning('Data types not compatible: fsdata+tsdata: hoorah!');
+end
+% tsdata + xydata
+r7 = minus(t1, x2);
+% tsdata + cdata
+r8 = minus(t1, c2);
+
+%---------- FSDATA
+% fsdata + tsdata
+try
+  r9 = minus(f1, t2);
+  error('Data types should be incompatible: booo!');
+catch
+  lasterr
+  warning('Data types not compatible: fsdata+tsdata: hoorah!');
+end
+% fsdata + fsdata
+r10 = minus(f1, f2);
+% fsdata + xydata
+r11 = minus(f1, x2);
+% fsdata + cdata
+r12 = minus(f1, c2);
+
+%---------- XYDATA
+% xydata + tsdata
+r13 = minus(x1, t2);
+% xydata + fsdata
+r14 = minus(x1, f2);
+% xydata + xydata
+r15 = minus(x1, x2);
+% xydata + cdata
+r16 = minus(x1, c2);
+
+
+
+
+%% Rule 2: vector + single
+
+v = [t1 t2];
+s = c2;
+
+% one output
+r1 = minus(v,s)
+r2 = minus(s,v)
+
+% multiple outputs
+[r1 r2] = minus(v,s)
+
+
+%% Rule 3: V_N + V_M
+
+v1 = [t1 t2];
+v2 = [c1 c2 c3];
+
+try
+  r = minus(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 = minus(v1,v2);
+[r1 r2 r3] = minus(v1,v2);
+
+%% Rule 5: matrix + single
+
+m = [t1 t2 t3; x1 x2 x3];
+s = c2;
+
+r1 = minus(m,s);
+r2 = minus(s,m);
+
+%% Rule 6: matrix + column vector
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1; c2];
+
+r1 = minus(m,v)
+r2 = minus(v,m)
+
+%% Rule 7: matrix + row vector
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1 c2 c3];
+
+r1 = minus(m,v)
+r2 = minus(v,m)
+
+%% Rule 8: M_NP + V_1Q (or V_Q1)
+
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1 c2];
+try
+  r = minus(m,v)
+  error('Rule 8 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 8 fails: hoorah!');
+end
+try
+  r = minus(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 = minus(m,v)
+  error('Rule 8 should fail: booo!');
+catch
+  lasterr
+  warning('Rule 8 fails: hoorah!');
+end
+try
+  r = minus(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 = minus(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 = minus(m,h)
+
+%% Rule 11: more than two inputs
+
+m = [t1 t2 t3; x1 x2 x3];
+v = [c1; c2];
+
+r = minus(m,v, c3, v, m)
+
+%% Check units
+
+
+% empty + m = m
+t1.setYunits('m');
+t2.setYunits('');
+r = minus(t1,t2)
+
+% V + m = error
+t2.setYunits('V');
+try
+  r = minus(t1,t2)
+  error('This unit check should fail: booo!');
+catch
+  lasterr
+  warning('The unit check failed: hoorah!');
+end
+