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

APPENDIX B. PROGRAM LISTINGS 316
/*
Woodwind.c
By Paul Dickens, 2005, 2006.
Listing B.31: Woodwind.c
A woodwind instrument. See Woodwind.h for interface details.
*/
#include
#include
#include
#include
#include
#include "Woodwind.h"
#include "Acoustics.h"
#include "Map.h"
#define P_ATM 101325 // 1 atm pressure
BoreSegment createBoreSegment(double radius1, double radius2,
double length) {
BoreSegment s = (BoreSegment)malloc(sizeof(*s));
s->radius1 = radius1;
s->radius2 = radius2;
s->length = length;
s->c = 0.0;
s->rho = 0.0;
return s;
}
Hole createHole(double radius, double length, double boreRadius,
Key key) {
Hole h = (Hole)malloc(sizeof(*h));
h->radius = radius;
h->length = length;
h->boreRadius = boreRadius;
h->key = key;
h->c = 0.0;
h->rho = 0.0;
return h;
}
Key createKey(double radius, double holeRadius, double height,
double thickness, double wallThickness, double chimneyHeight) {
Key k = (Key)malloc(sizeof(*k));
k->radius = radius;
k->holeRadius = holeRadius;
k->height = height;
k->thickness = thickness;
k->wallThickness = wallThickness;
k->chimneyHeight = chimneyHeight;
return k;
}
EmbouchureHole createEmbouchureHole(double radiusin, double radiusout,