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

APPENDIX B. PROGRAM LISTINGS 329
double innerRadiationLengthCorrection(double a, double b) {
double delta = b / a;
/* Dalmont et al. (2002), eqn. (4) */
/* Nederveen et al. (1998), eqn. (40) */
return b * (0.82 - 1.4 * pow(delta, 2) + 0.75 * pow(delta, 2.7));
}
double closedHoleSeriesLengthCorrection(double a, double b,
double t) {
double delta = b / a;
if (t > 0)
return -b * pow(delta, 2) / (1.78 / tanh(1.84 * t / b) + 0.940 +
0.540 * delta + 0.285 * pow(delta, 2));
else if (t < 0)
return -b * pow(delta, 2) / (1.78 / (1.84 * t / b) + 0.940 +
0.540 * delta + 0.285 * pow(delta, 2));
else
return 0;
}
double openHoleSeriesLengthCorrection(double a, double b) {
double delta = b / a;
/* Dubos et al. (1999), eqn. (74), modified to remove dependence on
t */
return -b * pow(delta, 2) / (1.78 + 0.940 + 0.540 * delta + 0.285
* pow(delta, 2));
}