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

int main(int argc,char *argv[])
{
	FILE *fpi1,*fpi2,*fpo;
	float *points;
	int *polygons;
	int data1;
	int data2;
	int i,j,k,l;
	int id;
	float xyz[3];

	if(argc!=6)
	{
		fprintf(stderr,"Usage: %s 5m_mesh_data.lem,5m_mesh_data.raw\n",argv[0]);
		exit(1);
	}
	printf("OPEN FILE NAME:%s\n",argv[1]);

	if((fpi1=fopen(argv[1],"rb"))==NULL)
	{
		fprintf(stderr,"Can not open\n");
		exit(1);
	}

	if((fpi2=fopen(argv[2],"rb"))==NULL)
	{
		fprintf(stderr,"Can not open\n");
		exit(1);
	}

	if((fpo=fopen(argv[3],"wb"))==NULL)
	{
		fprintf(stderr,"Can not open\n");
		exit(1);
	}

	data1=atoi(argv[4]);
	data2=atoi(argv[5]);

	points=(float *)malloc(data1*3*sizeof(float));
	fread((float *)points,sizeof(float),1,fpi1);
	polygons=(int *)malloc(data2*3*sizeof(int));
	fread((int *)polygons,sizeof(int),1,fpi2);

	for(k=0;k<data2;k++)
	{
		printf("1\n");
		for (l=0;l<3;l++)
		{
				printf("2\n");
			id=polygons[l+3*k];
				printf("%d\n",id);
			xyz[0]=points[3*id];
				printf("%f\n",xyz[0]);
			xyz[1]=points[3*id+1];
			printf("%f\n",xyz[1]);
			xyz[2]=points[3*id+2];
			printf("%f\n",xyz[2]);
			//fseek(fpi1,5*k+j*17+i*68,SEEK_SET);

			//fgets(line, 6, fpi1);
			//data=atof(line);
			//printf("%f\n",data);
			//fwrite(&x,sizeof(float),1,fpo);
			printf("7\n");
			//fwrite(&y,sizeof(float),1,fpo);
			printf("8\n");
			//fwrite(&z,sizeof(float),1,fpo);
			printf("9\n");
			printf("%f %f %f\n",&xyz[0],&xyz[1],&xyz[2]);
			fwrite(&xyz,sizeof(float),3,fpo);
		}
	}

	fclose(fpi1);
	fclose(fpi2);
	fclose(fpo);

	return 0;
}