1 /*
2  *  Copyright (c) 2012 The WebM project authors. All Rights Reserved.
3  *
4  *  Use of this source code is governed by a BSD-style license
5  *  that can be found in the LICENSE file in the root of the source
6  *  tree. An additional intellectual property rights grant can be found
7  *  in the file PATENTS.  All contributing project authors may
8  *  be found in the AUTHORS file in the root of the source tree.
9  */
10 
11 #ifndef TEST_UTIL_H_
12 #define TEST_UTIL_H_
13 
14 #include <stdio.h>
15 #include <math.h>
16 #include "third_party/googletest/src/include/gtest/gtest.h"
17 #include "vpx/vpx_image.h"
18 
19 // Macros
20 #define GET_PARAM(k) std::tr1::get< k >(GetParam())
21 
compute_psnr(const vpx_image_t * img1,const vpx_image_t * img2)22 static double compute_psnr(const vpx_image_t *img1,
23                            const vpx_image_t *img2) {
24   assert((img1->fmt == img2->fmt) &&
25          (img1->d_w == img2->d_w) &&
26          (img1->d_h == img2->d_h));
27 
28   const unsigned int width_y  = img1->d_w;
29   const unsigned int height_y = img1->d_h;
30   unsigned int i, j;
31 
32   int64_t sqrerr = 0;
33   for (i = 0; i < height_y; ++i)
34     for (j = 0; j < width_y; ++j) {
35       int64_t d = img1->planes[VPX_PLANE_Y][i * img1->stride[VPX_PLANE_Y] + j] -
36                   img2->planes[VPX_PLANE_Y][i * img2->stride[VPX_PLANE_Y] + j];
37       sqrerr += d * d;
38     }
39   double mse = static_cast<double>(sqrerr) / (width_y * height_y);
40   double psnr = 100.0;
41   if (mse > 0.0) {
42     psnr = 10 * log10(255.0 * 255.0 / mse);
43   }
44   return psnr;
45 }
46 
47 #endif  // TEST_UTIL_H_
48