mirror of
https://sourceware.org/git/glibc.git
synced 2024-11-30 08:40:07 +00:00
ec0ce0d3be
Similar to various other bugs in this area, some asin implementations do not raise the underflow exception for subnormal arguments, when the result is tiny and inexact. This patch forces the exception in a similar way to previous fixes. Tested for x86_64, x86, powerpc and mips64. [BZ #16351] * sysdeps/i386/fpu/e_asin.S (dbl_min): New object. (MO): New macro. (__ieee754_asin): Force underflow exception for results with small absolute value. * sysdeps/i386/fpu/e_asinf.S (flt_min): New object. (MO): New macro. (__ieee754_asinf): Force underflow exception for results with small absolute value. * sysdeps/ieee754/dbl-64/e_asin.c: Include <float.h> and <math.h>. (__ieee754_asin): Force underflow exception for results with small absolute value. * sysdeps/ieee754/flt-32/e_asinf.c: Include <float.h>. (__ieee754_asinf): Force underflow exception for results with small absolute value. * sysdeps/ieee754/ldbl-128/e_asinl.c: Include <float.h>. (__ieee754_asinl): Force underflow exception for results with small absolute value. * sysdeps/ieee754/ldbl-128ibm/e_asinl.c: Include <float.h>. (__ieee754_asinl): Force underflow exception for results with small absolute value. * sysdeps/ieee754/ldbl-96/e_asinl.c: Include <float.h>. (__ieee754_asinl): Force underflow exception for results with small absolute value. * sysdeps/x86_64/fpu/multiarch/e_asin.c [HAVE_FMA4_SUPPORT]: Include <math.h>. * math/auto-libm-test-in: Do not mark underflow exceptions as possibly missing for bug 16351. * math/auto-libm-test-out: Regenerated.
56 lines
975 B
ArmAsm
56 lines
975 B
ArmAsm
/*
|
|
* Written by J.T. Conklin <jtc@netbsd.org>.
|
|
* Public domain.
|
|
*/
|
|
|
|
#include <machine/asm.h>
|
|
|
|
RCSID("$NetBSD: e_asin.S,v 1.4 1995/05/08 23:45:40 jtc Exp $")
|
|
|
|
.section .rodata.cst8,"aM",@progbits,8
|
|
|
|
.p2align 3
|
|
.type dbl_min,@object
|
|
dbl_min: .byte 0, 0, 0, 0, 0, 0, 0x10, 0
|
|
ASM_SIZE_DIRECTIVE(dbl_min)
|
|
|
|
#ifdef PIC
|
|
# define MO(op) op##@GOTOFF(%ecx)
|
|
#else
|
|
# define MO(op) op
|
|
#endif
|
|
|
|
.text
|
|
|
|
/* asin = atan (x / sqrt((1-x) (1+x))) */
|
|
ENTRY(__ieee754_asin)
|
|
#ifdef PIC
|
|
LOAD_PIC_REG (cx)
|
|
#endif
|
|
fldl 4(%esp) /* x */
|
|
fld %st
|
|
fld1 /* 1 : x : x */
|
|
fsubp /* 1 - x : x */
|
|
fld1 /* 1 : 1 - x : x */
|
|
fadd %st(2) /* 1 + x : 1 - x : x */
|
|
fmulp /* 1 - x^2 */
|
|
fsqrt /* sqrt (1 - x^2) */
|
|
fpatan
|
|
fldl MO(dbl_min)
|
|
fld %st(1)
|
|
fabs
|
|
fucompp
|
|
fnstsw
|
|
sahf
|
|
jnc 1f
|
|
subl $8, %esp
|
|
cfi_adjust_cfa_offset (8)
|
|
fld %st(0)
|
|
fmul %st(0)
|
|
fstpl (%esp)
|
|
addl $8, %esp
|
|
cfi_adjust_cfa_offset (-8)
|
|
1: ret
|
|
END (__ieee754_asin)
|
|
strong_alias (__ieee754_asin, __asin_finite)
|