linux-sgx-trts-modified/external/sgx_libm/intel64/libm_scalbl_k80.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 "libm_scalbl_k80.c"
	.text
..TXTST0:
# -- Begin  __libm_scalbl_k80
	.text
       .align    16,0x90
	.globl __libm_scalbl_k80
__libm_scalbl_k80:
# parameter 1: 8 + %rsp
# parameter 2: 24 + %rsp
# parameter 3: %edi
..B1.1:
	.cfi_startproc
..___tag_value___libm_scalbl_k80.1:
..L2:

        fldt      24(%rsp)
        fldt      8(%rsp)
        fld       %st(0)
        cmpq      $0, 8(%rsp)
        fadd      %st(2), %st
        jne       ..B1.3
..B1.2:
        fstp      %st(2)
        fstp      %st(1)
        fstpt     -24(%rsp)
        fldt      -24(%rsp)
        ret       
..B1.3:
        movzwl    16(%rsp), %eax
        andl      $32767, %eax
        cmpl      $32767, %eax
        je        ..B1.2
..B1.4:
        fstpt     -24(%rsp)
        fldt      -24(%rsp)
        movzwl    -16(%rsp), %ecx
        movl      %ecx, %eax
        andl      $32767, %eax
        lea       -1(%rdi,%rax), %edx
        cmpl      $32765, %edx
        lea       (%rdi,%rax), %esi
        ja        ..B1.6
..B1.5:
        fstp      %st(1)
        fstp      %st(1)
        fstp      %st(0)
        andl      $-32768, %ecx
        andl      $32767, %esi
        orl       %esi, %ecx
        movw      %cx, -16(%rsp)
        fldt      -24(%rsp)
        ret       
..B1.6:
        cmpl      $32766, %esi
        jle       ..B1.8
..B1.7:
        fstp      %st(2)
        fstp      %st(1)
        lea       _CONSTANTS(%rip), %rax
        fldt      (%rax)
        fmul      %st(0), %st
        fmulp     %st, %st(1)
        fstpt     -24(%rsp)
        fldt      -24(%rsp)
        ret       
..B1.8:
        lea       16+_CONSTANTS(%rip), %rax
        cmpl      $-63, %esi
        fldt      (%rax)
        jge       ..B1.10
..B1.9:
        fstp      %st(3)
        fstp      %st(0)
        fxch      %st(1)
        fmul      %st(0), %st
        fmulp     %st, %st(1)
        fstpt     -24(%rsp)
        fldt      -24(%rsp)
        ret       
..B1.10:
        cmpl      $-19, %esi
        jge       ..B1.12
..B1.11:
        fstp      %st(2)
        fstp      %st(2)
        fldt      .L_2il0floatpacket.0(%rip)
        movzwl    8+.L_2il0floatpacket.0(%rip), %eax
        andl      $-32768, %eax
        fstpt     -56(%rsp)
        fmul      %st(0), %st
        faddp     %st, %st(1)
        fstpt     -24(%rsp)
        movzwl    -16(%rsp), %edx
        lea       256(%rdi,%rdx), %edi
        andl      $-32768, %edx
        andl      $32767, %edi
        orl       $-49409, %edx
        orl       %edi, %eax
        movw      %ax, -48(%rsp)
        fldt      -56(%rsp)
        movw      %dx, -16(%rsp)
        fldt      -24(%rsp)
        fmulp     %st, %st(1)
        fstpt     -24(%rsp)
        fldt      -24(%rsp)
        ret       
..B1.12:
        fstp      %st(1)
        fmul      %st(0), %st
        lea       32+_CONSTANTS(%rip), %rcx
        movzwl    8+.L_2il0floatpacket.0(%rip), %r8d
        faddp     %st, %st(1)
        fldt      .L_2il0floatpacket.1(%rip)
        andl      $-32768, %r8d
        fmul      %st(1), %st
        fld       %st(0)
        fsub      %st(2), %st
        fsubrp    %st, %st(1)
        fstpt     -40(%rsp)
        fldt      -40(%rsp)
        fsubrp    %st, %st(1)
        movzwl    -32(%rsp), %esi
        faddp     %st, %st(1)
        fldt      .L_2il0floatpacket.0(%rip)
        lea       -16383(%rsi), %eax
        negl      %eax
        movl      %esi, %edx
        addl      $16383, %eax
        lea       -32512(%rdi,%rsi), %edi
        andl      $32767, %eax
        andl      $-32768, %edx
        fstpt     -56(%rsp)
        orl       %eax, %r8d
        fldt      (%rcx)
        andl      $32767, %edi
        fmulp     %st, %st(1)
        orl       $-49409, %edx
        movw      %r8w, -48(%rsp)
        andl      $-32768, %r8d
        fldt      -56(%rsp)
        orl       %edi, %r8d
        fmulp     %st, %st(1)
        movw      %dx, -32(%rsp)
        fldt      -40(%rsp)
        movw      %r8w, -48(%rsp)
        fldt      -56(%rsp)
        fmul      %st, %st(1)
        fxch      %st(1)
        fstpt     -40(%rsp)
        fldt      -40(%rsp)
        fxch      %st(2)
        fmulp     %st, %st(1)
        faddp     %st, %st(1)
        fstpt     -24(%rsp)
        fldt      -24(%rsp)
..B1.13:
        ret       
        .align    16,0x90
	.cfi_endproc
	.type	__libm_scalbl_k80,@function
	.size	__libm_scalbl_k80,.-__libm_scalbl_k80
	.data
# -- End  __libm_scalbl_k80
	.section .rodata, "a"
	.align 16
	.align 16
.L_2il0floatpacket.0:
	.byte	0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0xff,0x3f,0x00,0x00,0x00,0x00,0x00,0x00
	.type	.L_2il0floatpacket.0,@object
	.size	.L_2il0floatpacket.0,16
	.align 16
.L_2il0floatpacket.1:
	.byte	0x00,0x00,0x00,0x40,0x00,0x00,0x00,0x80,0x20,0x40,0x00,0x00,0x00,0x00,0x00,0x00
	.type	.L_2il0floatpacket.1,@object
	.size	.L_2il0floatpacket.1,16
	.align 2
_CONSTANTS:
	.word	0
	.word	0
	.word	0
	.word	32768
	.word	32766
	.word	0
	.word	0
	.word	0
	.word	0
	.word	0
	.word	0
	.word	32768
	.word	1
	.word	0
	.word	0
	.word	0
	.word	0
	.word	0
	.word	0
	.word	32768
	.word	16127
	.word	0
	.word	0
	.word	0
	.type	_CONSTANTS,@object
	.size	_CONSTANTS,48
	.data
	.section .note.GNU-stack, ""
// -- Begin DWARF2 SEGMENT .eh_frame
	.section .eh_frame,"a",@progbits
.eh_frame_seg:
	.align 1
# End