1/**
  2 * This file has no copyright assigned and is placed in the Public Domain.
  3 * This file is part of the mingw-w64 runtime package.
  4 * No warranty is given; refer to the file DISCLAIMER.PD within this package.
  5 */
  6#include <_mingw_mac.h>
  7
  8	.file	"log1pl.S"
  9	.text
 10	/* The fyl2xp1 can only be used for values in
 11	   -1 + sqrt(2) / 2 <= x <= 1 - sqrt(2) / 2
 12	   0.29 is a safe value.
 13	 */
 14
 15	/* Only gcc understands the .tfloat type
 16	   The series of .long below represents
 17	   limit:	.tfloat 0.29
 18	 */
 19	.align 16
 20limit:
 21	.long 2920577761
 22	.long 2491081031
 23	.long 16381
 24#ifdef __x86_64__
 25	.align 8
 26#else
 27	.align 4
 28#endif
 29	/* Please note:	 we use a double value here.  Since 1.0 has
 30	   an exact representation this does not effect the accuracy
 31	   but it helps to optimize the code.  */
 32one:	.double 1.0
 33
 34/*
 35 * Use the fyl2xp1 function when the argument is in the range -0.29 to 0.29,
 36 * otherwise fyl2x with the needed extra computation.
 37 */
 38.globl __MINGW_USYMBOL(log1pl)
 39	.def	__MINGW_USYMBOL(log1pl);	.scl	2;	.type	32;	.endef
 40__MINGW_USYMBOL(log1pl):
 41#ifdef __x86_64__
 42	fldln2
 43	fldt	(%rdx)
 44	fxam
 45	fnstsw
 46	fld	%st
 47	sahf
 48	jc	3f		// in case x is NaN or �Inf
 494:
 50	fabs
 51	fldt	limit(%rip)
 52	fcompp
 53	fnstsw
 54	sahf
 55	jnc	2f
 56	faddl	one(%rip)
 57	fyl2x
 58	movq	%rcx,%rax
 59	movq	$0,8(%rcx)
 60	fstpt	(%rcx)
 61	ret
 62
 632:	fyl2xp1
 64	movq	%rcx,%rax
 65	movq	$0,8(%rcx)
 66	fstpt	(%rcx)
 67	ret
 68
 693:	jp	4b		// in case x is �Inf
 70	fstp	%st(1)
 71	fstp	%st(1)
 72	movq	%rcx,%rax
 73	movq	$0,8(%rcx)
 74	fstpt	(%rcx)
 75	ret
 76#else
 77	fldln2
 78	fldt	4(%esp)
 79	fxam
 80	fnstsw
 81	fld	%st
 82	sahf
 83	jc	3f		// in case x is NaN or �Inf
 844:
 85	fabs
 86	fldt	limit
 87	fcompp
 88	fnstsw
 89	sahf
 90	jnc	2f
 91	faddl	one
 92	fyl2x
 93	ret
 94
 952:	fyl2xp1
 96	ret
 97
 983:	jp	4b		// in case x is �Inf
 99	fstp	%st(1)
100	fstp	%st(1)
101	ret
102#endif