mc

%% Make test AOs

nsecs = 10000;
fs    = 10;

pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'sin(2*pi*7.433*t) + randn(size(t))');

a1 = ao(pl);
a1.setName('a1');
a1.setYunits('m');

% This should output N
iunits1 = unit('m');
ounits1 = unit('V');
iunits2 = unit('V');
ounits2 = unit('N');

lp = miir(plist('type', 'lowpass', 'fs', fs, 'fc', 0.1, 'order', 5,'iunits', iunits1, 'ounits', ounits1));
hp = miir(plist('type', 'highpass', 'fs', fs, 'fc', 1, 'order', 2,'iunits', iunits2, 'ounits', ounits2));

a1f = filter(a1, [lp hp], plist('bank', 'serial'));

%% Simplify

b = a1f.simplifyYunits()

%% Units with prefixes


a = ao(1:10);
a.setYunits('mm m^-1')

a.simplifyYunits


%% Units with prefixes


a = ao(1:10);
a.setYunits('mm m Hz um^-2')

a.simplifyYunits(plist('prefixes', true))

%% More tests on prefixes
N = 8;

test = true(N,1);

% Case 1
jj = 1;
% A unit is canceled out but there is a remaining prefix AND a remaining unit.
% Example: 'fN mN^-1 Hz'
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('fN mN^-1 Hz');

% The unit Newton will be canceled out and the remaining prefix (here 'p') will be added to the Hz
b = simplifyYunits(a, plist('prefixes', true));
c = simplifyYunits(a, plist('prefixes', false));

% Result: 'pHz'
if ne(b.yunits, unit('Hz')) || ne(c.yunits, unit('pHz'))
  test(jj) = false;
end

% Case 2
jj = 2;
%   A unit is canceled out but there is a remaining prefix AND NO remaining unit.
%   Example: 'fN mN^-1'
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('fN mN^-1');

%   In this case it is not possible to keep the remaining prefix
%   --> Don't simplify the unit (it is only possible if we compute the prefix to the values)
b = simplifyYunits(a, plist('prefixes', true));
c = simplifyYunits(a, plist('prefixes', false));

% Result:  'fN mN^-1'
if ne(b.yunits, unit('')) || ne(c.yunits, unit('fN mN^-1'))
  test(jj) = false;
end

% Case 3
jj = 3;
%   The 'normal' case
% Example: 'fN m^-2 mN^-2 Hz mm'
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('fN m^-2 mN^-2 Hz mm');

b = simplifyYunits(a, plist('prefixes', true));
c = simplifyYunits(a, plist('prefixes', false));

%   Result: 'GN^-1 km^-1 Hz'
if ne(b.yunits, unit('N^-1 m^-1 Hz')) || ne(c.yunits, unit('GN^-1 km^-1 Hz'))
  test(jj) = false;
end

% Case 4
jj = 4;
%   Power
% Example: ('mm')^2
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('mm');

b = simplifyYunits(a.^2, plist('prefixes', true));
c = simplifyYunits(a.^2, plist('prefixes', false));

%   Result: 'mm^2'
if ne(b.yunits, unit('m^2')) || ne(c.yunits, unit('mm^2'))
  test(jj) = false;
end

% Case 5
jj = 5;
%   Power
n = utils.math.randelement([-4:-1 1:4],1);
% Example: ('mm^-1 nm^-1')^n
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('mm^-1 nm^-1');

b = simplifyYunits(a.^n, plist('prefixes', true));
c = simplifyYunits(a.^n, plist('prefixes', false));

%   Result: 'um^-n'
if ne(b.yunits, unit(['m^' num2str(n*-2)])) || ne(c.yunits, unit(['um^' num2str(n*-2)]))
  test(jj) = false;
end

% Case 6
jj = 6;
%   A unit is canceled out but there is a remaining prefix AND a
%   remaining unit.
% Example: ('mm um ks^-1 s')
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('mm um ks^-1 s');

b = simplifyYunits(a, plist('prefixes', true))
c = simplifyYunits(a, plist('prefixes', false))

%   Result: 'um^2'
if ne(b.yunits, unit('m^2')) || ne(c.yunits, unit('um^2'))
  test(jj) = false;
end

% Case 7
jj = 7;
%   A unit is canceled out but there is a remaining prefix AND a
%   remaining unit.
% Example: ('mm um ks^-1 s')
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('mm um ks^-1 s');

b = simplifyYunits(sqrt(a), plist('prefixes', true))
c = simplifyYunits(sqrt(a), plist('prefixes', false))

%   Result: 'um^2'
if ne(b.yunits, unit('m')) || ne(c.yunits, unit('um'))
  test(jj) = false;
end

% Case 8
jj = 8;
%   Sqrt of Power
n = utils.math.randelement([-4:-2 2:4],1);
% Example: ('mm^-1 nm^-1')^2n
a = ao([1:10]',[2:2:20]',plist('xunits','s'));
a.setYunits('mm^-1 nm^-1');

b = simplifyYunits(sqrt(a.^(2*n)), plist('prefixes', true))
c = simplifyYunits(sqrt(a.^(2*n)), plist('prefixes', false))

%   Result: 'um^-n'
if ne(b.yunits, unit(['m^' num2str(n*-2)])) || ne(c.yunits, unit(['um^' num2str(n*-2)]))
  test(jj) = false;
end