mirrors/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2024-11-23 17:53:37 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

LoongArch: Add support for scalb[f]

Browse Source 
Add inline assembler for the scalb functions. Passes GLIBC regression.
This commit is contained in:
Xiaolin Tang 2022-11-23 11:49:25 +08:00 committed by caiyinyu
parent e3d9ca4e94
commit d3366b6b87
2 changed files with 120 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

60
sysdeps/loongarch/fpu/e_scalb.c Normal file
View File

@ -0,0 +1,60 @@
/* scalb(). LoongArch version.
Copyright (C) 2022 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/>. */
#define NO_MATH_REDIRECT
#include <math.h>
#include <libm-alias-finite.h>
#include <fpu_control.h>
#include <float.h>
double
__ieee754_scalb (double x, double fn)
{
int x_cond;
int fn_cond;
asm volatile ("fclass.d \t%0, %1" : "=f" (x_cond) : "f" (x));
asm volatile ("fclass.d \t%0, %1" : "=f" (fn_cond) : "f" (fn));
if (__glibc_unlikely(( x_cond | fn_cond) & _FCLASS_NAN))
return x * fn;
else if (__glibc_unlikely(fn_cond & _FCLASS_INF))
{
if (!(fn_cond & _FCLASS_MINF))
return x * fn;
else
return x / -fn;
}
else if (__glibc_likely(-DBL_MAX < fn && fn < DBL_MAX))
{
double rint_fn, tmp;
/* rint_fn = rint(fn) */
asm volatile ("frint.d \t%0, %1" : "=f" (rint_fn) : "f" (fn));
if (rint_fn != fn )
return (x - x) / (x - x);
asm volatile ("ftintrz.l.d \t%0, %1" : "=f" (tmp) : "f" (rint_fn));
asm volatile ("fscaleb.d \t%0, %1, %2" : "=f" (x) : "f" (x), "f" (tmp));
}
else
asm volatile ("fscaleb.d \t%0, %1, %2" : "=f" (x) : "f" (x), "f" (fn));
return x;
}
libm_alias_finite (__ieee754_scalb, __scalb)

60
sysdeps/loongarch/fpu/e_scalbf.c Normal file
View File

@ -0,0 +1,60 @@
/* scalbf(). LoongArch version.
Copyright (C) 2022 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/>. */
#define NO_MATH_REDIRECT
#include <math.h>
#include <libm-alias-finite.h>
#include <fpu_control.h>
#include <float.h>
float
__ieee754_scalbf (float x, float fn)
{
int x_cond;
int fn_cond;
asm volatile ("fclass.s \t%0, %1" : "=f" (x_cond) : "f" (x));
asm volatile ("fclass.s \t%0, %1" : "=f" (fn_cond) : "f" (fn));
if (__glibc_unlikely(( x_cond | fn_cond) & _FCLASS_NAN))
return x * fn;
else if (__glibc_unlikely(fn_cond & _FCLASS_INF))
{
if (!(fn_cond & _FCLASS_MINF))
return x * fn;
else
return x / -fn;
}
else if (__glibc_likely(-FLT_MAX < fn && fn < FLT_MAX))
{
float rintf_fn, tmp;
/* rintf_fn = rintf(fn) */
asm volatile ("frint.s \t%0, %1" : "=f" (rintf_fn) : "f" (fn));
if (rintf_fn != fn )
return (x - x) / (x - x);
asm volatile ("ftintrz.w.s \t%0, %1" : "=f" (tmp) : "f" (rintf_fn));
asm volatile ("fscaleb.s \t%0, %1, %2" : "=f" (x) : "f" (x), "f" (tmp));
}
else
asm volatile ("fscaleb.s \t%0, %1, %2" : "=f" (x) : "f" (x), "f" (fn));
return x;
}
libm_alias_finite (__ieee754_scalb, __scalb)