Next: Function: graph()
Up: GRASP Routines: General purpose
Previous: Function: ratio()
void reciprocal(float *b, float *a, int ncomplex)
This routine computes the array for an input array containing
complex numbers. The arguments are:
Note that the arrays a and b can be the same.
- a: Input. The array of complex numbers
a[0..2N-1] with a[2j] and a[2j+1] respectively
containing the real and imaginary parts.
- b: Output. The array of complex numbers
b[0..2N-1] with b[2j] and b[2j+1] respectively
containing the real and imaginary parts of .
- ncomplex: Input. The number of complex numbers in the
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:
- Author: Jolien Creighton, email@example.com
- Comments: None.