linux-sgx-trts-modified/external/sgx_libm/intel64/dpml_ux_sqrt.S

/*
* Math library
*
* Copyright (C) 2016 Intel Corporation. All rights reserved.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
* are met:
*
*   * Redistributions of source code must retain the above copyright
*     notice, this list of conditions and the following disclaimer.
*   * Redistributions in binary form must reproduce the above copyright
*     notice, this list of conditions and the following disclaimer in
*     the documentation and/or other materials provided with the
*     distribution.
*   * Neither the name of Intel Corporation nor the names of its
*     contributors may be used to endorse or promote products derived
*     from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
*
* Author Name <jingwei.zhang@intel.com>
*   History:
*   03-14-2016 Initial version. numerics svn rev. 12864
*/
	.file "dpml_ux_sqrt.c"
	.text
..TXTST0:
# -- Begin  __dpml_ux_sqrt_evaluation__
	.text
       .align    16,0x90
	.globl __dpml_ux_sqrt_evaluation__
__dpml_ux_sqrt_evaluation__:
# parameter 1: %rdi
# parameter 2: %rsi
# parameter 3: %rdx
..B1.1:
	.cfi_startproc
..___tag_value___dpml_ux_sqrt_evaluation__.1:
..L2:

        pushq     %r12
	.cfi_def_cfa_offset 16
	.cfi_offset 12, -16
        pushq     %r13
	.cfi_def_cfa_offset 24
	.cfi_offset 13, -24
        pushq     %r14
	.cfi_def_cfa_offset 32
	.cfi_offset 14, -32
        pushq     %r15
	.cfi_def_cfa_offset 40
	.cfi_offset 15, -40
        pushq     %rbx
	.cfi_def_cfa_offset 48
	.cfi_offset 3, -48
        pushq     %rbp
	.cfi_def_cfa_offset 56
	.cfi_offset 6, -56
        subq      $104, %rsp
	.cfi_def_cfa_offset 160
        movq      %rdi, %rbx
        movq      %fs:40, %rax
        lea       (%rsp), %rdi
        xorq      %rsp, %rax
        movq      %rdx, %r12
        movq      %rax, 88(%rdi)
        movq      %rsi, %rbp
..___tag_value___dpml_ux_sqrt_evaluation__.16:
        call      fegetenv@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.17:
..B1.2:
        movq      8(%rbx), %r8
        movq      %r8, %r10
        shrq      $11, %r10
        movq      $0x3fd0000000000000, %r9
        movl      4(%rbx), %eax
        addq      %r9, %r10
        movl      %eax, %r9d
        movq      %r8, %r13
        andl      $1, %r9d
        movq      %r8, %r15
        addl      %r9d, %eax
        pxor      %xmm2, %xmm2
        shrq      $56, %r8
        pxor      %xmm4, %xmm4
        movq      %r10, 80(%rsp)
        pxor      %xmm7, %xmm7
        movq      __sqrt_t_table@GOTPCREL(%rip), %r10
        lea       40(%r9), %r14d
        shll      $7, %r9d
        movl      %r14d, %ecx
        xorq      %r9, %r8
        negl      %ecx
        shlq      $4, %r8
        pxor      %xmm6, %xmm6
        shlq      %cl, %r13
        movl      %r14d, %ecx
        movsd     80(%rsp), %xmm1
        lea       56(%rsp), %rdx
        movq      16(%rbx), %r11
        movaps    %xmm1, %xmm0
        cvtss2sd  4(%r10,%r8), %xmm2
        cvtss2sd  (%r10,%r8), %xmm4
        mulsd     %xmm1, %xmm0
        mulsd     %xmm1, %xmm2
        mulsd     %xmm0, %xmm4
        addsd     8(%r10,%r8), %xmm2
        shrq      %cl, %r11
        movl      %r14d, %ecx
        shrq      %cl, %r15
        orq       %r11, %r13
        shrq      $11, %r13
        lea       128+__sqrt_x_table(%rip), %r11
        cvtsi2sdq %r15, %xmm7
        cvtsi2sdq %r13, %xmm6
        addsd     %xmm2, %xmm4
        mulsd     (%r11), %xmm7
        lea       136+__sqrt_x_table(%rip), %r13
        lea       96+__sqrt_x_table(%rip), %r8
        movaps    %xmm7, %xmm8
        lea       152+__sqrt_x_table(%rip), %r15
        lea       104+__sqrt_x_table(%rip), %r14
        lea       32(%rsp), %rdi
        movsd     .L_2il0floatpacket.21(%rip), %xmm12
        movq      %rbx, %rsi
        mulsd     (%r13), %xmm6
        mulsd     (%r8), %xmm4
        addsd     %xmm6, %xmm8
        movaps    %xmm8, %xmm3
        lea       144+__sqrt_x_table(%rip), %r8
        mulsd     %xmm4, %xmm3
        cvtsd2ss  %xmm4, %xmm4
        cvtsd2ss  %xmm3, %xmm3
        cvtss2sd  %xmm4, %xmm4
        cvtss2sd  %xmm3, %xmm3
        mulsd     %xmm4, %xmm7
        mulsd     %xmm4, %xmm6
        subsd     %xmm3, %xmm7
        movsd     (%r15), %xmm10
        movaps    %xmm4, %xmm9
        mulsd     %xmm8, %xmm10
        movaps    %xmm3, %xmm5
        addsd     %xmm6, %xmm7
        mulsd     %xmm4, %xmm9
        mulsd     %xmm4, %xmm5
        mulsd     %xmm4, %xmm7
        mulsd     %xmm9, %xmm10
        movsd     (%r14), %xmm11
        movsd     (%r8), %xmm15
        lea       112+__sqrt_x_table(%rip), %r8
        sarl      $1, %eax
        subsd     %xmm5, %xmm11
        subsd     %xmm10, %xmm15
        subsd     %xmm7, %xmm11
        mulsd     %xmm4, %xmm15
        movsd     (%r8), %xmm13
        lea       120+__sqrt_x_table(%rip), %r8
        mulsd     %xmm4, %xmm13
        negl      %eax
        mulsd     %xmm11, %xmm15
        comisd    %xmm12, %xmm13
        mulsd     (%r8), %xmm15
        movaps    %xmm13, %xmm14
        cvttsd2si %xmm15, %r9
        subsd     %xmm12, %xmm14
        jae       ..L18
        movaps    %xmm13, %xmm14
..L18:
        cvttsd2si %xmm14, %r10
        movq      %r9, %r8
        incl      %eax
        shlq      $39, %r10
        sarq      $12, %r8
        sarq      $11, %r9
        addq      %r8, %r10
        andq      $1, %r9
        movq      $0x4000000000000000, %r8
        addq      %r9, %r10
        andq      %r10, %r8
        testq     %r10, %r10
        movl      $0, -24(%rdx)
        movl      %eax, -20(%rdx)
        movq      $0, -8(%rdx)
        lea       -1(,%r8,2), %r8
        cmovns    %r8, %r10
        movq      %r10, -16(%rdx)
..___tag_value___dpml_ux_sqrt_evaluation__.19:
        call      __dpml_multiply__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.20:
..B1.3:
        movq      %r12, %rdx
        lea       32(%rsp), %rdi
        lea       56(%rsp), %rsi
..___tag_value___dpml_ux_sqrt_evaluation__.21:
        call      __dpml_multiply__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.22:
..B1.4:
        lea       160+__sqrt_x_table(%rip), %rdi
        movq      %r12, %rsi
        movl      $9, %edx
        movq      %r12, %rcx
..___tag_value___dpml_ux_sqrt_evaluation__.23:
        call      __dpml_addsub__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.24:
..B1.5:
        testq     $1, %rbp
        lea       32(%rsp), %rsi
        movq      %r12, %rdi
        lea       56(%rsp), %rax
        cmovne    %rax, %rsi
        movq      %r12, %rdx
..___tag_value___dpml_ux_sqrt_evaluation__.25:
        call      __dpml_multiply__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.26:
..B1.6:
        decl      4(%r12)
        lea       (%rsp), %rdi
..___tag_value___dpml_ux_sqrt_evaluation__.27:
        call      fesetenv@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.28:
..B1.7:
        xorl      %r15d, %r15d
        testq     $2, %rbp
        jne       ..B1.11
..B1.8:
        movq      88(%rsp), %rdx
        xorq      %rsp, %rdx
        cmpq      %fs:40, %rdx
        jne       ..B1.10
..B1.9:
        movq      %r15, %rax
        addq      $104, %rsp
	.cfi_def_cfa_offset 56
	.cfi_restore 6
        popq      %rbp
	.cfi_def_cfa_offset 48
	.cfi_restore 3
        popq      %rbx
	.cfi_def_cfa_offset 40
	.cfi_restore 15
        popq      %r15
	.cfi_def_cfa_offset 32
	.cfi_restore 14
        popq      %r14
	.cfi_def_cfa_offset 24
	.cfi_restore 13
        popq      %r13
	.cfi_def_cfa_offset 16
	.cfi_restore 12
        popq      %r12
	.cfi_def_cfa_offset 8
        ret       
	.cfi_def_cfa_offset 160
	.cfi_offset 3, -48
	.cfi_offset 6, -56
	.cfi_offset 12, -16
	.cfi_offset 13, -24
	.cfi_offset 14, -32
	.cfi_offset 15, -40
..B1.10:
        call      __stack_chk_fail@PLT
..B1.11:
        movq      %r12, %rdi
        xorl      %esi, %esi
..___tag_value___dpml_ux_sqrt_evaluation__.49:
        call      __dpml_ffs_and_shift__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.50:
..B1.12:
        movq      16(%r12), %r14
        lea       8(%r14), %r8
        testq     $16368, %r8
        jne       ..B1.8
..B1.13:
        andq      $-32768, %r14
        movq      %r12, %rdi
        movq      %r12, %rsi
        lea       56(%rsp), %rdx
        lea       32(%rsp), %rcx
        lea       16384(%r14), %r13
        movq      %r13, 16(%r12)
..___tag_value___dpml_ux_sqrt_evaluation__.51:
        call      __dpml_extended_multiply__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.52:
..B1.14:
        movq      %rbx, %rdi
        lea       56(%rsp), %rsi
        movl      $1, %edx
        movq      %rsi, %rcx
..___tag_value___dpml_ux_sqrt_evaluation__.53:
        call      __dpml_addsub__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.54:
..B1.15:
        movl      $1, %edx
        lea       56(%rsp), %rdi
        movq      %rdi, %rcx
        lea       32(%rsp), %rsi
..___tag_value___dpml_ux_sqrt_evaluation__.55:
        call      __dpml_addsub__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.56:
..B1.16:
        movl      56(%rsp), %r8d
        testl     %r8d, %r8d
        cmove     %r13, %r14
        testq     $4, %rbp
        movq      %r14, 16(%r12)
        je        ..B1.8
..B1.17:
        movl      4(%r12), %ebx
        xorl      %edx, %edx
        addl      $-113, %ebx
        lea       56(%rsp), %rdi
        subl      %ebx, 4(%rdi)
        movl      $1, %ebx
        testl     %r8d, %r8d
        movq      %r12, %rsi
        cmove     %ebx, %edx
        movq      %rdi, %rcx
        orl       $8, %edx
..___tag_value___dpml_ux_sqrt_evaluation__.57:
        call      __dpml_addsub__@PLT
..___tag_value___dpml_ux_sqrt_evaluation__.58:
..B1.18:
        xorl      %r15d, %r15d
        cmpq      $8192, 72(%rsp)
        cmove     %ebx, %r15d
        jmp       ..B1.8
        .align    16,0x90
	.cfi_endproc
	.type	__dpml_ux_sqrt_evaluation__,@function
	.size	__dpml_ux_sqrt_evaluation__,.-__dpml_ux_sqrt_evaluation__
	.data
# -- End  __dpml_ux_sqrt_evaluation__
	.text
# -- Begin  __sqrtq
	.text
       .align    16,0x90
	.globl __sqrtq
__sqrtq:
# parameter 1: %xmm0
..B2.1:
	.cfi_startproc
..___tag_value___sqrtq.60:
..L61:

        subq      $136, %rsp
	.cfi_def_cfa_offset 144
        xorl      %esi, %esi
        lea       __sqrt_x_table(%rip), %rcx
        lea       48(%rsp), %rdi
        movq      %fs:40, %rax
        lea       88(%rsp), %rdx
        xorq      %rsp, %rax
        lea       112(%rsp), %r8
        movaps    %xmm0, -40(%rdx)
        lea       (%rsp), %r9
        movq      %rax, 40(%rdx)
        movq      $1, -80(%rdx)
        movq      $0, (%r9)
..___tag_value___sqrtq.63:
        call      __dpml_unpack_x_or_y__@PLT
..___tag_value___sqrtq.64:
..B2.2:
        testq     %rax, %rax
        jl        ..B2.5
..B2.3:
        movl      $7, %esi
        lea       88(%rsp), %rdi
        lea       64(%rsp), %rdx
..___tag_value___sqrtq.65:
        call      __dpml_ux_sqrt_evaluation__@PLT
..___tag_value___sqrtq.66:
..B2.4:
        movl      $1, %edx
        xorl      %ecx, %ecx
        testq     %rax, %rax
        lea       64(%rsp), %rdi
        cmove     %edx, %ecx
        xorl      %edx, %edx
        lea       112(%rsp), %rsi
        movq      %rcx, -104(%rsi)
        xorl      %ecx, %ecx
        lea       (%rsp), %r8
..___tag_value___sqrtq.67:
        call      __dpml_pack__@PLT
..___tag_value___sqrtq.68:
..B2.5:
        movq      128(%rsp), %rax
        xorq      %rsp, %rax
        movaps    112(%rsp), %xmm0
        cmpq      %fs:40, %rax
        jne       ..B2.7
..B2.6:
        addq      $136, %rsp
	.cfi_def_cfa_offset 8
        ret       
	.cfi_def_cfa_offset 144
..B2.7:
        call      __stack_chk_fail@PLT
        .align    16,0x90
	.cfi_endproc
	.type	__sqrtq,@function
	.size	__sqrtq,.-__sqrtq
	.data
# -- End  __sqrtq
	.text
# -- Begin  __rsqrtq
	.text
       .align    16,0x90
	.globl __rsqrtq
__rsqrtq:
# parameter 1: %xmm0
..B3.1:
	.cfi_startproc
..___tag_value___rsqrtq.72:
..L73:

        subq      $136, %rsp
	.cfi_def_cfa_offset 144
        xorl      %esi, %esi
        lea       8+__sqrt_x_table(%rip), %rcx
        lea       48(%rsp), %rdi
        movq      %fs:40, %rax
        lea       88(%rsp), %rdx
        xorq      %rsp, %rax
        lea       112(%rsp), %r8
        movaps    %xmm0, -40(%rdx)
        lea       (%rsp), %r9
        movq      %rax, 40(%rdx)
        movq      $1, -80(%rdx)
        movq      $0, (%r9)
..___tag_value___rsqrtq.75:
        call      __dpml_unpack_x_or_y__@PLT
..___tag_value___rsqrtq.76:
..B3.2:
        testq     %rax, %rax
        jl        ..B3.5
..B3.3:
        xorl      %esi, %esi
        lea       88(%rsp), %rdi
        lea       64(%rsp), %rdx
..___tag_value___rsqrtq.77:
        call      __dpml_ux_sqrt_evaluation__@PLT
..___tag_value___rsqrtq.78:
..B3.4:
        movl      $1, %edx
        xorl      %ecx, %ecx
        testq     %rax, %rax
        lea       64(%rsp), %rdi
        cmove     %edx, %ecx
        xorl      %edx, %edx
        lea       112(%rsp), %rsi
        movq      %rcx, -104(%rsi)
        xorl      %ecx, %ecx
        lea       (%rsp), %r8
..___tag_value___rsqrtq.79:
        call      __dpml_pack__@PLT
..___tag_value___rsqrtq.80:
..B3.5:
        movq      128(%rsp), %rax
        xorq      %rsp, %rax
        movaps    112(%rsp), %xmm0
        cmpq      %fs:40, %rax
        jne       ..B3.7
..B3.6:
        addq      $136, %rsp
	.cfi_def_cfa_offset 8
        ret       
	.cfi_def_cfa_offset 144
..B3.7:
        call      __stack_chk_fail@PLT
        .align    16,0x90
	.cfi_endproc
	.type	__rsqrtq,@function
	.size	__rsqrtq,.-__rsqrtq
	.data
# -- End  __rsqrtq
	.text
# -- Begin  __dpml_ux_hypot__
	.text
       .align    16,0x90
	.globl __dpml_ux_hypot__
__dpml_ux_hypot__:
# parameter 1: %rdi
# parameter 2: %rsi
# parameter 3: %rdx
..B4.1:
	.cfi_startproc
..___tag_value___dpml_ux_hypot__.84:
..L85:

        pushq     %r12
	.cfi_def_cfa_offset 16
	.cfi_offset 12, -16
        pushq     %rbp
	.cfi_def_cfa_offset 24
	.cfi_offset 6, -24
        subq      $88, %rsp
	.cfi_def_cfa_offset 112
        movq      %rdx, %r12
        movq      %rsi, %rbp
        movq      %rdi, %rsi
        movq      %fs:40, %rax
        lea       (%rsp), %rdx
        xorq      %rsp, %rax
        movq      %rax, 72(%rdx)
..___tag_value___dpml_ux_hypot__.91:
        call      __dpml_multiply__@PLT
..___tag_value___dpml_ux_hypot__.92:
..B4.2:
        movq      %rbp, %rdi
        movq      %rbp, %rsi
        lea       24(%rsp), %rdx
..___tag_value___dpml_ux_hypot__.93:
        call      __dpml_multiply__@PLT
..___tag_value___dpml_ux_hypot__.94:
..B4.3:
        xorl      %edx, %edx
        lea       (%rsp), %rdi
        lea       24(%rsp), %rsi
        lea       48(%rsp), %rcx
..___tag_value___dpml_ux_hypot__.95:
        call      __dpml_addsub__@PLT
..___tag_value___dpml_ux_hypot__.96:
..B4.4:
        xorl      %esi, %esi
        lea       48(%rsp), %rdi
..___tag_value___dpml_ux_hypot__.97:
        call      __dpml_ffs_and_shift__@PLT
..___tag_value___dpml_ux_hypot__.98:
..B4.5:
        movl      $7, %esi
        lea       48(%rsp), %rdi
        movq      %r12, %rdx
..___tag_value___dpml_ux_hypot__.99:
        call      __dpml_ux_sqrt_evaluation__@PLT
..___tag_value___dpml_ux_hypot__.100:
..B4.6:
        testq     %rax, %rax
        je        ..B4.18
..B4.7:
        movl      4(%rsp), %ecx
        lea       (%rsp), %rbp
        subl      28(%rbp), %ecx
        movslq    %ecx, %rcx
        testq     %rcx, %rcx
        jns       ..B4.9
..B4.8:
        negq      %rcx
        lea       24(%rsp), %rbp
..B4.9:
        cmpq      $128, %rcx
        jle       ..B4.11
..B4.10:
        xorl      %eax, %eax
        jmp       ..B4.18
..B4.11:
        movl      $1, %edx
..B4.12:
        movq      8(%rbp,%rdx,8), %rax
        testq     %rax, %rax
        jne       ..B4.10
..B4.13:
        cmpq      $64, %rcx
        jl        ..B4.16
..B4.14:
        addq      $-64, %rcx
        decq      %rdx
        jns       ..B4.12
..B4.16:
        movq      %rax, %rdx
        movl      $1, %ebp
        shrq      %cl, %rdx
        xorl      %esi, %esi
        shlq      %cl, %rdx
        cmpq      %rdx, %rax
        cmove     %ebp, %esi
        movl      %esi, %eax
..B4.18:
        movq      72(%rsp), %rdx
        xorq      %rsp, %rdx
        cmpq      %fs:40, %rdx
        jne       ..B4.20
..B4.19:
        addq      $88, %rsp
	.cfi_def_cfa_offset 24
	.cfi_restore 6
        popq      %rbp
	.cfi_def_cfa_offset 16
	.cfi_restore 12
        popq      %r12
	.cfi_def_cfa_offset 8
        ret       
	.cfi_def_cfa_offset 112
	.cfi_offset 6, -24
	.cfi_offset 12, -16
..B4.20:
        call      __stack_chk_fail@PLT
        .align    16,0x90
	.cfi_endproc
	.type	__dpml_ux_hypot__,@function
	.size	__dpml_ux_hypot__,.-__dpml_ux_hypot__
	.data
# -- End  __dpml_ux_hypot__
	.section .rodata, "a"
	.align 16
	.align 16
__sqrt_x_table:
	.long	1088750600
	.long	605086734
	.long	1088766984
	.long	519758862
	.long	112
	.long	0
	.long	113
	.long	0
	.long	114
	.long	0
	.long	0
	.long	0
	.long	0
	.long	268435456
	.long	858997845
	.long	34
	.long	1636176969
	.long	409044504
	.long	1099306057
	.long	404751376
	.long	545392672
	.long	404783624
	.long	142938632
	.long	302522498
	.long	1719614413
	.long	1073127582
	.long	0
	.long	1072693248
	.long	0
	.long	1097859072
	.long	0
	.long	1151336448
	.long	0
	.long	1047527424
	.long	0
	.long	991952896
	.long	0
	.long	1072431104
	.long	0
	.long	1071120384
	.long	0
	.long	2
	.long	0
	.long	3221225472
	.long	0
	.long	0
	.long	0
	.long	0
	.type	__sqrt_x_table,@object
	.size	__sqrt_x_table,192
	.align 8
.L_2il0floatpacket.21:
	.long	0x00000000,0x43e00000
	.type	.L_2il0floatpacket.21,@object
	.size	.L_2il0floatpacket.21,8
	.align 4
.L_2il0floatpacket.20:
	.long	0x5f000000
	.type	.L_2il0floatpacket.20,@object
	.size	.L_2il0floatpacket.20,4
	.data
	.section .note.GNU-stack, ""
// -- Begin DWARF2 SEGMENT .eh_frame
	.section .eh_frame,"a",@progbits
.eh_frame_seg:
	.align 1
# End