[thirdparty/glibc.git] / sysdeps / ieee754 / ldbl-128ibm / s_getpayloadl.c

/* Get NaN payload.  ldbl-128ibm version.
   Copyright (C) 2016-2019 Free Software Foundation, Inc.
   This file is part of the GNU C Library.

   The GNU C Library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.

   The GNU C Library 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
   Lesser General Public License for more details.

   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/>.  */

#include <fix-int-fp-convert-zero.h>
#include <math.h>
#include <math_private.h>
#include <stdint.h>

long double
__getpayloadl (const long double *x)
{
  double xhi = ldbl_high (*x);
  uint64_t ix;
  EXTRACT_WORDS64 (ix, xhi);
  ix &= 0x7ffffffffffffULL;
  if (FIX_INT_FP_CONVERT_ZERO && ix == 0)
    return 0.0L;
  return (long double) ix;
}
weak_alias (__getpayloadl, getpayloadl)
Commit	Line	Data
f8e8b8ed	1	/* Get NaN payload. ldbl-128ibm version.
04277e02	2	Copyright (C) 2016-2019 Free Software Foundation, Inc.
f8e8b8ed JM	3	This file is part of the GNU C Library.
	4
	5	The GNU C Library is free software; you can redistribute it and/or
	6	modify it under the terms of the GNU Lesser General Public
	7	License as published by the Free Software Foundation; either
	8	version 2.1 of the License, or (at your option) any later version.
	9
	10	The GNU C Library is distributed in the hope that it will be useful,
	11	but WITHOUT ANY WARRANTY; without even the implied warranty of
	12	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
	13	Lesser General Public License for more details.
	14
	15	You should have received a copy of the GNU Lesser General Public
	16	License along with the GNU C Library; if not, see
	17	<http://www.gnu.org/licenses/>. */
	18
	19	#include <fix-int-fp-convert-zero.h>
	20	#include <math.h>
	21	#include <math_private.h>
	22	#include <stdint.h>
	23
	24	long double
01f28812	25	__getpayloadl (const long double *x)
f8e8b8ed JM	26	{
	27	double xhi = ldbl_high (*x);
	28	uint64_t ix;
	29	EXTRACT_WORDS64 (ix, xhi);
	30	ix &= 0x7ffffffffffffULL;
	31	if (FIX_INT_FP_CONVERT_ZERO && ix == 0)
	32	return 0.0L;
	33	return (long double) ix;
	34	}
01f28812	35	weak_alias (__getpayloadl, getpayloadl)