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

APPENDIX B. PROGRAM LISTINGS 252
}
if(descent && ascent) {
/* run parabola least squares fit on points, add minimum to
vector */
for(j = 0; j < (NUM_POINTS/2 - TWEAK); j++) {
progressPoints(points, allpoints, i);
i++;
}
e = parabolaExt(points, MINIMUM, WEIGHT);
if(e != NULL)
addElement(extv, e);
/* reset flags and remember we're ascending */
seq_inc = 0;
seq_dec = 0;
descent = 0;
ascent = 1;
}
/* if graph decreases... */
if(((Point)elementAt(points, sizeVector(points) - 1))->y
< ((Point)elementAt(points, sizeVector(points) - 2))->y) {
}
/* tally number of points increased ... set flag when
we have moved up far enough to expect minimum/maximum */
seq_dec+=2;
if(seq_dec/2 >= TWEAK)
descent = 1;
if(seq_inc/2 > 0)
seq_inc-=2;
if(ascent && descent) {
/* run parabola least squares fit on points, add minimum to
vector */
for(j = 0; j < (NUM_POINTS/2 - TWEAK); j++) {
progressPoints(points, allpoints, i);
i++;
}
e = parabolaExt(points, MAXIMUM, WEIGHT);
if(e != NULL)
addElement(extv, e);
/* reset flags and remember we're descending*/
seq_inc = 0;
seq_dec = 0;
descent = 1;
ascent = 0;
}
/* if graph neither increases nor decreases... */
if(((Point)elementAt(points, sizeVector(points) - 1))->y
== ((Point)elementAt(points, sizeVector(points) - 2))->y) {
if(descent) {
seq_inc++;
seq_dec--;