libc/test/src/math/exhaustive/fmabf16_test.cpp

//===-- Exhaustive tests for fmabf16 --------------------------------------===//
//
// 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
//
//===----------------------------------------------------------------------===//

#include "exhaustive_test.h"
#include "src/__support/FPUtil/FPBits.h"
#include "src/__support/FPUtil/bfloat16.h"
#include "src/math/fmabf16.h"
#include "test/UnitTest/FPMatcher.h"
#include "utils/MPFRWrapper/MPCommon.h"
#include "utils/MPFRWrapper/MPFRUtils.h"

namespace mpfr = LIBC_NAMESPACE::testing::mpfr;
using LIBC_NAMESPACE::fputil::BFloat16;

struct FmaBf16Checker : public virtual LIBC_NAMESPACE::testing::Test {

  using FloatType = BFloat16;
  using FPBits = LIBC_NAMESPACE::fputil::FPBits<bfloat16>;
  using StorageType = typename FPBits::StorageType;

  uint64_t check(uint16_t x_start, uint16_t x_stop, uint16_t y_start,
                 uint16_t y_stop, mpfr::RoundingMode rounding) {

    mpfr::ForceRoundingMode r(rounding);
    if (!r.success)
      return true;
    uint16_t xbits = x_start;
    uint64_t failed = 0;
    do {
      BFloat16 x = FPBits(xbits).get_val();
      uint16_t ybits = y_start;
      do {
        BFloat16 y = FPBits(ybits).get_val();
        BFloat16 z = FPBits(uint16_t(0x03E1)).get_val();
        mpfr::TernaryInput<BFloat16> input{x, y, z};
        bool correct = TEST_MPFR_MATCH_ROUNDING_SILENTLY(
            mpfr::Operation::Fma, input, LIBC_NAMESPACE::fmabf16(x, y, z), 0.5,
            rounding);
        failed += (!correct);

      } while (ybits++ < y_stop);
    } while (xbits++ < x_stop);
    return failed;
  }
};

using LlvmLibcBfloat16ExhaustiveFmaTest =
    LlvmLibcExhaustiveMathTest<FmaBf16Checker, 1 << 2>;

// range: [0, inf]
static constexpr uint16_t POS_START = 0x0000U;
static constexpr uint16_t POS_STOP = 0x7f80U;

// range: [-0, -inf]
static constexpr uint16_t NEG_START = 0x8000U;
static constexpr uint16_t NEG_STOP = 0xff80U;

TEST_F(LlvmLibcBfloat16ExhaustiveFmaTest, PositiveRange) {
  test_full_range_all_roundings(POS_START, POS_STOP, POS_START, POS_STOP);
}

TEST_F(LlvmLibcBfloat16ExhaustiveFmaTest, PositiveNegative) {
  test_full_range_all_roundings(POS_START, POS_STOP, NEG_START, NEG_STOP);
}

TEST_F(LlvmLibcBfloat16ExhaustiveFmaTest, NegativePositive) {
  test_full_range_all_roundings(NEG_START, NEG_STOP, POS_START, POS_STOP);
}

TEST_F(LlvmLibcBfloat16ExhaustiveFmaTest, NegativeRange) {
  test_full_range_all_roundings(NEG_START, NEG_STOP, NEG_START, NEG_STOP);
}