1 /*
2 * Copyright (c) 2023 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 /*
12 * Sum the square of the difference between every corresponding element of the
13 * buffers.
14 */
15
16 #include <stdlib.h>
17
18 #include "./vpx_config.h"
19 #include "./vpx_dsp_rtcd.h"
20
21 #include "vpx/vpx_integer.h"
22 #include "vpx_ports/mem.h"
23
vpx_sse_c(const uint8_t * a,int a_stride,const uint8_t * b,int b_stride,int width,int height)24 int64_t vpx_sse_c(const uint8_t *a, int a_stride, const uint8_t *b,
25 int b_stride, int width, int height) {
26 int y, x;
27 int64_t sse = 0;
28
29 for (y = 0; y < height; y++) {
30 for (x = 0; x < width; x++) {
31 const int32_t diff = abs(a[x] - b[x]);
32 sse += diff * diff;
33 }
34
35 a += a_stride;
36 b += b_stride;
37 }
38 return sse;
39 }
40
41 #if CONFIG_VP9_HIGHBITDEPTH
vpx_highbd_sse_c(const uint8_t * a8,int a_stride,const uint8_t * b8,int b_stride,int width,int height)42 int64_t vpx_highbd_sse_c(const uint8_t *a8, int a_stride, const uint8_t *b8,
43 int b_stride, int width, int height) {
44 int y, x;
45 int64_t sse = 0;
46 uint16_t *a = CONVERT_TO_SHORTPTR(a8);
47 uint16_t *b = CONVERT_TO_SHORTPTR(b8);
48 for (y = 0; y < height; y++) {
49 for (x = 0; x < width; x++) {
50 const int32_t diff = (int32_t)(a[x]) - (int32_t)(b[x]);
51 sse += diff * diff;
52 }
53
54 a += a_stride;
55 b += b_stride;
56 }
57 return sse;
58 }
59 #endif
60