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

/* Return classification value corresponding to argument.
   Copyright (C) 1997-2025 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
   <https://www.gnu.org/licenses/>.  */

#include <math.h>

#include <math_private.h>


int
__fpclassifyl (_Float128 x)
{
  uint64_t hx, lx;
  int retval = FP_NORMAL;

  GET_LDOUBLE_WORDS64 (hx, lx, x);
  lx |= (hx & 0x0000ffffffffffffLL);
  hx &= 0x7fff000000000000LL;
  if ((hx | lx) == 0)
    retval = FP_ZERO;
  else if (hx == 0)
    retval = FP_SUBNORMAL;
  else if (hx == 0x7fff000000000000LL)
    retval = lx != 0 ? FP_NAN : FP_INFINITE;

  return retval;
}
libm_hidden_def (__fpclassifyl)
Commit	Line	Data
	1	/* Return classification value corresponding to argument.
	2	Copyright (C) 1997-2025 Free Software Foundation, Inc.
	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	<https://www.gnu.org/licenses/>. */
	18
	19	#include <math.h>
	20
	21	#include <math_private.h>
	22
	23
	24	int
	25	__fpclassifyl (_Float128 x)
	26	{
	27	uint64_t hx, lx;
	28	int retval = FP_NORMAL;
	29
	30	GET_LDOUBLE_WORDS64 (hx, lx, x);
	31	lx \|= (hx & 0x0000ffffffffffffLL);
	32	hx &= 0x7fff000000000000LL;
	33	if ((hx \| lx) == 0)
	34	retval = FP_ZERO;
	35	else if (hx == 0)
	36	retval = FP_SUBNORMAL;
	37	else if (hx == 0x7fff000000000000LL)
	38	retval = lx != 0 ? FP_NAN : FP_INFINITE;
	39
	40	return retval;
	41	}
	42	libm_hidden_def (__fpclassifyl)