Flute acoustics: measurement, modelling and design - School of ...

APPENDIX B. PROGRAM LISTINGS 233
/*
Complex.c
By Andrew Botros, 2000-2004
Listing B.9: Complex.c
Hyperbolic trig fns and expz added by Paul Dickens, 2005
Complex arithmetic library.
Refer to Complex.h for interface details.
Refer to the following for algebraic expressions:
Brown, J.W., Churchill, R.W., 1996.
Complex Variables and Applications.
6th ed. McGraw-Hill, New York,
*/
#include
#include
#include
#include "Complex.h"
/* zero, one, j and inf defined for convenience */
complex zero = {0.0, 0.0};
complex j = {0.0, 1.0};
complex one = {1.0, 0.0};
complex inf = {DBL_MAX, 0.0};
double modz(complex z) {
return sqrt(z.Re*z.Re + z.Im*z.Im);
}
double argz(complex z) {
return atan2(z.Im, z.Re);
}
complex real(double x) {
complex z = {x, 0.0};
return z;
}
complex imaginary(double y) {
complex z = {0.0, y};
return z;
}
int equalz(complex z1, complex z2) {
return (z1.Re == z2.Re) && (z1.Im == z2.Im);
}
complex addz(complex z1, complex z2) {
complex z;
z.Re = z1.Re + z2.Re;
z.Im = z1.Im + z2.Im;
return z;
}