/* Test of sqrtl() function. Copyright (C) 2010-2023 Free Software Foundation, Inc. This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program. If not, see . */ /* Written by Bruno Haible , 2010. */ #include #include #include "signature.h" SIGNATURE_CHECK (sqrtl, long double, (long double)); #include #include "fpucw.h" #include "macros.h" #define DOUBLE long double #define L_(literal) literal##L #define MANT_DIG DBL_MANT_DIG #define SQRT sqrtl #define RANDOM randoml #include "test-sqrt.h" static long double my_ldexpl (long double x, int d) { for (; d > 0; d--) x *= 2.0L; for (; d < 0; d++) x *= 0.5L; return x; } int main () { DECL_LONG_DOUBLE_ROUNDING BEGIN_LONG_DOUBLE_ROUNDING (); /* A particular value. */ x = 0.6L; y = sqrtl (x); ASSERT (y >= 0.7745966692L && y <= 0.7745966693L); /* Another particular value. */ { long double z; long double err; x = 8.1974099812331540680810141969554806865L; y = sqrtl (x); z = y * y - x; err = my_ldexpl (z, LDBL_MANT_DIG); if (err < 0) err = - err; ASSERT (err <= 100.0L); } test_function (); return 0; }