/*
 * 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.
 */

/******************************************************************************
 * @file WedgeUtilTest.cc
 *
 * @brief Unit test for util functions in wedge prediction:
 * - svt_av1_wedge_sign_from_residuals_avx2
 * - svt_av1_wedge_compute_delta_squares_avx2
 * - svt_av1_wedge_sse_from_residuals_avx2
 * - svt_aom_sum_squares_i16_sse2
 *
 * @author Cidana-Wenyao
 *
 ******************************************************************************/
#include "gtest/gtest.h"
#include "definitions.h"
#include "utility.h"
#include "aom_dsp_rtcd.h"
#include "random.h"
#include "util.h"

using svt_av1_test_tool::SVTRandom;
namespace {
// test svt_av1_wedge_sign_from_residuals_avx2
// Choose the mask sign for a compound predictor.
class WedgeUtilTest : public ::testing::Test {
  public:
    void SetUp() override {
        r0 = (int16_t *)svt_aom_memalign(32, MAX_SB_SQUARE * sizeof(int16_t));
        r1 = (int16_t *)svt_aom_memalign(32, MAX_SB_SQUARE * sizeof(int16_t));
        m = (uint8_t *)svt_aom_memalign(32, MAX_SB_SQUARE);
    }

    void TearDown() override {
        svt_aom_free(r0);
        svt_aom_free(r1);
        svt_aom_free(m);
        aom_clear_system_state();
    }

    void wedge_sign_test(int N, int k) {
        DECLARE_ALIGNED(32, int16_t, ds[MAX_SB_SQUARE]);
        // pre-compute limit
        // MAX_MASK_VALUE/2 * (sum(r0**2) - sum(r1**2))
        int64_t limit;
        limit = (int64_t)svt_aom_sum_squares_i16_c(r0, N);
        limit -= (int64_t)svt_aom_sum_squares_i16_c(r1, N);
        limit *= (1 << WEDGE_WEIGHT_BITS) / 2;

        // calculate ds: r0**2 - r1**2
        for (int i = 0; i < N; ++i)
            ds[i] = clamp(r0[i] * r0[i] - r1[i] * r1[i], INT16_MIN, INT16_MAX);
        const int8_t ref_sign =
            svt_av1_wedge_sign_from_residuals_c(ds, m, N, limit);
        const int8_t sse2_tst_sign =
            svt_av1_wedge_sign_from_residuals_sse2(ds, m, N, limit);
        const int8_t avx2_tst_sign =
            svt_av1_wedge_sign_from_residuals_avx2(ds, m, N, limit);
        ASSERT_EQ(ref_sign, sse2_tst_sign)
            << "unit test for svt_av1_wedge_sign_from_residuals_sse2 fail at "
               "iteration "
            << k;
        ASSERT_EQ(ref_sign, avx2_tst_sign)
            << "unit test for svt_av1_wedge_sign_from_residuals_avx2 fail at "
               "iteration "
            << k;
    }

  protected:
    uint8_t *m;       /* mask */
    int16_t *r0, *r1; /* two predicted residual */
};

extern "C" uint64_t svt_aom_sum_squares_i16_c(const int16_t *src, uint32_t n);
#define MAX_MASK_VALUE (1 << WEDGE_WEIGHT_BITS)
TEST_F(WedgeUtilTest, MaskSignRandomTest) {
    const int iterations = 10000;
    SVTRandom rnd(13, true);             // max residual is 13-bit
    SVTRandom m_rnd(0, MAX_MASK_VALUE);  // [0, MAX_MASK_VALUE]
    SVTRandom n_rnd(1, 8191 / 64);       // required by assembly implementation

    for (int k = 0; k < iterations; ++k) {
        // populate the residual buffer randomly
        for (int i = 0; i < MAX_SB_SQUARE; ++i) {
            r0[i] = rnd.random();
            r1[i] = rnd.random();
            m[i] = m_rnd.random();
        }

        // N should be multiple of 64, required by
        // svt_av1_wedge_sign_from_residuals_avx2
        const int N = 64 * n_rnd.random();
        wedge_sign_test(N, k);
    }
}

TEST_F(WedgeUtilTest, MaskSignExtremeTest) {
    const int int_13bit_max = (1 << 12) - 1;
    SVTRandom rnd(13, true);             // max residual is 13-bit
    SVTRandom m_rnd(0, MAX_MASK_VALUE);  // [0, MAX_MASK_VALUE]
    SVTRandom n_rnd(1, 8191 / 64);       // required by assembly implementation

    for (int k = 0; k < 4; ++k) {
        switch (k) {
        case 0:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = 0;
                r1[i] = int_13bit_max;
            }
            break;
        case 1:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = int_13bit_max;
                r1[i] = 0;
            }
            break;
        case 2:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = -int_13bit_max;
                r1[i] = 0;
            }
            break;
        case 3:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = 0;
                r1[i] = -int_13bit_max;
            }
            break;
        default: assert(0); break;
        }

        for (int i = 0; i < MAX_SB_SQUARE; ++i)
            m[i] = MAX_MASK_VALUE;

        // N should be multiple of 64, required by
        // svt_av1_wedge_sign_from_residuals_avx2
        const int N = 64 * n_rnd.random();

        wedge_sign_test(N, k);
    }
}

// test svt_av1_wedge_compute_delta_squares_avx2
// element-by-element calculate the difference of square
TEST_F(WedgeUtilTest, ComputeDeltaSquareTest) {
    const int iterations = 10000;
    SVTRandom rnd(13, true);  // max residual is 13-bit
    SVTRandom n_rnd(1, MAX_SB_SQUARE / 64);
    DECLARE_ALIGNED(32, int16_t, ref_diff[MAX_SB_SQUARE]);
    DECLARE_ALIGNED(32, int16_t, tst_diff[MAX_SB_SQUARE]);

    for (int k = 0; k < iterations; ++k) {
        // populate the residual buffer randomly
        for (int i = 0; i < MAX_SB_SQUARE; ++i) {
            r0[i] = rnd.random();
            r1[i] = rnd.random();
        }
        memset(ref_diff, 0, sizeof(ref_diff));
        memset(tst_diff, 0, sizeof(ref_diff));

        // N should be multiple of 64, required by
        // svt_av1_wedge_compute_delta_squares_avx2
        const int N = 64 * n_rnd.random();

        svt_av1_wedge_compute_delta_squares_c(ref_diff, r0, r1, N);
        svt_av1_wedge_compute_delta_squares_avx2(tst_diff, r0, r1, N);

        // check the output
        for (int i = 0; i < N; ++i) {
            ASSERT_EQ(ref_diff[i], tst_diff[i])
                << "unit test for svt_av1_wedge_compute_delta_squares_avx2 "
                   "fail at "
                   "iteration "
                << k;
        }
    }
}

// test svt_av1_wedge_sse_from_residuals_avx2
// calculate the sse of two prediction combined with mask m
TEST_F(WedgeUtilTest, SseFromResidualRandomTest) {
    const int iterations = 10000;
    SVTRandom rnd(13, true);             // max residual is 13-bit
    SVTRandom m_rnd(0, MAX_MASK_VALUE);  // [0, MAX_MASK_VALUE]
    SVTRandom n_rnd(1, MAX_SB_SQUARE / 64);

    for (int k = 0; k < iterations; ++k) {
        // populate the residual buffer randomly
        for (int i = 0; i < MAX_SB_SQUARE; ++i) {
            r0[i] = rnd.random();
            r1[i] = rnd.random();
            m[i] = m_rnd.random();
        }

        // N should be multiple of 64, required by
        // svt_av1_wedge_sse_from_residuals_avx2
        const int N = 64 * n_rnd.random();

        uint64_t ref_sse = svt_av1_wedge_sse_from_residuals_c(r0, r1, m, N);
        uint64_t sse2_tst_sse =
            svt_av1_wedge_sse_from_residuals_sse2(r0, r1, m, N);
        uint64_t avx2_tst_sse =
            svt_av1_wedge_sse_from_residuals_avx2(r0, r1, m, N);

        // check output
        ASSERT_EQ(ref_sse, sse2_tst_sse)
            << "unit test for svt_av1_wedge_sse_from_residuals_sse2 fail at "
               "iteration "
            << k;
        ASSERT_EQ(ref_sse, avx2_tst_sse)
            << "unit test for svt_av1_wedge_sse_from_residuals_avx2 fail at "
               "iteration "
            << k;
    }
}

TEST_F(WedgeUtilTest, SseFromResidualExtremeTest) {
    const int int_13bit_max = (1 << 12) - 1;
    SVTRandom rnd(13, true);             // max residual is 13-bit
    SVTRandom m_rnd(0, MAX_MASK_VALUE);  // [0, MAX_MASK_VALUE]
    SVTRandom n_rnd(1, MAX_SB_SQUARE / 64);

    for (int k = 0; k < 4; ++k) {
        switch (k) {
        case 0:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = 0;
                r1[i] = int_13bit_max;
            }
            break;
        case 1:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = int_13bit_max;
                r1[i] = 0;
            }
            break;
        case 2:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = -int_13bit_max;
                r1[i] = 0;
            }
            break;
        case 3:
            for (int i = 0; i < MAX_SB_SQUARE; ++i) {
                r0[i] = 0;
                r1[i] = -int_13bit_max;
            }
            break;
        default: assert(0); break;
        }

        for (int i = 0; i < MAX_SB_SQUARE; ++i)
            m[i] = MAX_MASK_VALUE;

        // N should be multiple of 64, required by
        // svt_av1_wedge_sse_from_residuals_avx2
        const int N = 64 * n_rnd.random();

        uint64_t ref_sse = svt_av1_wedge_sse_from_residuals_c(r0, r1, m, N);
        uint64_t sse2_tst_sse =
            svt_av1_wedge_sse_from_residuals_sse2(r0, r1, m, N);
        uint64_t avx2_tst_sse =
            svt_av1_wedge_sse_from_residuals_avx2(r0, r1, m, N);

        // check output
        ASSERT_EQ(ref_sse, sse2_tst_sse)
            << "unit test for svt_av1_wedge_sse_from_residuals_sse2 fail at "
               "iteration "
            << k;
        ASSERT_EQ(ref_sse, avx2_tst_sse)
            << "unit test for svt_av1_wedge_sse_from_residuals_avx2 fail at "
               "iteration "
            << k;
    }
}

typedef uint64_t (*AomSumSquaresI16Func)(const int16_t *, uint32_t);
typedef ::testing::tuple<BlockSize, AomSumSquaresI16Func> AomHSumSquaresParam;

class AomSumSquaresTest : public ::testing::TestWithParam<AomHSumSquaresParam> {
  public:
    AomSumSquaresTest() : rnd_(0, 255){};
    virtual ~AomSumSquaresTest() {
    }

    void TearDown() override {
        aom_clear_system_state();
    }

    void run_test() {
        const int block_size = TEST_GET_PARAM(0);
        AomSumSquaresI16Func test_impl = TEST_GET_PARAM(1);
        const int width = block_size_wide[block_size];
        const int height = block_size_high[block_size];
        DECLARE_ALIGNED(16, uint16_t, src_[MAX_SB_SQUARE]);
        const int run_times = 100;
        for (int i = 0; i < run_times; ++i) {
            memset(src_, 0, sizeof(src_));
            for (int j = 0; j < width * height; j++) {
                src_[j] = rnd_.random();
            }

            uint64_t res_ref_ = svt_aom_sum_squares_i16_c((const int16_t *)src_,
                                                          width * height);

            uint64_t res_tst_ =
                test_impl((const int16_t *)src_, width * height);

            ASSERT_EQ(res_ref_, res_tst_);
        }
    }

  private:
    SVTRandom rnd_;
};

TEST_P(AomSumSquaresTest, MatchTest) {
    run_test();
}

INSTANTIATE_TEST_SUITE_P(
    SUM_SQUARES_TEST, AomSumSquaresTest,
    ::testing::Combine(::testing::Range(BLOCK_4X4, BlockSizeS_ALL),
                       ::testing::Values(svt_aom_sum_squares_i16_sse2)));

}  // namespace
