☆display_BMP_averave完成2_malloc 追加パフォーマンスチェック☆

#include <stdio.h>
#include <stdlib.h>
#include <GLUT/glut.h>
#include <time.h>

/* keep the picture data*/
#define  IMGFILE  "depth_combine48_cunbus.bmp"
#define  CHANNEL  3  // color(R, G, B)
FILE *fp;
int i,j;
int width;           // width of picture
int height;          // height of picture
float *file1,*file2,*file3;
/*
 float get_wx ;
 float get_wy ;
 float get_ww ;
 float get_wh ;
 */
//float Xshift,Yshift;

GLuint texture;      // texture number
GLubyte* image;      // keep the BMP data

void loadImage(const char* filepath);  // read BMP

float shiftX,shiftY;
float mouseX,mouseY;
float x,y;
//float a;
float g_pressed = 0;

//float *file1;
int WIDTH,HEIGHT;
int WIDTH1,HEIGHT1;
float dn1,dn2,dn3;
int b;
float total1,total2,total3;
float  average1,average2,average3;
float start, end;

void display(void)
{
    /*diaplay the picture */
    glEnable(GL_TEXTURE_2D);
    glClear(GL_COLOR_BUFFER_BIT);
    glBindTexture(GL_TEXTURE_2D, texture);

    glTranslated(-0.5,-0.5,0);

    /*glBegin(GL_QUADS);
     glTexCoord2f(0 , 0); glVertex2f(-0.9 , -0.9);
     glTexCoord2f(0 , 1); glVertex2f(-0.9 , 0.9);
     glTexCoord2f(1 , 1); glVertex2f(0.9 , 0.9);
     glTexCoord2f(1 , 0); glVertex2f(0.9 , -0.9);
     glEnd();
     */
    glBegin(GL_QUADS);

    glTexCoord2f(1.0, 0.0);  glVertex2f(1.0, 0.0);//position
    glTexCoord2f(1.0, 1.0);  glVertex2f(1.0, 1.0);
    glTexCoord2f(0.0, 1.0);  glVertex2f(0.0, 1.0);
    glTexCoord2f(0.0, 0.0);  glVertex2f(0.0, 0.0);

    glDisable(GL_TEXTURE_2D);

    glEnd();


    //glClear(GL_COLOR_BUFFER_BIT);


    //glColor3d(0.0, 0.0, 0.0);//color


    glTranslated(0.5,0.5,0);
    //glTranslatef(shiftX,shiftY,0);

    x=0.1;
    y=0.1;

    glBegin(GL_LINE_LOOP);
    /*
     get_wx=glutGet( GLUT_WINDOW_X );
     get_wx=glutGet( GLUT_WINDOW_Y );
     get_wx=glutGet( GLUT_WINDOW_WIDTH );
     get_wx=glutGet( GLUT_WINDOW_HEIGHT );
     */
    glVertex2d(-x+shiftX, -y+shiftY);//左下
    glVertex2d(x+shiftX, -y+shiftY);//右下
    glVertex2d(x+shiftX, y+shiftY);//右上
    glVertex2d(-x+shiftX, y+shiftY);//カーソルの左上
    printf("ShiftX1 = %f : ShiftY1 = %f\n" , -x+shiftX , -y+shiftY);//左下
    printf("ShiftX2 = %f : ShiftY2 = %f\n" , x+shiftX , -y+shiftY);//右下
    printf("ShiftX3 = %f : ShiftY3 = %f\n" , x+shiftX , y+shiftY);//右上
    printf("ShiftX4 = %f : ShiftY4 = %f\n" , -x+shiftX , y+shiftY);//左上
    printf("b=%d\n",b);

    printf("HEIGHT = %d : WIDTH = %d : HEIGHT1 = %d : WIDTH1 = %d : y+shiftY = %f : -x+shiftX = %f : b=%d  \n" , HEIGHT,WIDTH,HEIGHT1,WIDTH1,y+shiftY,-x+shiftX,b);

    total1=0;
    total2=0;
    total3=0;


    for(i= -(int)((float)HEIGHT*(y+shiftY))+b ;i<(int)(-(float)HEIGHT*(y+shiftY))+b+HEIGHT1;i++)
    {
        //printf("i=%d\n",i);
        //printf("i2=%d\n",-HEIGHT*(int)(y+shiftY)+b);
        for(j=(int)((float)WIDTH*(-x+shiftX))+b;j<(int)((float)WIDTH*(-x+shiftX))+b+WIDTH1;j++)
        {
            //   printf("j=%d\n",j);
            //  printf("check-7\n");
            // printf("position=%d\n",(i-(WIDTH-HEIGHT)/2)*WIDTH+j);
            //printf("position=%d¥n",(i-WIDTH1)*HEIGHT1+(j-HEIGHT1+((WIDTH-HEIGHT)/2)));
            dn1=file1[(i-(WIDTH-HEIGHT)/2)*WIDTH+j];
            total1=total1+dn1;
            dn2=file2[(i-(WIDTH-HEIGHT)/2)*WIDTH+j];
            total2=total2+dn2;
            dn3=file3[(i-(WIDTH-HEIGHT)/2)*WIDTH+j];
            total3=total3+dn3;

            // printf("check-8\n");
        }
        //printf("check-9\n");
    }

    printf("total1=%f\n",total1);
    printf("total2=%f\n",total2);
    printf("total3=%f\n",total3);
    average1=total1/(float)(WIDTH1*HEIGHT1);
    average2=total2/(float)(WIDTH1*HEIGHT1);
    average3=total3/(float)(WIDTH1*HEIGHT1);

    printf("the average water depth is %f m\n " ,average1 );
    if(average1>0 && average1<400.0)
    {
        printf("the averave depth water is so deep\n");
    }
    else if(average1>=400.0 && average1<1000.0)
    {
        printf("the averave depth water is too deep\n");
    }
    else if(average1>1000.0)
    {
        printf("the averave depth water is fucking deep \n");
        }

    printf("the average water depth is %f m\n " ,average2 );
    if(average2>0 && average2<400.0)
    {
        printf("the averave depth water is so deep\n");
    }
    else if(average2>=400.0 && average2<1000.0)
    {
        printf("the averave depth water is too deep\n");
    }
    else if(average2>1000.0)
    {
        printf("the averave depth water is fucking deep \n");
    }
    printf("the average water depth is %f m\n " ,average3 );
    if(average3>0 && average3<400.0)
    {
        printf("the averave depth water is so deep\n");
    }
    else if(average3>=400.0 && average3<1000.0)
    {
        printf("the averave depth water is too deep\n");
    }
    else if(average3>1000.0)
    {
        printf("the averave depth water is fucking deep \n");
    }
    end=clock();

    printf("start  %f /sec \n " ,start/1000000);
    printf("end  %f /sec \n " ,end/1000000);
    printf("star to end  %f /sec \n " ,(end-start)/1000000);
    //get_wx=glutGet( GLUT_WINDOW_X );
    //get_wy=glutGet( GLUT_WINDOW_Y );
    //get_ww=glutGet( GLUT_WINDOW_WIDTH );
    //get_wh=glutGet( GLUT_WINDOW_HEIGHT );

    // Xshift=x+shiftX;
    // Yshift=y+shiftY;



    //glPushMatrix();
    //glTranslatef(shiftX,shiftY,0);
    // glTranslatef(0.5,0.5,0);
    //glPopMatrix();
    glEnd();


    glutSwapBuffers();

    //glFlush();

}



void myInit(void)
{

    /* necessary setting for before create window*/
    //window1//
    glutInitWindowSize(600, 600);                  // setting for window size
    glutInitWindowPosition(0, 0);              // window position
    glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE);  // setting display mode
    glutCreateWindow("OpenGL");// create window

    glClearColor (0.0, 0.0, 0.0, 1.0);             // window back color

    glutDisplayFunc(display);
    //create texture
    loadImage(IMGFILE);
    glPixelStorei(GL_UNPACK_ALIGNMENT, 8);
    glGenTextures(1, &texture);
    glBindTexture(GL_TEXTURE_2D, texture);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
    glTexImage2D(GL_TEXTURE_2D, 0, GL_RGB,
                 width, height, 0, GL_RGB,
                 GL_UNSIGNED_BYTE, image);


    //gluBuild2DMipmaps(GL_TEXTURE_2D, GL_RGB, WIDTH, TEXHEIGHT,
    //        GL_RGB, GL_UNSIGNED_BYTE, texture);

    free(image);
}

void keyboard(unsigned char key, int x, int y)  //keyboard setting
{

    switch (key)
    {
        case 'q':
        case 'Q':
        case '\033':
            exit(0);
        default:
            break;

    }
}


void mouse(int btn, int state, int x, int y)  //mouze setting
{


    if(btn == GLUT_LEFT_BUTTON == state == GLUT_DOWN)  //when press the left button
    {
        //Xshift=x+shiftX;
        // Yshift=y+shiftY;
        mouseX = x;
        mouseY = y;
        g_pressed = 1;

    }
    else if(btn == GLUT_LEFT_BUTTON == state ==GLUT_UP)  //when take off the left button
        g_pressed = 0;

}


void motion(int x, int y)  //controll the motion by the mouse
{
    {

        shiftX = (float)((x - mouseX)/320);
        shiftY = (float)(-(y - mouseY)/320);
        start =clock();
        glutPostRedisplay();

    }

    /* printf("X = %d : Y = %d\n" , x , y);  //diplay the position of coordinate
     printf("ShiftX = %f : ShiftY = %f\n" , shiftX , shiftY);  //display the amount of the shift
     printf("mouseX = %f : mouseY = %f \n" , -x+shiftX , mouseY);
     //printf("mouseX = %f : mouseY = %f \n" , mouseX , mouseY);  //display the mouseX and mouseY
     printf("wx = %f : wy = %f \n" , get_wx , get_wy);
     // printf("windowX = %f : windowY = %f \n" , get_wx , get_wy);
     printf("windowW = %f : windowH = %f \n" , get_ww , get_wh);
     */
}

int main(int argc, char* argv[])
{
    //float *file1;
    printf("check-1\n");
     //float start, end;

    if(argc!=4)
    {
        fprintf(stderr,"Usage: %s\n(1)input file name\n(2)WIDTH\n(3)HEIGHT\n",argv[0]);
        exit(1);
    }
    //printf("check-2\n");
    WIDTH = atoi(argv[2]); printf("W = %d,", WIDTH);
    HEIGHT = atoi(argv[3]); printf("H = %d\n", HEIGHT);
    WIDTH1=WIDTH/5;
    HEIGHT1=WIDTH/5;
    b=WIDTH/2;
    // printf("check-3\n");
    //float *file1;
    file1= (float *)malloc(WIDTH*HEIGHT*sizeof(float));
    file2= (float *)malloc(WIDTH*HEIGHT*sizeof(float));
    file3= (float *)malloc(WIDTH*HEIGHT*sizeof(float));
    // printf("check-4\n");
    if((fp=fopen(argv[1],"rb"))==NULL)
    {
        fprintf(stderr,"input file open error\n");
        exit(1);
    }
    //printf("check-5\n");
    fread((float *)file1,sizeof(float),WIDTH*HEIGHT,fp);
    fread((float *)file2,sizeof(float),WIDTH*HEIGHT,fp);
    fread((float *)file3,sizeof(float),WIDTH*HEIGHT,fp);
    /*
     for(i=0;i<WIDTH*HEIGHT;i++)
     {
     printf("file_depth=%f\n",file1[i]);
     }
     */
    //total=0;

    //printf("HEIGHT = %d : WIDTH = %d : HEIGHT1 = %d : WIDTH1 = %d : y+shiftY = %f : -x+shiftX = %f : b=%d  \n" , HEIGHT,WIDTH,HEIGHT1,WIDTH1,y+shiftY,-x+shiftX,b);

   // printf("i=%d\n",-HEIGHT*(y+shiftY)+b);
   // printf("i_p=%d\n",-HEIGHT*(y+shiftY)+b+HEIGHT1);
    /*
    for(i= -HEIGHT*(y+shiftY)+b ;i<-HEIGHT*(y+shiftY)+b+HEIGHT1;i++)
    {
        //printf("i=%d\n",i);
        //printf("check-6\n");
        for(j=WIDTH*(-x+shiftX)+b;j<WIDTH*(-x+shiftX)+b+WIDTH1;j++)
        {
            //  printf("j=%d\n",j);
            //printf("check-7\n");
            //printf("position=%d¥n",(i-WIDTH1)*HEIGHT1+(j-HEIGHT1+((WIDTH-HEIGHT)/2)));
            //dn=file1[(i-(WIDTH1-HEIGHT1))*WIDTH1+j];
            dn=file1[(i-(WIDTH-HEIGHT)/2)*WIDTH+j];
            total=total+dn;
            // printf("check-8\n");
        }
        //printf("check-9\n");

    }
    averave=total/(WIDTH1*HEIGHT1);
    printf("the average water depth is %f m\n " ,averave );
    */
    printf("check-10\n");

    //free(file1);
    //fclose(fp);

    glutInit(&argc, argv);
    //glutInitDisplayMode(GLUT_RGBA);
    //glutCreateWindow(argv[0]);
    myInit();



    glutKeyboardFunc(keyboard);  //when need keyboard function
    glutMouseFunc(mouse);  //when need mouse function
    glutMotionFunc(motion);  //when need mouse motion function
    glutDisplayFunc(display);  //when need to set about display function
    //glutIdleFunc(idle_func);
    glutMainLoop();

    return 0;

    free(file1);
    free(file2);
    free(file3);
    fclose(fp);

}


void loadImage(const char* filepath)
{
    int i;
    int size;  // size of picture (width , height, color)
    FILE* fp;

    /* read file */
    fp = fopen(filepath, "rb");  //read binary
    if (fp == NULL)
    {
        fprintf(stderr, "%sfail loading\n", filepath);
        exit(EXIT_FAILURE);
    }

    /* get width and height  */
    fseek(fp, 18, SEEK_SET);  //move fp to place where has been put width of picture
    fread(&width, 4, 1, fp);
    fread(&height, 4, 1, fp);

    /* read data of picture */
    fseek(fp, 54, SEEK_SET);  //move fp to place where has been put dataof the picture
    size = width * height * CHANNEL;
    image = (GLubyte* ) malloc(size);
    if (image == NULL)
    {
        fprintf(stderr, "fail to keep the memory\n");
        exit(EXIT_FAILURE);
    }

    for (i = 0; i < size; ++i)
        fread(&image[i], sizeof(GLubyte), 1, fp);

    /* exit */
    fclose(fp);
}