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


Function: splitup_freq2()

0 float splitup_freq2(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 more computationally-efficient version, which does not filter $\tilde h$ through each of the $p$ independent time domain filters. The arguments are identical to those of splitup_freq().

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_freq3() Up: GRASP Routines: Gravitational Radiation Previous: Function: splitup_freq()   Contents
Bruce Allen 2000-11-19