// clang-format off
/*
* Copyright(c) 2019 Intel Corporation
* Copyright(c) 2019 Netflix, Inc.
*
* This source code is subject to the terms of the BSD 2 Clause License and
* the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
* was not distributed with this source code in the LICENSE file, you can
* obtain it at https://www.aomedia.org/license/software-license. If the Alliance for Open
* Media Patent License 1.0 was not distributed with this source code in the
* PATENTS file, you can obtain it at https://www.aomedia.org/license/patent-license.
*/

#include <stdio.h>
#include <inttypes.h>

#include "aom_dsp_rtcd.h"
#include "EbPictureControlSet.h"
#include "EbSequenceControlSet.h"
#include "EbMotionEstimation.h"
#include "EbUtility.h"

#include "EbComputeSAD.h"
#include "EbReferenceObject.h"
#include "EbAvcStyleMcp.h"

#include "EbEncIntraPrediction.h"
#include "EbLambdaRateTables.h"
#include "EbTransforms.h"

#include "EbLog.h"
#include "EbResize.h"

/********************************************
 * Constants
 ********************************************/
#define REFERENCE_PIC_LIST_0 0
#define REFERENCE_PIC_LIST_1 1
/*******************************************
 * Compute8x4SAD_Default
 *   Unoptimized 8x4 SAD
 *******************************************/
uint32_t svt_aom_compute8x4_sad_kernel_c(uint8_t *src, // input parameter, source samples Ptr
                                 uint32_t src_stride, // input parameter, source stride
                                 uint8_t *ref, // input parameter, reference samples Ptr
                                 uint32_t ref_stride) // input parameter, reference stride
{
    uint32_t row_number_in_blocks_8x4;
    uint32_t sad_block_8x4 = 0;

    for (row_number_in_blocks_8x4 = 0; row_number_in_blocks_8x4 < 4; ++row_number_in_blocks_8x4) {
        sad_block_8x4 += EB_ABS_DIFF(src[0x00], ref[0x00]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x01], ref[0x01]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x02], ref[0x02]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x03], ref[0x03]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x04], ref[0x04]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x05], ref[0x05]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x06], ref[0x06]);
        sad_block_8x4 += EB_ABS_DIFF(src[0x07], ref[0x07]);
        src += src_stride;
        ref += ref_stride;
    }

    return sad_block_8x4;
}
/*******************************************
 * Compute8x8SAD_Default
 *   Unoptimized 8x8 SAD
 *******************************************/
static uint32_t compute8x8_sad_kernel_c(uint8_t *src, // input parameter, source samples Ptr
                                 uint32_t src_stride, // input parameter, source stride
                                 uint8_t *ref, // input parameter, reference samples Ptr
                                 uint32_t ref_stride) // input parameter, reference stride
{
    uint32_t row_number_in_blocks_8x8;
    uint32_t sad_block_8x8 = 0;

    for (row_number_in_blocks_8x8 = 0; row_number_in_blocks_8x8 < 8; ++row_number_in_blocks_8x8) {
        sad_block_8x8 += EB_ABS_DIFF(src[0x00], ref[0x00]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x01], ref[0x01]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x02], ref[0x02]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x03], ref[0x03]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x04], ref[0x04]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x05], ref[0x05]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x06], ref[0x06]);
        sad_block_8x8 += EB_ABS_DIFF(src[0x07], ref[0x07]);
        src += src_stride;
        ref += ref_stride;
    }

    return sad_block_8x8;
}

/*******************************************
Calculate SAD for 16x16 and its 8x8 sublcoks
and check if there is improvment, if yes keep
the best SAD+MV
*******************************************/
void svt_ext_sad_calculation_8x8_16x16_c(uint8_t *src, uint32_t src_stride, uint8_t *ref,
                                         uint32_t ref_stride, uint32_t *p_best_sad_8x8,
                                         uint32_t *p_best_sad_16x16, uint32_t *p_best_mv8x8,
                                         uint32_t *p_best_mv16x16, uint32_t mv,
                                         uint32_t *p_sad16x16, uint32_t *p_sad8x8, Bool sub_sad) {
    uint32_t sad16x16;

    if (sub_sad) {
        p_sad8x8[0] = (svt_aom_compute8x4_sad_kernel_c(src + 0 * src_stride + 0,
                                               2 * src_stride,
                                               ref + 0 * ref_stride + 0,
                                               2 * ref_stride))
            << 1;
        p_sad8x8[1] = (svt_aom_compute8x4_sad_kernel_c(src + 0 * src_stride + 8,
                                               2 * src_stride,
                                               ref + 0 * ref_stride + 8,
                                               2 * ref_stride))
            << 1;
        p_sad8x8[2] = (svt_aom_compute8x4_sad_kernel_c(src + 8 * src_stride + 0,
                                               2 * src_stride,
                                               ref + 8 * ref_stride + 0,
                                               2 * ref_stride))
            << 1;
        p_sad8x8[3] = (svt_aom_compute8x4_sad_kernel_c(src + 8 * src_stride + 8,
                                               2 * src_stride,
                                               ref + 8 * ref_stride + 8,
                                               2 * ref_stride))
            << 1;
    } else {
        p_sad8x8[0] = compute8x8_sad_kernel_c(
            src + 0 * src_stride + 0, src_stride, ref + 0 * ref_stride + 0, ref_stride);
        p_sad8x8[1] = compute8x8_sad_kernel_c(
            src + 0 * src_stride + 8, src_stride, ref + 0 * ref_stride + 8, ref_stride);
        p_sad8x8[2] = compute8x8_sad_kernel_c(
            src + 8 * src_stride + 0, src_stride, ref + 8 * ref_stride + 0, ref_stride);
        p_sad8x8[3] = compute8x8_sad_kernel_c(
            src + 8 * src_stride + 8, src_stride, ref + 8 * ref_stride + 8, ref_stride);
    }

    if (p_sad8x8[0] < p_best_sad_8x8[0]) {
        p_best_sad_8x8[0] = (uint32_t)p_sad8x8[0];
        p_best_mv8x8[0]   = mv;
    }

    if (p_sad8x8[1] < p_best_sad_8x8[1]) {
        p_best_sad_8x8[1] = (uint32_t)p_sad8x8[1];
        p_best_mv8x8[1]   = mv;
    }

    if (p_sad8x8[2] < p_best_sad_8x8[2]) {
        p_best_sad_8x8[2] = (uint32_t)p_sad8x8[2];
        p_best_mv8x8[2]   = mv;
    }

    if (p_sad8x8[3] < p_best_sad_8x8[3]) {
        p_best_sad_8x8[3] = (uint32_t)p_sad8x8[3];
        p_best_mv8x8[3]   = mv;
    }

    sad16x16 = p_sad8x8[0] + p_sad8x8[1] + p_sad8x8[2] + p_sad8x8[3];
    if (sad16x16 < p_best_sad_16x16[0]) {
        p_best_sad_16x16[0] = (uint32_t)sad16x16;
        p_best_mv16x16[0]   = mv;
    }

    *p_sad16x16 = (uint32_t)sad16x16;
}

/*******************************************
Calculate SAD for 32x32,64x64 from 16x16
and check if there is improvment, if yes keep
the best SAD+MV
*******************************************/
void svt_ext_sad_calculation_32x32_64x64_c(uint32_t *p_sad16x16, uint32_t *p_best_sad_32x32,
                                           uint32_t *p_best_sad_64x64, uint32_t *p_best_mv32x32,
                                           uint32_t *p_best_mv64x64, uint32_t mv,
                                           uint32_t *p_sad32x32) {
    uint32_t sad32x32_0, sad32x32_1, sad32x32_2, sad32x32_3, sad64x64;

    p_sad32x32[0] = sad32x32_0 = p_sad16x16[0] + p_sad16x16[1] + p_sad16x16[2] + p_sad16x16[3];
    if (sad32x32_0 < p_best_sad_32x32[0]) {
        p_best_sad_32x32[0] = sad32x32_0;
        p_best_mv32x32[0]   = mv;
    }

    p_sad32x32[1] = sad32x32_1 = p_sad16x16[4] + p_sad16x16[5] + p_sad16x16[6] + p_sad16x16[7];
    if (sad32x32_1 < p_best_sad_32x32[1]) {
        p_best_sad_32x32[1] = sad32x32_1;
        p_best_mv32x32[1]   = mv;
    }

    p_sad32x32[2] = sad32x32_2 = p_sad16x16[8] + p_sad16x16[9] + p_sad16x16[10] + p_sad16x16[11];
    if (sad32x32_2 < p_best_sad_32x32[2]) {
        p_best_sad_32x32[2] = sad32x32_2;
        p_best_mv32x32[2]   = mv;
    }

    p_sad32x32[3] = sad32x32_3 = p_sad16x16[12] + p_sad16x16[13] + p_sad16x16[14] + p_sad16x16[15];
    if (sad32x32_3 < p_best_sad_32x32[3]) {
        p_best_sad_32x32[3] = sad32x32_3;
        p_best_mv32x32[3]   = mv;
    }
    sad64x64 = sad32x32_0 + sad32x32_1 + sad32x32_2 + sad32x32_3;
    if (sad64x64 < p_best_sad_64x64[0]) {
        p_best_sad_64x64[0] = sad64x64;
        p_best_mv64x64[0]   = mv;
    }
}

/*******************************************
 * svt_ext_eight_sad_calculation_8x8_16x16
 *******************************************/
static void svt_ext_eight_sad_calculation_8x8_16x16(
    uint8_t *src, uint32_t src_stride, uint8_t *ref, uint32_t ref_stride, uint32_t mv,
    uint32_t start_16x16_pos, uint32_t *p_best_sad_8x8, uint32_t *p_best_sad_16x16,
    uint32_t *p_best_mv8x8, uint32_t *p_best_mv16x16, uint32_t p_eight_sad16x16[16][8],
    uint32_t p_eight_sad8x8[64][8], Bool sub_sad) {
    const uint32_t start_8x8_pos = 4 * start_16x16_pos;
    int16_t        x_mv, y_mv;

    (void)p_eight_sad8x8;

    p_best_sad_8x8 += start_8x8_pos;
    p_best_mv8x8 += start_8x8_pos;
    p_best_sad_16x16 += start_16x16_pos;
    p_best_mv16x16 += start_16x16_pos;
    if (sub_sad) {
        uint32_t src_stride_sub = (src_stride << 1);
        uint32_t ref_stride_sub = (ref_stride << 1);
        for (int search_index = 0; search_index < 8; search_index++) {
            uint32_t sad8x8_0 =
                (svt_aom_compute8x4_sad_kernel_c(src, src_stride_sub, ref + search_index, ref_stride_sub))
                << 1;
            if (sad8x8_0 < p_best_sad_8x8[0]) {
                p_best_sad_8x8[0] = (uint32_t)sad8x8_0;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[0]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }

            uint32_t sad8x8_1 = (svt_aom_compute8x4_sad_kernel_c(src + 8,
                                                         src_stride_sub,
                                                         ref + 8 + search_index,
                                                         ref_stride_sub))
                << 1;
            if (sad8x8_1 < p_best_sad_8x8[1]) {
                p_best_sad_8x8[1] = (uint32_t)sad8x8_1;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[1]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }

            uint32_t sad8x8_2 = (svt_aom_compute8x4_sad_kernel_c(src + (src_stride << 3),
                                                         src_stride_sub,
                                                         ref + (ref_stride << 3) + search_index,
                                                         ref_stride_sub))
                << 1;
            if (sad8x8_2 < p_best_sad_8x8[2]) {
                p_best_sad_8x8[2] = (uint32_t)sad8x8_2;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[2]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }

            uint32_t sad8x8_3 = (svt_aom_compute8x4_sad_kernel_c(src + (src_stride << 3) + 8,
                                                         src_stride_sub,
                                                         ref + (ref_stride << 3) + 8 + search_index,
                                                         ref_stride_sub))
                << 1;
            if (sad8x8_3 < p_best_sad_8x8[3]) {
                p_best_sad_8x8[3] = (uint32_t)sad8x8_3;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[3]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }
            uint32_t sad16x16 = p_eight_sad16x16[start_16x16_pos][search_index] = sad8x8_0 +
                sad8x8_1 + sad8x8_2 + sad8x8_3;
            if (sad16x16 < p_best_sad_16x16[0]) {
                p_best_sad_16x16[0] = (uint32_t)sad16x16;
                x_mv                = _MVXT(mv) + (int16_t)search_index;
                y_mv                = _MVYT(mv);
                p_best_mv16x16[0]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }
        }
    } else {
        for (int search_index = 0; search_index < 8; search_index++) {
            uint32_t sad8x8_0 = compute8x8_sad_kernel_c(
                src, src_stride, ref + search_index, ref_stride);
            if (sad8x8_0 < p_best_sad_8x8[0]) {
                p_best_sad_8x8[0] = (uint32_t)sad8x8_0;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[0]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }

            uint32_t sad8x8_1 = (compute8x8_sad_kernel_c(
                src + 8, src_stride, ref + 8 + search_index, ref_stride));
            if (sad8x8_1 < p_best_sad_8x8[1]) {
                p_best_sad_8x8[1] = (uint32_t)sad8x8_1;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[1]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }

            uint32_t sad8x8_2 = (compute8x8_sad_kernel_c(src + (src_stride << 3),
                                                         src_stride,
                                                         ref + (ref_stride << 3) + search_index,
                                                         ref_stride));
            if (sad8x8_2 < p_best_sad_8x8[2]) {
                p_best_sad_8x8[2] = (uint32_t)sad8x8_2;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[2]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }

            uint32_t sad8x8_3 = (compute8x8_sad_kernel_c(src + (src_stride << 3) + 8,
                                                         src_stride,
                                                         ref + (ref_stride << 3) + 8 + search_index,
                                                         ref_stride));
            if (sad8x8_3 < p_best_sad_8x8[3]) {
                p_best_sad_8x8[3] = (uint32_t)sad8x8_3;
                x_mv              = _MVXT(mv) + (int16_t)search_index;
                y_mv              = _MVYT(mv);
                p_best_mv8x8[3]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }
            uint32_t sad16x16 = p_eight_sad16x16[start_16x16_pos][search_index] = sad8x8_0 +
                sad8x8_1 + sad8x8_2 + sad8x8_3;
            if (sad16x16 < p_best_sad_16x16[0]) {
                p_best_sad_16x16[0] = (uint32_t)sad16x16;
                x_mv                = _MVXT(mv) + (int16_t)search_index;
                y_mv                = _MVYT(mv);
                p_best_mv16x16[0]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
            }
        }
    }
}

void svt_ext_all_sad_calculation_8x8_16x16_c(uint8_t *src, uint32_t src_stride, uint8_t *ref,
                                             uint32_t ref_stride, uint32_t mv,
                                             uint32_t *p_best_sad_8x8, uint32_t *p_best_sad_16x16,
                                             uint32_t *p_best_mv8x8, uint32_t *p_best_mv16x16,
                                             uint32_t p_eight_sad16x16[16][8],
                                             uint32_t p_eight_sad8x8[64][8], Bool sub_sad) {
    static const char offsets[16] = {0, 1, 4, 5, 2, 3, 6, 7, 8, 9, 12, 13, 10, 11, 14, 15};
    //---- 16x16 : 0, 1, 4, 5, 2, 3, 6, 7, 8, 9, 12, 13, 10, 11, 14, 15
    for (int y = 0; y < 4; y++) {
        for (int x = 0; x < 4; x++) {
            const uint32_t block_index           = 16 * y * src_stride + 16 * x;
            const uint32_t search_position_index = 16 * y * ref_stride + 16 * x;
            svt_ext_eight_sad_calculation_8x8_16x16(src + block_index,
                                                    src_stride,
                                                    ref + search_position_index,
                                                    ref_stride,
                                                    mv,
                                                    offsets[4 * y + x],
                                                    p_best_sad_8x8,
                                                    p_best_sad_16x16,
                                                    p_best_mv8x8,
                                                    p_best_mv16x16,
                                                    p_eight_sad16x16,
                                                    p_eight_sad8x8,
                                                    sub_sad);
        }
    }
}

/*******************************************
Calculate SAD for 32x32,64x64 from 16x16
and check if there is improvment, if yes keep
the best SAD+MV
*******************************************/
void svt_ext_eight_sad_calculation_32x32_64x64_c(uint32_t  p_sad16x16[16][8],
                                                 uint32_t *p_best_sad_32x32,
                                                 uint32_t *p_best_sad_64x64,
                                                 uint32_t *p_best_mv32x32, uint32_t *p_best_mv64x64,
                                                 uint32_t mv, uint32_t p_sad32x32[4][8]) {
    uint32_t search_index;
    int16_t  x_mv, y_mv;
    for (search_index = 0; search_index < 8; search_index++) {
        uint32_t sad32x32_0, sad32x32_1, sad32x32_2, sad32x32_3, sad64x64;

        p_sad32x32[0][search_index] = sad32x32_0 = p_sad16x16[0][search_index] +
            p_sad16x16[1][search_index] + p_sad16x16[2][search_index] + p_sad16x16[3][search_index];
        if (sad32x32_0 < p_best_sad_32x32[0]) {
            p_best_sad_32x32[0] = sad32x32_0;
            x_mv                = _MVXT(mv) + (int16_t)search_index;
            y_mv                = _MVYT(mv);
            p_best_mv32x32[0]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
        }

        p_sad32x32[1][search_index] = sad32x32_1 = p_sad16x16[4][search_index] +
            p_sad16x16[5][search_index] + p_sad16x16[6][search_index] + p_sad16x16[7][search_index];
        if (sad32x32_1 < p_best_sad_32x32[1]) {
            p_best_sad_32x32[1] = sad32x32_1;
            x_mv                = _MVXT(mv) + (int16_t)search_index;
            y_mv                = _MVYT(mv);
            p_best_mv32x32[1]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
        }

        p_sad32x32[2][search_index] = sad32x32_2 = p_sad16x16[8][search_index] +
            p_sad16x16[9][search_index] + p_sad16x16[10][search_index] +
            p_sad16x16[11][search_index];
        if (sad32x32_2 < p_best_sad_32x32[2]) {
            p_best_sad_32x32[2] = sad32x32_2;
            x_mv                = _MVXT(mv) + (int16_t)search_index;
            y_mv                = _MVYT(mv);
            p_best_mv32x32[2]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
        }

        p_sad32x32[3][search_index] = sad32x32_3 = p_sad16x16[12][search_index] +
            p_sad16x16[13][search_index] + p_sad16x16[14][search_index] +
            p_sad16x16[15][search_index];
        if (sad32x32_3 < p_best_sad_32x32[3]) {
            p_best_sad_32x32[3] = sad32x32_3;
            x_mv                = _MVXT(mv) + (int16_t)search_index;
            y_mv                = _MVYT(mv);
            p_best_mv32x32[3]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
        }

        sad64x64 = sad32x32_0 + sad32x32_1 + sad32x32_2 + sad32x32_3;
        if (sad64x64 < p_best_sad_64x64[0]) {
            p_best_sad_64x64[0] = sad64x64;
            x_mv                = _MVXT(mv) + (int16_t)search_index;
            y_mv                = _MVYT(mv);
            p_best_mv64x64[0]   = ((uint16_t)y_mv << 16) | ((uint16_t)x_mv);
        }
    }
}
/*******************************************
 * open_loop_me_get_search_point_results_block
 *******************************************/
static void open_loop_me_get_eight_search_point_results_block(
    MeContext *me_ctx, // input parameter, ME context Ptr, used to get SB Ptr
    uint32_t   list_index, // input parameter, reference list index
    uint32_t   ref_pic_index,
    uint32_t   search_region_index, // input parameter, search area origin, used to
    // point to reference samples
    int32_t x_search_index, // input parameter, search region position in the
    // horizontal direction, used to derive xMV
    int32_t y_search_index // input parameter, search region position in the
    // vertical direction, used to derive yMV
) {
    // uint32_t ref_luma_stride = ref_pic_ptr->stride_y; // NADER
    // uint8_t  *ref_ptr = ref_pic_ptr->buffer_y; // NADER
    const Bool sub_sad         = (me_ctx->me_search_method == SUB_SAD_SEARCH);
    uint32_t     ref_luma_stride = me_ctx->interpolated_full_stride[list_index][ref_pic_index];
    uint8_t *    ref_ptr         = me_ctx->integer_buffer_ptr[list_index][ref_pic_index] +
        ((ME_FILTER_TAP >> 1) * me_ctx->interpolated_full_stride[list_index][ref_pic_index]) +
        (ME_FILTER_TAP >> 1) + search_region_index;

    uint32_t curr_mv_1 = (((uint16_t)y_search_index) << 16);
    uint16_t curr_mv_2 = ((uint16_t)x_search_index);
    uint32_t curr_mv   = curr_mv_1 | curr_mv_2;

    svt_ext_all_sad_calculation_8x8_16x16(me_ctx->b64_src_ptr,
                                          me_ctx->b64_src_stride,
                                          ref_ptr,
                                          ref_luma_stride,
                                          curr_mv,
                                          me_ctx->p_best_sad_8x8,
                                          me_ctx->p_best_sad_16x16,
                                          me_ctx->p_best_mv8x8,
                                          me_ctx->p_best_mv16x16,
                                          me_ctx->p_eight_sad16x16,
                                          me_ctx->p_eight_sad8x8,
                                          sub_sad);

    svt_ext_eight_sad_calculation_32x32_64x64(me_ctx->p_eight_sad16x16,
                                              me_ctx->p_best_sad_32x32,
                                              me_ctx->p_best_sad_64x64,
                                              me_ctx->p_best_mv32x32,
                                              me_ctx->p_best_mv64x64,
                                              curr_mv,
                                              me_ctx->p_eight_sad32x32);
}
/*******************************************
 * open_loop_me_get_search_point_results_block
 *******************************************/
static void open_loop_me_get_search_point_results_block(
    MeContext *me_ctx, // input parameter, ME context Ptr, used to get SB Ptr
    uint32_t   list_index, // input parameter, reference list index
    uint32_t   ref_pic_index,
    uint32_t   search_region_index, // input parameter, search area origin, used to
    // point to reference samples
    int32_t x_search_index, // input parameter, search region position in the
    // horizontal direction, used to derive xMV
    int32_t y_search_index) // input parameter, search region position in the
// vertical direction, used to derive yMV
{
    const Bool sub_sad = (me_ctx->me_search_method == SUB_SAD_SEARCH);
    uint8_t *    src_ptr = me_ctx->b64_src_ptr;

    // uint8_t  *ref_ptr = ref_pic_ptr->buffer_y; // NADER
    uint8_t *ref_ptr = me_ctx->integer_buffer_ptr[list_index][ref_pic_index] +
        (ME_FILTER_TAP >> 1) +
        ((ME_FILTER_TAP >> 1) * me_ctx->interpolated_full_stride[list_index][ref_pic_index]);
    // uint32_t ref_luma_stride = ref_pic_ptr->stride_y; // NADER
    uint32_t ref_luma_stride = me_ctx->interpolated_full_stride[list_index][ref_pic_index];
    uint32_t search_position_tl_index = search_region_index;
    uint32_t search_position_index;
    uint32_t block_index;
    uint32_t src_next_16x16_offset;
    // uint32_t ref_next_16x16_offset = (ref_pic_ptr->stride_y << 4); // NADER
    uint32_t  ref_next_16x16_offset = (ref_luma_stride << 4);
    uint32_t  curr_mv_1             = (((uint16_t)y_search_index) << 16);
    uint16_t  curr_mv_2             = ((uint16_t)x_search_index);
    uint32_t  curr_mv               = curr_mv_1 | curr_mv_2;
    uint32_t *p_best_sad_8x8        = me_ctx->p_best_sad_8x8;
    uint32_t *p_best_sad_16x16      = me_ctx->p_best_sad_16x16;
    uint32_t *p_best_sad_32x32      = me_ctx->p_best_sad_32x32;
    uint32_t *p_best_sad_64x64      = me_ctx->p_best_sad_64x64;
    uint32_t *p_best_mv8x8          = me_ctx->p_best_mv8x8;
    uint32_t *p_best_mv16x16        = me_ctx->p_best_mv16x16;
    uint32_t *p_best_mv32x32        = me_ctx->p_best_mv32x32;
    uint32_t *p_best_mv64x64        = me_ctx->p_best_mv64x64;
    uint32_t *p_sad32x32            = me_ctx->p_sad32x32;
    uint32_t *p_sad16x16            = me_ctx->p_sad16x16;
    uint32_t *p_sad8x8              = me_ctx->p_sad8x8;

    // TODO: block_index search_position_index could be removed
    const uint32_t src_stride = me_ctx->b64_src_stride;
    src_next_16x16_offset     = src_stride << 4;

    //---- 16x16 : 0
    block_index           = 0;
    search_position_index = search_position_tl_index;

    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[0],
                                      &p_best_sad_16x16[0],
                                      &p_best_mv8x8[0],
                                      &p_best_mv16x16[0],
                                      curr_mv,
                                      &p_sad16x16[0],
                                      &p_sad8x8[0],
                                      sub_sad);

    //---- 16x16 : 1
    block_index           = block_index + 16;
    search_position_index = search_position_tl_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[4],
                                      &p_best_sad_16x16[1],
                                      &p_best_mv8x8[4],
                                      &p_best_mv16x16[1],
                                      curr_mv,
                                      &p_sad16x16[1],
                                      &p_sad8x8[4],
                                      sub_sad);
    //---- 16x16 : 4
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;

    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[16],
                                      &p_best_sad_16x16[4],
                                      &p_best_mv8x8[16],
                                      &p_best_mv16x16[4],
                                      curr_mv,
                                      &p_sad16x16[4],
                                      &p_sad8x8[16],
                                      sub_sad);

    //---- 16x16 : 5
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[20],
                                      &p_best_sad_16x16[5],
                                      &p_best_mv8x8[20],
                                      &p_best_mv16x16[5],
                                      curr_mv,
                                      &p_sad16x16[5],
                                      &p_sad8x8[20],
                                      sub_sad);

    //---- 16x16 : 2
    block_index           = src_next_16x16_offset;
    search_position_index = search_position_tl_index + ref_next_16x16_offset;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[8],
                                      &p_best_sad_16x16[2],
                                      &p_best_mv8x8[8],
                                      &p_best_mv16x16[2],
                                      curr_mv,
                                      &p_sad16x16[2],
                                      &p_sad8x8[8],
                                      sub_sad);
    //---- 16x16 : 3
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[12],
                                      &p_best_sad_16x16[3],
                                      &p_best_mv8x8[12],
                                      &p_best_mv16x16[3],
                                      curr_mv,
                                      &p_sad16x16[3],
                                      &p_sad8x8[12],
                                      sub_sad);
    //---- 16x16 : 6
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[24],
                                      &p_best_sad_16x16[6],
                                      &p_best_mv8x8[24],
                                      &p_best_mv16x16[6],
                                      curr_mv,
                                      &p_sad16x16[6],
                                      &p_sad8x8[24],
                                      sub_sad);
    //---- 16x16 : 7
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[28],
                                      &p_best_sad_16x16[7],
                                      &p_best_mv8x8[28],
                                      &p_best_mv16x16[7],
                                      curr_mv,
                                      &p_sad16x16[7],
                                      &p_sad8x8[28],
                                      sub_sad);

    //---- 16x16 : 8
    block_index           = (src_next_16x16_offset << 1);
    search_position_index = search_position_tl_index + (ref_next_16x16_offset << 1);
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[32],
                                      &p_best_sad_16x16[8],
                                      &p_best_mv8x8[32],
                                      &p_best_mv16x16[8],
                                      curr_mv,
                                      &p_sad16x16[8],
                                      &p_sad8x8[32],
                                      sub_sad);
    //---- 16x16 : 9
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[36],
                                      &p_best_sad_16x16[9],
                                      &p_best_mv8x8[36],
                                      &p_best_mv16x16[9],
                                      curr_mv,
                                      &p_sad16x16[9],
                                      &p_sad8x8[36],
                                      sub_sad);
    //---- 16x16 : 12
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[48],
                                      &p_best_sad_16x16[12],
                                      &p_best_mv8x8[48],
                                      &p_best_mv16x16[12],
                                      curr_mv,
                                      &p_sad16x16[12],
                                      &p_sad8x8[48],
                                      sub_sad);
    //---- 16x16 : 13
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[52],
                                      &p_best_sad_16x16[13],
                                      &p_best_mv8x8[52],
                                      &p_best_mv16x16[13],
                                      curr_mv,
                                      &p_sad16x16[13],
                                      &p_sad8x8[52],
                                      sub_sad);

    //---- 16x16 : 10
    block_index           = (src_next_16x16_offset * 3);
    search_position_index = search_position_tl_index + (ref_next_16x16_offset * 3);
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[40],
                                      &p_best_sad_16x16[10],
                                      &p_best_mv8x8[40],
                                      &p_best_mv16x16[10],
                                      curr_mv,
                                      &p_sad16x16[10],
                                      &p_sad8x8[40],
                                      sub_sad);
    //---- 16x16 : 11
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[44],
                                      &p_best_sad_16x16[11],
                                      &p_best_mv8x8[44],
                                      &p_best_mv16x16[11],
                                      curr_mv,
                                      &p_sad16x16[11],
                                      &p_sad8x8[44],
                                      sub_sad);
    //---- 16x16 : 14
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[56],
                                      &p_best_sad_16x16[14],
                                      &p_best_mv8x8[56],
                                      &p_best_mv16x16[14],
                                      curr_mv,
                                      &p_sad16x16[14],
                                      &p_sad8x8[56],
                                      sub_sad);
    //---- 16x16 : 15
    block_index           = block_index + 16;
    search_position_index = search_position_index + 16;
    svt_ext_sad_calculation_8x8_16x16(src_ptr + block_index,
                                      src_stride,
                                      ref_ptr + search_position_index,
                                      ref_luma_stride,
                                      &p_best_sad_8x8[60],
                                      &p_best_sad_16x16[15],
                                      &p_best_mv8x8[60],
                                      &p_best_mv16x16[15],
                                      curr_mv,
                                      &p_sad16x16[15],
                                      &p_sad8x8[60],
                                      sub_sad);

    svt_ext_sad_calculation_32x32_64x64(p_sad16x16,
                                        p_best_sad_32x32,
                                        p_best_sad_64x64,
                                        p_best_mv32x32,
                                        p_best_mv64x64,
                                        curr_mv,
                                        &p_sad32x32[0]);
}

/*******************************************
 * open_loop_me_fullpel_search_sblock
 *******************************************/
static void open_loop_me_fullpel_search_sblock(MeContext *me_ctx, uint32_t list_index,
                                               uint32_t ref_pic_index, int16_t x_search_area_origin,
                                               int16_t  y_search_area_origin,
                                               uint32_t search_area_width,
                                               uint32_t search_area_height) {
    uint32_t x_search_index, y_search_index;
    uint32_t search_area_width_rest_8 = search_area_width & 7;
    uint32_t search_area_width_mult_8 = search_area_width - search_area_width_rest_8;

    for (y_search_index = 0; y_search_index < search_area_height; y_search_index++) {
        for (x_search_index = 0; x_search_index < search_area_width_mult_8; x_search_index += 8) {
            // this function will do:  x_search_index, +1, +2, ..., +7
            open_loop_me_get_eight_search_point_results_block(
                me_ctx,
                list_index,
                ref_pic_index,
                x_search_index +
                    y_search_index *
                        me_ctx->interpolated_full_stride[list_index][ref_pic_index],
                (int32_t)x_search_index + x_search_area_origin,
                (int32_t)y_search_index + y_search_area_origin);
        }

        for (x_search_index = search_area_width_mult_8; x_search_index < search_area_width;
             x_search_index++) {
            open_loop_me_get_search_point_results_block(
                me_ctx,
                list_index,
                ref_pic_index,
                x_search_index +
                    y_search_index *
                        me_ctx->interpolated_full_stride[list_index][ref_pic_index],
                (int32_t)x_search_index + x_search_area_origin,
                (int32_t)y_search_index + y_search_area_origin);
        }
    }
}

// Perform HME Level 0 for one 64x64 block on the given picture
static void hme_level_0(
    MeContext *me_ctx, // ME context Ptr, used to get/update ME results
    int16_t    org_x, // Block position in the horizontal direction- sixteenth resolution
    int16_t    org_y, // Block position in the vertical direction- sixteenth resolution
    uint32_t   block_width, // Block width - sixteenth resolution
    uint32_t   block_height, // Block height - sixteenth resolution
    int16_t    sa_width, // search area width
    int16_t    sa_height, // search area height
    EbPictureBufferDesc *sixteenth_ref_pic_ptr, // sixteenth-downsampled reference picture
    uint32_t             sr_w, // current search region index in the horizontal direction
    uint32_t             sr_h, // current search region index in the vertical direction
    uint64_t *           best_sad, // output: Level0 SAD at (sr_w, sr_h)
    int16_t *            hme_l0_sc_x, // output: Level0 xMV at (sr_w, sr_h)
    int16_t *            hme_l0_sc_y // output: Level0 yMV at (sr_w, sr_h)
) {
    // round up the search region width to nearest multiple of 8 because the SAD calculation performance (for
    // intrinsic functions) is the same for search region width from 1 to 8
    sa_width = (int16_t)((sa_width + 7) & ~0x07);
    int16_t pad_width  = (int16_t)(sixteenth_ref_pic_ptr->org_x) - 1;
    int16_t pad_height = (int16_t)(sixteenth_ref_pic_ptr->org_y) - 1;

    int16_t x_search_region_distance = sa_width * sr_w;
    int16_t y_search_region_distance = sa_height * sr_h;
    int16_t sa_origin_x = -(int16_t)((sa_width * me_ctx->num_hme_sa_w) >> 1) +
        x_search_region_distance;
    int16_t sa_origin_y = -(int16_t)((sa_height * me_ctx->num_hme_sa_h) >> 1) +
        y_search_region_distance;
    // Correct the left edge of the Search Area if it is not on the reference picture
    if (((org_x + sa_origin_x) < -pad_width)) {
        sa_origin_x = -pad_width - org_x;
        sa_width    = sa_width - (-pad_width - (org_x + sa_origin_x));
    }

    // Correct the right edge of the Search Area if its not on the reference picture
    if (((org_x + sa_origin_x) > (int16_t)sixteenth_ref_pic_ptr->width - 1))
        sa_origin_x = sa_origin_x -
            ((org_x + sa_origin_x) - ((int16_t)sixteenth_ref_pic_ptr->width - 1));

    if (((org_x + sa_origin_x + sa_width) > (int16_t)sixteenth_ref_pic_ptr->width))
        sa_width = MAX(
            1,
            sa_width -
                ((org_x + sa_origin_x + sa_width) - (int16_t)sixteenth_ref_pic_ptr->width));
    // Constrain x_HME_L1 to be a multiple of 8 (round down as cropping alrea performed)
    sa_width = (sa_width < 8) ? sa_width : sa_width & ~0x07;
    // Correct the top edge of the Search Area if it is not on the reference picture
    if (((org_y + sa_origin_y) < -pad_height)) {
        sa_origin_y = -pad_height - org_y;
        sa_height   = sa_height - (-pad_height - (org_y + sa_origin_y));
    }

    // Correct the bottom edge of the Search Area if its not on the reference picture
    if (((org_y + sa_origin_y) > (int16_t)sixteenth_ref_pic_ptr->height - 1))
        sa_origin_y = sa_origin_y -
            ((org_y + sa_origin_y) - ((int16_t)sixteenth_ref_pic_ptr->height - 1));

    if ((org_y + sa_origin_y + sa_height > (int16_t)sixteenth_ref_pic_ptr->height))
        sa_height = MAX(
            1,
            sa_height -
                ((org_y + sa_origin_y + sa_height) - (int16_t)sixteenth_ref_pic_ptr->height));

    // Move to the top left of the search region
    int16_t x_top_left_search_region = ((int16_t)sixteenth_ref_pic_ptr->org_x + org_x) +
        sa_origin_x;
    int16_t y_top_left_search_region = ((int16_t)sixteenth_ref_pic_ptr->org_y + org_y) +
        sa_origin_y;
    uint32_t search_region_index = x_top_left_search_region +
        y_top_left_search_region * sixteenth_ref_pic_ptr->stride_y;

    // Put the first search location into level0 results
    svt_sad_loop_kernel(
        &me_ctx->sixteenth_b64_buffer[0],
        (me_ctx->hme_search_method == FULL_SAD_SEARCH)
        ? me_ctx->sixteenth_b64_buffer_stride
        : me_ctx->sixteenth_b64_buffer_stride * 2,
        &sixteenth_ref_pic_ptr->buffer_y[search_region_index],
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? sixteenth_ref_pic_ptr->stride_y
        : sixteenth_ref_pic_ptr->stride_y * 2,
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? block_height : block_height >> 1,
        block_width,
        /* results */
        best_sad,
        hme_l0_sc_x,
        hme_l0_sc_y,
        /* range */
        sixteenth_ref_pic_ptr->stride_y,
        0, // skip search line
        sa_width,
        sa_height);

    *best_sad = (me_ctx->hme_search_method == FULL_SAD_SEARCH)
        ? *best_sad
        : *best_sad * 2; // Multiply by 2 because considered only ever other line
    *hme_l0_sc_x += sa_origin_x;
    *hme_l0_sc_x *= 4; // Multiply by 4 because operating on 1/4 resolution
    *hme_l0_sc_y += sa_origin_y;
    *hme_l0_sc_y *= 4; // Multiply by 4 because operating on 1/4 resolution

    return;
}

// Perform HME Level 1 for one 64x64 block on the given picture
static void hme_level_1(
    MeContext *me_ctx, // ME context Ptr, used to get/update ME results
    int16_t    org_x, // Block position in the horizontal direction - quarter resolution
    int16_t    org_y, // Block position in the vertical direction - quarter resolution
    uint32_t   block_width, // Block width - quarter resolution
    uint32_t   block_height, // Block height - quarter resolution
    EbPictureBufferDesc *quarter_ref_pic_ptr, // quarter reference picture
    int16_t              sa_width, // hme level 1 search area in width
    int16_t              sa_height, // hme level 1 search area in height
    int16_t              hme_l0_sc_x, // input parameter, best Level0 xMV at (sr_w, sr_h)
    int16_t              hme_l0_sc_y, // input parameter, best Level0 yMV at (sr_w, sr_h)
    uint64_t *           best_sad, // output parameter, Level1 SAD at (sr_w, sr_h)
    int16_t *            hme_l1_sc_x, // output parameter, Level1 xMV at (sr_w, sr_h)
    int16_t *            hme_l1_sc_y // output parameter, Level1 yMV at (sr_w, sr_h)
) {
    // round up the search region width to nearest multiple of 8 because the SAD calculation performance (for
    // intrinsic functions) is the same for search region width from 1 to 8
    sa_width = (int16_t)((sa_width + 7) & ~0x07);

    int16_t pad_width  = (int16_t)(quarter_ref_pic_ptr->org_x) - 1;
    int16_t pad_height = (int16_t)(quarter_ref_pic_ptr->org_y) - 1;

    int16_t sa_origin_x = -(sa_width >> 1) + hme_l0_sc_x;
    int16_t sa_origin_y = -(sa_height >> 1) + hme_l0_sc_y;

    // Correct the left edge of the Search Area if it is not on the reference picture
    if (((org_x + sa_origin_x) < -pad_width)) {
        sa_origin_x = -pad_width - org_x;
        sa_width    = sa_width - (-pad_width - (org_x + sa_origin_x));
    }

    // Correct the right edge of the Search Area if its not on the reference picture
    if (((org_x + sa_origin_x) > (int16_t)quarter_ref_pic_ptr->width - 1))
        sa_origin_x = sa_origin_x -
            ((org_x + sa_origin_x) - ((int16_t)quarter_ref_pic_ptr->width - 1));

    if (((org_x + sa_origin_x + sa_width) > (int16_t)quarter_ref_pic_ptr->width))
        sa_width = MAX(
            1,
            sa_width - ((org_x + sa_origin_x + sa_width) - (int16_t)quarter_ref_pic_ptr->width));

    // Constrain x_HME_L1 to be a multiple of 8 (round down as cropping alrea performed)
    sa_width = (sa_width < 8) ? sa_width : sa_width & ~0x07;

    // Correct the top edge of the Search Area if it is not on the reference picture
    if (((org_y + sa_origin_y) < -pad_height)) {
        sa_origin_y = -pad_height - org_y;
        sa_height   = sa_height - (-pad_height - (org_y + sa_origin_y));
    }

    // Correct the bottom edge of the Search Area if its not on the reference picture
    if (((org_y + sa_origin_y) > (int16_t)quarter_ref_pic_ptr->height - 1))
        sa_origin_y = sa_origin_y -
            ((org_y + sa_origin_y) - ((int16_t)quarter_ref_pic_ptr->height - 1));

    if ((org_y + sa_origin_y + sa_height > (int16_t)quarter_ref_pic_ptr->height))
        sa_height = MAX(
            1,
            sa_height -
                ((org_y + sa_origin_y + sa_height) - (int16_t)quarter_ref_pic_ptr->height));

    // Move to the top left of the search region
    int16_t x_top_left_search_region = ((int16_t)quarter_ref_pic_ptr->org_x + org_x) +
        sa_origin_x;
    int16_t y_top_left_search_region = ((int16_t)quarter_ref_pic_ptr->org_y + org_y) +
        sa_origin_y;
    uint32_t search_region_index = x_top_left_search_region +
        y_top_left_search_region * quarter_ref_pic_ptr->stride_y;

    // Put the first search location into level1 results
    svt_sad_loop_kernel(
        &me_ctx->quarter_b64_buffer[0],
        (me_ctx->hme_search_method == FULL_SAD_SEARCH)
        ? me_ctx->quarter_b64_buffer_stride
        : me_ctx->quarter_b64_buffer_stride * 2,
        &quarter_ref_pic_ptr->buffer_y[search_region_index],
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? quarter_ref_pic_ptr->stride_y
                                                            : quarter_ref_pic_ptr->stride_y * 2,
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? block_height : block_height >> 1,
        block_width,
        /* results */
        best_sad,
        hme_l1_sc_x,
        hme_l1_sc_y,
        /* range */
        quarter_ref_pic_ptr->stride_y,
        0, // skip search line
        sa_width,
        sa_height);

    *best_sad = (me_ctx->hme_search_method == FULL_SAD_SEARCH)
        ? *best_sad
        : *best_sad * 2; // Multiply by 2 because considered only ever other line
    *hme_l1_sc_x += sa_origin_x;
    *hme_l1_sc_x *= 2; // Multiply by 2 because operating on 1/2 resolution
    *hme_l1_sc_y += sa_origin_y;
    *hme_l1_sc_y *= 2; // Multiply by 2 because operating on 1/2 resolution

    return;
}

// Perform HME Level 2 for one 64x64 block on the given picture
void hme_level_2(MeContext *          me_ctx, // ME context Ptr, used to get/update ME results
                 int16_t              org_x, // Block position in the horizontal direction
                 int16_t              org_y, // Block position in the vertical direction
                 uint32_t             block_width, // Block pwidth - full resolution
                 uint32_t             block_height, // Block height - full resolution
                 EbPictureBufferDesc *ref_pic_ptr, // reference picture
                 int16_t              sa_width, // hme level 1 search area in width
                 int16_t              sa_height, // hme level 1 search area in height
                 int16_t              hme_l1_sc_x, // best Level1 xMV at (sr_w, sr_h)
                 int16_t              hme_l1_sc_y, // best Level1 yMV at (sr_w, sr_h)
                 uint64_t *           best_sad, // Level2 SAD at (sr_w, sr_h)
                 int16_t *            hme_l2_sc_x, // Level2 xMV at (sr_w, sr_h)
                 int16_t *            hme_l2_sc_y // Level2 yMV at (sr_w, sr_h)
) {
    // round up the search region width to nearest multiple of 8 because the SAD calculation performance (for
    // intrinsic functions) is the same for search region width from 1 to 8
    sa_width = (int16_t)((sa_width + 7) & ~0x07);

    int16_t pad_width  = (int16_t)BLOCK_SIZE_64 - 1;
    int16_t pad_height = (int16_t)BLOCK_SIZE_64 - 1;

    int16_t sa_origin_x = -(sa_width >> 1) + hme_l1_sc_x;
    int16_t sa_origin_y = -(sa_height >> 1) + hme_l1_sc_y;

    // Correct the left edge of the Search Area if it is not on the reference picture
    if (((org_x + sa_origin_x) < -pad_width)) {
        sa_origin_x = -pad_width - org_x;
        sa_width    = sa_width - (-pad_width - (org_x + sa_origin_x));
    }

    // Correct the right edge of the Search Area if its not on the reference picture
    if (((org_x + sa_origin_x) > (int16_t)ref_pic_ptr->width - 1))
        sa_origin_x = sa_origin_x - ((org_x + sa_origin_x) - ((int16_t)ref_pic_ptr->width - 1));

    if (((org_x + sa_origin_x + sa_width) > (int16_t)ref_pic_ptr->width))
        sa_width = MAX(
            1, sa_width - ((org_x + sa_origin_x + sa_width) - (int16_t)ref_pic_ptr->width));

    // Constrain x_HME_L1 to be a multiple of 8 (round down as cropping already performed)
    sa_width = (sa_width < 8) ? sa_width : sa_width & ~0x07;

    // Correct the top edge of the Search Area if it is not on the reference picture
    if (((org_y + sa_origin_y) < -pad_height)) {
        sa_origin_y = -pad_height - org_y;
        sa_height   = sa_height - (-pad_height - (org_y + sa_origin_y));
    }

    // Correct the bottom edge of the Search Area if its not on the reference picture
    if (((org_y + sa_origin_y) > (int16_t)ref_pic_ptr->height - 1))
        sa_origin_y = sa_origin_y - ((org_y + sa_origin_y) - ((int16_t)ref_pic_ptr->height - 1));

    if ((org_y + sa_origin_y + sa_height > (int16_t)ref_pic_ptr->height))
        sa_height = MAX(
            1, sa_height - ((org_y + sa_origin_y + sa_height) - (int16_t)ref_pic_ptr->height));

    // Move to the top left of the search region
    int16_t  x_top_left_search_region = ((int16_t)ref_pic_ptr->org_x + org_x) + sa_origin_x;
    int16_t  y_top_left_search_region = ((int16_t)ref_pic_ptr->org_y + org_y) + sa_origin_y;
    uint32_t search_region_index      = x_top_left_search_region +
        y_top_left_search_region * ref_pic_ptr->stride_y;

    // Put the first search location into level2 results
    svt_sad_loop_kernel(
        me_ctx->b64_src_ptr,
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? me_ctx->b64_src_stride
        : me_ctx->b64_src_stride * 2,
        &ref_pic_ptr->buffer_y[search_region_index],
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? ref_pic_ptr->stride_y
                                                            : ref_pic_ptr->stride_y * 2,
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? block_height : block_height >> 1,
        block_width,
        /* results */
        best_sad,
        hme_l2_sc_x,
        hme_l2_sc_y,
        /* range */
        ref_pic_ptr->stride_y,
        0, // skip search line
        sa_width,
        sa_height);

    *best_sad = (me_ctx->hme_search_method == FULL_SAD_SEARCH)
        ? *best_sad
        : *best_sad * 2; // Multiply by 2 because considered only ever other line
    *hme_l2_sc_x += sa_origin_x;
    *hme_l2_sc_y += sa_origin_y;

    return;
}
// Nader - to be replaced by loock-up table
/*******************************************
 * svt_aom_get_me_info_index
 *   search the correct index of the motion
 *   info that corresponds to the input
 *   md candidate
 *******************************************/
uint32_t svt_aom_get_me_info_index(uint32_t max_me_block, const BlockGeom *blk_geom, uint32_t geom_offset_x,
                           uint32_t geom_offset_y) {
    // search for motion info
    uint32_t block_index;
    uint32_t me_info_index = 0xFFFFFFF;

    for (block_index = 0; block_index < max_me_block; block_index++) {
        if ((blk_geom->bwidth == partition_width[block_index]) &&
            (blk_geom->bheight == partition_height[block_index]) &&
            ((blk_geom->org_x - geom_offset_x) == pu_search_index_map[block_index][0]) &&
            ((blk_geom->org_y - geom_offset_y) == pu_search_index_map[block_index][1])) {
            me_info_index = block_index;
            break;
        }
    }
    return me_info_index;
}

uint32_t check_00_center(EbPictureBufferDesc *ref_pic_ptr, MeContext *me_ctx,
                         uint32_t sb_origin_x, uint32_t sb_origin_y, uint32_t sb_width,
                         uint32_t sb_height, int16_t *x_search_center, int16_t *y_search_center,
                         uint32_t zz_sad)

{
    uint32_t search_region_index, zero_mv_sad, hme_mv_sad;
    uint64_t hme_mv_cost, zero_mv_cost, search_center_cost;
    int16_t  org_x      = (int16_t)sb_origin_x;
    int16_t  org_y      = (int16_t)sb_origin_y;
    uint32_t subsample_sad = 1;
    int16_t  pad_width     = (int16_t)BLOCK_SIZE_64 - 1;
    int16_t  pad_height    = (int16_t)BLOCK_SIZE_64 - 1;

    search_region_index = (int16_t)ref_pic_ptr->org_x + org_x +
        ((int16_t)ref_pic_ptr->org_y + org_y) * ref_pic_ptr->stride_y;
    if (me_ctx->me_early_exit_th)
        zero_mv_sad = zz_sad;
    else
        zero_mv_sad = svt_nxm_sad_kernel(me_ctx->b64_src_ptr,
            me_ctx->b64_src_stride << subsample_sad,
            &(ref_pic_ptr->buffer_y[search_region_index]),
            ref_pic_ptr->stride_y << subsample_sad,
            sb_height >> subsample_sad,
            sb_width);

    zero_mv_sad = zero_mv_sad << subsample_sad;

    // FIX
    // Correct the left edge of the Search Area if it is not on the reference
    // Picture
    *x_search_center = ((org_x + *x_search_center) < -pad_width) ? -pad_width - org_x
                                                                    : *x_search_center;
    // Correct the right edge of the Search Area if its not on the reference
    // Picture
    *x_search_center = ((org_x + *x_search_center) > (int16_t)ref_pic_ptr->width - 1)
        ? *x_search_center - ((org_x + *x_search_center) - ((int16_t)ref_pic_ptr->width - 1))
        : *x_search_center;
    // Correct the top edge of the Search Area if it is not on the reference
    // Picture
    *y_search_center = ((org_y + *y_search_center) < -pad_height) ? -pad_height - org_y
                                                                     : *y_search_center;
    // Correct the bottom edge of the Search Area if its not on the reference
    // Picture
    *y_search_center = ((org_y + *y_search_center) > (int16_t)ref_pic_ptr->height - 1)
        ? *y_search_center - ((org_y + *y_search_center) - ((int16_t)ref_pic_ptr->height - 1))
        : *y_search_center;
    ///

    zero_mv_cost        = zero_mv_sad << COST_PRECISION;
    search_region_index = (int16_t)(ref_pic_ptr->org_x + org_x) + *x_search_center +
        ((int16_t)(ref_pic_ptr->org_y + org_y) + *y_search_center) * ref_pic_ptr->stride_y;

    hme_mv_sad = svt_nxm_sad_kernel(me_ctx->b64_src_ptr,
                                    me_ctx->b64_src_stride << subsample_sad,
                                    &(ref_pic_ptr->buffer_y[search_region_index]),
                                    ref_pic_ptr->stride_y << subsample_sad,
                                    sb_height >> subsample_sad,
                                    sb_width);

    hme_mv_sad = hme_mv_sad << subsample_sad;
    hme_mv_cost = hme_mv_sad << COST_PRECISION;
    search_center_cost = MIN(zero_mv_cost, hme_mv_cost);

    *x_search_center = (search_center_cost == zero_mv_cost) ? 0 : *x_search_center;
    *y_search_center = (search_center_cost == zero_mv_cost) ? 0 : *y_search_center;
    return hme_mv_sad;
}
// get ME references based on level:
// level: 0 => sixteenth, 1 => quarter, 2 => original

static EbPictureBufferDesc *get_me_reference(PictureParentControlSet *pcs,
                                             MeContext *me_ctx, uint8_t list_index,
                                             uint8_t ref_pic_index, uint8_t level, uint16_t *dist,
                                             uint16_t input_width, uint16_t input_height) {
    EbPictureBufferDesc *ref_pic_ptr;
    ref_pic_ptr = level == 0
        ? me_ctx->me_ds_ref_array[list_index][ref_pic_index].sixteenth_picture_ptr
        : level == 1 ? me_ctx->me_ds_ref_array[list_index][ref_pic_index].quarter_picture_ptr
                     : me_ctx->me_ds_ref_array[list_index][ref_pic_index].picture_ptr;

    if ((input_width >> (2 - level)) != ref_pic_ptr->width ||
        (input_height >> (2 - level)) != ref_pic_ptr->height) {
        SVT_WARN(
            "picture %3llu: HME level%d resolution mismatch! input (%dx%d) != (%dx%d) pa ref. \n",
            pcs->picture_number,
            level,
            input_width >> (2 - level),
            input_height >> (2 - level),
            ref_pic_ptr->width,
            ref_pic_ptr->height);
    }

    *dist = (int16_t)ABS(
        (int64_t)pcs->picture_number -
        (int64_t)me_ctx->me_ds_ref_array[list_index][ref_pic_index].picture_number);
    return ref_pic_ptr;
}

// factor to slowdown the ME search region growth to MAX
uint16_t svt_aom_get_scaled_picture_distance(uint16_t dist) {

    uint8_t round_up = ((dist % 8) == 0) ? 0 : 1;
    return ((dist * 5) / 8) + round_up;
}

/*******************************************
 *   performs integer search motion estimation for
 all avaiable references frames
 *******************************************/
static void integer_search_b64(PictureParentControlSet *pcs, uint32_t b64_index, uint32_t b64_origin_x,
    uint32_t b64_origin_y, MeContext *me_ctx,
    EbPictureBufferDesc *input_ptr) {
    SequenceControlSet *scs        = pcs->scs;
    int16_t             picture_width  = pcs->aligned_width;
    int16_t             picture_height = pcs->aligned_height;
    uint32_t            b64_width      = me_ctx->b64_width;
    uint32_t            b64_height     = me_ctx->b64_height;
    int16_t             pad_width      = (int16_t)BLOCK_SIZE_64 - 1;
    int16_t             pad_height     = (int16_t)BLOCK_SIZE_64 - 1;
    int16_t             org_x       = (int16_t)b64_origin_x;
    int16_t             org_y       = (int16_t)b64_origin_y;
    int16_t             search_area_width;
    int16_t             search_area_height;
    int16_t             x_search_area_origin;
    int16_t             y_search_area_origin;
    int16_t             x_top_left_search_region;
    int16_t             y_top_left_search_region;
    uint32_t            search_region_index;
    uint32_t            num_of_list_to_search;
    uint32_t            list_index;
    uint8_t             ref_pic_index;
    // Final ME Search Center
    int16_t              x_search_center = 0;
    int16_t              y_search_center = 0;
    EbPictureBufferDesc *ref_pic_ptr;
    num_of_list_to_search = me_ctx->num_of_list_to_search;

    // Uni-Prediction motion estimation loop
    // List Loop
    for (list_index = REF_LIST_0; list_index < num_of_list_to_search; ++list_index) {
        uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];

        // Ref Picture Loop
        for (ref_pic_index = 0; ref_pic_index < num_of_ref_pic_to_search; ++ref_pic_index) {
            uint16_t dist = 0;
            ref_pic_ptr   = get_me_reference(pcs,
                                           me_ctx,
                                           list_index,
                                           ref_pic_index,
                                           2,
                                           &dist,
                                           input_ptr->width,
                                           input_ptr->height);
            // Get hme results
            if (me_ctx->search_results[list_index][ref_pic_index].do_ref == 0)
                continue; //so will not get ME results for those references.
            x_search_center = me_ctx->search_results[list_index][ref_pic_index].hme_sc_x;
            y_search_center = me_ctx->search_results[list_index][ref_pic_index].hme_sc_y;
            search_area_width  = me_ctx->me_sa.sa_min.width;
            search_area_height = me_ctx->me_sa.sa_min.height;

            // factor to slowdown the ME search region growth to MAX
            if (me_ctx->me_type != ME_MCTF) {
                dist = svt_aom_get_scaled_picture_distance(dist);
            }
            search_area_width  = MIN((search_area_width * dist), me_ctx->me_sa.sa_max.width);
            search_area_height = MIN((search_area_height * dist), me_ctx->me_sa.sa_max.height);

            // Constrain x_ME to be a multiple of 8 (round up)
            // Update ME search reagion size based on hme-data
            search_area_width = (MAX(1, (search_area_width / me_ctx->reduce_me_sr_divisor[list_index][ref_pic_index])) + 7) & ~0x07;
            search_area_height = MAX(
                3,
                (search_area_height /
                 me_ctx->reduce_me_sr_divisor[list_index][ref_pic_index]));
            int16_t search_area_height_before_sr_reduction = search_area_height;
            int16_t search_area_width_before_sr_reduction  = search_area_width;
            uint64_t best_hme_sad = (uint64_t)~0;
            if (me_ctx->me_early_exit_th) {
                if (me_ctx->zz_sad[list_index][ref_pic_index] <
                    (me_ctx->me_early_exit_th / 6)) {
                    search_area_width  = 1;
                    search_area_height = 1;
                }
            } else {
                uint8_t hme_is_accuarte = 1;
                    if ((x_search_center != 0 || y_search_center != 0) &&
                        (me_ctx->is_ref == TRUE)) {

                        best_hme_sad = check_00_center(
                            ref_pic_ptr,
                            me_ctx,
                            b64_origin_x,
                            b64_origin_y,
                            b64_width,
                            b64_height,
                            &x_search_center,
                            &y_search_center
                            ,
                            me_ctx->zz_sad[list_index][ref_pic_index]
                        );

                        if (x_search_center == 0 && y_search_center == 0)
                            hme_is_accuarte = 0;
                    }
                if (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment == 2) {
                    if ((hme_is_accuarte && (best_hme_sad < (24 * 24))) ||
                        (me_ctx->is_ref && me_ctx->search_results[list_index][ref_pic_index].hme_sad < (24 * 24))) {
                        search_area_height = search_area_height / 2;
                    }
                }
                if (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment == 2) {
                    if (list_index || ref_pic_index) {
                        if (me_ctx->p_sb_best_sad[0][0][0] < 5000)
                            if (search_area_height == search_area_height_before_sr_reduction &&
                                search_area_width == search_area_width_before_sr_reduction) {
                                search_area_height = search_area_height >> 1;
                                search_area_width  = search_area_width >> 1;
                            }
                    }
                }
            }
            svt_initialize_buffer_32bits(
                me_ctx->p_sb_best_sad[list_index][ref_pic_index], 21, 1, MAX_SAD_VALUE);
            me_ctx->p_best_sad_64x64 = &(
                me_ctx->p_sb_best_sad[list_index][ref_pic_index][ME_TIER_ZERO_PU_64x64]);
            me_ctx->p_best_sad_32x32 = &(
                me_ctx->p_sb_best_sad[list_index][ref_pic_index][ME_TIER_ZERO_PU_32x32_0]);
            me_ctx->p_best_sad_16x16 = &(
                me_ctx->p_sb_best_sad[list_index][ref_pic_index][ME_TIER_ZERO_PU_16x16_0]);
            me_ctx->p_best_sad_8x8 = &(
                me_ctx->p_sb_best_sad[list_index][ref_pic_index][ME_TIER_ZERO_PU_8x8_0]);

            me_ctx->p_best_mv64x64 = &(
                me_ctx->p_sb_best_mv[list_index][ref_pic_index][ME_TIER_ZERO_PU_64x64]);
            me_ctx->p_best_mv32x32 = &(
                me_ctx->p_sb_best_mv[list_index][ref_pic_index][ME_TIER_ZERO_PU_32x32_0]);
            me_ctx->p_best_mv16x16 = &(
                me_ctx->p_sb_best_mv[list_index][ref_pic_index][ME_TIER_ZERO_PU_16x16_0]);
            me_ctx->p_best_mv8x8 = &(
                me_ctx->p_sb_best_mv[list_index][ref_pic_index][ME_TIER_ZERO_PU_8x8_0]);

            /* If search area is large enough, check the ME 8x8 SAD variance, and if low, reduce search area
            * (as the 64x64 MVs are likely good for all the 8x8 blocks that make it up).  If the search area
            * is already low, the overhead of searching one additional point will be high (and fruitless, since
            * the minimum search size that will be set by the 8x8 SAD variance algorithm is 8x3.
            */
            if (me_ctx->me_8x8_var_ctrls.enabled && (search_area_width * search_area_height > 24)) {

                x_search_area_origin = x_search_center;
                y_search_area_origin = y_search_center;
                x_top_left_search_region = (int16_t)(ref_pic_ptr->org_x + b64_origin_x) -
                    (ME_FILTER_TAP >> 1) + x_search_area_origin;
                y_top_left_search_region = (int16_t)(ref_pic_ptr->org_y + b64_origin_y) -
                    (ME_FILTER_TAP >> 1) + y_search_area_origin;
                search_region_index = (x_top_left_search_region)+
                    (y_top_left_search_region)*ref_pic_ptr->stride_y;
                me_ctx->integer_buffer_ptr[list_index][ref_pic_index] = &(
                    ref_pic_ptr->buffer_y[search_region_index]);
                me_ctx->interpolated_full_stride[list_index][ref_pic_index] =
                    ref_pic_ptr->stride_y;

                open_loop_me_fullpel_search_sblock(me_ctx,
                    list_index,
                    ref_pic_index,
                    x_search_center,
                    y_search_center,
                    1,
                    1);

                // Since only one point was searched, the 64x64 SAD will be the same as the sum of the 8x8 SADs
                const uint32_t mean_dist_8x8 = me_ctx->p_best_sad_64x64[0] / 64;
                uint32_t sum_ofsq_dist_8x8 = 0;
                for (unsigned i = 0; i < 64; i++) {
                    const int32_t diff = ((int32_t)me_ctx->p_best_sad_8x8[i] - (int32_t)mean_dist_8x8);
                    sum_ofsq_dist_8x8 += diff * diff;
                }

                uint32_t me_8x8_cost_var = (uint32_t)(sum_ofsq_dist_8x8 / 64);

                if (me_8x8_cost_var > me_ctx->me_8x8_var_ctrls.me_sr_mult2_th) {
                    search_area_width = (MAX(1, search_area_width * 3 / 2) + 7) & ~0x7;
                    search_area_height = MAX(1, search_area_height * 3 / 2);
                }

                if (me_8x8_cost_var < me_ctx->me_8x8_var_ctrls.me_sr_div4_th) {
                    search_area_width = (MAX(1, search_area_width >> 2) + 7) & ~0x7;
                    search_area_height = MAX(1, search_area_height >> 2);
                    search_area_height = MAX(3, search_area_height);
                }
                else if (me_8x8_cost_var < me_ctx->me_8x8_var_ctrls.me_sr_div2_th) {
                    search_area_width = (MIN(search_area_width, search_area_width >> 1) + 7) & ~0x7;
                    search_area_height = MIN(search_area_height, search_area_height >> 1);
                    search_area_height = MAX(3, search_area_height);
                }
            }
            x_search_area_origin = x_search_center - (search_area_width >> 1);
            y_search_area_origin = y_search_center - (search_area_height >> 1);

            if (scs->static_config.restricted_motion_vector) {
                // b64_geom in scs and ppcs are different when super-res is enabled
                // ME_OPEN_LOOP is performed on downscaled frames while others (ME_MCTF and ME_FIRST_PASS) are performed on unscaled frames
                B64Geom *b64_geom = me_ctx->me_type != ME_OPEN_LOOP
                    ? scs->b64_geom
                    : pcs->b64_geom;

                int       tile_start_x    = b64_geom[b64_index].tile_start_x;
                int       tile_end_x      = b64_geom[b64_index].tile_end_x;

                // Correct the left edge of the Search Area if it is not on the
                // reference Picture
                x_search_area_origin = ((org_x + x_search_area_origin) < tile_start_x)
                    ? tile_start_x - org_x
                    : x_search_area_origin;
                search_area_width    = ((org_x + x_search_area_origin) < tile_start_x)
                       ? search_area_width - (tile_start_x - (org_x + x_search_area_origin))
                       : search_area_width;
                // Correct the right edge of the Search Area if its not on the
                // reference Picture
                x_search_area_origin = ((org_x + x_search_area_origin) > tile_end_x - 1)
                    ? x_search_area_origin - ((org_x + x_search_area_origin) - (tile_end_x - 1))
                    : x_search_area_origin;
                search_area_width    = ((org_x + x_search_area_origin + search_area_width) >
                                     tile_end_x)
                       ? MAX(1,
                          search_area_width -
                              ((org_x + x_search_area_origin + search_area_width) - tile_end_x))
                       : search_area_width;
                // Constrain x_ME to be a multiple of 8 (round down as cropping
                // already performed)
                search_area_width = (search_area_width < 8) ? search_area_width
                                                            : search_area_width & ~0x07;
            } else {
                // Correct the left edge of the Search Area if it is not on the
                // reference Picture
                x_search_area_origin = ((org_x + x_search_area_origin) < -pad_width)
                    ? -pad_width - org_x
                    : x_search_area_origin;
                search_area_width    = ((org_x + x_search_area_origin) < -pad_width)
                       ? search_area_width - (-pad_width - (org_x + x_search_area_origin))
                       : search_area_width;
                // Correct the right edge of the Search Area if its not on the
                // reference Picture
                x_search_area_origin = ((org_x + x_search_area_origin) > picture_width - 1)
                    ? x_search_area_origin -
                        ((org_x + x_search_area_origin) - (picture_width - 1))
                    : x_search_area_origin;

                search_area_width = ((org_x + x_search_area_origin + search_area_width) >
                                     picture_width)
                    ? MAX(1,
                          search_area_width -
                              ((org_x + x_search_area_origin + search_area_width) -
                               picture_width))
                    : search_area_width;

                // Constrain x_ME to be a multiple of 8 (round down as cropping
                // already performed)
                search_area_width = (search_area_width < 8) ? search_area_width
                                                            : search_area_width & ~0x07;
            }
            if (scs->static_config.restricted_motion_vector) {
                // b64_geom in scs and ppcs are different when super-res is enabled
                // ME_OPEN_LOOP is performed on downscaled frames while others (ME_MCTF and ME_FIRST_PASS) are performed on unscaled frames
                B64Geom *b64_geom = me_ctx->me_type != ME_OPEN_LOOP
                    ? scs->b64_geom
                    : pcs->b64_geom;

                int       tile_start_y    = b64_geom[b64_index].tile_start_y;
                int       tile_end_y      = b64_geom[b64_index].tile_end_y;

                // Correct the top edge of the Search Area if it is not on the
                // reference Picture
                y_search_area_origin = ((org_y + y_search_area_origin) < tile_start_y)
                    ? tile_start_y - org_y
                    : y_search_area_origin;

                search_area_height = ((org_y + y_search_area_origin) < tile_start_y)
                    ? search_area_height - (tile_start_y - (org_y + y_search_area_origin))
                    : search_area_height;

                // Correct the bottom edge of the Search Area if its not on the
                // reference Picture
                y_search_area_origin = ((org_y + y_search_area_origin) > tile_end_y - 1)
                    ? y_search_area_origin - ((org_y + y_search_area_origin) - (tile_end_y - 1))
                    : y_search_area_origin;

                search_area_height = (org_y + y_search_area_origin + search_area_height >
                                      tile_end_y)
                    ? MAX(1,
                          search_area_height -
                              ((org_y + y_search_area_origin + search_area_height) - tile_end_y))
                    : search_area_height;
            } else {
                // Correct the top edge of the Search Area if it is not on the
                // reference Picture
                y_search_area_origin = ((org_y + y_search_area_origin) < -pad_height)
                    ? -pad_height - org_y
                    : y_search_area_origin;
                search_area_height   = ((org_y + y_search_area_origin) < -pad_height)
                      ? search_area_height - (-pad_height - (org_y + y_search_area_origin))
                      : search_area_height;
                // Correct the bottom edge of the Search Area if its not on the
                // reference Picture
                y_search_area_origin = ((org_y + y_search_area_origin) > picture_height - 1)
                    ? y_search_area_origin -
                        ((org_y + y_search_area_origin) - (picture_height - 1))
                    : y_search_area_origin;
                search_area_height   = (org_y + y_search_area_origin + search_area_height >
                                      picture_height)
                      ? MAX(1,
                          search_area_height -
                              ((org_y + y_search_area_origin + search_area_height) -
                               picture_height))
                      : search_area_height;
            }

            x_top_left_search_region            = (int16_t)(ref_pic_ptr->org_x + b64_origin_x) -
                (ME_FILTER_TAP >> 1) + x_search_area_origin;
            y_top_left_search_region = (int16_t)(ref_pic_ptr->org_y + b64_origin_y) -
                (ME_FILTER_TAP >> 1) + y_search_area_origin;
            search_region_index = (x_top_left_search_region) +
                (y_top_left_search_region)*ref_pic_ptr->stride_y;
            me_ctx->integer_buffer_ptr[list_index][ref_pic_index] = &(
                ref_pic_ptr->buffer_y[search_region_index]);
            me_ctx->interpolated_full_stride[list_index][ref_pic_index] =
                ref_pic_ptr->stride_y;

            // Move to the top left of the search region
            x_top_left_search_region = (int16_t)(ref_pic_ptr->org_x + b64_origin_x) + x_search_area_origin;
            y_top_left_search_region = (int16_t)(ref_pic_ptr->org_y + b64_origin_y) + y_search_area_origin;
            open_loop_me_fullpel_search_sblock(me_ctx,
                                               list_index,
                                               ref_pic_index,
                                               x_search_area_origin,
                                               y_search_area_origin,
                                               search_area_width,
                                               search_area_height);
        }
    }
}

/*
  using previous stage ME results (Integer Search) for each reference
  frame. keep only the references that are close to the best reference.
*/
static void me_prune_ref(MeContext *me_ctx) {
    uint8_t num_of_list_to_search = me_ctx->num_of_list_to_search;
    for (uint8_t list_index = REF_LIST_0; list_index < num_of_list_to_search; ++list_index) {
        uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];
        // Ref Picture Loop
        for (uint8_t ref_pic_index = 0; ref_pic_index < num_of_ref_pic_to_search; ++ref_pic_index) {

            me_ctx->search_results[list_index][ref_pic_index].hme_sad = 0;
            // Get hme results
            if (me_ctx->search_results[list_index][ref_pic_index].do_ref == 0) {
                me_ctx->search_results[list_index][ref_pic_index].hme_sad = MAX_SAD_VALUE * 64;
                continue;
            }
            me_ctx->p_best_sad_8x8 = &(
                me_ctx->p_sb_best_sad[list_index][ref_pic_index][ME_TIER_ZERO_PU_8x8_0]);
            // 8x8   [64 partitions]
            for (uint32_t pu_index = 0; pu_index < 64; ++pu_index) {
                uint32_t idx = tab8x8[pu_index];
                me_ctx->search_results[list_index][ref_pic_index].hme_sad += me_ctx->p_best_sad_8x8[idx];
            }
        }
    }
    uint64_t best = (uint64_t)~0;
    for (int i = 0; i < MAX_NUM_OF_REF_PIC_LIST; ++i) {
        for (int j = 0; j < REF_LIST_MAX_DEPTH; ++j) {
            if (me_ctx->search_results[i][j].hme_sad < best) {
                best = me_ctx->search_results[i][j].hme_sad;
            }
        }
    }
    for (uint32_t li = 0; li < MAX_NUM_OF_REF_PIC_LIST; li++) {
        for (uint32_t ri = 0; ri < REF_LIST_MAX_DEPTH; ri++) {
            // Prune references based on ME sad
            uint16_t prune_ref_th =
                me_ctx->me_hme_prune_ctrls.prune_ref_if_me_sad_dev_bigger_than_th;
            if (me_ctx->me_hme_prune_ctrls.enable_me_hme_ref_pruning &&
                (!me_ctx->me_hme_prune_ctrls.protect_closest_refs || ri > 0) &&
                (prune_ref_th != (uint16_t)~0) && (me_ctx->search_results[li][ri].hme_sad - best) * 100 > (prune_ref_th * best)) {
                me_ctx->search_results[li][ri].do_ref = 0;
            }
        }
    }
}

/* perform  motion search over a given search area*/
static void prehme_core(MeContext *me_ctx, int16_t org_x, int16_t org_y, uint32_t sb_width,
                 uint32_t sb_height, EbPictureBufferDesc *sixteenth_ref_pic_ptr,
                 SearchInfo *prehme_data) {
    int16_t  x_top_left_search_region;
    int16_t  y_top_left_search_region;
    uint32_t search_region_index;

    int16_t pad_width  = (int16_t)(sixteenth_ref_pic_ptr->org_x) - 1;
    int16_t pad_height = (int16_t)(sixteenth_ref_pic_ptr->org_y) - 1;

    int16_t search_area_width  = prehme_data->sa.width;
    int16_t search_area_height = prehme_data->sa.height;

    int16_t x_search_area_origin = -(int16_t)(search_area_width >> 1);
    int16_t y_search_area_origin = -(int16_t)(search_area_height >> 1);

    // Correct the left edge of the Search Area if it is not on the reference Picture
    x_search_area_origin = ((org_x + x_search_area_origin) < -pad_width) ? -pad_width - org_x
                                                                            : x_search_area_origin;

    search_area_width = ((org_x + x_search_area_origin) < -pad_width)
        ? search_area_width - (-pad_width - (org_x + x_search_area_origin))
        : search_area_width;

    // Correct the right edge of the Search Area if its not on the reference Picture
    x_search_area_origin = ((org_x + x_search_area_origin) >
                            (int16_t)sixteenth_ref_pic_ptr->width - 1)
        ? x_search_area_origin -
            ((org_x + x_search_area_origin) - ((int16_t)sixteenth_ref_pic_ptr->width - 1))
        : x_search_area_origin;

    search_area_width = ((org_x + x_search_area_origin + search_area_width) >
                         (int16_t)sixteenth_ref_pic_ptr->width)
        ? MAX(1,
              search_area_width -
                  ((org_x + x_search_area_origin + search_area_width) -
                   (int16_t)sixteenth_ref_pic_ptr->width))
        : search_area_width;

    // Correct the top edge of the Search Area if it is not on the reference Picture
    y_search_area_origin = ((org_y + y_search_area_origin) < -pad_height)
        ? -pad_height - org_y
        : y_search_area_origin;

    search_area_height = ((org_y + y_search_area_origin) < -pad_height)
        ? search_area_height - (-pad_height - (org_y + y_search_area_origin))
        : search_area_height;

    // Correct the bottom edge of the Search Area if its not on the reference Picture
    y_search_area_origin = ((org_y + y_search_area_origin) >
                            (int16_t)sixteenth_ref_pic_ptr->height - 1)
        ? y_search_area_origin -
            ((org_y + y_search_area_origin) - ((int16_t)sixteenth_ref_pic_ptr->height - 1))
        : y_search_area_origin;

    search_area_height = (org_y + y_search_area_origin + search_area_height >
                          (int16_t)sixteenth_ref_pic_ptr->height)
        ? MAX(1,
              search_area_height -
                  ((org_y + y_search_area_origin + search_area_height) -
                   (int16_t)sixteenth_ref_pic_ptr->height))
        : search_area_height;

    x_top_left_search_region = ((int16_t)sixteenth_ref_pic_ptr->org_x + org_x) +
        x_search_area_origin;
    y_top_left_search_region = ((int16_t)sixteenth_ref_pic_ptr->org_y + org_y) +
        y_search_area_origin;
    search_region_index = x_top_left_search_region +
        y_top_left_search_region * sixteenth_ref_pic_ptr->stride_y;

    svt_sad_loop_kernel(
        &me_ctx->sixteenth_b64_buffer[0],
        me_ctx->hme_search_method == FULL_SAD_SEARCH
        ? me_ctx->sixteenth_b64_buffer_stride
        : me_ctx->sixteenth_b64_buffer_stride * 2,
        &sixteenth_ref_pic_ptr->buffer_y[search_region_index],
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? sixteenth_ref_pic_ptr->stride_y
                                                            : sixteenth_ref_pic_ptr->stride_y * 2,
        (me_ctx->hme_search_method == FULL_SAD_SEARCH) ? sb_height : sb_height >> 1,
        sb_width,
        /* results */
        &prehme_data->sad,
        &prehme_data->best_mv.as_mv.col,
        &prehme_data->best_mv.as_mv.row,
        sixteenth_ref_pic_ptr->stride_y,
        me_ctx->prehme_ctrl.skip_search_line,
        search_area_width,
        search_area_height);

    prehme_data->sad = (me_ctx->hme_search_method == FULL_SAD_SEARCH)
        ? prehme_data->sad
        : prehme_data->sad * 2; // Multiply by 2 because considered only ever other line
    prehme_data->best_mv.as_mv.col += x_search_area_origin;
    prehme_data->best_mv.as_mv.col *= 4; // Multiply by 4 because operating on 1/4 resolution
    prehme_data->best_mv.as_mv.row += y_search_area_origin;
    prehme_data->best_mv.as_mv.row *= 4; // Multiply by 4 because operating on 1/4 resolution
    prehme_data->valid = 1;
    return;
}
static uint32_t get_zz_sad(EbPictureBufferDesc *ref_pic_ptr, MeContext *me_ctx, uint32_t sb_origin_x,
                    uint32_t sb_origin_y, uint32_t sb_width, uint32_t sb_height)

{
    uint32_t search_region_index, zero_mv_sad;
    int16_t  org_x      = (int16_t)sb_origin_x;
    int16_t  org_y      = (int16_t)sb_origin_y;
    uint32_t subsample_sad = 1;

    search_region_index = (int16_t)ref_pic_ptr->org_x + org_x +
        ((int16_t)ref_pic_ptr->org_y + org_y) * ref_pic_ptr->stride_y;

    zero_mv_sad = svt_nxm_sad_kernel(me_ctx->b64_src_ptr,
                                     me_ctx->b64_src_stride << subsample_sad,
                                     &(ref_pic_ptr->buffer_y[search_region_index]),
                                     ref_pic_ptr->stride_y << subsample_sad,
                                     sb_height >> subsample_sad,
                                     sb_width);

    zero_mv_sad = zero_mv_sad << subsample_sad;

    return zero_mv_sad;
}
// Determine if pre-HME for the current picture and search region should be skipped.
// Return 1 if can early exit (i.e. skip pre-hme for current frame and search region)
// Return 0 if can't skip
static Bool check_prehme_early_exit(MeContext *me_ctx, uint8_t list_i, uint8_t ref_i, uint8_t sr_i) {
    SearchInfo *prehme_data = &me_ctx->prehme_data[list_i][ref_i][sr_i];

    if (me_ctx->me_early_exit_th) {
        if (me_ctx->zz_sad[list_i][ref_i] < me_ctx->me_early_exit_th) {
            prehme_data->best_mv.as_mv.col = 0;
            prehme_data->best_mv.as_mv.row = 0;
            prehme_data->sad               = 0;
            prehme_data->valid             = 1;
            return 1;
        }
    }

    if (me_ctx->prehme_ctrl.l1_early_exit) {
        if (list_i == 1 && me_ctx->prehme_data[0][ref_i][sr_i].valid  &&
            ((me_ctx->prehme_data[0][ref_i][sr_i].sad < (32 * 32)) ||
             ((ABS(me_ctx->prehme_data[0][ref_i][sr_i].best_mv.as_mv.col) < 16) &&
              (ABS(me_ctx->prehme_data[0][ref_i][sr_i].best_mv.as_mv.row) < 16)))) {
            prehme_data->best_mv.as_mv.col = -me_ctx->prehme_data[0][ref_i][sr_i].best_mv.as_mv.col;
            prehme_data->best_mv.as_mv.row = -me_ctx->prehme_data[0][ref_i][sr_i].best_mv.as_mv.row;
            prehme_data->sad               = me_ctx->prehme_data[0][ref_i][sr_i].sad;
            prehme_data->valid             = 1;
            return 1;
        }
    }
    return 0;
}

/* Perform Pre-HME for one Block 64x64*/
static void prehme_b64(PictureParentControlSet *pcs, uint32_t org_x, uint32_t org_y,
                       MeContext *me_ctx, EbPictureBufferDesc *input_ptr) {
    const uint32_t block_width  = me_ctx->b64_width;
    const uint32_t block_height = me_ctx->b64_height;
    uint32_t best_sad = MAX_U32;
    // List Loop
    for (int list_i = REF_LIST_0; list_i < me_ctx->num_of_list_to_search; ++list_i) {
        // Ref Picture Loop
        const uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_i];
        for (uint8_t ref_i = 0; ref_i < num_of_ref_pic_to_search; ++ref_i) {
            uint16_t             dist              = 0;
            EbPictureBufferDesc *sixteenth_ref_pic = get_me_reference(
                pcs, me_ctx, list_i, ref_i, 0, &dist, input_ptr->width, input_ptr->height);

            if (me_ctx->temporal_layer_index > 0 || list_i == 0) {

                uint32_t hme_sr_factor = svt_aom_get_scaled_picture_distance(dist);

                for (uint8_t sr_i = 0; sr_i < SEARCH_REGION_COUNT; sr_i++) {
                    if (check_prehme_early_exit(me_ctx, list_i, ref_i, sr_i))
                        continue;

                    SearchInfo *prehme_data = &me_ctx->prehme_data[list_i][ref_i][sr_i];
                    if (!me_ctx->search_results[list_i][ref_i].do_ref) {
                        prehme_data->best_mv.as_mv.col = 0;
                        prehme_data->best_mv.as_mv.row = 0;
                        prehme_data->sad = MAX_U32;
                        continue;
                    }
                    prehme_data->sa.width = MIN(
                        (me_ctx->prehme_ctrl.prehme_sa_cfg[sr_i].sa_min.width * hme_sr_factor),
                        me_ctx->prehme_ctrl.prehme_sa_cfg[sr_i].sa_max.width);
                    prehme_data->sa.height = MIN(
                        (me_ctx->prehme_ctrl.prehme_sa_cfg[sr_i].sa_min.height * hme_sr_factor),
                        me_ctx->prehme_ctrl.prehme_sa_cfg[sr_i].sa_max.height);

                    prehme_core(me_ctx,
                                ((int16_t)org_x) >> 2,
                                ((int16_t)org_y) >> 2,
                                block_width >> 2,
                                block_height >> 2,
                                sixteenth_ref_pic,
                                prehme_data);
                    me_ctx->performed_phme[list_i][ref_i][sr_i] = 1;
                }
                uint32_t min_sad = (uint32_t)MIN(me_ctx->prehme_data[list_i][ref_i][0].sad, me_ctx->prehme_data[list_i][ref_i][1].sad);
                best_sad = MIN(best_sad, min_sad);
            } else {
                // PW: Does this account for base pictures
                for (uint8_t sr_i = 0; sr_i < SEARCH_REGION_COUNT; sr_i++) {
                    me_ctx->prehme_data[1][ref_i][sr_i].best_mv.as_mv.col =
                        -me_ctx->prehme_data[0][ref_i][sr_i].best_mv.as_mv.col;
                    me_ctx->prehme_data[1][ref_i][sr_i].best_mv.as_mv.row =
                        -me_ctx->prehme_data[0][ref_i][sr_i].best_mv.as_mv.row;
                    me_ctx->prehme_data[1][ref_i][sr_i].sad = me_ctx->prehme_data[0][ref_i][sr_i].sad;
                }
            }
        } // End ref pic loop
    } // End list loop
    if (me_ctx->temporal_layer_index > 0 && best_sad < me_ctx->me_hme_prune_ctrls.phme_sad_th) {
        for (int list_i = REF_LIST_0; list_i < me_ctx->num_of_list_to_search; ++list_i) {
            for (uint8_t ref_i = 0; ref_i < me_ctx->num_of_ref_pic_to_search[list_i]; ++ref_i) {
                if (!me_ctx->search_results[list_i][ref_i].do_ref) continue;
                if (ref_i == 0) continue;

                const uint32_t prhme_th = me_ctx->me_hme_prune_ctrls.phme_sad_pct;
                uint32_t prehme_sad = (uint32_t)MIN(me_ctx->prehme_data[list_i][ref_i][0].sad,
                                                    me_ctx->prehme_data[list_i][ref_i][1].sad);
                if ((prehme_sad - best_sad) * 100 > (prhme_th * best_sad)) {
                    me_ctx->search_results[list_i][ref_i].do_ref = 0;
                }
            }
        }
    }
}

// Set the HME L0 search area.  Perform scaling based on list index and ref index.
// HME L0 search area should be the same for each search region
static void get_hme_l0_search_area(MeContext *me_ctx, uint8_t list_index,
                                   uint8_t ref_pic_index, uint16_t dist, int16_t *sa_width,
                                   int16_t *sa_height) {
    // Reduce HME search area for higher ref indices
    if (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment &&
        me_ctx->me_sr_adjustment_ctrls.distance_based_hme_resizing) {
        uint8_t is_hor   = 1;
        uint8_t is_ver   = 1;
        uint8_t is_still = 0;

        if (me_ctx->reduce_hme_l0_sr_th_min && me_ctx->reduce_hme_l0_sr_th_max) {
            if (list_index || ref_pic_index) {
                int16_t l0_mvx = me_ctx->x_hme_level0_search_center[0][0][0 /*quadrant-x*/]
                                                                        [0 /*quadrant-y*/];
                int16_t l0_mvy = me_ctx->y_hme_level0_search_center[0][0][0 /*quadrant-x*/]
                                                                        [0 /*quadrant-y*/];

                // Determine whether the computed motion from list0/ref_index0 is in vertical or horizintal direction
                is_ver   = ((ABS(l0_mvx) < me_ctx->reduce_hme_l0_sr_th_min) &&
                          (ABS(l0_mvy) > me_ctx->reduce_hme_l0_sr_th_max));
                is_hor   = ((ABS(l0_mvx) > me_ctx->reduce_hme_l0_sr_th_max) &&
                          (ABS(l0_mvy) < me_ctx->reduce_hme_l0_sr_th_min));
                is_still = ((ABS(l0_mvx) < (me_ctx->reduce_hme_l0_sr_th_min * 3)) &&
                            (ABS(l0_mvy) < (me_ctx->reduce_hme_l0_sr_th_min * 3)));
            }
        }

        uint8_t x_offset = 1;
        uint8_t y_offset = 1;
        if (!is_ver) {
            y_offset = 2;
        }
        if (!is_hor) {
            x_offset = 2;
        }

        if (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment == 2) {
            if (is_still) {
                x_offset = 4;
                y_offset = 4;
            }
        }

        me_ctx->hme_l0_sa.sa_min.width = me_ctx->hme_l0_sa.sa_min.width /
            (x_offset + ref_pic_index);
        me_ctx->hme_l0_sa.sa_min.height = me_ctx->hme_l0_sa.sa_min.height /
            (y_offset + ref_pic_index);
        me_ctx->hme_l0_sa.sa_max.width = me_ctx->hme_l0_sa.sa_max.width /
            (x_offset + ref_pic_index);
        me_ctx->hme_l0_sa.sa_max.height = me_ctx->hme_l0_sa.sa_max.height /
            (y_offset + ref_pic_index);
    }

    int32_t hme_sr_factor = svt_aom_get_scaled_picture_distance(dist);

    // Derive the search area width and height, rounding the width up to the nearest sixteenth
    int16_t search_area_width = me_ctx->hme_l0_sa.sa_min.width / me_ctx->num_hme_sa_w;
    search_area_width         = (int16_t)MIN(
        (((search_area_width * hme_sr_factor) + 15) & ~0x0F),
        (((me_ctx->hme_l0_sa.sa_max.width / me_ctx->num_hme_sa_w) + 15) & ~0x0F));
    int16_t search_area_height = me_ctx->hme_l0_sa.sa_min.height / me_ctx->num_hme_sa_h;
    search_area_height         = (int16_t)MIN(
        (search_area_height * hme_sr_factor),
        me_ctx->hme_l0_sa.sa_max.height / me_ctx->num_hme_sa_h);

    *sa_width  = search_area_width;
    *sa_height = search_area_height;
}

//this functions returns the worst quadrant in terms of sad.
//it is implemented w/o for loops to get away from a VS2022 compiler issue.
//it then assumes a fixed quadrant sizes of 2 each direction.
static void get_worst_quadrant( MeContext *me_ctx,uint32_t list_index,uint32_t ref_pic_index ,uint8_t *best_w, uint8_t *best_h){


    if( me_ctx->num_hme_sa_w!=2 || me_ctx->num_hme_sa_h!=2) {
        svt_aom_assert_err(0,"update other quadrant sizes");
        return;
    }
    uint64_t max_sad  = 0;

    if (me_ctx->hme_level0_sad[list_index][ref_pic_index][0][0] > max_sad) {
         max_sad  =  me_ctx->hme_level0_sad[list_index][ref_pic_index][0][0];
         *best_w = 0;
         *best_h = 0;
    }
    if (me_ctx->hme_level0_sad[list_index][ref_pic_index][1][0] > max_sad) {
         max_sad  =  me_ctx->hme_level0_sad[list_index][ref_pic_index][1][0];
         *best_w = 1;
         *best_h = 0;
    }
    if (me_ctx->hme_level0_sad[list_index][ref_pic_index][0][1] > max_sad) {
         max_sad  =  me_ctx->hme_level0_sad[list_index][ref_pic_index][0][1];
         *best_w = 0;
         *best_h = 1;
    }
    if (me_ctx->hme_level0_sad[list_index][ref_pic_index][1][1] > max_sad) {
         *best_w = 1;
         *best_h = 1;
    }

}

/*******************************************
 * performs hierarchical ME level 0 for one 64x64 block (uni-prediction only)
 *******************************************/
static void hme_level0_b64(PictureParentControlSet *pcs, uint32_t org_x, uint32_t org_y,
                           MeContext *me_ctx, EbPictureBufferDesc *input_ptr) {
    const uint32_t block_width  = me_ctx->b64_width;
    const uint32_t block_height = me_ctx->b64_height;

    // store base HME sizes, to be used if using ref-index based HME resizing
    SearchAreaMinMax base_hme_sa;
    base_hme_sa.sa_min = (SearchArea){me_ctx->hme_l0_sa.sa_min.width, me_ctx->hme_l0_sa.sa_min.height};
    base_hme_sa.sa_max = (SearchArea){me_ctx->hme_l0_sa.sa_max.width, me_ctx->hme_l0_sa.sa_max.height};

    // List Loop
    const uint8_t num_of_list_to_search = me_ctx->num_of_list_to_search;
    for (uint8_t list_index = REF_LIST_0; list_index < num_of_list_to_search; ++list_index) {
        // Ref Picture Loop
        const uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];
        for (uint8_t ref_pic_index = 0; ref_pic_index < num_of_ref_pic_to_search; ++ref_pic_index) {
            // If me_early_exit_th is enabled, skip HME L0 for the current block if the zero-zero SAD is low
            if (me_ctx->me_early_exit_th) {
                if (me_ctx->zz_sad[list_index][ref_pic_index] < (me_ctx->me_early_exit_th >> 2)) {
                    for (uint32_t sr_idx_y = 0; sr_idx_y < me_ctx->num_hme_sa_h; sr_idx_y++) {
                        for (uint32_t sr_idx_x = 0; sr_idx_x < me_ctx->num_hme_sa_w; sr_idx_x++) {
                            me_ctx->x_hme_level0_search_center[list_index][ref_pic_index][sr_idx_x]
                                                           [sr_idx_y] = 0;
                            me_ctx->y_hme_level0_search_center[list_index][ref_pic_index][sr_idx_x]
                                                           [sr_idx_y]                          = 0;
                            me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = 0;
                        }
                    }
                    continue;
                }
            }
            if (me_ctx->prev_me_stage_based_exit_th) {
                uint8_t sr_i = me_ctx->prehme_data[list_index][ref_pic_index][0].sad <=
                    me_ctx->prehme_data[list_index][ref_pic_index][1].sad
                    ? 0
                    : 1;
                if (me_ctx->performed_phme[list_index][ref_pic_index][sr_i]) {
                    if (me_ctx->prehme_data[list_index][ref_pic_index][sr_i].sad < (me_ctx->prev_me_stage_based_exit_th >> 4)) {
                        for (uint32_t sr_idx_y = 0; sr_idx_y < me_ctx->num_hme_sa_h; sr_idx_y++) {
                            for (uint32_t sr_idx_x = 0; sr_idx_x < me_ctx->num_hme_sa_w; sr_idx_x++) {
                                me_ctx->x_hme_level0_search_center[list_index][ref_pic_index][sr_idx_x]
                                    [sr_idx_y] = me_ctx->prehme_data[list_index][ref_pic_index][sr_i].best_mv.as_mv.col;
                                me_ctx->y_hme_level0_search_center[list_index][ref_pic_index][sr_idx_x]
                                    [sr_idx_y] = me_ctx->prehme_data[list_index][ref_pic_index][sr_i].best_mv.as_mv.row;
                                me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = me_ctx->prehme_data[list_index][ref_pic_index][sr_i].sad;
                            }
                        }
                        continue;
                    }
                }
            }

            if (!me_ctx->search_results[list_index][ref_pic_index].do_ref) {

                for (uint32_t sr_idx_y = 0; sr_idx_y < me_ctx->num_hme_sa_h; sr_idx_y++) {
                    for (uint32_t sr_idx_x = 0; sr_idx_x < me_ctx->num_hme_sa_w; sr_idx_x++) {
                        me_ctx->x_hme_level0_search_center[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = 0;
                        me_ctx->y_hme_level0_search_center[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = 0;
                        me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = MAX_U32;
                    }
                }
                continue;
            }
            // Get the sixteenth downsampled reference picture
            uint16_t             dist              = 0;
            EbPictureBufferDesc *sixteenth_ref_pic = get_me_reference(
                pcs, me_ctx, list_index, ref_pic_index, 0, &dist, input_ptr->width, input_ptr->height);

            if (me_ctx->temporal_layer_index > 0 || list_index == 0) {
                // Get the HME L0 search dimensions for the current frame
                int16_t sa_width = 0, sa_height = 0;
                get_hme_l0_search_area(me_ctx, list_index, ref_pic_index, dist, &sa_width, &sa_height);
                for (uint8_t sr_h = 0; sr_h < me_ctx->num_hme_sa_h; sr_h++) {
                    for (uint8_t sr_w = 0; sr_w < me_ctx->num_hme_sa_w; sr_w++) {
                        hme_level_0(me_ctx,
                                    ((int16_t)org_x) >> 2,
                                    ((int16_t)org_y) >> 2,
                                    block_width >> 2,
                                    block_height >> 2,
                                    sa_width,
                                    sa_height,
                                    sixteenth_ref_pic,
                                    sr_w,
                                    sr_h,
                                    &(me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_w][sr_h]),
                                    &(me_ctx->x_hme_level0_search_center[list_index][ref_pic_index]
                                                                     [sr_w][sr_h]),
                                    &(me_ctx->y_hme_level0_search_center[list_index][ref_pic_index]
                                                                     [sr_w][sr_h]));
                    }
                }

                // reset base HME area
                if (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment &&
                    me_ctx->me_sr_adjustment_ctrls.distance_based_hme_resizing) {
                    me_ctx->hme_l0_sa.sa_min = base_hme_sa.sa_min;
                    me_ctx->hme_l0_sa.sa_max = base_hme_sa.sa_max;
                }

                if (me_ctx->prehme_ctrl.enable) {
                    //get the worst quadrant
                    uint8_t  sr_h_max = 0, sr_w_max = 0;
                    get_worst_quadrant(
                               me_ctx,
                               list_index,
                               ref_pic_index,
                               &sr_w_max, &sr_h_max);

                    uint8_t sr_i = me_ctx->prehme_data[list_index][ref_pic_index][0].sad <=
                            me_ctx->prehme_data[list_index][ref_pic_index][1].sad
                        ? 0
                        : 1;
                    //replace worst with pre-hme
                    if (me_ctx->prehme_data[list_index][ref_pic_index][sr_i].sad <
                        me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_w_max][sr_h_max]) {
                        me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_w_max][sr_h_max] =
                            me_ctx->prehme_data[list_index][ref_pic_index][sr_i].sad;

                        me_ctx->x_hme_level0_search_center[list_index][ref_pic_index][sr_w_max]
                                                       [sr_h_max] =
                            me_ctx->prehme_data[list_index][ref_pic_index][sr_i].best_mv.as_mv.col;

                        me_ctx->y_hme_level0_search_center[list_index][ref_pic_index][sr_w_max]
                                                       [sr_h_max] =
                            me_ctx->prehme_data[list_index][ref_pic_index][sr_i].best_mv.as_mv.row;
                    }
                }
            }
        } // End ref pic loop
    } // End list loop
}

/*******************************************
 * performs hierarchical ME level 1 for one 64x64 block (uni-prediction only)
 *******************************************/
static void hme_level1_b64(PictureParentControlSet *pcs, uint32_t org_x, uint32_t org_y,
                    MeContext *me_ctx, EbPictureBufferDesc *input_ptr) {
    const uint32_t block_width  = me_ctx->b64_width;
    const uint32_t block_height = me_ctx->b64_height;

    // List Loop
    const uint8_t num_of_list_to_search = me_ctx->num_of_list_to_search;
    for (uint32_t list_index = REF_LIST_0; list_index < num_of_list_to_search; ++list_index) {
        // Ref Picture Loop
        const uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];
        for (uint8_t ref_pic_index = 0; ref_pic_index < num_of_ref_pic_to_search; ++ref_pic_index) {
            uint16_t             dist            = 0;
            EbPictureBufferDesc *quarter_ref_pic = get_me_reference(
                pcs, me_ctx, list_index, ref_pic_index, 1, &dist, input_ptr->width, input_ptr->height);

            if (me_ctx->temporal_layer_index > 0 || list_index == 0) {
                // If me_early_exit_th is enabled, skip HME L0 for the current block if the zero-zero SAD is low
                if (me_ctx->me_early_exit_th) {
                    if (me_ctx->zz_sad[list_index][ref_pic_index] < (me_ctx->me_early_exit_th >> 2)) {
                        for (uint32_t sr_idx_y = 0; sr_idx_y < me_ctx->num_hme_sa_h; sr_idx_y++) {
                            for (uint32_t sr_idx_x = 0; sr_idx_x < me_ctx->num_hme_sa_w; sr_idx_x++) {
                                me_ctx->x_hme_level1_search_center[list_index][ref_pic_index][sr_idx_x]
                                                               [sr_idx_y] = 0;
                                me_ctx->y_hme_level1_search_center[list_index][ref_pic_index][sr_idx_x]
                                                               [sr_idx_y] = 0;
                                me_ctx->hme_level1_sad[list_index][ref_pic_index][sr_idx_x][sr_idx_y] =
                                    0;
                            }
                        }
                        continue;
                    }
                }
                if (!me_ctx->search_results[list_index][ref_pic_index].do_ref) {

                    for (uint32_t sr_idx_y = 0; sr_idx_y < me_ctx->num_hme_sa_h; sr_idx_y++) {
                        for (uint32_t sr_idx_x = 0; sr_idx_x < me_ctx->num_hme_sa_w; sr_idx_x++) {
                            me_ctx->x_hme_level1_search_center[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = 0;
                            me_ctx->y_hme_level1_search_center[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = 0;
                            me_ctx->hme_level1_sad[list_index][ref_pic_index][sr_idx_x][sr_idx_y] = MAX_U32;
                        }
                    }
                    continue;
                }
                for (uint8_t sr_h = 0; sr_h < me_ctx->num_hme_sa_h; sr_h++) {
                    for (uint8_t sr_w = 0; sr_w < me_ctx->num_hme_sa_w; sr_w++) {
                        if (me_ctx->prev_me_stage_based_exit_th) {
                            if (me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_w][sr_h] < (me_ctx->prev_me_stage_based_exit_th >> 5)) {
                                me_ctx->x_hme_level1_search_center[list_index][ref_pic_index][sr_w][sr_h] =
                                    me_ctx->x_hme_level0_search_center[list_index][ref_pic_index][sr_w][sr_h];
                                me_ctx->y_hme_level1_search_center[list_index][ref_pic_index][sr_w][sr_h] =
                                    me_ctx->y_hme_level0_search_center[list_index][ref_pic_index][sr_w][sr_h];
                                me_ctx->hme_level1_sad[list_index][ref_pic_index][sr_w][sr_h] =
                                    me_ctx->hme_level0_sad[list_index][ref_pic_index][sr_w][sr_h];
                                continue;
                            }
                        }

                        hme_level_1(me_ctx,
                                    ((int16_t)org_x) >> 1,
                                    ((int16_t)org_y) >> 1,
                                    block_width >> 1,
                                    block_height >> 1,
                                    quarter_ref_pic,
                                    (int16_t)me_ctx->hme_l1_sa.width,
                                    (int16_t)me_ctx->hme_l1_sa.height,
                                    me_ctx->x_hme_level0_search_center[list_index][ref_pic_index][sr_w]
                                                                   [sr_h] >>
                                        1,
                                    me_ctx->y_hme_level0_search_center[list_index][ref_pic_index][sr_w]
                                                                   [sr_h] >>
                                        1,
                                    &(me_ctx->hme_level1_sad[list_index][ref_pic_index][sr_w][sr_h]),
                                    &(me_ctx->x_hme_level1_search_center[list_index][ref_pic_index]
                                                                     [sr_w][sr_h]),
                                    &(me_ctx->y_hme_level1_search_center[list_index][ref_pic_index]
                                                                     [sr_w][sr_h]));
                    }
                }
            }
        } // End ref pic loop
    } // End list loop
}

/*******************************************
 * performs hierarchical ME level 2 for one 64x64 block (uni-prediction only)
 *******************************************/
static void hme_level2_b64(PictureParentControlSet *pcs, uint32_t org_x, uint32_t org_y,
                           MeContext *me_ctx, EbPictureBufferDesc *input_ptr) {
    const uint32_t block_width  = me_ctx->b64_width;
    const uint32_t block_height = me_ctx->b64_height;
    // List Loop
    const uint8_t num_of_list_to_search = me_ctx->num_of_list_to_search;
    for (int list_index = REF_LIST_0; list_index < num_of_list_to_search; ++list_index) {
        // Ref Picture Loop
        const uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];
        for (uint8_t ref_pic_index = 0; ref_pic_index < num_of_ref_pic_to_search; ++ref_pic_index) {
            uint16_t             dist    = 0;
            EbPictureBufferDesc *ref_pic = get_me_reference(
                pcs, me_ctx, list_index, ref_pic_index, 2, &dist, input_ptr->width, input_ptr->height);

            if (me_ctx->temporal_layer_index > 0 || list_index == 0) {
                for (uint8_t sr_h = 0; sr_h < me_ctx->num_hme_sa_h; sr_h++) {
                    for (uint8_t sr_w = 0; sr_w < me_ctx->num_hme_sa_w; sr_w++) {
                        if (me_ctx->prev_me_stage_based_exit_th) {
                            if (me_ctx->hme_level1_sad[list_index][ref_pic_index][sr_w][sr_h] < (me_ctx->prev_me_stage_based_exit_th >> 2)) {
                                me_ctx->x_hme_level2_search_center[list_index][ref_pic_index][sr_w][sr_h] =
                                    me_ctx->x_hme_level1_search_center[list_index][ref_pic_index][sr_w][sr_h];
                                me_ctx->y_hme_level2_search_center[list_index][ref_pic_index][sr_w][sr_h] =
                                    me_ctx->y_hme_level1_search_center[list_index][ref_pic_index][sr_w][sr_h];
                                me_ctx->hme_level2_sad[list_index][ref_pic_index][sr_w][sr_h] =
                                    me_ctx->hme_level1_sad[list_index][ref_pic_index][sr_w][sr_h];
                                continue;
                            }
                        }

                        hme_level_2(
                            me_ctx,
                            (int16_t)org_x,
                            (int16_t)org_y,
                            block_width,
                            block_height,
                            ref_pic,
                            (int16_t)me_ctx->hme_l2_sa.width,
                            (int16_t)me_ctx->hme_l2_sa.height,
                            me_ctx->x_hme_level1_search_center[list_index][ref_pic_index][sr_w][sr_h],
                            me_ctx->y_hme_level1_search_center[list_index][ref_pic_index][sr_w][sr_h],
                            &(me_ctx->hme_level2_sad[list_index][ref_pic_index][sr_w][sr_h]),
                            &(me_ctx->x_hme_level2_search_center[list_index][ref_pic_index][sr_w]
                                                             [sr_h]),
                            &(me_ctx->y_hme_level2_search_center[list_index][ref_pic_index][sr_w]
                                                             [sr_h]));
                    }
                }
            }
        } // End ref pic loop
    } // End list loop
}
/*******************************************
 *   Set the final search centre
 *******************************************/

void set_final_seach_centre_sb(PictureParentControlSet *pcs, MeContext *me_ctx) {
    UNUSED(pcs);
    // Hierarchical ME Search Center
    int16_t xHmeSearchCenter = 0;
    int16_t yHmeSearchCenter = 0;

    // Final ME Search Center
    int16_t x_search_center = 0;
    int16_t y_search_center = 0;

    // Search Center SADs
    uint64_t hmeMvSad = 0;
    uint32_t num_of_list_to_search;
    uint32_t list_index;
    uint8_t  ref_pic_index;
    // Configure HME level 0, level 1 and level 2 from static config parameters
    Bool enable_hme_level0_flag = me_ctx->enable_hme_level0_flag;
    Bool enable_hme_level1_flag = me_ctx->enable_hme_level1_flag;
    Bool enable_hme_level2_flag = me_ctx->enable_hme_level2_flag;

    uint64_t best_cost         = (uint64_t)~0;
    me_ctx->best_list_idx = 0;
    me_ctx->best_ref_idx  = 0;
    num_of_list_to_search      = me_ctx->num_of_list_to_search;

    // Uni-Prediction motion estimation loop
    // List Loop
    for (list_index = REF_LIST_0; list_index < num_of_list_to_search; ++list_index) {
        uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];
        // Ref Picture Loop
        for (ref_pic_index = 0; ref_pic_index < num_of_ref_pic_to_search; ++ref_pic_index) {
            if (me_ctx->temporal_layer_index > 0 || list_index == 0) {
                if (me_ctx->enable_hme_flag) {
                    // Hierarchical ME - Search Center
                    if (enable_hme_level0_flag && !enable_hme_level1_flag &&
                        !enable_hme_level2_flag) {
                        xHmeSearchCenter =
                            me_ctx
                                ->x_hme_level0_search_center[list_index][ref_pic_index][0][0];
                        yHmeSearchCenter =
                            me_ctx
                                ->y_hme_level0_search_center[list_index][ref_pic_index][0][0];
                        hmeMvSad = me_ctx->hme_level0_sad[list_index][ref_pic_index][0][0];

                        uint32_t search_region_number_in_width  = 1;
                        uint32_t search_region_number_in_height = 0;
                        while (search_region_number_in_height < me_ctx->num_hme_sa_h) {
                            while (search_region_number_in_width < me_ctx->num_hme_sa_w) {
                                xHmeSearchCenter =
                                    (me_ctx->hme_level0_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->x_hme_level0_search_center
                                          [list_index][ref_pic_index][search_region_number_in_width]
                                          [search_region_number_in_height]
                                    : xHmeSearchCenter;
                                yHmeSearchCenter =
                                    (me_ctx->hme_level0_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->y_hme_level0_search_center
                                          [list_index][ref_pic_index][search_region_number_in_width]
                                          [search_region_number_in_height]
                                    : yHmeSearchCenter;
                                hmeMvSad =
                                    (me_ctx->hme_level0_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->hme_level0_sad[list_index][ref_pic_index]
                                                                 [search_region_number_in_width]
                                                                 [search_region_number_in_height]
                                    : hmeMvSad;
                                search_region_number_in_width++;
                            }
                            search_region_number_in_width = 0;
                            search_region_number_in_height++;
                        }
                    }

                    if (enable_hme_level1_flag && !enable_hme_level2_flag) {
                        xHmeSearchCenter =
                            me_ctx
                                ->x_hme_level1_search_center[list_index][ref_pic_index][0][0];
                        yHmeSearchCenter =
                            me_ctx
                                ->y_hme_level1_search_center[list_index][ref_pic_index][0][0];
                        hmeMvSad = me_ctx->hme_level1_sad[list_index][ref_pic_index][0][0];

                        uint32_t search_region_number_in_width  = 1;
                        uint32_t search_region_number_in_height = 0;
                        while (search_region_number_in_height < me_ctx->num_hme_sa_h) {
                            while (search_region_number_in_width < me_ctx->num_hme_sa_w) {
                                xHmeSearchCenter =
                                    (me_ctx->hme_level1_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->x_hme_level1_search_center
                                          [list_index][ref_pic_index][search_region_number_in_width]
                                          [search_region_number_in_height]
                                    : xHmeSearchCenter;
                                yHmeSearchCenter =
                                    (me_ctx->hme_level1_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->y_hme_level1_search_center
                                          [list_index][ref_pic_index][search_region_number_in_width]
                                          [search_region_number_in_height]
                                    : yHmeSearchCenter;
                                hmeMvSad =
                                    (me_ctx->hme_level1_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->hme_level1_sad[list_index][ref_pic_index]
                                                                 [search_region_number_in_width]
                                                                 [search_region_number_in_height]
                                    : hmeMvSad;
                                search_region_number_in_width++;
                            }
                            search_region_number_in_width = 0;
                            search_region_number_in_height++;
                        }
                    }

                    if (enable_hme_level2_flag) {
                        xHmeSearchCenter =
                            me_ctx
                                ->x_hme_level2_search_center[list_index][ref_pic_index][0][0];
                        yHmeSearchCenter =
                            me_ctx
                                ->y_hme_level2_search_center[list_index][ref_pic_index][0][0];
                        hmeMvSad = me_ctx->hme_level2_sad[list_index][ref_pic_index][0][0];

                        uint32_t search_region_number_in_width  = 1;
                        uint32_t search_region_number_in_height = 0;
                        while (search_region_number_in_height < me_ctx->num_hme_sa_h) {
                            while (search_region_number_in_width < me_ctx->num_hme_sa_w) {
                                xHmeSearchCenter =
                                    (me_ctx->hme_level2_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->x_hme_level2_search_center
                                          [list_index][ref_pic_index][search_region_number_in_width]
                                          [search_region_number_in_height]
                                    : xHmeSearchCenter;
                                yHmeSearchCenter =
                                    (me_ctx->hme_level2_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->y_hme_level2_search_center
                                          [list_index][ref_pic_index][search_region_number_in_width]
                                          [search_region_number_in_height]
                                    : yHmeSearchCenter;
                                hmeMvSad =
                                    (me_ctx->hme_level2_sad[list_index][ref_pic_index]
                                                                [search_region_number_in_width]
                                                                [search_region_number_in_height] <
                                     hmeMvSad)
                                    ? me_ctx->hme_level2_sad[list_index][ref_pic_index]
                                                                 [search_region_number_in_width]
                                                                 [search_region_number_in_height]
                                    : hmeMvSad;
                                search_region_number_in_width++;
                            }
                            search_region_number_in_width = 0;
                            search_region_number_in_height++;
                        }
                    }

                    x_search_center = xHmeSearchCenter;
                    y_search_center = yHmeSearchCenter;
                }
            } else {
                x_search_center = 0;
                y_search_center = 0;
            }

            //sc valid for all cases. 0,0 if hme not done.
            me_ctx->search_results[list_index][ref_pic_index].hme_sc_x = x_search_center;
            me_ctx->search_results[list_index][ref_pic_index].hme_sc_y = y_search_center;

            me_ctx->search_results[list_index][ref_pic_index].hme_sad =
                hmeMvSad; //this is not valid in all cases. only when HME is done, and when HMELevel2 is done
                          //also for base layer some references are redundant!!
            if (hmeMvSad < best_cost) {
                best_cost                  = hmeMvSad;
                me_ctx->best_list_idx = list_index;
                me_ctx->best_ref_idx  = ref_pic_index;
            }
        }
    }
}
// Initialize zz SAD array
static void init_zz_sad(PictureParentControlSet *pcs, MeContext *me_ctx, uint32_t org_x, uint32_t org_y) {
    const uint32_t block_width  = me_ctx->b64_width;
    const uint32_t block_height = me_ctx->b64_height;
    uint32_t best_zz_sad = MAX_U32;
    // List Loop
    for (int list_i = REF_LIST_0; list_i < me_ctx->num_of_list_to_search; ++list_i) {
        // Ref Picture Loop
        for (uint8_t ref_i = 0; ref_i < me_ctx->num_of_ref_pic_to_search[list_i]; ++ref_i) {
            if (me_ctx->temporal_layer_index > 0 || list_i == 0) {
                EbPictureBufferDesc *ref_pic = me_ctx->me_ds_ref_array[list_i][ref_i].picture_ptr;
                uint32_t             zz_sad  = get_zz_sad(
                        ref_pic, me_ctx, org_x, org_y, block_width, block_height);
                //normalize for incomplete b64
                zz_sad = (zz_sad * 64 * 64) / (block_width * block_height);
                me_ctx->zz_sad[list_i][ref_i] = zz_sad;
                best_zz_sad = MIN(best_zz_sad, zz_sad);
            }
        }
    }
    const uint32_t zz_th = me_ctx->me_hme_prune_ctrls.zz_sad_th;
    if (me_ctx->temporal_layer_index > 0 && best_zz_sad < zz_th) {
        for (int list_i = REF_LIST_0; list_i < me_ctx->num_of_list_to_search; ++list_i) {
            for (uint8_t ref_i = 0; ref_i < me_ctx->num_of_ref_pic_to_search[list_i]; ++ref_i) {

                if (ref_i == 0) continue;

                const uint32_t zz_sad_pct = me_ctx->me_hme_prune_ctrls.zz_sad_pct;
                if ((me_ctx->zz_sad[list_i][ref_i] - best_zz_sad) * 100 > (zz_sad_pct * best_zz_sad)) {
                       me_ctx->search_results[list_i][ref_i].do_ref = 0;
                }
            }
        }
    }

    const uint32_t safe_limit_zz_th = me_ctx->me_safe_limit_zz_th;
    if (safe_limit_zz_th) {
        bool me_safe_limit_refs = false;
        if (pcs->hierarchical_levels > 0 && me_ctx->num_of_list_to_search == 2 &&
            pcs->temporal_layer_index >= pcs->hierarchical_levels && pcs->similar_brightness_refs &&
            me_ctx->zz_sad[0][0] < safe_limit_zz_th &&  me_ctx->zz_sad[1][0] < safe_limit_zz_th
            ) {
            me_safe_limit_refs = true;
        }

        for (int list_i = REF_LIST_0; list_i < me_ctx->num_of_list_to_search; ++list_i) {
            for (uint8_t ref_i = 0; ref_i < me_ctx->num_of_ref_pic_to_search[list_i]; ++ref_i) {

                if (me_safe_limit_refs && ref_i > 0) {
                    me_ctx->search_results[list_i][ref_i].do_ref = 0;
                }
            }
        }
    }


}
/*******************************************
 * performs hierarchical ME for a 64x64 block for every ref frame
 *******************************************/
static void hme_b64(PictureParentControlSet *pcs, uint32_t org_x, uint32_t org_y,
             MeContext *me_ctx, EbPictureBufferDesc *input_ptr) {
    // If needed, initialize the zz sad array
    if (me_ctx->me_early_exit_th || me_ctx->me_safe_limit_zz_th  )
        init_zz_sad(pcs, me_ctx, org_x, org_y);

    if (me_ctx->prehme_ctrl.enable) {
        // perform pre-HME
        prehme_b64(pcs, org_x, org_y, me_ctx, input_ptr);
    }

    if (me_ctx->enable_hme_flag) {
        // perform hierarchical ME level 0
        if (me_ctx->enable_hme_level0_flag)
            hme_level0_b64(pcs, org_x, org_y, me_ctx, input_ptr);

        // perform hierarchical ME level 1
        if (me_ctx->enable_hme_level1_flag)
            hme_level1_b64(pcs, org_x, org_y, me_ctx, input_ptr);

        // perform hierarchical ME level 2
        if (me_ctx->enable_hme_level2_flag)
            hme_level2_b64(pcs, org_x, org_y, me_ctx, input_ptr);
    }

    // Set final MV centre
    set_final_seach_centre_sb(pcs, me_ctx);

    if (me_ctx->me_type == ME_MCTF) {
        if (ABS(me_ctx->search_results[0][0].hme_sc_x) > ABS(me_ctx->search_results[0][0].hme_sc_y))
            me_ctx->tf_tot_horz_blks++;
        else
            me_ctx->tf_tot_vert_blks++;
    }
}

static void hme_prune_ref_and_adjust_sr(MeContext *me_ctx) {
    uint64_t best = (uint64_t)~0;
    for (int i = 0; i < MAX_NUM_OF_REF_PIC_LIST; ++i) {
        for (int j = 0; j < REF_LIST_MAX_DEPTH; ++j) {
            if (me_ctx->search_results[i][j].hme_sad < best) {
                best = me_ctx->search_results[i][j].hme_sad;
            }
        }
    }
    uint16_t prune_ref_th = me_ctx->me_hme_prune_ctrls.prune_ref_if_hme_sad_dev_bigger_than_th;
    uint16_t mv_length_th = me_ctx->me_sr_adjustment_ctrls.reduce_me_sr_based_on_mv_length_th;
    uint16_t stationary_hme_sad_abs_th =
        me_ctx->me_sr_adjustment_ctrls.stationary_hme_sad_abs_th;
    uint16_t reduce_me_sr_based_on_hme_sad_abs_th =
        me_ctx->me_sr_adjustment_ctrls.reduce_me_sr_based_on_hme_sad_abs_th;
    for (uint32_t li = 0; li < MAX_NUM_OF_REF_PIC_LIST; li++) {
        for (uint32_t ri = 0; ri < REF_LIST_MAX_DEPTH; ri++) {
            // Prune references based on HME sad
            if (me_ctx->me_hme_prune_ctrls.enable_me_hme_ref_pruning &&
                (!me_ctx->me_hme_prune_ctrls.protect_closest_refs || ri > 0) &&
                (prune_ref_th != (uint16_t)~0) && ((me_ctx->search_results[li][ri].hme_sad - best) * 100 > (prune_ref_th * best))) {
                me_ctx->search_results[li][ri].do_ref = 0;
            }

            // Reduce the ME search region if the hme sad is low
            if (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment) {
                if (ABS(me_ctx->search_results[li][ri].hme_sc_x) <= mv_length_th &&
                    ABS(me_ctx->search_results[li][ri].hme_sc_y) <= mv_length_th &&
                    me_ctx->search_results[li][ri].hme_sad < stationary_hme_sad_abs_th) {
                    me_ctx->reduce_me_sr_divisor[li][ri] = me_ctx->me_sr_adjustment_ctrls.stationary_me_sr_divisor;
                }
                else if (me_ctx->search_results[li][ri].hme_sad < reduce_me_sr_based_on_hme_sad_abs_th) {
                    me_ctx->reduce_me_sr_divisor[li][ri] = me_ctx->me_sr_adjustment_ctrls.me_sr_divisor_for_low_hme_sad;
                }
            }
        }
    }
}

static const uint8_t z_to_raster[85] = {0,  1,  2,  3,  4,
                                 5,  6,  9,  10, 7,  8,  11, 12, 13,
                                 14, 17, 18, 15, 16, 19, 20,
                                 21, 22, 29, 30, 23, 24, 31, 32, 37,
                                 38, 45, 46, 39, 40, 47, 48,
                                 25, 26, 33, 34, 27, 28, 35, 36, 41,
                                 42, 49, 50, 43, 44, 51, 52,
                                 53, 54, 61, 62, 55, 56, 63, 64, 69,
                                 70, 77, 78, 71, 72, 79, 80,
                                 57, 58, 65, 66, 59, 60, 67, 68, 73,
                                 74, 81, 82, 75, 76, 83, 84
};
static void construct_me_candidate_array_mrp_off(PictureParentControlSet *pcs, MeContext *me_ctx,
                                          uint32_t num_of_list_to_search, uint32_t sb_index) {
    // This function should only be called if there is one ref frame in each list
    assert(me_ctx->num_of_ref_pic_to_search[0] == 1);
    assert(me_ctx->num_of_ref_pic_to_search[1] == 1);
    const uint8_t ref_pic_idx = 0;

    // Set whether the reference from each list is allowed
    uint8_t blk_do_ref_org[MAX_NUM_OF_REF_PIC_LIST];
    blk_do_ref_org[REF_LIST_0] = me_ctx->search_results[REF_LIST_0][0].do_ref;
    blk_do_ref_org[REF_LIST_1] = (num_of_list_to_search == 1) ? 0 : me_ctx->search_results[REF_LIST_1][0].do_ref;

    if (num_of_list_to_search < 2 || !me_ctx->search_results[REF_LIST_1][0].do_ref)
        num_of_list_to_search = 1;
    const uint32_t me_prune_th = (blk_do_ref_org[0] && blk_do_ref_org[1])
        ? me_ctx->prune_me_candidates_th
        : 0;

    // Set the count to 1 for all PUs using memset, which is faster than setting at the end of each loop.  The count will only need
    // to be updated if both reference frames are allowed.
    uint8_t number_of_pus = pcs->enable_me_16x16
        ? pcs->enable_me_8x8 ? pcs->max_number_of_pus_per_sb : MAX_SB64_PU_COUNT_NO_8X8
        : MAX_SB64_PU_COUNT_WO_16X16;
    memset(pcs->pa_me_data->me_results[sb_index]->total_me_candidate_index, 1, number_of_pus);

    for (uint8_t n_idx = 0; n_idx < pcs->max_number_of_pus_per_sb; ++n_idx) {
        const uint8_t pu_index       = z_to_raster[n_idx];
        uint8_t       me_cand_offset = 0;

        uint8_t      use_me_pu          = pcs->enable_me_16x16
                          ? pcs->enable_me_8x8 || n_idx < MAX_SB64_PU_COUNT_NO_8X8
                          : n_idx < MAX_SB64_PU_COUNT_WO_16X16;
        MeCandidate *me_candidate_array = NULL;
        if (use_me_pu)
            me_candidate_array =
                &pcs->pa_me_data->me_results[sb_index]
                     ->me_candidate_array[pu_index * pcs->pa_me_data->max_cand];
        uint8_t  blk_do_ref[MAX_NUM_OF_REF_PIC_LIST] = {blk_do_ref_org[REF_LIST_0],
                                                       blk_do_ref_org[REF_LIST_1]};
        const uint32_t best_me_dist =
            blk_do_ref_org[REF_LIST_0] && blk_do_ref_org[REF_LIST_1]
            ? MIN(me_ctx->p_sb_best_sad[REF_LIST_0][ref_pic_idx][n_idx], me_ctx->p_sb_best_sad[REF_LIST_1][ref_pic_idx][n_idx])
            : blk_do_ref_org[REF_LIST_0]
            ? me_ctx->p_sb_best_sad[REF_LIST_0][ref_pic_idx][n_idx]
            : me_ctx->p_sb_best_sad[REF_LIST_1][ref_pic_idx][n_idx];

        me_ctx->me_distortion[pu_index] = best_me_dist;
        int8_t min_dist_list = -1;
        // If both refs have a candidate, use only the best one for unipred
        if (me_ctx->use_best_unipred_cand_only && blk_do_ref[REF_LIST_0] &&
            blk_do_ref[REF_LIST_1])
            min_dist_list = me_ctx->p_sb_best_sad[REF_LIST_0][ref_pic_idx][n_idx] <
                    me_ctx->p_sb_best_sad[REF_LIST_1][ref_pic_idx][n_idx]
                ? 0
                : 1;
            // Unipred candidates
        for (int list_index = REF_LIST_0;
             (uint32_t)list_index < num_of_list_to_search && (use_me_pu || me_cand_offset == 0);
             ++list_index) {

            //ME was skipped, so do not add this Unipred candidate
            if (blk_do_ref[list_index] == 0)
                continue;

            if (me_prune_th > 0) {
                uint32_t current_to_best_dist_distance =
                    (me_ctx->p_sb_best_sad[list_index][ref_pic_idx][n_idx] - best_me_dist) *
                    100;
                if (current_to_best_dist_distance > (best_me_dist * me_prune_th)) {
                    blk_do_ref[list_index] = 0;
                    continue;
                }
            }
            if (min_dist_list != -1 && min_dist_list != list_index) {
                // Need to save the MV in case bipred is injected
                if (use_me_pu)
                    pcs->pa_me_data->me_results[sb_index]
                        ->me_mv_array[pu_index * pcs->pa_me_data->max_refs +
                                      (list_index ? pcs->pa_me_data->max_l0 : 0) + ref_pic_idx]
                        .as_int = me_ctx->p_sb_best_mv[list_index][ref_pic_idx][n_idx];
                continue;
            }
            if (use_me_pu) {
                me_candidate_array[me_cand_offset].direction  = list_index;
                me_candidate_array[me_cand_offset].ref_idx_l0 = ref_pic_idx;
                me_candidate_array[me_cand_offset].ref_idx_l1 = ref_pic_idx;
                me_candidate_array[me_cand_offset].ref0_list  = list_index == 0 ? list_index : 24;
                me_candidate_array[me_cand_offset].ref1_list  = list_index == 1 ? list_index : 24;

                pcs->pa_me_data->me_results[sb_index]
                    ->me_mv_array[pu_index * pcs->pa_me_data->max_refs +
                                  (list_index ? pcs->pa_me_data->max_l0 : 0) + ref_pic_idx]
                    .as_int = me_ctx->p_sb_best_mv[list_index][ref_pic_idx][n_idx];
            }

            me_cand_offset++;
        }

        // Can have up to one bipred cand (LAST ,BWD)
        if (blk_do_ref[REF_LIST_0] && blk_do_ref[REF_LIST_1] && use_me_pu) {
            // If get here, will have 3 candidates, since both unipred directions are valid
            assert(num_of_list_to_search == 2);
            me_candidate_array[me_cand_offset].direction  = BI_PRED;
            me_candidate_array[me_cand_offset].ref_idx_l0 = ref_pic_idx;
            me_candidate_array[me_cand_offset].ref_idx_l1 = ref_pic_idx;
            me_candidate_array[me_cand_offset].ref0_list  = REFERENCE_PIC_LIST_0;
            me_candidate_array[me_cand_offset].ref1_list  = REFERENCE_PIC_LIST_1;

            // store total me candidate count
            pcs->pa_me_data->me_results[sb_index]->total_me_candidate_index[pu_index] =
                me_cand_offset + 1;
        }
    }
}
static void construct_me_candidate_array_single_ref(PictureParentControlSet *pcs, MeContext *ctx,
    uint32_t num_of_list_to_search, uint32_t sb_index) {
    // This function should only be called if there is one ref frame in list 0
    assert(ctx->num_of_ref_pic_to_search[0] == 1);
    assert(ctx->num_of_ref_pic_to_search[1] == 0);
    const uint8_t ref_pic_idx = 0;

    // Set whether the reference from each list is allowed
    uint8_t blk_do_ref = ctx->search_results[REF_LIST_0][0].do_ref;

    if (num_of_list_to_search < 2 || !ctx->search_results[REF_LIST_1][0].do_ref)
        num_of_list_to_search = 1;

    // Set the count to 1 for all PUs using memset, which is faster than setting at the end of each loop.  The count will only need
    // to be updated if both reference frames are allowed.
    uint8_t number_of_pus = pcs->enable_me_16x16
        ? pcs->enable_me_8x8 ? pcs->max_number_of_pus_per_sb : MAX_SB64_PU_COUNT_NO_8X8
        : MAX_SB64_PU_COUNT_WO_16X16;
    memset(pcs->pa_me_data->me_results[sb_index]->total_me_candidate_index, 1, number_of_pus);

    for (uint8_t n_idx = 0; n_idx < pcs->max_number_of_pus_per_sb; ++n_idx) {
        const uint8_t pu_index       = z_to_raster[n_idx];

        uint8_t      use_me_pu          = pcs->enable_me_16x16
            ? pcs->enable_me_8x8 || n_idx < MAX_SB64_PU_COUNT_NO_8X8
            : n_idx < MAX_SB64_PU_COUNT_WO_16X16;
        MeCandidate *me_candidate_array = NULL;
        if (use_me_pu)
            me_candidate_array =
            &pcs->pa_me_data->me_results[sb_index]
            ->me_candidate_array[pu_index * pcs->pa_me_data->max_cand];
        ctx->me_distortion[pu_index] = ctx->p_sb_best_sad[REF_LIST_0][ref_pic_idx][n_idx];;

        //ME was skipped, so do not add this Unipred candidate
        if (blk_do_ref == 0)
            continue;


        if (use_me_pu) {
            me_candidate_array[0].direction  = REF_LIST_0;
            me_candidate_array[0].ref_idx_l0 = ref_pic_idx;
            me_candidate_array[0].ref_idx_l1 = ref_pic_idx;
            me_candidate_array[0].ref0_list  = 0;
            me_candidate_array[0].ref1_list  = 0;

            pcs->pa_me_data->me_results[sb_index]
                ->me_mv_array[pu_index * pcs->pa_me_data->max_refs + ref_pic_idx]
                .as_int = ctx->p_sb_best_mv[0][ref_pic_idx][n_idx];
        }

    }
}
static void construct_me_candidate_array(PictureParentControlSet *pcs, MeContext *me_ctx,
                                  uint32_t num_of_list_to_search, uint32_t sb_index) {
    for (uint32_t n_idx = 0; n_idx < pcs->max_number_of_pus_per_sb; ++n_idx) {
        uint8_t pu_index       = (n_idx > 4) ? z_to_raster[n_idx] : n_idx;
        uint8_t me_cand_offset = 0;

        uint8_t      use_me_pu          = pcs->enable_me_16x16
                          ? pcs->enable_me_8x8 || n_idx < MAX_SB64_PU_COUNT_NO_8X8
                          : n_idx < MAX_SB64_PU_COUNT_WO_16X16;
        MeCandidate *me_candidate_array = NULL;
        if (use_me_pu)
            me_candidate_array =
                &pcs->pa_me_data->me_results[sb_index]
                     ->me_candidate_array[pu_index * pcs->pa_me_data->max_cand];
        uint8_t  blk_do_ref[MAX_NUM_OF_REF_PIC_LIST][MAX_REF_IDX];
        uint32_t current_to_best_dist_distance;
        const uint32_t me_prune_th = me_ctx->prune_me_candidates_th; //to change to 32bit
        uint32_t best_me_dist = (uint32_t)~0;

        // Determine the best ME distortion
        for (uint32_t list_index = REF_LIST_0; list_index < num_of_list_to_search; list_index++) {
            const uint8_t num_of_ref_pic_to_search = me_ctx->num_of_ref_pic_to_search[list_index];
            for (uint32_t ref_pic = 0; ref_pic < num_of_ref_pic_to_search; ref_pic++) {
                blk_do_ref[list_index][ref_pic] = me_ctx->search_results[list_index][ref_pic].do_ref;
                if (blk_do_ref[list_index][ref_pic] == 0)
                    continue;

                best_me_dist = me_ctx->p_sb_best_sad[list_index][ref_pic][n_idx] < best_me_dist
                    ? me_ctx->p_sb_best_sad[list_index][ref_pic][n_idx] : best_me_dist;
            }
        }

        me_ctx->me_distortion[pu_index] = best_me_dist;
        // Unipred candidates
        for (uint32_t list_index = REF_LIST_0;
             list_index < num_of_list_to_search && (use_me_pu || me_cand_offset == 0);
             ++list_index) {
            const uint8_t num_of_ref_pic_to_search =
                me_ctx->num_of_ref_pic_to_search[list_index];

            for (uint32_t ref_pic_index = 0; (ref_pic_index < num_of_ref_pic_to_search) && (use_me_pu || (me_cand_offset == 0)); ++ref_pic_index) {
                //ME was skipped, so do not add this Unipred candidate
                if (blk_do_ref[list_index][ref_pic_index] == 0)
                    continue;

                if (me_prune_th > 0) {
                    current_to_best_dist_distance =
                        (me_ctx->p_sb_best_sad[list_index][ref_pic_index][n_idx] -
                         best_me_dist) *
                        100;
                    if (current_to_best_dist_distance > (best_me_dist * me_prune_th)) {
                        blk_do_ref[list_index][ref_pic_index] = 0;
                        continue;
                    }
                }
                if (use_me_pu) {
                    me_candidate_array[me_cand_offset].direction  = list_index;
                    me_candidate_array[me_cand_offset].ref_idx_l0 = ref_pic_index;
                    me_candidate_array[me_cand_offset].ref_idx_l1 = ref_pic_index;
                    me_candidate_array[me_cand_offset].ref0_list  = list_index == 0 ? list_index
                                                                                    : 24;
                    me_candidate_array[me_cand_offset].ref1_list  = list_index == 1 ? list_index
                                                                                    : 24;

                    pcs->pa_me_data->me_results[sb_index]
                        ->me_mv_array[pu_index * pcs->pa_me_data->max_refs +
                                      (list_index ? pcs->pa_me_data->max_l0 : 0) +
                                      ref_pic_index]
                        .as_int = me_ctx->p_sb_best_mv[list_index][ref_pic_index][n_idx];
                }
                me_cand_offset++;
            }
        }
        if (num_of_list_to_search == 2 && use_me_pu) {
            // 1st set of BIPRED cand
            // (LAST ,BWD), (LAST,ALT ), (LAST,ALT2 )
            // (LAST2,BWD), (LAST2,ALT), (LAST2,ALT2)
            // (LAST3,BWD), (LAST3,ALT), (LAST3,ALT2)
            // (GOLD ,BWD), (GOLD,ALT ), (GOLD,ALT2 )
            for (uint32_t first_list_ref_pict_idx = 0;
                 first_list_ref_pict_idx < me_ctx->num_of_ref_pic_to_search[REF_LIST_0];
                 first_list_ref_pict_idx++) {
                for (uint32_t second_list_ref_pict_idx = 0;
                     second_list_ref_pict_idx < me_ctx->num_of_ref_pic_to_search[REF_LIST_1];
                     second_list_ref_pict_idx++) {
                    if (pcs->scs->mrp_ctrls.only_l_bwd && (first_list_ref_pict_idx > 0 || second_list_ref_pict_idx > 0))
                        continue;
                    if (blk_do_ref[REF_LIST_0][first_list_ref_pict_idx] &&
                        blk_do_ref[REF_LIST_1][second_list_ref_pict_idx]) {
                        me_candidate_array[me_cand_offset].direction  = BI_PRED;
                        me_candidate_array[me_cand_offset].ref_idx_l0 = first_list_ref_pict_idx;
                        me_candidate_array[me_cand_offset].ref_idx_l1 = second_list_ref_pict_idx;
                        me_candidate_array[me_cand_offset].ref0_list  = REFERENCE_PIC_LIST_0;
                        me_candidate_array[me_cand_offset].ref1_list  = REFERENCE_PIC_LIST_1;
                        me_cand_offset++;
                    }
                }
            }
            if (!pcs->scs->mrp_ctrls.only_l_bwd) {
            // 2nd set of BIPRED cand: (LAST,LAST2) (LAST,LAST3) (LAST,GOLD)
            for (uint32_t first_list_ref_pict_idx = 1;
                 first_list_ref_pict_idx < me_ctx->num_of_ref_pic_to_search[REF_LIST_0];
                 first_list_ref_pict_idx++) {
                if (blk_do_ref[REF_LIST_0][0] && blk_do_ref[REF_LIST_0][first_list_ref_pict_idx]) {
                    me_candidate_array[me_cand_offset].direction  = BI_PRED;
                    me_candidate_array[me_cand_offset].ref_idx_l0 = 0;
                    me_candidate_array[me_cand_offset].ref_idx_l1 = first_list_ref_pict_idx;
                    me_candidate_array[me_cand_offset].ref0_list  = REFERENCE_PIC_LIST_0;
                    me_candidate_array[me_cand_offset].ref1_list  = REFERENCE_PIC_LIST_0;
                    me_cand_offset++;
                }
            }
            }

            // 3rd set of BIPRED cand: (BWD, ALT)
            if (!pcs->scs->mrp_ctrls.only_l_bwd) {

            if (me_ctx->num_of_ref_pic_to_search[REF_LIST_1] == 3 &&
                blk_do_ref[REF_LIST_1][0] && blk_do_ref[REF_LIST_1][2]) {
                {
                    me_candidate_array[me_cand_offset].direction  = BI_PRED;
                    me_candidate_array[me_cand_offset].ref_idx_l0 = 0;
                    me_candidate_array[me_cand_offset].ref_idx_l1 = 2;
                    me_candidate_array[me_cand_offset].ref0_list  = REFERENCE_PIC_LIST_1;
                    me_candidate_array[me_cand_offset].ref1_list  = REFERENCE_PIC_LIST_1;
                    me_cand_offset++;
                }
            }
            }
        }


        // store total me candidate count
        if (use_me_pu)
            pcs->pa_me_data->me_results[sb_index]->total_me_candidate_index[pu_index] =
                me_cand_offset;
    }
}

// Active and stationary detection for global motion
static void perform_gm_detection(
    PictureParentControlSet *pcs, // input parameter, Picture Control Set Ptr
    uint32_t                 sb_index, // input parameter, SB Index
    MeContext
        *me_ctx // input parameter, ME Context Ptr, used to store decimated/interpolated SB/SR
) {
    SequenceControlSet *scs        = pcs->scs;
    uint64_t            stationary_cnt = 0;
    uint64_t            per_sig_cnt[MAX_NUM_OF_REF_PIC_LIST][REF_LIST_MAX_DEPTH][NUM_MV_COMPONENTS]
                        [NUM_MV_HIST];
    uint64_t tot_cnt = 0;
    memset(per_sig_cnt, 0, sizeof(uint64_t) * MAX_MV_HIST_SIZE);

    if (scs->input_resolution <= INPUT_SIZE_480p_RANGE) {
        for (unsigned i = 0; i < 64; i++) {
            uint8_t n_idx = 21 + i;
            if (!pcs->enable_me_8x8) {
                if (n_idx >= MAX_SB64_PU_COUNT_NO_8X8)
                    n_idx = me_idx_85_8x8_to_16x16_conversion[n_idx - MAX_SB64_PU_COUNT_NO_8X8];
                if (!pcs->enable_me_16x16)
                    if (n_idx >= MAX_SB64_PU_COUNT_WO_16X16)
                        n_idx = me_idx_16x16_to_parent_32x32_conversion[n_idx -
                                                                        MAX_SB64_PU_COUNT_WO_16X16];
            }
            MeCandidate *me_candidate = &(
                pcs->pa_me_data->me_results[sb_index]
                    ->me_candidate_array[n_idx * pcs->pa_me_data->max_cand]);

            uint32_t list_index    = (me_candidate->direction == 0 || me_candidate->direction == 2)
                   ? me_candidate->ref0_list
                   : me_candidate->ref1_list;
            uint32_t ref_pic_index = (me_candidate->direction == 0 || me_candidate->direction == 2)
                ? me_candidate->ref_idx_l0
                : me_candidate->ref_idx_l1;

            // Active block detection
            uint64_t pcs_pic_num = pcs->picture_number;
            uint64_t ref_pic_num =
                me_ctx->me_ds_ref_array[list_index][ref_pic_index].picture_number;
            uint16_t dist = ABS(
                (int16_t)(MAX(pcs_pic_num, ref_pic_num) - MIN(pcs_pic_num, ref_pic_num)));
            int active_th = (pcs->gm_ctrls.use_distance_based_active_th) ? MAX(dist >> 1, 4)
                                                                             : 4;

            int mx = _MVXT(me_ctx->p_sb_best_mv[list_index][ref_pic_index][n_idx]) << 2;
            if (mx < -active_th)
                per_sig_cnt[list_index][ref_pic_index][0][0]++;
            else if (mx > active_th)
                per_sig_cnt[list_index][ref_pic_index][0][1]++;
            int my = _MVYT(me_ctx->p_sb_best_mv[list_index][ref_pic_index][n_idx]) << 2;
            if (my < -active_th)
                per_sig_cnt[list_index][ref_pic_index][1][0]++;
            else if (my > active_th)
                per_sig_cnt[list_index][ref_pic_index][1][1]++;

            // Stationary block detection
            int stationary_th = 0;
            if (abs(mx) <= stationary_th && abs(my) <= stationary_th)
                stationary_cnt++;

            tot_cnt++;
        }
    } else {
        for (unsigned i = 0; i < 16; i++) {
            uint8_t n_idx = 5 + i;
            if (!pcs->enable_me_16x16)
                if (n_idx >= MAX_SB64_PU_COUNT_WO_16X16)
                    n_idx =
                        me_idx_16x16_to_parent_32x32_conversion[n_idx - MAX_SB64_PU_COUNT_WO_16X16];
            MeCandidate *me_candidate = &(
                pcs->pa_me_data->me_results[sb_index]
                    ->me_candidate_array[n_idx * pcs->pa_me_data->max_cand]);

            uint32_t list_index    = (me_candidate->direction == 0 || me_candidate->direction == 2)
                   ? me_candidate->ref0_list
                   : me_candidate->ref1_list;
            uint32_t ref_pic_index = (me_candidate->direction == 0 || me_candidate->direction == 2)
                ? me_candidate->ref_idx_l0
                : me_candidate->ref_idx_l1;

            // Active block detection
            uint16_t dist = ABS(
                (int16_t)(pcs->picture_number -
                          me_ctx->me_ds_ref_array[list_index][ref_pic_index].picture_number));
            int active_th = (pcs->gm_ctrls.use_distance_based_active_th) ? MAX(dist * 16, 32)
                                                                             : 32;

            int mx = _MVXT(me_ctx->p_sb_best_mv[list_index][ref_pic_index][n_idx]) << 2;
            if (mx < -active_th)
                per_sig_cnt[list_index][ref_pic_index][0][0]++;
            else if (mx > active_th)
                per_sig_cnt[list_index][ref_pic_index][0][1]++;
            int my = _MVYT(me_ctx->p_sb_best_mv[list_index][ref_pic_index][n_idx]) << 2;
            if (my < -active_th)
                per_sig_cnt[list_index][ref_pic_index][1][0]++;
            else if (my > active_th)
                per_sig_cnt[list_index][ref_pic_index][1][1]++;

            // Stationary block detection
            int stationary_th = 4;
            if (abs(mx) <= stationary_th && abs(my) <= stationary_th)
                stationary_cnt++;

            tot_cnt++;
        }
    }

    // Set stationary_block_present_sb to 1 if stationary_cnt is higher than 5%
    if (stationary_cnt > ((tot_cnt * 5) / 100))
        pcs->stationary_block_present_sb[sb_index] = 1;

    for (int l = 0; l < MAX_NUM_OF_REF_PIC_LIST; l++) {
        for (int r = 0; r < REF_LIST_MAX_DEPTH; r++) {
            for (int c = 0; c < NUM_MV_COMPONENTS; c++) {
                for (int s = 0; s < NUM_MV_HIST; s++) {
                    if (per_sig_cnt[l][r][c][s] > (tot_cnt / 2)) {
                        pcs->rc_me_allow_gm[sb_index] = 1;
                        break;
                    }
                }
            }
        }
    }
}

// Compute the distortion per block size based on the ME results
static void compute_distortion(
    PictureParentControlSet *pcs, // input parameter, Picture Control Set Ptr
    uint32_t                 b64_index, // input parameter, B64 Index
    MeContext
        *me_ctx // input parameter, ME Context Ptr, used to store decimated/interpolated SB/SR
) {
    SequenceControlSet *scs = pcs->scs;
    // Determine sb_64x64_me_class
    B64Geom *b64_geom  = &pcs->b64_geom[b64_index];
    uint32_t  b64_size    = 64 * 64;
    uint32_t  dist_64x64 = 0, dist_32x32 = 0, dist_16x16 = 0, dist_8x8 = 0;

    // 64x64
    { dist_64x64 = me_ctx->me_distortion[0]; }

    // 32x32
    for (unsigned i = 0; i < 4; i++) { dist_32x32 += me_ctx->me_distortion[1 + i]; }

    // 16x16
    for (unsigned i = 0; i < 16; i++) { dist_16x16 += me_ctx->me_distortion[5 + i]; }

    // 8x8
    for (unsigned i = 0; i < 64; i++) { dist_8x8 += me_ctx->me_distortion[21 + i]; }

    uint64_t mean_dist_8x8     = dist_8x8 / 64;
    uint64_t sum_ofsq_dist_8x8 = 0;
    for (unsigned i = 0; i < 64; i++) {
        const int64_t diff = ((int64_t)me_ctx->me_distortion[21 + i] - (int64_t)mean_dist_8x8);
        sum_ofsq_dist_8x8 += diff * diff;
    }

    pcs->me_8x8_cost_variance[b64_index] = (uint32_t)(sum_ofsq_dist_8x8 / 64);
    // Compute the sum of the distortion of all 16 16x16 (720 and above) and
    // 64 8x8 (for lower resolutions) blocks in the SB
    pcs->rc_me_distortion[b64_index] = (scs->input_resolution <= INPUT_SIZE_480p_RANGE)
        ? dist_8x8
        : dist_16x16;
    const uint32_t pix_num              = b64_geom->width * b64_geom->height;
    // Normalize
    pcs->me_64x64_distortion[b64_index] = (dist_64x64 * b64_size) / (pix_num);
    pcs->me_32x32_distortion[b64_index] = (dist_32x32 * b64_size) / (pix_num);
    pcs->me_16x16_distortion[b64_index] = (dist_16x16 * b64_size) / (pix_num);
    pcs->me_8x8_distortion[b64_index]   = (dist_8x8 * b64_size) / (pix_num);
}

// Initalize data used in ME/HME
static INLINE void init_me_hme_data(MeContext *me_ctx) {
    // Initialize HME search centres to 0
    if (me_ctx->enable_hme_flag) {
        memset(me_ctx->x_hme_level0_search_center,
               0,
               sizeof(me_ctx->x_hme_level0_search_center[0][0][0][0]) *
                   MAX_NUM_OF_REF_PIC_LIST * MAX_REF_IDX * EB_HME_SEARCH_AREA_COLUMN_MAX_COUNT *
                   EB_HME_SEARCH_AREA_ROW_MAX_COUNT);
        memset(me_ctx->y_hme_level0_search_center,
               0,
               sizeof(me_ctx->y_hme_level0_search_center[0][0][0][0]) *
                   MAX_NUM_OF_REF_PIC_LIST * MAX_REF_IDX * EB_HME_SEARCH_AREA_COLUMN_MAX_COUNT *
                   EB_HME_SEARCH_AREA_ROW_MAX_COUNT);

        memset(me_ctx->x_hme_level1_search_center,
               0,
               sizeof(me_ctx->x_hme_level1_search_center[0][0][0][0]) *
                   MAX_NUM_OF_REF_PIC_LIST * MAX_REF_IDX * EB_HME_SEARCH_AREA_COLUMN_MAX_COUNT *
                   EB_HME_SEARCH_AREA_ROW_MAX_COUNT);
        memset(me_ctx->y_hme_level1_search_center,
               0,
               sizeof(me_ctx->y_hme_level1_search_center[0][0][0][0]) *
                   MAX_NUM_OF_REF_PIC_LIST * MAX_REF_IDX * EB_HME_SEARCH_AREA_COLUMN_MAX_COUNT *
                   EB_HME_SEARCH_AREA_ROW_MAX_COUNT);

        memset(me_ctx->x_hme_level2_search_center,
               0,
               sizeof(me_ctx->x_hme_level2_search_center[0][0][0][0]) *
                   MAX_NUM_OF_REF_PIC_LIST * MAX_REF_IDX * EB_HME_SEARCH_AREA_COLUMN_MAX_COUNT *
                   EB_HME_SEARCH_AREA_ROW_MAX_COUNT);
        memset(me_ctx->y_hme_level2_search_center,
               0,
               sizeof(me_ctx->y_hme_level2_search_center[0][0][0][0]) *
                   MAX_NUM_OF_REF_PIC_LIST * MAX_REF_IDX * EB_HME_SEARCH_AREA_COLUMN_MAX_COUNT *
                   EB_HME_SEARCH_AREA_ROW_MAX_COUNT);
    }

    // R2R FIX: no winner integer MV is set in special case like initial p_sb_best_mv for overlay case,
    // then it sends dirty p_sb_best_mv to MD, initializing it is necessary
    memset(me_ctx->p_sb_best_mv,
           0,
           sizeof(me_ctx->p_sb_best_mv[0][0][0]) * MAX_NUM_OF_REF_PIC_LIST *
               REF_LIST_MAX_DEPTH * SQUARE_PU_COUNT);

    //init hme results buffer
    for (uint32_t li = 0; li < MAX_NUM_OF_REF_PIC_LIST; li++) {
        for (uint32_t ri = 0; ri < REF_LIST_MAX_DEPTH; ri++) {
            if (me_ctx->me_type != ME_MCTF)
                me_ctx->search_results[li][ri].list_i = li;
            me_ctx->search_results[li][ri].ref_i    = ri;
            me_ctx->search_results[li][ri].do_ref   = 1;
            me_ctx->search_results[li][ri].hme_sad  = MAX_U32;
            me_ctx->reduce_me_sr_divisor[li][ri] = 1;
            me_ctx->zz_sad[li][ri] = (uint32_t)~0;
            me_ctx->prehme_data[li][ri][0].valid = 0;
            me_ctx->prehme_data[li][ri][1].valid = 0;
        }
    }
    memset(me_ctx->performed_phme,0,
        sizeof(uint8_t)*MAX_NUM_OF_REF_PIC_LIST*REF_LIST_MAX_DEPTH*SEARCH_REGION_COUNT);
}
/*******************************************
* motion_estimation
*   performs ME on 64x64 blocks
*******************************************/

EbErrorType svt_aom_motion_estimation_b64(
    PictureParentControlSet *pcs, // input parameter, Picture Control Set Ptr
    uint32_t                 b64_index, // input parameter, SB Index
    uint32_t                 b64_origin_x, // input parameter, SB Origin X
    uint32_t                 b64_origin_y, // input parameter, SB Origin X
    MeContext
        *me_ctx, // input parameter, ME Context Ptr, used to store decimated/interpolated SB/SR
    EbPictureBufferDesc *input_ptr) // input parameter, source Picture Ptr

{
    EbErrorType return_error = EB_ErrorNone;

    uint32_t num_of_list_to_search = me_ctx->num_of_list_to_search;

    // input picture width and height might be disaligned after resizing
    // we use aligned width and height to avoid disalignment of calculation
    // of block size
    uint16_t aligned_width  = (uint16_t)ALIGN_POWER_OF_TWO(input_ptr->width, 3);
    uint16_t aligned_height = (uint16_t)ALIGN_POWER_OF_TWO(input_ptr->height, 3);
    me_ctx->b64_width  = (aligned_width - b64_origin_x) < BLOCK_SIZE_64
                              ? aligned_width - b64_origin_x
                              : BLOCK_SIZE_64;
    me_ctx->b64_height = (aligned_height - b64_origin_y) < BLOCK_SIZE_64
                              ? aligned_height - b64_origin_y
                              : BLOCK_SIZE_64;

    //pruning of the references is not done for alt-ref / when HMeLevel2 not done
    uint8_t prune_ref = me_ctx->enable_hme_flag && me_ctx->me_type != ME_MCTF;
    // Initialize ME/HME buffers
    init_me_hme_data(me_ctx);
    // HME: Perform Hierachical Motion Estimation for all refrence frames for the current 64x64 block.
    hme_b64(pcs, b64_origin_x, b64_origin_y, me_ctx, input_ptr);

    if (me_ctx->me_type == ME_MCTF &&
        me_ctx->search_results[0][0].hme_sad < me_ctx->tf_me_exit_th) {
        me_ctx->tf_use_pred_64x64_only_th = (uint8_t)~0;
        return return_error;
    }
    // prune the refrence frames based on the HME outputs.
    if (prune_ref &&
        (me_ctx->me_sr_adjustment_ctrls.enable_me_sr_adjustment ||
         me_ctx->me_hme_prune_ctrls.enable_me_hme_ref_pruning)) {
        hme_prune_ref_and_adjust_sr(me_ctx);
    }
    // Full pel: Perform the Integer Motion Estimation on the allowed refrence frames.
    integer_search_b64(pcs, b64_index, b64_origin_x, b64_origin_y, me_ctx, input_ptr);

    // prune the refrence frames
    if (prune_ref && me_ctx->me_hme_prune_ctrls.enable_me_hme_ref_pruning) {
        me_prune_ref(me_ctx);
    }

    if (me_ctx->me_type != ME_MCTF) {
        {
            if (me_ctx->num_of_ref_pic_to_search[REF_LIST_0] == 1 &&
                me_ctx->num_of_ref_pic_to_search[REF_LIST_1] == 0)
                construct_me_candidate_array_single_ref(
                    pcs, me_ctx, num_of_list_to_search, b64_index);
            else if (me_ctx->num_of_ref_pic_to_search[REF_LIST_0] == 1 &&
                me_ctx->num_of_ref_pic_to_search[REF_LIST_1] == 1)
                construct_me_candidate_array_mrp_off(
                    pcs, me_ctx, num_of_list_to_search, b64_index);
            else
                construct_me_candidate_array(pcs,
                    me_ctx,
                    num_of_list_to_search,
                    b64_index);
        }
            // Save the distortion per block size
            compute_distortion(pcs, b64_index, me_ctx);

        // Perform GM detection if GM is enabled
        pcs->stationary_block_present_sb[b64_index] = 0;
        pcs->rc_me_allow_gm[b64_index]              = 0;

        if (pcs->gm_ctrls.enabled)
            perform_gm_detection(pcs, b64_index, me_ctx);
    }
    return return_error;
}

EbErrorType svt_aom_open_loop_intra_search_mb(PictureParentControlSet *pcs, uint32_t b64_index,
                                      EbPictureBufferDesc *input_ptr) {
    EbErrorType         return_error = EB_ErrorNone;
    SequenceControlSet *scs      = pcs->scs;

    uint32_t      cu_origin_x;
    uint32_t      cu_origin_y;
    uint32_t      pa_blk_index = 0;
    B64Geom *    b64_geom    = &scs->b64_geom[b64_index];
    OisMbResults *ois_mb_results_ptr;
    uint8_t *     above_row;
    uint8_t *     left_col;
    uint8_t *     above0_row;
    uint8_t *     left0_col;
    uint32_t      mb_stride = (scs->max_input_luma_width + 15) / 16;

    DECLARE_ALIGNED(16, uint8_t, left0_data[MAX_TX_SIZE * 2 + 32]);
    DECLARE_ALIGNED(16, uint8_t, above0_data[MAX_TX_SIZE * 2 + 32]);
    DECLARE_ALIGNED(16, uint8_t, left_data[MAX_TX_SIZE * 2 + 32]);
    DECLARE_ALIGNED(16, uint8_t, above_data[MAX_TX_SIZE * 2 + 32]);

    DECLARE_ALIGNED(32, uint8_t, predictor8[256 * 2]);
    DECLARE_ALIGNED(32, int16_t, src_diff[256]);
    DECLARE_ALIGNED(32, int32_t, coeff[256]);
    uint8_t *predictor = predictor8;

    while (pa_blk_index < CU_MAX_COUNT) {
        const CodedBlockStats *blk_stats_ptr;
        blk_stats_ptr              = svt_aom_get_coded_blk_stats(pa_blk_index);
        uint8_t bsize              = blk_stats_ptr->size;
        Bool  small_boundary_blk = FALSE;

        //if(sb_params->raster_scan_blk_validity[md_scan_to_raster_scan[pa_blk_index]])
        {
            cu_origin_x = b64_geom->org_x + blk_stats_ptr->org_x;
            cu_origin_y = b64_geom->org_y + blk_stats_ptr->org_y;
            if ((blk_stats_ptr->org_x % 16) == 0 && (blk_stats_ptr->org_y % 16) == 0 &&
                ((pcs->enhanced_pic->width - cu_origin_x) < 16 ||
                 (pcs->enhanced_pic->height - cu_origin_y) < 16))
                small_boundary_blk = TRUE;
        }

        if (bsize != 16 && !small_boundary_blk) {
            pa_blk_index++;
            continue;
        }
        if (b64_geom->raster_scan_blk_validity[md_scan_to_raster_scan[pa_blk_index]]) {
            // always process as block16x16 even bsize or tx_size is 8x8
            TxSize tx_size = TX_16X16;
            bsize          = 16;
            cu_origin_x    = b64_geom->org_x + blk_stats_ptr->org_x;
            cu_origin_y    = b64_geom->org_y + blk_stats_ptr->org_y;
            above0_row     = above0_data + 16;
            left0_col      = left0_data + 16;
            above_row      = above_data + 16;
            left_col       = left_data + 16;
            ois_mb_results_ptr =
                pcs->pa_me_data
                    ->ois_mb_results[(cu_origin_y >> 4) * mb_stride + (cu_origin_x >> 4)];
            memset(ois_mb_results_ptr, 0, sizeof(*ois_mb_results_ptr));
            uint8_t *src = input_ptr->buffer_y + pcs->enhanced_pic->org_x +
                cu_origin_x +
                (pcs->enhanced_pic->org_y + cu_origin_y) * input_ptr->stride_y;

            // Fill Neighbor Arrays
            svt_aom_update_neighbor_samples_array_open_loop_mb(1, // use_top_righ_bottom_left
                                                       1, // update_top_neighbor
                                                       above0_row - 1,
                                                       left0_col - 1,
                                                       input_ptr,
                                                       input_ptr->stride_y,
                                                       cu_origin_x,
                                                       cu_origin_y,
                                                       bsize,
                                                       bsize);
            uint8_t        ois_intra_mode;
            uint8_t        intra_mode_start = DC_PRED;
            uint8_t intra_mode_end = pcs->tpl_ctrls.intra_mode_end;
            PredictionMode best_mode      = DC_PRED;
            int64_t        best_intra_cost = INT64_MAX;

            for (ois_intra_mode = intra_mode_start; ois_intra_mode <= intra_mode_end;
                 ++ois_intra_mode) {
                int32_t p_angle = av1_is_directional_mode((PredictionMode)ois_intra_mode)
                    ? mode_to_angle_map[(PredictionMode)ois_intra_mode]
                    : 0;
                // Edge filter
                if (av1_is_directional_mode((PredictionMode)ois_intra_mode) &&
                    1 /*scs->seq_header.enable_intra_edge_filter*/) {
                    EB_MEMCPY(left_data, left0_data, sizeof(uint8_t) * (MAX_TX_SIZE * 2 + 32));
                    EB_MEMCPY(above_data, above0_data, sizeof(uint8_t) * (MAX_TX_SIZE * 2 + 32));
                    above_row = above_data + 16;
                    left_col  = left_data + 16;
                    filter_intra_edge(ois_mb_results_ptr,
                                      ois_intra_mode,
                                      scs->max_input_luma_width,
                                      scs->max_input_luma_height,
                                      p_angle,
                                      (int32_t)cu_origin_x,
                                      (int32_t)cu_origin_y,
                                      above_row,
                                      left_col);
                } else {
                    above_row = above0_row;
                    left_col  = left0_col;
                }
                // PRED
                svt_aom_intra_prediction_open_loop_mb(p_angle,
                                              ois_intra_mode,
                                              cu_origin_x,
                                              cu_origin_y,
                                              tx_size,
                                              above_row,
                                              left_col,
                                              predictor,
                                              16);
                // Distortion
                int64_t intra_cost;
                if (pcs->tpl_ctrls.use_sad_in_src_search) {
                    intra_cost = svt_nxm_sad_kernel_sub_sampled(
                        src, input_ptr->stride_y, predictor, 16, 16, 16);
                } else {
                    TxSize cost_tx_size = pcs->tpl_ctrls.subsample_tx ? TX_16X8 : TX_16X16;
                    svt_aom_subtract_block(16 >> pcs->tpl_ctrls.subsample_tx,
                                           16,
                                           src_diff,
                                           16 << pcs->tpl_ctrls.subsample_tx,
                                           src,
                                           input_ptr->stride_y << pcs->tpl_ctrls.subsample_tx,
                                           predictor,
                                           16 << pcs->tpl_ctrls.subsample_tx);

                    EB_TRANS_COEFF_SHAPE pf_shape = pcs->tpl_ctrls.pf_shape;
                    svt_av1_wht_fwd_txfm(src_diff,
                                         16 << pcs->tpl_ctrls.subsample_tx,
                                         coeff,
                                         cost_tx_size,
                                         pf_shape,
                                         8,
                                         0);

                    intra_cost = svt_aom_satd(coeff, 256 >> pcs->tpl_ctrls.subsample_tx)
                        << pcs->tpl_ctrls.subsample_tx;
                }
                // printf("svt_aom_open_loop_intra_search_mb aom_satd mbxy %d %d, mode=%d, satd=%d, dst[0~4]=0x%d,%d,%d,%d\n", cu_origin_x, cu_origin_y, ois_intra_mode, intra_cost, predictor[0], predictor[1], predictor[2], predictor[3]);
                if (intra_cost < best_intra_cost) {
                    best_intra_cost = intra_cost;
                    best_mode       = ois_intra_mode;
                }
            }
            // store intra_cost to pcs
            ois_mb_results_ptr->intra_mode = best_mode;
            ois_mb_results_ptr->intra_cost = best_intra_cost;
            //if(pcs->picture_number == 16 && cu_origin_x <= 15 && cu_origin_y == 0)
            //    printf("svt_aom_open_loop_intra_search_mb cost0 poc%d b64_index=%d, mb_origin_xy=%d %d, best_mode=%d, best_intra_cost=%d, offset=%d, src[0~3]= %d %d %d %d\n", pcs->picture_number, b64_index, cu_origin_x, cu_origin_y, best_mode, best_intra_cost, (cu_origin_y >> 4) * mb_stride + (cu_origin_x >> 4), src[0], src[1], src[2], src[3]);
        }
        pa_blk_index++;
    }
    return return_error;
}
// clang-format on
