[libc-commits] [libc] [libc] asin fuzz testing (PR #147786)

[via libc-commits]

https://github.com/sribee8 created https://github.com/llvm/llvm-project/pull/147786

Added fuzz test for asin

>From fb5c77064bf9905e093083acc9092a1e920b55cf Mon Sep 17 00:00:00 2001
From: Sriya Pratipati <sriyap at google.com>
Date: Wed, 9 Jul 2025 17:18:12 +0000
Subject: [PATCH] [libc] asin fuzz testing

Added fuzz test for asin
---
 libc/fuzzing/math/CMakeLists.txt |  9 ++++++++
 libc/fuzzing/math/asin_fuzz.cpp  | 38 ++++++++++++++++++++++++++++++++
 2 files changed, 47 insertions(+)
 create mode 100644 libc/fuzzing/math/asin_fuzz.cpp

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..20e7a58d4816f
--- /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 <math.h>
+
+extern "C" int LLVMFuzzerTestOneInput(const 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;
+}