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

308 APPENDIX B. PROGRAM LISTINGS
/*
Waves.c
By Paul Dickens, 2005, 2006
Listing B.29: Waves.c
Calculates the pressure and flow distribution along a flute.
*/
#include
#include
#include
#include "Acoustics.h"
#include "Woodwind.h"
#include "ParseXML.h"
int parseCommandLine(int argc, char** argv,
char** holestring, double* xres, int* midi,
double* f, char** xml_filename);
/* Default x resolution */
#define XRES 2.0e-3;
int main(int argc, char** argv) {
double f;
double x, xres, xmin, xmax;
int midi;
char* holestring;
char* xml_filename;
Woodwind instrument;
complex Z0, Zin, pin, Uin;
double c;
TransferMatrix m;
complex p, U;
double entryradius = WW_EMB_RADIUS;
/* check correct usage */
if(!parseCommandLine(argc, argv, &holestring, &xres, &midi,
&f, &xml_filename)) {
fprintf(stderr,
"Usage: Waves [OPTIONS] \n\n");
fprintf(stderr, " Options:\n");
fprintf(stderr, "\t-s \n");
fprintf(stderr, "\t-r (default 2.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,
"\t- Length must be equal to the number of defined holes.\n");
fprintf(stderr, "\t- e.g. \"XXOOOOOOXOOOOXOOO\"\n\n");
return -1;
}
/* retrieve data structures from XML file */
if(!parseXMLFile(xml_filename, &instrument)) {