1 /******************************************************************************
2 *
3 * Copyright (C) 2022 The Android Open Source Project
4 *
5 * Licensed under the Apache License, Version 2.0 (the "License");
6 * you may not use this file except in compliance with the License.
7 * You may obtain a copy of the License at:
8 *
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Unless required by applicable law or agreed to in writing, software
12 * distributed under the License is distributed on an "AS IS" BASIS,
13 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 * See the License for the specific language governing permissions and
15 * limitations under the License.
16 *
17 *****************************************************************************/
18
19 /* System Include Files */
20 #include <math.h>
21
22 #include "ih264_macros.h"
23 #include "ih264_typedefs.h"
24 #include "psnr.h"
25
get_sse(UWORD8 * pu1_src_luma,UWORD8 * pu1_rec_luma,UWORD8 * pu1_src_chroma,UWORD8 * pu1_rec_chroma,WORD32 src_strd,WORD32 rec_strd,WORD32 width,WORD32 height,DOUBLE pd_sse[3])26 void get_sse(UWORD8 *pu1_src_luma, UWORD8 *pu1_rec_luma, UWORD8 *pu1_src_chroma,
27 UWORD8 *pu1_rec_chroma, WORD32 src_strd, WORD32 rec_strd, WORD32 width, WORD32 height,
28 DOUBLE pd_sse[3])
29 {
30 WORD32 i, j;
31
32 for(j = 0; j < height; j++)
33 {
34 for(i = 0; i < width; i++)
35 {
36 WORD32 diff = pu1_src_luma[i] - pu1_rec_luma[i];
37 pd_sse[0] += diff * diff;
38 }
39 pu1_src_luma += src_strd;
40 pu1_rec_luma += rec_strd;
41 }
42
43 for(j = 0; j < height / 2; j++)
44 {
45 for(i = 0; i < width / 2; i++)
46 {
47 WORD32 diff = pu1_src_chroma[i * 2] - pu1_rec_chroma[i * 2];
48 pd_sse[1] += diff * diff;
49 diff = pu1_src_chroma[i * 2 + 1] - pu1_rec_chroma[i * 2 + 1];
50 pd_sse[2] += diff * diff;
51 }
52 pu1_src_chroma += src_strd;
53 pu1_rec_chroma += rec_strd;
54 }
55 }
56
sse_to_psnr(DOUBLE samples,DOUBLE sse)57 DOUBLE sse_to_psnr(DOUBLE samples, DOUBLE sse)
58 {
59 DOUBLE psnr;
60 if(samples <= 0) return -1;
61 if (sse<=0) return MAX_PSNR;
62 psnr = 10.0 * (log10(samples) + 2*log10(255) - log10(sse));
63 psnr = MIN(MAX_PSNR, psnr);
64 return psnr;
65 }
66