//==---------------- fma.cpp - DPC++ ESIMD on-device test -------------==//
//
// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
// See https://llvm.org/LICENSE.txt for license information.
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
//
//===----------------------------------------------------------------------===//
// REQUIRES-INTEL-DRIVER: lin: 29138, win: 101.5499
// RUN: %{build} -fsycl-device-code-split=per_kernel -o %t.out
// RUN: %{run} %t.out

#include "esimd_test_utils.hpp"

#include <sycl/ext/intel/experimental/esimd/math.hpp>

using namespace sycl;

template <typename T> bool test(queue q) {
  std::cout << "Running case: T=" << esimd_test::type_name<T>() << std::endl;
  constexpr unsigned Size = 16;
  constexpr unsigned VL = 16;

  T *A = new T[Size];
  T *B = new T[Size];
  T *C = new T[Size];

  for (unsigned i = 0; i < Size; ++i) {
    A[i] = B[i] = i;
    C[i] = 0.0f;
  }

  try {
    buffer<T, 1> bufa(A, range<1>(Size));
    buffer<T, 1> bufb(B, range<1>(Size));
    buffer<T, 1> bufc(C, range<1>(Size));

    auto e = q.submit([&](handler &cgh) {
      auto PA = bufa.template get_access<access::mode::read>(cgh);
      auto PB = bufb.template get_access<access::mode::read>(cgh);
      auto PC = bufc.template get_access<access::mode::write>(cgh);
      cgh.single_task([=]() SYCL_ESIMD_KERNEL {
        using namespace sycl::ext::intel::esimd;
        unsigned int offset = 0;
        simd<T, VL> va;
        va.copy_from(PA, offset);
        simd<T, VL> vb;
        vb.copy_from(PB, offset);
        simd<T, VL> vc = 5.0f;
        simd<T, VL> vd = ext::intel::experimental::esimd::fma(va, vb, vc);
        vd.copy_to(PC, offset);
      });
    });
    e.wait();
  } catch (sycl::exception const &e) {
    std::cout << "SYCL exception caught: " << e.what() << '\n';

    delete[] A;
    delete[] B;
    delete[] C;
    return 0;
  }

  int err_cnt = 0;

  for (unsigned i = 0; i < Size; ++i) {
    if ((A[i] * B[i]) + 5.0f != C[i]) {
      if (++err_cnt < 10) {
        std::cout << "failed at index " << i << ", " << C[i] << " != " << A[i]
                  << " + " << B[i] << "\n";
      }
    }
  }
  if (err_cnt > 0) {
    std::cout << "  pass rate: "
              << ((float)(Size - err_cnt) / (float)Size) * 100.0f << "% ("
              << (Size - err_cnt) << "/" << Size << ")\n";
  }

  delete[] A;
  delete[] B;
  delete[] C;

  std::cout << (err_cnt > 0 ? "FAILED\n" : "Passed\n");
  return err_cnt == 0;
}

int main() {
  auto q = queue{gpu_selector_v};
  esimd_test::printTestLabel(q);
  bool passed = true;

  passed &= test<float>(q);
#ifdef TEST_PVC
  passed &= test<sycl::ext::oneapi::bfloat16>(q);
  passed &= test<ext::intel::experimental::esimd::tfloat32>(q);
#endif
  if (q.get_device().has(sycl::aspect::fp16))
    passed &= test<sycl::half>(q);
  if (q.get_device().has(sycl::aspect::fp64))
    passed &= test<double>(q);

  return passed ? 0 : 1;
}