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

242 APPENDIX B. PROGRAM LISTINGS
/*
Impedance.c
By Paul Dickens, 2005
Listing B.10: Impedance.c
Physical model of the acoustic impedance of a flute.
*/
#include
#include
#include
#include
#include "Woodwind.h"
#include "ParseXML.h"
int parseCommandLine(int argc, char** argv, char** holestring,
double* temp, double* humid, double* flo, double* fhi,
double* fres, double* entryratio, char** xml_filename);
/* Default parameter values */
#define TEMP 25.0
#define HUMID 0.5
#define X_CO2 0.0004
#define FLO 200.0
#define FHI 4000.0
#define FRES 2.0
#define ENTRYRATIO 1.0
int main(int argc, char** argv) {
char* holestring;
double temp, humid;
double f, flo, fhi, fres, entryratio;
char* xml_filename;
Woodwind instrument;
complex Z;
/* check correct usage */
if(!parseCommandLine(argc, argv, &holestring, &temp, &humid,
&flo, &fhi, &fres, &entryratio, &xml_filename)) {
fprintf(stderr, "Usage: Impedance [OPTIONS] \n\n");
fprintf(stderr, " Options:\n");
fprintf(stderr, "\t-s \n");
fprintf(stderr, "\t-t (default 25.0 degC)\n");
fprintf(stderr, "\t-u (default 0.5)\n");
fprintf(stderr, "\t-l (default 100.0)\n");
fprintf(stderr, "\t-h (default 4000.0)\n");
fprintf(stderr, "\t-r (default 2.0)\n");
fprintf(stderr, "\t-e (default 1.0)\n\n");
fprintf(stderr, " :\n");
fprintf(stderr,
"\t- Optional if no holes are defined in XML file.\n");
fprintf(stderr, "\t- Must be a sequence of 'O' (open hole) ");
fprintf(stderr, "and 'X' (closed hole) characters.\n");
fprintf(stderr,