54 lines · cpp
1//===-- asin_fuzz.cpp -----------------------------------------------------===//2//3// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.4// See https://llvm.org/LICENSE.txt for license information.5// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception6//7//===----------------------------------------------------------------------===//8///9/// Fuzzing test for llvm-libc asin implementation.10///11//===----------------------------------------------------------------------===//12 13#include "src/math/asin.h"14#include "utils/MPFRWrapper/mpfr_inc.h"15#include <cstdint>16#include <cstring>17#include <iostream>18#include <math.h>19 20extern "C" int LLVMFuzzerTestOneInput(const uint8_t *data, size_t size) {21 mpfr_t input;22 mpfr_init2(input, 53);23 for (size_t i = 0; i < size / sizeof(double); ++i) {24 double x;25 std::memcpy(&x, data, sizeof(double));26 data += sizeof(double);27 28 // remove NaN and inf and values outside accepted range29 if (isnan(x) || isinf(x) || x > 1 || x < -1)30 continue;31 32 // signed zeros already tested in unit tests33 if (signbit(x) && x == 0.0)34 continue;35 36 mpfr_set_d(input, x, MPFR_RNDN);37 int output = mpfr_asin(input, input, MPFR_RNDN);38 mpfr_subnormalize(input, output, MPFR_RNDN);39 double to_compare = mpfr_get_d(input, MPFR_RNDN);40 41 double result = LIBC_NAMESPACE::asin(x);42 43 if (result != to_compare) {44 std::cout << std::hexfloat << "Failing input: " << x << std::endl;45 std::cout << std::hexfloat << "Failing output: " << result << std::endl;46 std::cout << std::hexfloat << "Expected: " << to_compare << std::endl;47 __builtin_trap();48 }49 }50 51 mpfr_clear(input);52 return 0;53}54