xref: /aosp_15_r20/external/libavc/encoder/psnr.c (revision 495ae853bb871d1e5a258cb02c2cc13cde8ddb9a)

/******************************************************************************
 *
 * Copyright (C) 2022 The Android Open Source Project
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at:
 *
 * http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *
 *****************************************************************************/

/* System Include Files */
#include <math.h>

#include "ih264_macros.h"
#include "ih264_typedefs.h"
#include "psnr.h"

void 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])
{
    WORD32 i, j;

    for(j = 0; j < height; j++)
    {
        for(i = 0; i < width; i++)
        {
            WORD32 diff = pu1_src_luma[i] - pu1_rec_luma[i];
            pd_sse[0] += diff * diff;
        }
        pu1_src_luma += src_strd;
        pu1_rec_luma += rec_strd;
    }

    for(j = 0; j < height / 2; j++)
    {
        for(i = 0; i < width / 2; i++)
        {
            WORD32 diff = pu1_src_chroma[i * 2] - pu1_rec_chroma[i * 2];
            pd_sse[1] += diff * diff;
            diff = pu1_src_chroma[i * 2 + 1] - pu1_rec_chroma[i * 2 + 1];
            pd_sse[2] += diff * diff;
        }
        pu1_src_chroma += src_strd;
        pu1_rec_chroma += rec_strd;
    }
}

DOUBLE sse_to_psnr(DOUBLE samples, DOUBLE sse)
{
    DOUBLE psnr;
    if(samples <= 0) return -1;
    if (sse<=0) return MAX_PSNR;
    psnr = 10.0 * (log10(samples) + 2*log10(255) - log10(sse));
    psnr = MIN(MAX_PSNR, psnr);
    return psnr;
}