#include<stdio.h>
#include<stddef.h>
#include <string.h>
#include <stdlib.h>


// simple reader for identification problem format
// Dag Wedelin 2008


int
main(int argc,char *argv[]) 
{
	char X[200],fil1[30];
	FILE *filp;
	char outfile[50];
	int i;
	
	if(argc==1) {
        printf("Please give the name of your file!\n");
        scanf("%s",fil1);
		readFile(fil1);
	}
	else {
		for(i=1;i<argc;i++) readFile(argv[i]);
	}
}


int
endOfLine(FILE *infile) // checks if we have eol in front of us
{
	char c;
	while((c=getc(infile))==' ');  // eat blanks
	ungetc(c,infile); // put back non-blank character
	if(c=='\n') return 1;
	else return 0;
}


int
readFile(char *filename) /* this routine does not change */
{
	FILE *infile; int lines;
	char rowbuffer[1000];   /* not so good really */
	int rowlength;
	char keyToken[100], *keyword, *indexString;
	int accept=0,i,keyIndex,end=0;
	char c;
	char formatVersion[100], type[100],date[100],url[100],rName[100],lvName[100],name[100];
	char lvName2[100],lpName[100];
	char property[100],reactionType[100],varSetX1[100],rangeK1[100],varSetX2[100],object[100];
	char errorType[100],equation[100];
	int setIndex,varIndex,index,lvIndex,lvIndex2,lpIndex;
	double lb,ub,lambda,time,value;
	
	char problemName[40];
	int l,w;
	
	infile=fopen(filename,"r");
	if(!infile) {printf("Error: file ""%s"" does not exist!\n",filename); return;}

	while(1) {
		
		// find the first identifier = keyToken
		
		keyToken[0]='\0';
		if(fscanf(infile,"%s",&keyToken)==EOF) break; // find keyToken
		
		// split keyToken into keyword_indexString
		
		keyword=keyToken;
		indexString=keyToken;
		while(*indexString != '\0') {
			if(*indexString == '_') {*indexString='\0'; indexString++; break;}
			indexString++;
		}
		sscanf(indexString,"%d",&keyIndex); // translate indexString to integer

		// now check the keyword and read accordingly
		// in this code we choose to reprint the problem to stdout, this can be replaced as appropriate
			
		if(!strcmp(keyword,"begin")) {
			fscanf(infile," problem %s",problemName);
			printf("begin problem %s\n",problemName);
			end=0;
		}
		else if(!strcmp(keyword,"end")) {
			fscanf(infile," problem %s",problemName);
			printf("end problem %s\n",problemName);
			printf("\n\n\n"); // nice with blank lines if files are concatenated
			end=1;
			//break;
		}
		else if(!strcmp(keyword,"format")) {
			fscanf(infile,"format version %s",formatVersion);
			printf("format version 1.0\n");
		}
		else if(!strcmp(keyword,"type")) {
			fscanf(infile," = %s",type);
			printf("type = %s\n",type);
		}
		else if(!strcmp(keyword,"date")) {
			fscanf(infile," = %s",date);
			printf("date = %s\n",date);
		}
		else if(!strcmp(keyword,"url")) {
			fscanf(infile," = %s",url);
			printf("url = %s\n",url);
		}
		else if(!strcmp(keyword,"date")) {
			fscanf(infile," = %s",date);
			printf("date = %s\n",date);
		}
		else if(!strcmp(keyword,"reaction")) {   // borde detta vara med alls?
			fscanf(infile," has name = %s",rName);
			printf("reaction_%d\n",keyIndex);
			printf("has name = %s\n",rName);
			if(!strcmp(rName,"uniMolecularMassAction")) {
				fscanf(infile," has localVariableName_%d = %s",&lvIndex,lvName);
				fscanf(infile," has localParameterName_%d = %s",&lpIndex,lpName);
				fscanf(infile," has equation = %s",equation);
				printf("has localVariableName_%d = %s\n",lvIndex,lvName);
				printf("has localParameterName_%d = %s\n",lpIndex,lpName);
				printf("has equation = %s\n",equation);
			}
			else if(!strcmp(rName,"biMolecularMassAction")) {
				fscanf(infile," has localVariableName_%d = %s",&lvIndex,lvName);
				fscanf(infile," has localVariableName_%d = %s",&lvIndex2,lvName2);
				fscanf(infile," has localParameterName_%d = %s",&lpIndex,lpName);
				fscanf(infile," has equation = %s",equation);
				printf("has localVariableName_%d = %s\n",lvIndex,lvName);
				printf("has localVariableName_%d = %s\n",lvIndex2,lvName2);
				printf("has localParameterName_%d = %s\n",lpIndex,lpName);
				printf("has equation = %s\n",equation);
			}
		}
		else if(!strcmp(keyword,"variable")) {
			fscanf(infile," has name = %s is %s",name,property);
			printf("variable_%d has name = %s is %s\n",keyIndex,name,property);
		}
		else if(!strcmp(keyword,"range")) {
			fscanf(infile," has lowerBound = %lf has upperBound = %lf",&lb,&ub);
			printf("range_%d has lowerBound = %lf has upperBound = %lf\n",keyIndex,lb,ub);
		}
		else if(!strcmp(keyword,"memberOfSet")) {
			printf("memberOfSet_%d",keyIndex);
			while(1) {
				fscanf(infile," variable_%d",&varIndex);
				printf(" variable_%d",varIndex);
				if(endOfLine(infile)) break;
			}
			printf("\n");
		}
		else if(!strcmp(keyword,"possibleReaction")) {
			fscanf(infile," of variable_%d",&varIndex);
			fscanf(infile," has type = %s",reactionType);
			printf("possibleReaction_%d of variable_%d\n",keyIndex,varIndex);
			printf("has type = %s\n",reactionType);
			if(!strcmp(reactionType,"uniMolecularMassAction")) {
				fscanf(infile," has spaceOfVariable X1 = %s",varSetX1);
				fscanf(infile," has rangeOfParameter k1 = %s",rangeK1);
				
				printf("has spaceOfVariable X1 = %s\n",varSetX1);
				printf("has rangeOfParameter k1 = %s\n",rangeK1);
				//printf("\n");
			}
			else if(!strcmp(reactionType,"biMolecularMassAction")) {
				fscanf(infile," has spaceOfVariable X1 = %s",varSetX1);
				fscanf(infile," has spaceOfVariable X2 = %s",varSetX2);
				fscanf(infile," has rangeOfParameter k1 = %s",rangeK1);
				
				printf("has spaceOfVariable X1 = %s\n",varSetX1);
				printf("has spaceOfVariable X2 = %s\n",varSetX2);
				printf("has rangeOfParameter k1 = %s\n",rangeK1);
				//printf("\n");
			}
			//else discard rest of row...
		}
		else if(!strcmp(keyword,"errorFunction")) {
			fscanf(infile," has type = %s",errorType);
			if(!strcmp(errorType,"minusLogLikelihoodPlusLambdaK")) {
				fscanf(infile," has equation = %s",equation);
				fscanf(infile," has lambda = %lf",&lambda);
				printf("errorFunction\n");
				printf("has type = %s\n",errorType);
				printf("has equation = %s\n",equation);
				printf("has lambda = %lf\n",lambda);
				//printf("\n");
			}
			//else discard rest of row...
		}
		else if(!strcmp(keyword,"experiment")) {
			fscanf(infile," has name = %s has %s",name,object);
			printf("experiment_%d has name = %s has %s\n",keyIndex,name,object);
		}
		else if(!strcmp(keyword,"sample")) {
			fscanf(infile," of experiment_%d",&index);
			fscanf(infile," has time = %lf",&time);
			printf("sample_%d of experiment_%d\n",keyIndex,index);
			printf("has time = %lf\n",time);
			fscanf(infile," has variable_ =");
			printf("has variable_ =");
			while(1) {
				fscanf(infile," %lf",&value);
				printf(" %lf",value);
				if(endOfLine(infile)) break;
			}
			printf("\n");
			fscanf(infile," has sdev of variable_ =");
			printf("has sdev of variable_ =");
			while(1) {
				fscanf(infile," %lf",&value);
				printf(" %lf",value);
				if(endOfLine(infile)) break;
			}
			printf("\n");
		}
		
		if(!strcmp(keyword,"//")) {
			fgets(rowbuffer,999,infile);  // read rest of line
			printf("//%s",rowbuffer);
		}
		else while(getc(infile)!='\n' && !feof(infile)); // always read to eol including the eol
		
		if(endOfLine(infile) && !end) printf("\n");	// pass on a blank line after a statement
		
		/*
		if(!strcmp(keyword,"length")){
			fscanf(infile," of person = %d \n",&l);
			printf("length of person = %d\n",l);
		}
		else if(!strcmp(keyword,"weight")) {
			fscanf(infile," of person = %d \n",&w);
			printf("weight of person = %d\n",w);
		}
		else if(!strcmp(keyword,"person")) {
			fscanf(infile," is %s \n",property);
			printf("person_%d is %s\n",keyIndex,property);
		}
		 */


	}
	fclose(infile);
}