iX=i*X;
        iX1=(i+1)*X;
        iX2=(i+2)*X;
        iX3=(i+3)*X;
        iX4=(i+4)*X;
        for(j=0; j < Y; j=j+5){
           valor+=HIM[iX+j+sXY]*HIM[iX+j+(mXY)]/XYZ;               valor+=HIM[iX+j+1+sXY]*HIM[iX+j+1+(mXY)]/XYZ;
                    valor+=HIM[iX+j+2+sXY]*HIM[iX+j+2+(mXY)]/XYZ;
                    valor+=HIM[iX+j+3+sXY]*HIM[iX+j+3+(mXY)]/XYZ;
                    valor+=HIM[iX+j+4+sXY]*HIM[iX+j+4+(mXY)]/XYZ;
        }
        for(j=0; j < Y; j=j+5){
            valor+=HIM[iX1+j+sXY]*HIM[iX1+j+(mXY)]/XYZ;
                    valor+=HIM[iX1+j+1+sXY]*HIM[iX1+j+1+(mXY)]/XYZ;
                    valor+=HIM[iX1+j+2+sXY]*HIM[iX1+j+2+(mXY)]/XYZ;
                    valor+=HIM[iX1+j+3+sXY]*HIM[iX1+j+3+(mXY)]/XYZ;
                    valor+=HIM[iX1+j+4+sXY]*HIM[iX1+j+4+(mXY)]/XYZ;
            }
        for(j=0; j < Y; j=j+5){
           valor+=HIM[iX2+j+sXY]*HIM[iX2+j+(mXY)]/XYZ;
                   valor+=HIM[iX2+j+1+sXY]*HIM[iX2+j+1+(mXY)]/XYZ;
                   valor+=HIM[iX2+j+2+sXY]*HIM[iX2+j+2+(mXY)]/XYZ;
                   valor+=HIM[iX2+j+3+sXY]*HIM[iX2+j+3+(mXY)]/XYZ;
                   valor+=HIM[iX2+j+4+sXY]*HIM[iX2+j+4+(mXY)]/XYZ;
        }
        for(j=0; j < Y; j=j+5){
           valor+=HIM[iX3+j+sXY]*HIM[iX3+j+(mXY)]/XYZ;
                   valor+=HIM[iX3+j+1+sXY]*HIM[iX3+j+1+(mXY)]/XYZ;
                   valor+=HIM[iX3+j+2+sXY]*HIM[iX3+j+2+(mXY)]/XYZ;
                   valor+=HIM[iX3+j+3+sXY]*HIM[iX3+j+3+(mXY)]/XYZ;
                   valor+=HIM[iX3+j+4+sXY]*HIM[iX3+j+4+(mXY)]/XYZ;
        }
        for(j=0; j < Y; j=j+5){
           valor+=HIM[iX4+j+sXY]*HIM[iX4+j+(mXY)]/XYZ;
                   valor+=HIM[iX4+j+1+sXY]*HIM[iX4+j+1+(mXY)]/XYZ;
                   valor+=HIM[iX4+j+2+sXY]*HIM[iX4+j+2+(mXY)]/XYZ;
                   valor+=HIM[iX4+j+3+sXY]*HIM[iX4+j+3+(mXY)]/XYZ;
                   valor+=HIM[iX4+j+4+sXY]*HIM[iX4+j+4+(mXY)]/XYZ;
        }