1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
|
//===-- sincos_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 cos implementation.
///
//===----------------------------------------------------------------------===//
#include "src/math/sincos.h"
#include "utils/MPFRWrapper/mpfr_inc.h"
#include <math.h>
extern "C" int LLVMFuzzerTestOneInput(double x) {
// remove NaN and inf as preconditions
if (isnan(x) || isinf(x))
return 0;
// signed zeros already tested in unit tests
if (signbit(x) && x == 0.0)
return 0;
mpfr_t input;
mpfr_t sin_x;
mpfr_t cos_x;
mpfr_init2(input, 53);
mpfr_init2(sin_x, 53);
mpfr_init2(cos_x, 53);
mpfr_set_d(input, x, MPFR_RNDN);
int output = mpfr_sin_cos(sin_x, cos_x, input, MPFR_RNDN);
mpfr_subnormalize(sin_x, output, MPFR_RNDN);
mpfr_subnormalize(cos_x, output, MPFR_RNDN);
double to_compare_sin = mpfr_get_d(sin_x, MPFR_RNDN);
double to_compare_cos = mpfr_get_d(cos_x, MPFR_RNDN);
double sin_res, cos_res;
LIBC_NAMESPACE::sincos(x, &sin_res, &cos_res);
if (sin_res != to_compare_sin || cos_res != to_compare_cos)
__builtin_trap();
mpfr_clear(input);
mpfr_clear(sin_x);
mpfr_clear(cos_x);
return 0;
}
|