]> git.ipfire.org Git - thirdparty/glibc.git/blobdiff - sysdeps/ieee754/ldbl-96/k_sinl.c
Update copyright dates with scripts/update-copyrights
[thirdparty/glibc.git] / sysdeps / ieee754 / ldbl-96 / k_sinl.c
index b7b5ae359db721cfeba2b10b81a8f648dd8c0841..06a8f7c14af76494c82613fd15029bd9628e9706 100644 (file)
@@ -1,7 +1,6 @@
 /* Quad-precision floating point sine on <-pi/4,pi/4>.
-   Copyright (C) 1999-2015 Free Software Foundation, Inc.
+   Copyright (C) 1999-2024 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
-   Based on quad-precision sine by Jakub Jelinek <jj@ultra.linux.cz>
 
    The GNU C Library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
@@ -15,7 +14,7 @@
 
    You should have received a copy of the GNU Lesser General Public
    License along with the GNU C Library; if not, see
-   <http://www.gnu.org/licenses/>.  */
+   <https://www.gnu.org/licenses/>.  */
 
 /* The polynomials have not been optimized for extended-precision and
    may contain more terms than needed.  */
@@ -23,6 +22,7 @@
 #include <float.h>
 #include <math.h>
 #include <math_private.h>
+#include <math-underflow.h>
 
 /* The polynomials have not been optimized for extended-precision and
    may contain more terms than needed.  */
@@ -96,11 +96,7 @@ __kernel_sinl(long double x, long double y, int iy)
         polynomial of degree 17.  */
       if (absx < 0x1p-33L)
        {
-         if (fabsl (x) < LDBL_MIN)
-           {
-             long double force_underflow = x * x;
-             math_force_eval (force_underflow);
-           }
+         math_check_force_underflow (x);
          if (!((int)x)) return x;      /* generate inexact */
        }
       z = x * x;