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

APPENDIX B. PROGRAM LISTINGS 296
p2 = (equalz(Zload, inf)) ?
one :
divz(Zload, addz(Zload, one));
U2 = (equalz(Zload, inf)) ?
zero :
divz(one, addz(Zload, one));
p1 = addz(multz(m->A, p2), multz(m->B, U2));
U1 = addz(multz(m->C, p2), multz(m->D, U2));
Zin = divz(p1, U1);
return Zin;
}
void invertm(TransferMatrix m) {
complex det = subz(multz(m->A, m->D), multz(m->B, m->C));
complex A = m->A;
m->A = divz(m->D, det);
m->B = divz(multz(real(-1.0), m->B), det);
m->C = divz(multz(real(-1.0), m->C), det);
m->D = divz(A, det);
}