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

APPENDIX B. PROGRAM LISTINGS 318
s = (BoreSegment)elementAt(h->upstreamBore, segmentCount);
temp = t_0 + t_grad * (x + s->length / 2);
if (temp < t_amb) temp = t_amb;
s->c = speedSound(temp, P_ATM, humid, x_CO2);
s->rho = densityAirGiacomo(temp, P_ATM, humid, x_CO2);
x += s->length;
}
}
x = 0;
/* for each bore segment in downstream */
for (segmentCount = 0; segmentCount < sizeVector(h->downstreamBore);
segmentCount++) {
s = (BoreSegment)elementAt(h->downstreamBore, segmentCount);
temp = t_0 + t_grad * (x + s->length / 2);
if (temp < t_amb) temp = t_amb;
s->c = speedSound(temp, P_ATM, humid, x_CO2);
s->rho = densityAirGiacomo(temp, P_ATM, humid, x_CO2);
x += s->length;
}
/* for each unit cell */
for (cellCount = 0; cellCount < sizeVector(w->cells); cellCount++) {
cell = (UnitCell)elementAt(w->cells, cellCount);
temp = t_0 + t_grad * x;
if (temp < t_amb) temp = t_amb;
cell->hole->c = speedSound(temp, P_ATM, humid, x_CO2);
cell->hole->rho = densityAirGiacomo(temp, P_ATM, humid, x_CO2);
/* for each bore segment in unit cell */
for (segmentCount = 0; segmentCount < sizeVector(cell->bore);
segmentCount++) {
s = (BoreSegment)elementAt(cell->bore, segmentCount);
temp = t_0 + t_grad * (x + s->length / 2);
if (temp < t_amb) temp = t_amb;
s->c = speedSound(temp, P_ATM, humid, x_CO2);
s->rho = densityAirGiacomo(temp, P_ATM, humid, x_CO2);
x += s->length;
}
}
}
void discretiseWoodwind(Woodwind w, double maxLength) {
int cellCount;
UnitCell cell;
discretiseBore(w->head->upstreamBore, maxLength);
discretiseBore(w->head->downstreamBore, maxLength);
for (cellCount = 0; cellCount < sizeVector(w->cells); cellCount++) {
cell = (UnitCell)elementAt(w->cells, cellCount);
discretiseBore(cell->bore, maxLength);
}
}
void discretiseBore(Vector bore, double maxLength) {
int segmentCount, newSegmentCount, numSegments;
BoreSegment s, newSegment;
double radius1, radius2, length;
/* for each bore segment in bore */