next up previous contents
Next: Function: splitup_freq2() Up: GRASP Routines: Gravitational Radiation Previous: Function: splitup()   Contents


Function: splitup_freq()

0 float splitup_freq(float c0, float c90, float *chirp0, float *chirp90, float norm, float* twice_inv_noise, int n, int offset, int p, int* indices, float* stats, float* working, float* htilde)
This routine returns the value of the statistic $r^2=\sum_{i=1}^p (\Delta S_i)^2$. This is a less-efficient version, which internally constructs filters for each of the different frequency subintervals, and then filters the metric perturbation through those filters. It is useful to understand how the different frequency components behave in the time domain, after filtering.

The arguments are:

c0: Input. The coefficient of the 0-phase template.
c90: Input. The coefficient of the $90^\circ$-phase template. Note that $c_0^2 + c_{90}^2$ should be 1.
chirp0: Input. An array chirp0[0..n-1] containing the FFT of the 0-phase chirp.
chirp90: Input. An array chirp90[0..n-1] containing the FFT of the $90^\circ$-phase chirp.
norm: Input. The normalization of the 0-phase chirp.
twice_inv_noise: Input. The array twice_inv_noise[0..n/2] contains $2/S_h(f)$, as described previously. The array element twice_inv_noise[0] contains the DC value, and the array element twice_inv_noise[n/2] contains the value at the Nyquist frequency.
n: Input. Defines the lengths of the previous arrays.
offset: Input. The offset of the moment of maximum signal in the filter output.
p: Input. The number of frequency bands $p$ for the vetoing test.
indices: Output. An array indices[0..p-1] used for internal storage of the frequency subintervals (see splitup().
stats: Output. An array stats[0..p-1] containing the values of the $S_i$ for $i=1,\cdots,p$.
working: Output. An array working[0..n-1] used for internal storage.
htilde: Input. An array htilde[0..n-1] containing the positive frequency part of $\tilde h(f)$.
Author: Bruce Allen, ballen@dirac.phys.uwm.edu
Comments: None.


next up previous contents
Next: Function: splitup_freq2() Up: GRASP Routines: Gravitational Radiation Previous: Function: splitup()   Contents
Bruce Allen 2000-11-19