Next: Function: graph()
Up: GRASP Routines: General purpose
Previous: Function: ratio()
Contents
Function: reciprocal()
void reciprocal(float *b, float *a, int ncomplex)
This routine computes the array for an input array containing
complex numbers. The arguments are:
 a: Input. The array of complex numbers
a[0..2N1] with a[2j] and a[2j+1] respectively
containing the real and imaginary parts.
 b: Output. The array of complex numbers
b[0..2N1] with b[2j] and b[2j+1] respectively
containing the real and imaginary parts of .
 ncomplex: Input. The number of complex numbers in the
arrays.
Note that the arrays a[] and b[] can be the same.
In order to reduce the potential for overflows
(since floating point arithmetic was used), the reciprocals of the
complex numbers were computed according to the following formula:

(16.10.330) 
 Author: Jolien Creighton, jolien@tapir.caltech.edu
 Comments: None.
Bruce Allen
20001119