#include <stdio.h>
#include <string.h>
#define MAX(A, B) ((A) > (B) ? (A) : (B))


#ifdef _MSC_VER
#define ANSILIBS
#define casecmp stricmp
#else
#define casecmp strcasecmp
#endif


struct stressdef {char *name ;  /* rhythm designator */
                  char *meter;
		  int nseg;    /* number of segments; */
                  int nval;     /* number of values; */
                  int vel[8];	/* segment velocities */
		  float expcoef[8]; /* segment expander coefficient */
                  } stresspat[32];


/* most of these externals link to variables in genmidi.c */
extern int segnum,segden,nseg;
extern float fdursum[32],fdur[32];
extern int ngain[32];
extern float maxdur;
extern int time_num,time_denom;
extern int verbose;
extern int beatmodel,stressmodel;
extern char timesigstring[16];  /* from parseabc.c */

void reduce(int *, int *);

void init_stresspat () {
   stresspat[0].name = "Hornpipe";
   stresspat[0].meter = "4/4";
   stresspat[0].nseg = 8;
   stresspat[0].nval = 2;
   stresspat[0].vel[0] = 110;
   stresspat[0].vel[1] = 90;
   stresspat[0].expcoef[0] = (float) 1.4;
   stresspat[0].expcoef[1] = (float) 0.6;

   stresspat[1].name = "Hornpipe";
   stresspat[1].meter = "C|";
   stresspat[1].nseg = 8;
   stresspat[1].nval = 2;
   stresspat[1].vel[0] = 110;
   stresspat[1].vel[1] = 90;
   stresspat[1].expcoef[0] = (float) 1.4;
   stresspat[1].expcoef[1] = (float) 0.6;

   stresspat[2].name = "Hornpipe";
   stresspat[2].meter = "2/4";
   stresspat[2].nseg = 8;
   stresspat[2].nval = 2;
   stresspat[2].vel[0] = 110;
   stresspat[2].vel[1] = 90;
   stresspat[2].expcoef[0] = (float) 1.4;
   stresspat[2].expcoef[1] = (float) 0.6;

   stresspat[3].name = "Hornpipe";
   stresspat[3].meter = "9/4";
   stresspat[3].nseg = 9;
   stresspat[3].nval = 3;
   stresspat[3].vel[0] = 110;
   stresspat[3].vel[1] = 80;
   stresspat[3].vel[2] = 100;
   stresspat[3].expcoef[0] = (float) 1.4;
   stresspat[3].expcoef[1] = (float) 0.6;
   stresspat[3].expcoef[2] = (float) 1.0;

   stresspat[4].name = "Hornpipe";
   stresspat[4].meter = "3/2";
   stresspat[4].nseg = 12;
   stresspat[4].nval = 2;
   stresspat[4].vel[0] = 110;
   stresspat[4].vel[1] =  90;
   stresspat[4].expcoef[0] = (float) 1.4;
   stresspat[4].expcoef[1] = (float) 0.6;

   stresspat[5].name = "Hornpipe";
   stresspat[5].meter = "12/8";
   stresspat[5].nseg = 12;
   stresspat[5].nval = 3;
   stresspat[5].vel[0] = 110;
   stresspat[5].vel[1] = 80;
   stresspat[5].vel[2] = 100;
   stresspat[5].expcoef[0] = (float) 1.4;
   stresspat[5].expcoef[1] = (float) 0.6;
   stresspat[5].expcoef[2] = (float) 1.0;

   stresspat[6].name = "Double hornpipe";
   stresspat[6].meter = "6/2";
   stresspat[6].nseg = 12;
   stresspat[6].nval = 6;
   stresspat[6].vel[0] = 110;
   stresspat[6].vel[1] = 80;
   stresspat[6].vel[2] = 80;
   stresspat[6].vel[3] = 105;
   stresspat[6].vel[4] = 80;
   stresspat[6].vel[5] = 80;
   stresspat[6].expcoef[0] = (float) 1.2;
   stresspat[6].expcoef[1] = (float) 0.9;
   stresspat[6].expcoef[2] = (float) 0.9;
   stresspat[6].expcoef[3] = (float) 1.2;
   stresspat[6].expcoef[4] = (float) 0.9;
   stresspat[6].expcoef[5] = (float) 0.9;

   stresspat[7].name = "Reel";
   stresspat[7].meter = "4/4";
   stresspat[7].nseg = 8;
   stresspat[7].nval = 2;
   stresspat[7].vel[0] = 120;
   stresspat[7].vel[1] = 60;
   stresspat[7].expcoef[0] = (float) 1.1;
   stresspat[7].expcoef[1] = (float) 0.9;

   stresspat[8].name = "Reel";
   stresspat[8].meter = "C|";
   stresspat[8].nseg = 8;
   stresspat[8].nval = 8;
   stresspat[8].vel[0] = 80;
   stresspat[8].vel[1] = 60;
   stresspat[8].vel[2] = 120;
   stresspat[8].vel[3] = 60;
   stresspat[8].vel[4] = 80;
   stresspat[8].vel[5] = 60;
   stresspat[8].vel[6] = 120;
   stresspat[8].vel[7] = 60;
   stresspat[8].expcoef[0] = (float) 1.1;
   stresspat[8].expcoef[1] = (float) 0.9;
   stresspat[8].expcoef[2] = (float) 1.1;
   stresspat[8].expcoef[3] = (float) 0.9;
   stresspat[8].expcoef[4] = (float) 1.1;
   stresspat[8].expcoef[5] = (float) 0.9;
   stresspat[8].expcoef[6] = (float) 1.1;
   stresspat[8].expcoef[7] = (float) 0.9;

   stresspat[9].name = "Reel";
   stresspat[9].meter = "2/4";
   stresspat[9].nseg = 8;
   stresspat[9].nval = 2;
   stresspat[9].vel[0] = 120;
   stresspat[9].vel[1] = 60;
   stresspat[9].expcoef[0] = (float) 1.1;
   stresspat[9].expcoef[1] = (float) 0.9;

   stresspat[10].name = "Slip Jig";
   stresspat[10].meter = "9/8";
   stresspat[10].nseg = 9;
   stresspat[10].nval = 3;
   stresspat[10].vel[0] = 110;
   stresspat[10].vel[1] = 70;
   stresspat[10].vel[2] = 80;
   stresspat[10].expcoef[0] = (float) 1.4;
   stresspat[10].expcoef[1] = (float) 0.5;
   stresspat[10].expcoef[2] = (float) 1.1;

   stresspat[11].name = "Double Jig";
   stresspat[11].meter = "6/8";
   stresspat[11].nseg = 6;
   stresspat[11].nval = 3;
   stresspat[11].vel[0] = 110;
   stresspat[11].vel[1] = 70;
   stresspat[11].vel[2] = 80;
   stresspat[11].expcoef[0] = (float) 1.2;
   stresspat[11].expcoef[1] = (float) 0.7;
   stresspat[11].expcoef[2] = (float) 1.1;

   stresspat[12].name = "Single Jig";
   stresspat[12].meter = "6/8";
   stresspat[12].nseg = 6;
   stresspat[12].nval = 3;
   stresspat[12].vel[0] = 110;
   stresspat[12].vel[1] = 80;
   stresspat[12].vel[2] = 90;
   stresspat[12].expcoef[0] = (float) 1.2;
   stresspat[12].expcoef[1] = (float) 0.7;
   stresspat[12].expcoef[2] = (float) 1.1;

   stresspat[13].name = "Slide";
   stresspat[13].meter = "6/8";
   stresspat[13].nseg = 6;
   stresspat[13].nval = 3;
   stresspat[13].vel[0] = 110;
   stresspat[13].vel[1] = 80;
   stresspat[13].vel[2] = 90;
   stresspat[13].expcoef[0] = (float) 1.3;
   stresspat[13].expcoef[1] = (float) 0.8;
   stresspat[13].expcoef[2] = (float) 0.9;

   stresspat[14].name = "Jig";
   stresspat[14].meter = "6/8";
   stresspat[14].nseg = 6;
   stresspat[14].nval = 3;
   stresspat[14].vel[0] = 110;
   stresspat[14].vel[1] = 80;
   stresspat[14].vel[2] = 90;
   stresspat[14].expcoef[0] = (float) 1.2;
   stresspat[14].expcoef[1] = (float) 0.7;
   stresspat[14].expcoef[2] = (float) 1.1;

   stresspat[15].name = "Ragtime";
   stresspat[15].meter = "12/8";
   stresspat[15].nseg = 12;
   stresspat[15].nval = 3;
   stresspat[15].vel[0] = 110;
   stresspat[15].vel[1] = 70;
   stresspat[15].vel[2] = 90;
   stresspat[15].expcoef[0] = (float) 1.4;
   stresspat[15].expcoef[1] = (float) 0.6;
   stresspat[15].expcoef[2] = (float) 1.0;

   stresspat[16].name = "Strathspey";
   stresspat[16].meter = "C";
   stresspat[16].nseg = 8;
   stresspat[16].nval = 2;
   stresspat[16].vel[0] = 120;
   stresspat[16].vel[1] = 80;
   stresspat[16].expcoef[0] = (float) 1.0;
   stresspat[16].expcoef[1] = (float) 1.0;

   stresspat[17].name = "Fling";
   stresspat[17].meter = "C";
   stresspat[17].nseg = 8;
   stresspat[17].nval = 2;
   stresspat[17].vel[0] = 110;
   stresspat[17].vel[1] = 90;
   stresspat[17].expcoef[0] = (float) 1.4;
   stresspat[17].expcoef[1] = (float) 0.6;

   stresspat[18].name = "Set Dance";
   stresspat[18].meter = "4/4";
   stresspat[18].nseg = 8;
   stresspat[18].nval = 2;
   stresspat[18].vel[0] = 110;
   stresspat[18].vel[1] = 90;
   stresspat[18].expcoef[0] = (float) 1.4;
   stresspat[18].expcoef[1] = (float) 0.6;

   stresspat[19].name = "Set Dance";
   stresspat[19].meter = "C|";
   stresspat[19].nseg = 8;
   stresspat[19].nval = 2;
   stresspat[19].vel[0] = 110;
   stresspat[19].vel[1] =  90;
   stresspat[19].expcoef[0] = (float) 1.4;
   stresspat[19].expcoef[1] = (float) 0.6;

   stresspat[20].name = "Waltz";
   stresspat[20].meter = "3/4";
   stresspat[20].nseg = 3;
   stresspat[20].nval = 3;
   stresspat[20].vel[0] = 110;
   stresspat[20].vel[1] = 70;
   stresspat[20].vel[2] = 70;
   stresspat[20].expcoef[0] = (float) 1.04;
   stresspat[20].expcoef[1] = (float) 0.98;
   stresspat[20].expcoef[2] = (float) 0.98;

   stresspat[21].name = "Slow March";
   stresspat[21].meter = "C|";
   stresspat[21].nseg = 8;
   stresspat[21].nval = 3;
   stresspat[21].vel[0] = 115;
   stresspat[21].vel[1] =  85;
   stresspat[21].vel[2] = 100;
   stresspat[21].vel[3] =  85;
   stresspat[21].expcoef[0] = (float) 1.1;
   stresspat[21].expcoef[1] = (float) 0.9;
   stresspat[21].expcoef[2] = (float) 1.1;
   stresspat[21].expcoef[3] = (float) 0.9;

   stresspat[22].name = "Slow March";
   stresspat[22].meter = "C";
   stresspat[22].nseg = 8;
   stresspat[22].nval = 3;
   stresspat[22].vel[0] = 115;
   stresspat[22].vel[1] =  85;
   stresspat[22].vel[2] = 100;
   stresspat[22].vel[3] =  85;
   stresspat[22].expcoef[0] = (float) 1.1;
   stresspat[22].expcoef[1] = (float) 0.9;
   stresspat[22].expcoef[2] = (float) 1.1;
   stresspat[22].expcoef[3] = (float) 0.9;

   stresspat[23].name = "March";
   stresspat[23].meter = "C|";
   stresspat[23].nseg = 8;
   stresspat[23].nval = 2;
   stresspat[23].vel[0] = 115;
   stresspat[23].vel[1] =  85;
   stresspat[23].expcoef[0] = (float) 1.1;
   stresspat[23].expcoef[1] = (float) 0.9;

   stresspat[24].name = "March";
   stresspat[24].meter = "C";
   stresspat[24].nseg = 8;
   stresspat[24].nval = 4;
   stresspat[24].vel[0] = 115;
   stresspat[24].vel[1] =  85;
   stresspat[24].vel[2] = 100;
   stresspat[24].vel[3] =  85;
   stresspat[24].expcoef[0] = (float) 1.1;
   stresspat[24].expcoef[1] = (float) 0.9;
   stresspat[24].expcoef[2] = (float) 1.1;
   stresspat[24].expcoef[3] = (float) 0.9;

   stresspat[25].name = "March";
   stresspat[25].meter = "6/8";
   stresspat[25].nseg = 8;
   stresspat[25].nval = 3;
   stresspat[25].vel[0] = 110;
   stresspat[25].vel[1] =  70;
   stresspat[25].vel[2] =  80;
   stresspat[25].expcoef[0] = (float) 1.1;
   stresspat[25].expcoef[1] = (float) 0.95;
   stresspat[25].expcoef[2] = (float) 0.95;

   stresspat[26].name = "March";
   stresspat[26].meter = "2/4";
   stresspat[26].nseg = 8;
   stresspat[26].nval = 2;
   stresspat[26].vel[0] = 115;
   stresspat[26].vel[1] =  85;
   stresspat[26].expcoef[0] = (float) 1.1;
   stresspat[26].expcoef[1] = (float) 0.9;

   stresspat[27].name = "Polka k1";
   stresspat[27].meter = "3/4";
   stresspat[27].nseg = 3;
   stresspat[27].nval = 3;
   stresspat[27].vel[0] = 90;
   stresspat[27].vel[1] = 110;
   stresspat[27].vel[2] = 100;
   stresspat[27].expcoef[0] = (float) 0.75;
   stresspat[27].expcoef[1] = (float) 1.25;
   stresspat[27].expcoef[2] = (float) 1.00;

   stresspat[28].name = "Polka";
   stresspat[28].meter = "4/4";
   stresspat[28].nseg = 8;
   stresspat[28].nval = 2;
   stresspat[28].vel[0] = 110;
   stresspat[28].vel[1] = 90;
   stresspat[28].expcoef[0] = (float) 1.4;
   stresspat[28].expcoef[1] = (float) 0.6;

   stresspat[29].name = "saucy";
   stresspat[29].meter = "3/4";
   stresspat[29].nseg = 6;
   stresspat[29].nval = 6;
   stresspat[29].vel[0] = 115;
   stresspat[29].vel[1] =  85;
   stresspat[29].vel[2] = 120;
   stresspat[29].vel[3] =  85;
   stresspat[29].vel[4] = 115;
   stresspat[29].vel[5] =  85;
   stresspat[29].expcoef[0] = (float) 1.2;
   stresspat[29].expcoef[1] = (float) 0.8;
   stresspat[29].expcoef[2] = (float) 1.3;
   stresspat[29].expcoef[3] = (float) 0.7;
   stresspat[29].expcoef[4] = (float) 1.1;
   stresspat[29].expcoef[5] = (float) 0.9;

   stresspat[30].name = "Slip jig";
   stresspat[30].meter = "3/4";
   stresspat[30].nseg = 9;
   stresspat[30].nval = 3;
   stresspat[30].vel[0] = 110;
   stresspat[30].vel[1] =  80;
   stresspat[30].vel[2] = 100;
   stresspat[30].expcoef[0] = (float) 1.3;
   stresspat[30].expcoef[1] = (float) 0.7;
   stresspat[30].expcoef[2] = (float) 1.0;

   stresspat[31].name = "Tango";
   stresspat[31].meter = "2/4";
   stresspat[31].nseg = 8;
   stresspat[31].nval = 4;
   stresspat[31].vel[4] = 110;
   stresspat[31].vel[4] =  90;
   stresspat[31].vel[4] =  90;
   stresspat[31].vel[4] = 100;
   stresspat[31].expcoef[4] = (float) 1.6;
   stresspat[31].expcoef[4] = (float) 0.8;
   stresspat[31].expcoef[4] = (float) 0.8;
   stresspat[31].expcoef[4] = (float) 0.8;

  }



int i;

void load_stress_parameters (char * rhythmdesignator)
{
int n,i,index,nval;
int qnotenum,qnoteden;
maxdur = 0;
for (n=0;n<32;n++) {fdur[n]= 0.0; ngain[n] = 0;}
fdursum[0]=fdur[0];
index = -1;
if (strlen(rhythmdesignator) < 2) {
        beatmodel = 0;
        return;
        }
for (i=0;i<32;i++) {
  if (casecmp(rhythmdesignator,stresspat[i].name) == 0 &&
      casecmp(timesigstring,stresspat[i].meter) == 0) {
     index = i;
     break;
     }
   }
if (index == -1) {printf("**warning** rhythm designator %s %s is not one of\n",
   rhythmdesignator,timesigstring);
       for (i = 0;i<32;i++) {printf("%s %s ",stresspat[i].name,stresspat[i].meter);
                         if (i % 5 == 4) printf("\n");
                         }
       printf("\n");
       beatmodel = 0;
       return;
       }

if (stressmodel == 0) beatmodel = 2;
else beatmodel = stressmodel;

nseg = stresspat[index].nseg;
nval = stresspat[index].nval;

segnum = time_num;
segden = time_denom*nseg;
reduce(&segnum,&segden);
/* compute number of segments in quarter note */
qnotenum = segden;
qnoteden = segnum*4;
reduce(&qnotenum,&qnoteden);

for (n=0;n < nseg+1;n++) {
  i = n % nval;
  ngain[n] = stresspat[index].vel[i];
  fdur[n] = stresspat[index].expcoef[i];
  if (verbose) printf("%d %f\n",ngain[n],fdur[n]);
  maxdur = MAX(maxdur,fdur[n]);
  if (n > 0) fdursum[n] = fdursum[n-1]+fdur[n-1]*(float) qnoteden /(float) qnotenum;
/* num[],denom[] use quarter note units = 1/1, segment units are usually
   half that size, so we divide by 2.0
*/
  }
/*printf("maxdur = %f\n",maxdur);*/
}