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

APPENDIX B. PROGRAM LISTINGS 317
double length, double boreRadius) {
EmbouchureHole h = (EmbouchureHole)malloc(sizeof(*h));
h->radiusin = radiusin;
h->radiusout = radiusout;
h->length = length;
h->boreRadius = boreRadius;
h->c = 0.0;
h->rho = 0.0;
return h;
}
Head createHead(EmbouchureHole embouchureHole, Vector upstreamBore,
double upstreamFlange, Vector downstreamBore) {
Head h = (Head)malloc(sizeof(*h));
h->embouchureHole = embouchureHole;
h->upstreamBore = upstreamBore;
h->upstreamFlange = upstreamFlange;
h->downstreamBore = downstreamBore;
h->matrixMap = createMap();
return h;
}
UnitCell createUnitCell(Hole hole, Vector bore) {
UnitCell c = (UnitCell)malloc(sizeof(*c));
c->hole = hole;
c->bore = bore;
c->openMatrixMap = createMap();
c->closedMatrixMap = createMap();
return c;
}
Woodwind createWoodwind(Head head, Vector cells, double flange) {
Woodwind w = (Woodwind)malloc(sizeof(*w));
w->head = head;
w->cells = cells;
w->flange = flange;
return w;
}
void setAirProperties(Woodwind w, double t_0, double t_amb,
double t_grad, double humid, double x_CO2) {
Head h = w->head;
double temp, x;
int segmentCount, cellCount;
BoreSegment s;
UnitCell cell;
temp = t_0;
if (w->head->embouchureHole != NULL) {
/* set c and rho for the embouchure hole */
h->embouchureHole->c = speedSound(temp, P_ATM, humid, x_CO2);
h->embouchureHole->rho =
densityAirGiacomo(temp, P_ATM, humid, x_CO2);
x = 0;
/* for each bore segment in upstream */
for (segmentCount = 0; segmentCount < sizeVector(h->upstreamBore);
segmentCount++) {