for(i=0; i<length; i++) for(j=0; j<width; j++) { num += powf(image_mod[i][j] - g_mod[i][j] ,2); denum += powf(image_mod[i][j] - f_real[i][j] ,2); } return 10*log10f(num/denum);