This function calculates the high resolution multitaper estimate of the (complex-valued) cross-correlation spectrum of the two input time-series and .

The arguments are:

`o1:`Input.`o1[0..npoints-1]`is an array of floating point variables containing the values of the input time-series .`o1[i]`contains the value of evaluated at the discrete time , where .`o2:`Input.`o2[0..npoints-1]`is an array of floating point variables containing the values of the input time-series , in exactly the same format as the previous argument.`npoints:`Input. The total number of data points contained in the two input time-series.`padded_length:`Input. The padded length is an integer power of 2, greater than (or equal to) the total number of data points. The tapered data sets are zero-padded out to this length. The total number of frequency bins (including DC and Nyquist) in the output cross-correlation spectrum is .`delta_t:`Input. The sampling period (in sec).`nwin:`Input. The total number of data tapers used when forming the multitaper spectral estimate.`nwdt:`Input. The (total sample time) (frequency resolution bandwidth) product .`ReImSpec12:`Output.`ReImSpec12[0...padded_length+1]`is an array of double precision variables containing the values of the high resolution multitaper spectral estimate of the (complex-valued) cross-correlation spectrum .`ReImSpec12[2*j]`and`ReImSpec12[2*j+1]`contain, respectively, the values of the real and imaginary parts of

(16.22.333)

If you want to obtain the same normalization as that used in the `avg_spec()` routine described by equation () for
the case where then the output array `ReImSpec12`
should be multiplied by a factor of .

- Author:
Adapted from the original code (Lees and Park) by Bruce Allen (ballen@dirac.phys.uwm.edu), Adrian Ottewill (ottewill@relativity.ucd.ie), and Joseph Romano (romano@csd.uwm.edu).

- Comments: None.