Next: Function: splitup_freq() Up: GRASP Routines: Gravitational Radiation Previous: How does the test   Contents

## Function: splitup()

0 void splitup(float *working, float template, float *r, int n, float total, int p, int *indices)
This routine takes as inputs a template and a noise-power spectrum, and splits up the frequency spectrum into a set of sub-intervals to use with the vetoing technique just described.

The arguments are:

working: Input. An array working[0..n-1] used for working space.
template: Input. The array template[0..n-1] contains the positive frequency () part of the complex function . The packing of into this array follows the scheme used by the Numerical Recipes routine realft(), which is described between equations (12.3.5) and (12.3.6) of [1]. The DC component is real, and located in template[0]. The Nyquist-frequency component is also real, and is located in template[1]. The array elements template[2] and template[3] contain the real and imaginary parts, respectively, of where . Array elements template[2j] and template[2j+1] contain the real and imaginary parts of for .
r: Input. The array r[0..n/2] contains the values of the real function which is twice the inverse of the receiver noise, as in equation (), so that . The array elements are arranged in order of increasing frequency, from the DC value at subscript 0, to the Nyquist frequency at subscript n/2. Thus, the 'th array element r[j] contains the real value , for . Again it is assumed that .
n: Input. The total length of the complex arrays template and working, and the number of points in the output array s. Note that the array r contains points. n must be even.
total: Input. This is the total value of the integrated template squared over ; the frequency subintervals are choose so that each of the p subintervals contains of this total.
p: Input. The number of frequency bands into which you want to divide the range from DC to .
indices: Ouput. The frequency bins of the first frequency band are i=0..indices[0]. The next frequency band is i=indices[0]+1..indices[1]. The p'th frequency band is i=indices[p-2]+1..indices[p-1]. Note that indices[p-1]=n-1.
Author: Bruce Allen, ballen@dirac.phys.uwm.edu