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

APPENDIX B. PROGRAM LISTINGS 321
double rho = head->embouchureHole->rho;
double entryradius = WW_EMB_RADIUS;
double corr = 2.9370 * log(midi) - 11.6284;
return multz(flangedZ(f, c, rho, entryradius), real(corr));
}
TransferMatrix unitCellMatrix(double f, UnitCell c, double x) {
TransferMatrix m;
Map map = (strcmp(c->hole->fingering, "OPEN") == 0) ?
c->openMatrixMap :
c->closedMatrixMap;
if((x >= boreLength(c->bore)) && containsKey(map, f))
m = (TransferMatrix)get(map, f);
else {
m = traverseHoleMatrix(f, c->hole);
if((x > 0) && (c->bore != NULL))
rmultm(m, boreMatrix(f, c->bore, x));
if (x >= boreLength(c->bore)) {
clear(map);
put(map, f, m);
}
}
return m;
}
TransferMatrix woodwindMatrix(double f, Woodwind w, double entryratio,
double x) {
TransferMatrix m = identitym();
Head h;
complex branchZ;
UnitCell cell;
int cellCount = 0;
if (x >= 0) {
rmultm(m, headMatrix(f, w->head, entryratio, x));
x -= boreLength(w->head->downstreamBore);
while (x > 0 && cellCount < sizeVector(w->cells)) {
cell = (UnitCell)elementAt(w->cells, cellCount);
rmultm(m, unitCellMatrix(f, cell, x));
x -= boreLength(cell->bore);
cellCount++;
}
}
else {
h = w->head;
if (h->embouchureHole != NULL) {
branchZ = woodwindDownstreamZ(f, w);
rmultm(m, embouchureMatrix(f, h->embouchureHole, entryratio,
branchZ));
}
rmultm(m, boreMatrix(f, h->upstreamBore, -x));
}
return m;
}
int getZ0_c(Woodwind w, double x, complex* Z0, double* c) {
Vector bore;