linux-sgx-trts-modified/external/sgx_libm/intel64/scalb_gen.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 "scalb_gen.c"
	.text
..TXTST0:
# -- Begin  scalb
	.text
       .align    16,0x90
	.globl scalb
scalb:
# parameter 1: %xmm0
# parameter 2: %xmm1
..B1.1:
	.cfi_startproc
..___tag_value_scalb.1:
..L2:

        movaps    %xmm1, %xmm2
        movsd     %xmm0, -16(%rsp)
        movzwl    -10(%rsp), %eax
        movl      %eax, %edx
        movsd     %xmm2, -8(%rsp)
        andl      $32752, %edx
        movzwl    -2(%rsp), %ecx
        andl      $32752, %ecx
        shrl      $4, %edx
        shrl      $4, %ecx
        cmpl      $2047, %edx
        je        ..B1.56
..B1.2:
        cmpl      $2047, %ecx
        je        ..B1.37
..B1.3:
        cmpl      $1023, %ecx
        jl        ..B1.20
..B1.4:
        cmpl      $1044, %ecx
        jge       ..B1.7
..B1.5:
        lea       _TWO_52H(%rip), %rcx
        pxor      %xmm1, %xmm1
        movsd     (%rcx), %xmm0
        addsd     %xmm2, %xmm0
        movsd     %xmm0, -24(%rsp)
        movl      -24(%rsp), %edi
        cvtsi2sd  %edi, %xmm1
        ucomisd   %xmm2, %xmm1
        jp        ..B1.6
        je        ..B1.10
..B1.6:
        lea       _infs(%rip), %rax
        lea       _zeros(%rip), %rdx
        movsd     (%rax), %xmm0
        mulsd     (%rdx), %xmm0
        ret       
..B1.7:
        movb      -1(%rsp), %sil
        lea       large_integers(%rip), %rdi
        andb      $-128, %sil
        shrb      $7, %sil
        movzbl    %sil, %r8d
        cmpl      $1075, %ecx
        movl      (%rdi,%r8,4), %edi
        jge       ..B1.10
..B1.8:
        addl      $13, %ecx
        movl      -8(%rsp), %esi
        shll      %cl, %esi
        testl     %esi, %esi
        jne       ..B1.6
..B1.10:
        testl     %edx, %edx
        jne       ..B1.15
..B1.11:
        testl     $1048575, -12(%rsp)
        jne       ..B1.14
..B1.12:
        cmpl      $0, -16(%rsp)
        jne       ..B1.14
..B1.13:
        movsd     -16(%rsp), %xmm0
        ret       
..B1.14:
        lea       _TWO_55(%rip), %rax
        movsd     -16(%rsp), %xmm0
        mulsd     (%rax), %xmm0
        movsd     %xmm0, -16(%rsp)
        movzwl    -10(%rsp), %eax
        movl      %eax, %edx
        andl      $32752, %edx
        shrl      $4, %edx
        addl      $-55, %edx
..B1.15:
        movl      $65536, %ecx
        cmpl      $65536, %edi
        cmovg     %ecx, %edi
        movl      $-65536, %esi
        cmpl      $-65536, %edi
        cmovge    %edi, %esi
        lea       (%rdx,%rsi), %ecx
        testl     %ecx, %ecx
        jle       ..B1.33
..B1.16:
        cmpl      $2047, %ecx
        jl        ..B1.18
..B1.17:
        movb      -9(%rsp), %al
        lea       _large_value_64(%rip), %rcx
        andb      $-128, %al
        shrb      $7, %al
        movzbl    %al, %edx
        movsd     (%rcx,%rdx,8), %xmm0
        mulsd     (%rcx), %xmm0
        ret       
..B1.18:
        andl      $2047, %ecx
        andl      $-32753, %eax
        shll      $4, %ecx
        orl       %ecx, %eax
        movw      %ax, -10(%rsp)
        movsd     -16(%rsp), %xmm0
..B1.19:
        ret       
..B1.20:
        testl     %ecx, %ecx
        jne       ..B1.24
..B1.21:
        testl     $1048575, -4(%rsp)
        jne       ..B1.24
..B1.22:
        cmpl      $0, -8(%rsp)
        jne       ..B1.24
..B1.23:
        lea       _ones(%rip), %rax
        movsd     -16(%rsp), %xmm0
        mulsd     (%rax), %xmm0
        ret       
..B1.24:
        cmpl      $2047, %edx
        jne       ..B1.6
..B1.25:
        testl     %ecx, %ecx
        jle       ..B1.6
..B1.26:
        cmpl      $1, %ecx
        jne       ..B1.30
..B1.27:
        movl      -4(%rsp), %eax
        andl      $1048575, %eax
..B1.28:
        jne       ..B1.30
..B1.29:
        cmpl      $0, -8(%rsp)
        jbe       ..B1.6
..B1.30:
        movb      -9(%rsp), %al
        lea       _infs(%rip), %rdx
        andb      $-128, %al
        shrb      $7, %al
        movzbl    %al, %ecx
        movsd     (%rdx,%rcx,8), %xmm0
        ret       
..B1.33:
        cmpl      $-52, %ecx
        jge       ..B1.35
..B1.34:
        movb      -9(%rsp), %al
        lea       _small_value_64(%rip), %rcx
        andb      $-128, %al
        shrb      $7, %al
        movzbl    %al, %edx
        movsd     (%rcx,%rdx,8), %xmm0
        mulsd     (%rcx), %xmm0
        ret       
..B1.35:
        andl      $-32753, %eax
        lea       55(%rdx,%rsi), %edx
        andl      $2047, %edx
        shll      $4, %edx
        orl       %edx, %eax
        movw      %ax, -10(%rsp)
        lea       8+_TWO_55(%rip), %rax
        movsd     -16(%rsp), %xmm0
        mulsd     (%rax), %xmm0
..B1.36:
        ret       
..B1.37:
        testl     $1048575, -4(%rsp)
        jne       ..B1.51
..B1.38:
        cmpl      $0, -8(%rsp)
        jne       ..B1.51
..B1.39:
        testl     %edx, %edx
        jne       ..B1.46
..B1.40:
        testl     $1048575, -12(%rsp)
        jne       ..B1.47
..B1.41:
        cmpl      $0, -16(%rsp)
        jne       ..B1.47
..B1.42:
        testb     $-128, -1(%rsp)
        je        ..B1.44
..B1.43:
        movsd     -16(%rsp), %xmm0
        ret       
..B1.44:
        lea       _infs(%rip), %rax
        lea       _zeros(%rip), %rdx
        movsd     (%rax), %xmm0
        mulsd     (%rdx), %xmm0
..B1.45:
        ret       
..B1.46:
        cmpl      $2047, %edx
        je        ..B1.52
..B1.47:
        movb      -9(%rsp), %al
        andb      $-128, %al
        shrb      $7, %al
        movzbl    %al, %edx
        testb     $-128, -1(%rsp)
        je        ..B1.49
..B1.48:
        lea       _zeros(%rip), %rax
        movsd     (%rax,%rdx,8), %xmm0
        ret       
..B1.49:
        lea       _infs(%rip), %rax
        movsd     (%rax,%rdx,8), %xmm0
..B1.50:
        ret       
..B1.51:
        movsd     -16(%rsp), %xmm0
        addsd     -8(%rsp), %xmm0
        ret       
..B1.52:
        testb     $-128, -1(%rsp)
        jne       ..B1.54
..B1.53:
        movsd     -16(%rsp), %xmm0
        ret       
..B1.54:
        lea       _infs(%rip), %rax
        lea       _zeros(%rip), %rdx
        movsd     (%rax), %xmm0
        mulsd     (%rdx), %xmm0
..B1.55:
        ret       
..B1.56:
        testl     $1048575, -12(%rsp)
        jne       ..B1.51
..B1.57:
        cmpl      $0, -16(%rsp)
        jne       ..B1.51
..B1.58:
        cmpl      $2047, %ecx
        je        ..B1.37
        jmp       ..B1.20
        .align    16,0x90
	.cfi_endproc
	.type	scalb,@function
	.size	scalb,.-scalb
	.data
# -- End  scalb
	.section .rodata, "a"
	.align 4
	.align 4
_TWO_52H:
	.long	0
	.long	1127743488
	.type	_TWO_52H,@object
	.size	_TWO_52H,8
	.align 4
_infs:
	.long	0
	.long	2146435072
	.long	0
	.long	4293918720
	.type	_infs,@object
	.size	_infs,16
	.align 4
_zeros:
	.long	0
	.long	0
	.long	0
	.long	2147483648
	.type	_zeros,@object
	.size	_zeros,16
	.align 4
large_integers:
	.long	65536
	.long	-65536
	.type	large_integers,@object
	.size	large_integers,8
	.align 4
_TWO_55:
	.long	0
	.long	1130364928
	.long	0
	.long	1015021568
	.type	_TWO_55,@object
	.size	_TWO_55,16
	.align 4
_large_value_64:
	.long	0
	.long	2121269248
	.long	0
	.long	4268752896
	.type	_large_value_64,@object
	.size	_large_value_64,16
	.align 4
_ones:
	.long	0
	.long	1072693248
	.long	0
	.long	3220176896
	.type	_ones,@object
	.size	_ones,16
	.align 4
_small_value_64:
	.long	0
	.long	24117248
	.long	0
	.long	2171600896
	.type	_small_value_64,@object
	.size	_small_value_64,16
	.data
	.section .note.GNU-stack, ""
// -- Begin DWARF2 SEGMENT .eh_frame
	.section .eh_frame,"a",@progbits
.eh_frame_seg:
	.align 1
# End