/**
 * Copyright (c) 2025 Huawei Technologies Co., Ltd.
 * This program is free software, you can redistribute it and/or modify it under the terms and conditions of
 * CANN Open Software License Agreement Version 2.0 (the "License").
 * Please refer to the License for details. You may not use this file except in compliance with the License.
 * THIS SOFTWARE IS PROVIDED ON AN "AS IS" BASIS, WITHOUT WARRANTIES OF ANY KIND, EITHER EXPRESS OR IMPLIED,
 * INCLUDING BUT NOT LIMITED TO NON-INFRINGEMENT, MERCHANTABILITY, OR FITNESS FOR A PARTICULAR PURPOSE.
 * See LICENSE in the root of the software repository for the full text of the License.
 */

#include "aicpu_test_utils.h"

using namespace std;

uint64_t CalTotalElements(std::vector<std::vector<int64_t>>& shapes, uint32_t index)
{
    if (index < 0U) {
        return 0;
    }
    uint64_t nums = 1;
    for (auto shape : shapes[index]) {
        nums = nums * shape;
    }
    return nums;
}

template <>
bool CompareResult(float output[], float expect_output[], uint64_t num)
{
    bool result = true;
    for (uint64_t i = 0; i < num; ++i) {
        double absolute_error = std::fabs(output[i] - expect_output[i]);
        double relative_error = 0;
        if (expect_output[i] == 0) {
            relative_error = 2e-6;
        } else {
            relative_error = absolute_error / std::fabs(expect_output[i]);
        }
        if ((absolute_error > 1e-6) && (relative_error > 1e-6)) {
            std::cout << "output[" << i << "] = ";
            std::cout << output[i];
            std::cout << ", expect_output[" << i << "] = ";
            std::cout << expect_output[i] << std::endl;
            result = false;
        }
    }
    return result;
}

bool CompareResultAllClose(float output[], float expect_output[], uint64_t num)
{
    bool result = true;
    for (uint64_t i = 0; i < num; ++i) {
        double absolute_error = std::fabs(output[i] - expect_output[i]);
        float atol = 1e-4;
        float rtol = 1e-4;
        double relative_error = 0;
        if (!(absolute_error <= (atol + rtol * std::fabs(expect_output[i])))) {
            std::cout << "output[" << i << "] = ";
            std::cout << output[i];
            std::cout << ", expect_output[" << i << "] = ";
            std::cout << expect_output[i] << std::endl;
            result = false;
        }
    }

    return result;
}

template <>
bool CompareResult(double output[], double expect_output[], uint64_t num)
{
    bool result_code = true;
    for (uint64_t i = 0; i < num; ++i) {
        double absolute_error = std::fabs(output[i] - expect_output[i]);
        double relative_error = 0;
        if (expect_output[i] == 0) {
            relative_error = 2e-10;
        } else {
            relative_error = absolute_error / std::fabs(expect_output[i]);
        }
        if ((absolute_error > 1e-12) && (relative_error > 1e-10)) {
            std::cout << "output[" << i << "] = ";
            std::cout << output[i];
            std::cout << ", expect_output[" << i << "] = ";
            std::cout << expect_output[i] << std::endl;
            result_code = false;
        }
    }
    return result_code;
}

template <>
bool CompareResult(complex<float> output[], complex<float> expect_output[], uint64_t num)
{
    bool result = true;
    for (uint64_t i = 0; i < num; ++i) {
        double absolute_real_error = std::fabs(output[i].real() - expect_output[i].real());
        double relative_real_error = 0;
        if (expect_output[i].real() == 0) {
            relative_real_error = 2e-6;
        } else {
            relative_real_error = absolute_real_error / std::fabs(expect_output[i].real());
        }

        double absolute_imag_error = std::fabs(output[i].imag() - expect_output[i].imag());
        double relative_imag_error = 0;
        if (expect_output[i].imag() == 0) {
            relative_imag_error = 2e-6;
        } else {
            relative_imag_error = absolute_imag_error / std::fabs(expect_output[i].imag());
        }

        if ((absolute_real_error > 1e-6) && (relative_real_error > 1e-6) ||
            (absolute_imag_error > 1e-6) && (relative_imag_error > 1e-6)) {
            result = false;
            std::cout << "output[" << i << "] = ";
            std::cout << output[i];
            std::cout << ", expect_output[" << i << "] = ";
            std::cout << expect_output[i] << std::endl;
        }
    }

    return result;
}

template <>
bool CompareResult(complex<double> output[], complex<double> expect_output[], uint64_t num)
{
    bool result = true;
    for (uint64_t i = 0; i < num; ++i) {
        double absolute_real_error = std::fabs(output[i].real() - expect_output[i].real());
        double relative_real_error = 0;
        if (expect_output[i].real() == 0) {
            relative_real_error = 2e-10;
        } else {
            relative_real_error = absolute_real_error / std::fabs(expect_output[i].real());
        }

        double absolute_imag_error = std::fabs(output[i].imag() - expect_output[i].imag());
        double relative_imag_error = 0;
        if (expect_output[i].imag() == 0) {
            relative_imag_error = 2e-10;
        } else {
            relative_imag_error = absolute_imag_error / std::fabs(expect_output[i].imag());
        }

        if ((absolute_real_error > 10e-8) && (relative_real_error > 10e-8) ||
            (absolute_imag_error > 10e-8) && (relative_imag_error > 10e-8)) {
            std::cout << "output[" << i << "] = ";
            std::cout << output[i];
            std::cout << ", expect_output[" << i << "] = ";
            std::cout << expect_output[i] << std::endl;
            result = false;
        }
    }
    return result;
}

bool CompareResult(Eigen::half output[], Eigen::half expect_output[], uint64_t num)
{
    bool result = true;
    for (uint64_t i = 0; i < num; ++i) {
        Eigen::half absolute_error = (output[i] - expect_output[i]);
        absolute_error = absolute_error >= Eigen::half(0) ? absolute_error : -absolute_error;
        Eigen::half relative_error(0);
        if (expect_output[i] == Eigen::half(0)) {
            relative_error = Eigen::half(2e-3);
        } else {
            relative_error = absolute_error / expect_output[i];
            relative_error = relative_error >= Eigen::half(0) ? relative_error : -relative_error;
        }
        if ((absolute_error > Eigen::half(1e-3)) && (relative_error > Eigen::half(1e-3))) {
            std::cout << "output[" << i << "] = ";
            std::cout << output[i];
            std::cout << ", expect_output[" << i << "] = ";
            std::cout << expect_output[i] << std::endl;
            result = false;
        }
    }
    return result;
}