diff --git a/libc/fuzzing/math/CMakeLists.txt b/libc/fuzzing/math/CMakeLists.txt index 31336fa3e49f3..1ac72c1007ef2 100644 --- a/libc/fuzzing/math/CMakeLists.txt +++ b/libc/fuzzing/math/CMakeLists.txt @@ -62,6 +62,15 @@ add_libc_fuzzer( libc.src.math.nextafterl ) +add_libc_fuzzer( + asin_fuzz + NEED_MPFR + SRCS + asin_fuzz.cpp + DEPENDS + libc.src.math.asin +) + add_libc_fuzzer( sin_fuzz NEED_MPFR diff --git a/libc/fuzzing/math/asin_fuzz.cpp b/libc/fuzzing/math/asin_fuzz.cpp new file mode 100644 index 0000000000000..94ae5c7bfdeee --- /dev/null +++ b/libc/fuzzing/math/asin_fuzz.cpp @@ -0,0 +1,38 @@ +//===-- asin_fuzz.cpp -----------------------------------------------------===// +// +// 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 +// +//===----------------------------------------------------------------------===// +/// +/// Fuzzing test for llvm-libc asin implementation. +/// +//===----------------------------------------------------------------------===// + +#include "src/math/asin.h" +#include "utils/MPFRWrapper/mpfr_inc.h" +#include + +extern "C" int LLVMFuzzerTestOneInput(double x) { + // remove NaN and inf and values outside accepted range + if (isnan(x) || isinf(x) || x > 1 || x < -1) + return 0; + // signed zeros already tested in unit tests + if (signbit(x) && x == 0.0) + return 0; + mpfr_t input; + mpfr_init2(input, 53); + mpfr_set_d(input, x, MPFR_RNDN); + int output = mpfr_asin(input, input, MPFR_RNDN); + mpfr_subnormalize(input, output, MPFR_RNDN); + double to_compare = mpfr_get_d(input, MPFR_RNDN); + + double result = LIBC_NAMESPACE::asin(x); + + if (result != to_compare) + __builtin_trap(); + + mpfr_clear(input); + return 0; +}