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

APPENDIX B. PROGRAM LISTINGS 325
for (segmentCount = 0; segmentCount < sizeVector(bore);
segmentCount++) {
s = (BoreSegment)elementAt(bore, segmentCount);
if (s->length length / maxLength);
radius1 = s->radius1;
length = s->length / numSegments;
for (newSegmentCount = 1; newSegmentCount radius2
+ (numSegments - newSegmentCount) * s->radius1) / numSegments;
newSegment = createBoreSegment(radius1, radius2, length);
if (newSegmentCount == 1)
setAt(bore, newSegment, segmentCount);
else
insertAt(bore, newSegment, ++segmentCount);
radius1 = radius2;
}
}
}
int setFingering(Woodwind w, char* holestring) {
int numholes = sizeVector(w->cells);
int i;
UnitCell cell;
/* check for no holes */
if((holestring == NULL) || (numholes == 0)) {
if((holestring == NULL) && (numholes == 0))
return 1;
else
return 0;
}
/* check for wrong number of holes entered */
if(strlen(holestring) != numholes)
return 0;
/* set fingering array based on command line hole string */
for(i = 0; i < numholes; i++) {
cell = (UnitCell)elementAt(w->cells, i);
if(holestring[i] == 'O')
cell->hole->fingering = "OPEN";
if(holestring[i] == 'X')
cell->hole->fingering = "CLOSED";
if((holestring[i] != 'O') && (holestring[i] != 'X'))
return 0;
}
return 1;
}
TransferMatrix boreSegmentMatrix(double f, BoreSegment s, double x) {
TransferMatrix m;
double length, radius1, radius2;
radius1 = s->radius1;
if (x >= s->length) {
length = s->length;
radius2 = s->radius2;
}
else {