linux-sgx-trts-modified/external/sgx_libm/ia32/hypotf.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 "hypotf.c"
	.text
..TXTST0:
# -- Begin  hypotf
	.text
       .align    16,0x90
	.globl hypotf
hypotf:
# parameter 1: 8 + %ebp
# parameter 2: 12 + %ebp
..B1.1:
..L1:

        pushl     %ebp
        movl      %esp, %ebp
        andl      $-16, %esp
        subl      $16, %esp
        movzwl    10(%ebp), %ecx
        andl      $32640, %ecx
        movzwl    14(%ebp), %edx
        andl      $32640, %edx
        shrl      $7, %ecx
        shrl      $7, %edx
        cmpl      $255, %ecx
        je        ..B1.9
..B1.2:
        cmpl      $255, %edx
        je        ..B1.9
..B1.3:
        movl      %ecx, %eax
        subl      %edx, %eax
        addl      $28, %eax
        cmpl      $56, %eax
        ja        ..B1.5
..B1.4:
        flds      8(%ebp)
        addl      %edx, %ecx
        fmul      %st(0), %st
        cmpl      $378, %ecx
        flds      12(%ebp)
        fmul      %st(0), %st
        faddp     %st, %st(1)
        fsqrt     
        jg        ..B1.7
        jmp       ..B1.6
..B1.5:
        flds      8(%ebp)
        call      ..L2
..L2:
        popl      %eax
        lea       _GLOBAL_OFFSET_TABLE_+[. - ..L2](%eax), %eax
        fabs      
        flds      12(%ebp)
        fabs      
        fldl      .L_2il0floatpacket.0@GOTOFF(%eax)
        fmul      %st, %st(2)
        fmulp     %st, %st(1)
        faddp     %st, %st(1)
        fstps     8(%esp)
        flds      8(%esp)
        fmuls     .L_2il0floatpacket.1@GOTOFF(%eax)
        cmpl      $2130706432, 8(%esp)
..B1.6:
        movl      %ebp, %esp
        popl      %ebp
        ret       
..B1.7:
        fstl      (%esp)
        cmpl      $1206910976, 4(%esp)
        jl        ..B1.6
..B1.8:
        movl      %ebp, %esp
        popl      %ebp
        ret       
..B1.9:
        andb      $127, 11(%ebp)
        andb      $127, 15(%ebp)
        cmpl      $255, %ecx
        je        ..B1.21
..B1.10:
        cmpl      $255, %edx
        je        ..B1.19
..B1.11:
        flds      8(%ebp)
        fld       %st(0)
        fadds     12(%ebp)
        fstps     8(%esp)
..B1.12:
        cmpl      $255, %ecx
        je        ..B1.17
..B1.13:
        fstp      %st(0)
        cmpl      $255, %edx
        jne       ..B1.16
..B1.14:
        testl     $8388607, 12(%ebp)
        jne       ..B1.16
..B1.15:
        flds      12(%ebp)
        movl      %ebp, %esp
        popl      %ebp
        ret       
..B1.16:
        flds      8(%esp)
        movl      %ebp, %esp
        popl      %ebp
        ret       
..B1.17:
        testl     $8388607, 8(%ebp)
        jne       ..B1.13
..B1.18:
        movl      %ebp, %esp
        popl      %ebp
        ret       
..B1.19:
        testl     $8388607, 12(%ebp)
        je        ..B1.11
..B1.20:
        flds      8(%ebp)
        fld       %st(0)
        fmuls     12(%ebp)
        fstps     8(%esp)
        jmp       ..B1.12
..B1.21:
        testl     $8388607, 8(%ebp)
        jne       ..B1.20
        jmp       ..B1.10
        .align    16,0x90
	.type	hypotf,@function
	.size	hypotf,.-hypotf
	.data
# -- End  hypotf
	.section .rodata, "a"
	.align 8
	.align 8
.L_2il0floatpacket.0:
	.long	0x00000000,0x3fe00000
	.type	.L_2il0floatpacket.0,@object
	.size	.L_2il0floatpacket.0,8
	.align 4
.L_2il0floatpacket.1:
	.long	0x40000000
	.type	.L_2il0floatpacket.1,@object
	.size	.L_2il0floatpacket.1,4
	.data
	.section .note.GNU-stack, ""
# End