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

326 APPENDIX B. PROGRAM LISTINGS
length = x;
radius2 = (s->radius2 * length + s->radius1
* (s->length - length)) / (s->length);
}
m = (radius1 == radius2) ?
tubeMatrix(f, s->c, s->rho, length, radius1, 1) :
coneMatrix(f, s->c, s->rho, length, radius1, radius2, 1);
return m;
}
TransferMatrix boreMatrix(double f, Vector bore, double x) {
TransferMatrix m = identitym();
int n = 0;
BoreSegment s;
while ((x > 0) && (n < sizeVector(bore))) {
s = (BoreSegment)elementAt(bore, n);
rmultm(m, boreSegmentMatrix(f, s, x));
x -= s->length;
n++;
}
return m;
}
TransferMatrix headMatrix(double f, Head h, double entryratio,
double x) {
TransferMatrix m, branchMatrix;
complex branchZ, ZL;
BoreSegment lastSegment;
Map map = h->matrixMap;
if((x >= boreLength(h->downstreamBore)) && containsKey(map, f))
m = (TransferMatrix)get(map,f);
else {
m = identitym();
if (h->embouchureHole != NULL) {
branchMatrix =
boreMatrix(f, h->upstreamBore, boreLength(h->upstreamBore));
lastSegment = (BoreSegment)elementAt(h->upstreamBore,
sizeVector(h->upstreamBore) - 1);
ZL = radiationZ(f, lastSegment->c, lastSegment->rho,
lastSegment->radius2, h->upstreamFlange);
branchZ = calcZin(branchMatrix, ZL);
rmultm(m,
embouchureMatrix(f, h->embouchureHole, entryratio, branchZ));
}
if (x > 0)
rmultm(m, boreMatrix(f, h->downstreamBore, x));
if (x >= boreLength(h->downstreamBore)) {
clear(map);
put(map, f, m);
}
}
return m;
}
complex faceZ(double f, Head head, int midi) {
double c = head->embouchureHole->c;