#include<stdio.h>
#include<math.h>
#define MAX_INSTRUMENTS 20
#define MAX_MEASURES 1000
typedef struct{
	int ticker[64];
	int ender[64];
	int pitches[64];
}melody;
	


typedef struct{
	int total_thismel; /*number of times it will play current melody*/
	int count_thismel; /*number of times it HAS played it*/
	int octave;
	int currfrag;
	melody currentmelody;
}instrument;

int randpitch(int octave){
	int answer=(octave+1)*12;
	int rando=rand()%8;
	switch(rando){
	case 0:return answer;/*do*/
	case 1:return answer+2;/*re*/
	case 2:return answer+4;
	case 3:return answer+5;
	case 4:return answer+7;
	case 5:return answer+9;
	case 6:return answer+11;
	case 7:return answer+12;/*do*/
	}
}
		
		
melody makemelody(int tickpermeasure, int sloppiness, int density, int octave){
	melody currentmelody;
	int i;int tickercount=0;int endercount=0;
	for(i=0;i<tickpermeasure;i++){
		if (rand()%100<density){
			currentmelody.ticker[i]=1;
			tickercount++;
		}
		else currentmelody.ticker[i]=0;
	}
	     printf("dens=%d", density);
	for(i=0;i<tickpermeasure-1;i++){
		if (rand()%100<density){
			currentmelody.ender[i]=1;
			endercount++;
		}
		else currentmelody.ender[i]=0;
	}	
	currentmelody.ender[tickpermeasure-1]=1;

	for(i=0;i<tickpermeasure;i++){
		currentmelody.pitches[i]=randpitch(octave);
	}
	

	printf("\nmelody made:");
	for(i=0;i<tickpermeasure;i++){
		printf("%d ", currentmelody.ticker[i]);
	}
	printf("\nender  made:");
	for(i=0;i<tickpermeasure;i++){
		printf("%d ", currentmelody.ender[i]);
	}
	printf("\n    pitch:");
	for(i=0;i<tickpermeasure;i++){
		printf("%d ",currentmelody.pitches[i]);
	}

	return currentmelody;
}

/*takes number, returns rand within 0 and that number with bell curve distribution*/
int bellcurve(int input){
     return rand()%input+1;
}

melody makepulse(int tickpermeasure, int sloppiness, int octave, int resolution){
	melody currentmelody;
	int i;int pulseeveryx;
	switch(resolution){
	case 1: pulseeveryx=1;break;
	case 2: pulseeveryx=1;break;
	case 3: pulseeveryx=1;break;
	case 4: pulseeveryx=2;break;
	case 5: pulseeveryx=4;break;
	case 6: pulseeveryx=8;break;
	}


	for(i=0;i<tickpermeasure;i+=pulseeveryx){
		currentmelody.ticker[i]=1;
		
	}
	for(i=pulseeveryx-1;i<tickpermeasure;i+=pulseeveryx){
		
			currentmelody.ender[i]=1;
	}
	for(i=0;i<8;i++){
		currentmelody.pitches[i]=24;
	}			
	printf("\nmelody made:");
	for(i=0;i<tickpermeasure;i++){
		printf("%d ", currentmelody.ticker[i]);
	}
	printf("\nender  made:");
	for(i=0;i<tickpermeasure;i++){
		printf("%d ", currentmelody.ender[i]);
	}
	return currentmelody;
}

int printSKINI(instrument inst[],int instcount, float tick,int sloppiness, int tickpermeasure, FILE* output,int measure){
	int i,j=0;int last[20]={0};int pitchcount[20]={0};int justended=0;float zero=0;
	for(i=0;i<tickpermeasure;i++){
		
		
		
		for(j=0;j<instcount;j++){
			if(inst[j].currentmelody.ticker[i]==1){
				fprintf(output,"\nNoteOn %f %d %d %d",zero,j+1,inst[j].currentmelody.pitches[pitchcount[j]]+(inst[j].octave*12),100);
			}
		}
		justended=0;		
		for(j=0;j<instcount;j++){
			if(inst[j].currentmelody.ender[i]==1){
				if(justended==1){
					fprintf(output,"\nNoteOff %f %d %d %d",zero,j+1,inst[j].currentmelody.pitches[pitchcount[j]++]+(inst[j].octave*12),100);
				}
				else fprintf(output,"\nNoteOff %f %d %d %d",tick-.01,j+1,inst[j].currentmelody.pitches[pitchcount[j]++]+(inst[j].octave*12),100);
				justended=1;
			}
		}
		if(justended==0)fprintf(output,"\nNoteOff %f %d %d %d",tick,instcount+1,0,0);
		else fprintf(output,"\nNoteOff %f %d %d %d",.01,instcount+1,0,0);
	}
	return 0;
}


int main(int argc, char *argv[]){
	int resolution=4;/*6 or less.2=quart 3=eigth*/
	int tempo=bellcurve(90)+45;
	melody themelodies[MAX_MEASURES];
	
	int duration=100; /*measures*/
	int tickpermeasure=pow(2,resolution);/*64 or less*/
	int sloppiness=0;
	int density=40;
	instrument inst[MAX_INSTRUMENTS];
	int i,j,instcount=rand()%6+5;
	float tick;
	FILE *output;
/*	int instmeasurecount[10]={0};*/

	if(argc<2){
     printf("usage: riley output.ski flags\n");
     printf("       -t tempo: in bpm, default random bellcurve around 90\n");
     printf("       -n number of instruments: value 1(boring) to 20. default random 5-10\n");
     printf("       -r resolution: subdivision of beat, value 1-6: 1=halfnotes, 2=quarter,3=eighth,etc. default 4\n");
     printf("       -d duration, in 4/4 measures. default 100\n");
		exit(0);
	}
     for(i=2;i<argc;i++){
		if(argv[i][1]=='t'){
			tempo=atoi(argv[i+1]);
			printf("input:tempo=%d", tempo);
		}
		if(argv[i][1]=='n'){
		    instcount=atoi(argv[i+1]);
			if (instcount>MAX_INSTRUMENTS){
			 printf("\nI'm sorry, please keep instrument count under %d",MAX_INSTRUMENTS);exit(0);
		    }
		    printf("input:instcount=%d", instcount);
	  
		  }
		if(argv[i][1]=='r'){
		    resolution=atoi(argv[i+1]);
			tickpermeasure=pow(2,resolution);
		    printf("input:res=%d", resolution);
		}
		if(argv[i][1]=='d'){
		    duration=atoi(argv[i+1]);
		    if (duration>MAX_MEASURES){
			 printf("\nI'm sorry, please keep duration under %d measures",MAX_MEASURES);exit(0);
		    }
		    printf("input:dur=%d", duration);
		}
	 }
	output=fopen(argv[1],"w");

	fprintf(output,"//Skini file created by RILEY.C by Daniel Iglesia, 5/10/03\n");
	fprintf(output,"//original parameters:\n");
	fprintf(output,"//   number of instruments:%d\n",instcount);
	fprintf(output,"//   tempo in bpm:%d\n",tempo);
	fprintf(output,"//   resolution (or, minimum note value= 1/(2^resolution):%d\n",resolution);
	fprintf(output,"//   duration:%d measures\n",duration);
	






	tick=(float)240/(tempo*tickpermeasure);

	printf("tick=%f ", tick);

	inst[0].total_thismel=duration;
	inst[0].count_thismel=0;
	inst[0].octave=3;
	inst[0].currentmelody=makepulse(tickpermeasure, sloppiness, inst[0].octave,resolution);

	for(i=0;i<MAX_MEASURES;i++){
		themelodies[i]=makemelody(tickpermeasure, sloppiness, density, 0);
	}
	
	for(i=1;i<instcount;i++){
		inst[i].total_thismel=(bellcurve(6)+1)*i;
		inst[i].count_thismel=0;
		inst[i].octave=bellcurve(6);
		inst[i].currfrag=0;
	
	}

	for(i=0;i<duration;i++){
		printSKINI(inst, instcount,tick, sloppiness, tickpermeasure, output,i);
		
		for(j=0;j<instcount;j++){
			inst[j].count_thismel++;
			if(inst[j].total_thismel==inst[j].count_thismel){/*done with current melody*/
				inst[j].currentmelody=themelodies[inst[j].currfrag++];
				inst[j].total_thismel=bellcurve(10)+1;
				inst[j].count_thismel=0;
				printf("inst %d moved to new at measure %d...will do it for %d", j, i, inst[j].total_thismel);
			}
		}
	}

		
	


	return 0;
}