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 103