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

APPENDIX B. PROGRAM LISTINGS 323
int n;
BoreSegment s;
double length = 0;
for(n = 0; n < sizeVector(bore); n++) {
s = (BoreSegment)elementAt(bore, n);
length += s->length;
}
return length;
}
double woodwindLengthPos(Woodwind w) {
int cellCount;
UnitCell cell;
double length = 0;
length += boreLength(w->head->downstreamBore);
for (cellCount = 0; cellCount < sizeVector(w->cells); cellCount++) {
cell = (UnitCell)elementAt(w->cells, cellCount);
length += boreLength(cell->bore);
}
return length;
}
double woodwindLengthNeg(Woodwind w) {
return boreLength(w->head->upstreamBore);
}
complex impedance(double f, Woodwind w, double entryratio) {
TransferMatrix matrix =
woodwindMatrix(f, w, entryratio, woodwindLengthPos(w));
return calcZin(matrix, woodwindLoadZ(f, w));
}
complex playedImpedance(double f, Woodwind w, int midi) {
double entryradius = WW_EMB_RADIUS;
complex Z = impedance(f, w, entryradius / woodwindEntryRadius(w));
Z = addz(Z, faceZ(f, w->head, midi));
return Z;
}
double woodwindEntryRadius(Woodwind w) {
double a;
BoreSegment s;
if (w->head->embouchureHole != NULL)
a = w->head->embouchureHole->radiusout;
else {
s = (BoreSegment)elementAt(w->head->downstreamBore, 0);
a = s->radius1;
}
return a;
}
complex woodwindLoadZ(double f, Woodwind w) {
Vector lastBore;
UnitCell lastCell;
BoreSegment s;