mirror of
https://git.savannah.gnu.org/git/gnulib.git
synced 2025-08-20 11:01:20 +03:00
7fd1d13152
* lib/math.in.h (sqrtl): Replace it if REPLACE_SQRTL is 1. * m4/sqrtl.m4 (gl_FUNC_SQRTL_WORKS): New macro. (gl_FUNC_SQRTL): Invoke it. Set REPLACE_SQRTL to 1 if sqrtl() produces too big rounding errors. * m4/math_h.m4 (gl_MATH_H_DEFAULTS): Initialize REPLACE_SQRTL. * modules/math (Makefile.am): Substitute REPLACE_SQRTL. * modules/sqrtl (configure.ac): Consider REPLACE_SQRTL. (Depends-on): Update conditions. * tests/test-sqrtl.c (my_ldexpl): New function. (main): Add test of a particular value. * doc/posix-functions/sqrtl.texi: Mention the OpenBSD 5.1/SPARC bug.
78 lines
1.7 KiB
C
78 lines
1.7 KiB
C
/* Test of sqrtl() function.
|
|
Copyright (C) 2010-2012 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 <http://www.gnu.org/licenses/>. */
|
|
|
|
/* Written by Bruno Haible <bruno@clisp.org>, 2010. */
|
|
|
|
#include <config.h>
|
|
|
|
#include <math.h>
|
|
|
|
#include "signature.h"
|
|
SIGNATURE_CHECK (sqrtl, long double, (long double));
|
|
|
|
#include <float.h>
|
|
|
|
#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;
|
|
}
|