void qn_normalize(float *u, float *q, float *r, int n, float *norm)Given a filter, , and twice the inverse power spectrum, , this routine generates a normalized template for which . The arguments are:

`u`: Output. The array`u[0..n-1]`contains the positive frequency part of the complex template function , packed as described in the Numerical Recipes routine`realft()`.`q`: Input. The array`q[0..n-1]`contains the positive frequency part of the complex filter function , also packed as described in the Numerical Recipes routine`realft()`.`r`: Input. The array`r[0..n/2]`contains the values of the real function used as a weight in the normalization. The array elements are arranged in order of increasing frequency from the DC component at subscript 0 to the Nyquist frequency at subscript`n`/2.`n`: Input. The total length of the arrays`u`and`q`. Must be even.`norm`: Output. The normalization constant, , defined below.

Given a filter, , this routine computes a template,
,
which is normalized so that , where is the inner
product defined by equation (). Thus, the normalization
constant is given by

(8.11.184) |

- Author: Jolien Creighton, jolien@tapir.caltech.edu