points-polygons-x

#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;
}