% Test the computation of the in-loop acceleration terms 
% 
% M Hewitson 11-04-08
% 
% $Id: test_ltpda_mdc1_ifo2acc_inloop.m,v 1.1 2008/04/14 15:39:44 hewitson Exp $
% 


% Make a white-noise AO
fs    = 10;
nsecs = 100000;

o1  = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs));
o12 = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs));

[a1, a2] = ltpda_mdc1_ifo2acc_inloop(plist('o1', o1, 'o12', o12));

% Spectra
a1xx = ltpda_lpsd(a1);
a2xx = ltpda_lpsd(a2);

iplot(a1xx,a2xx)

% END