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

240 APPENDIX B. PROGRAM LISTINGS
complex subz(complex z1, complex z2) {
complex z;
z.Re = z1.Re - z2.Re;
z.Im = z1.Im - z2.Im;
return z;
}
complex multz(complex z1, complex z2) {
complex z;
z.Re = z1.Re*z2.Re - z1.Im*z2.Im;
z.Im = z1.Re*z2.Im + z1.Im*z2.Re;
return z;
}
complex divz(complex z1, complex z2) {
if(equalz(z2, zero)) return inf;
complex z;
z.Re = (z1.Re*z2.Re + z1.Im*z2.Im)/(z2.Re*z2.Re + z2.Im*z2.Im);
z.Im = (z1.Im*z2.Re - z1.Re*z2.Im)/(z2.Re*z2.Re + z2.Im*z2.Im);
return z;
}
complex expz(complex z1) {
complex z;
z.Re = exp(z1.Re)*cos(z1.Im);
z.Im = exp(z1.Re)*sin(z1.Im);
return z;
}
complex expjz(complex z1) {
return expz(multz(j, z1));
}
complex coshz(complex z1) {
complex minus_one = {-1.0, 0.0};
complex two = {2.0, 0.0};
complex minus_z1 = multz(minus_one, z1);
return divz(addz(expz(z1), expz(minus_z1)), two);
}
complex sinhz(complex z1) {
complex minus_one = {-1.0, 0.0};
complex two = {2.0, 0.0};
complex minus_z1 = multz(minus_one, z1);
return divz(subz(expz(z1), expz(minus_z1)), two);
}
complex cosz(complex z1) {
return coshz(multz(j, z1));
}
complex sinz(complex z1) {
return divz(sinhz(multz(j, z1)), j);
}
complex tanz(complex z1) {