** 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
- Comments:
None.

** Next:** Function: splitup_freq()
** Up:** GRASP Routines: Gravitational Radiation
** Previous:** How does the test
** Contents**
Bruce Allen
2000-11-19