/* Function sinf vectorized with SSE4.

   Copyright (C) 2014-2018 Free Software Foundation, Inc.

   This file is part of the GNU C Library.

   The GNU C Library is free software; you can redistribute it and/or

   modify it under the terms of the GNU Lesser General Public

   License as published by the Free Software Foundation; either

   version 2.1 of the License, or (at your option) any later version.

   The GNU C Library is distributed in the hope that it will be useful,

   but WITHOUT ANY WARRANTY; without even the implied warranty of

   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

   Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public

   License along with the GNU C Library; if not, see

   <http://www.gnu.org/licenses/>.  */

#include <sysdep.h>

#include "svml_s_trig_data.h"

	.text

ENTRY(_ZGVbN4v_sinf_sse4)

/*

   ALGORITHM DESCRIPTION:

   1) Range reduction to [-Pi/2; +Pi/2] interval

      a) Grab sign from source argument and save it.

      b) Remove sign using AND operation

      c) Getting octant Y by 1/Pi multiplication

      d) Add "Right Shifter" value

      e) Treat obtained value as integer for destination sign setting.

         Shift first bit of this value to the last (sign) position

      f) Change destination sign if source sign is negative

         using XOR operation.

      g) Subtract "Right Shifter" value

      h) Subtract Y*PI from X argument, where PI divided to 4 parts:

         X = X - Y*PI1 - Y*PI2 - Y*PI3 - Y*PI4;

   2) Polynomial (minimax for sin within [-Pi/2; +Pi/2] interval)

      a) Calculate X^2 = X * X

      b) Calculate polynomial:

         R = X + X * X^2 * (A3 + x^2 * (A5 + ......

   3) Destination sign setting

      a) Set shifted destination sign using XOR operation:

         R = XOR( R, S );

 */

        pushq     %rbp

        cfi_adjust_cfa_offset (8)

        cfi_rel_offset (%rbp, 0)

        movq      %rsp, %rbp

        cfi_def_cfa_register (%rbp)

        andq      $-64, %rsp

        subq      $320, %rsp

        movaps    %xmm0, %xmm5

        movq      __svml_s_trig_data@GOTPCREL(%rip), %rax

        movups    __sAbsMask(%rax), %xmm2

/* b) Remove sign using AND operation */

        movaps    %xmm2, %xmm4

/*

  f) Change destination sign if source sign is negative

  using XOR operation.

 */

        andnps    %xmm5, %xmm2

        movups    __sInvPI(%rax), %xmm1

        andps     %xmm5, %xmm4

/* c) Getting octant Y by 1/Pi multiplication

   d) Add "Right Shifter" value  */

        mulps     %xmm4, %xmm1

/* h) Subtract Y*PI from X argument, where PI divided to 4 parts:

   X = X - Y*PI1 - Y*PI2 - Y*PI3 - Y*PI4  */

        movaps    %xmm4, %xmm0

/* Check for large and special values */

        cmpnleps  __sRangeReductionVal(%rax), %xmm4

        movups    __sRShifter(%rax), %xmm6

        movups    __sPI1(%rax), %xmm7

        addps     %xmm6, %xmm1

        movmskps  %xmm4, %ecx

/* e) Treat obtained value as integer for destination sign setting.

   Shift first bit of this value to the last (sign) position  */

        movaps    %xmm1, %xmm3

/* g) Subtract "Right Shifter" value */

        subps     %xmm6, %xmm1

        mulps     %xmm1, %xmm7

        pslld     $31, %xmm3

        movups    __sPI2(%rax), %xmm6

        subps     %xmm7, %xmm0

        mulps     %xmm1, %xmm6

        movups    __sPI3(%rax), %xmm7

        subps     %xmm6, %xmm0

        mulps     %xmm1, %xmm7

        movups    __sPI4(%rax), %xmm6

        subps     %xmm7, %xmm0

        mulps     %xmm6, %xmm1

        subps     %xmm1, %xmm0

/* 2) Polynomial (minimax for sin within [-Pi/2; +Pi/2] interval)

   a) Calculate X^2 = X * X

   b) Calculate polynomial:

   R = X + X * X^2 * (A3 + x^2 * (A5 + ...... */

        movaps    %xmm0, %xmm1

        mulps     %xmm0, %xmm1

        xorps     %xmm3, %xmm0

        movups    __sA9(%rax), %xmm3

        mulps     %xmm1, %xmm3

        addps     __sA7(%rax), %xmm3

        mulps     %xmm1, %xmm3

        addps     __sA5(%rax), %xmm3

        mulps     %xmm1, %xmm3

        addps     __sA3(%rax), %xmm3

        mulps     %xmm3, %xmm1

        mulps     %xmm0, %xmm1

        addps     %xmm1, %xmm0

/* 3) Destination sign setting

   a) Set shifted destination sign using XOR operation:

   R = XOR( R, S ); */

        xorps     %xmm2, %xmm0

        testl     %ecx, %ecx

        jne       .LBL_1_3

.LBL_1_2:

        cfi_remember_state

        movq      %rbp, %rsp

        cfi_def_cfa_register (%rsp)

        popq      %rbp

        cfi_adjust_cfa_offset (-8)

        cfi_restore (%rbp)

        ret

.LBL_1_3:

        cfi_restore_state

        movups    %xmm5, 192(%rsp)

        movups    %xmm0, 256(%rsp)

        je        .LBL_1_2

        xorb      %dl, %dl

        xorl      %eax, %eax

        movups    %xmm8, 112(%rsp)

        movups    %xmm9, 96(%rsp)

        movups    %xmm10, 80(%rsp)

        movups    %xmm11, 64(%rsp)

        movups    %xmm12, 48(%rsp)

        movups    %xmm13, 32(%rsp)

        movups    %xmm14, 16(%rsp)

        movups    %xmm15, (%rsp)

        movq      %rsi, 136(%rsp)

        movq      %rdi, 128(%rsp)

        movq      %r12, 168(%rsp)

        cfi_offset_rel_rsp (12, 168)

        movb      %dl, %r12b

        movq      %r13, 160(%rsp)

        cfi_offset_rel_rsp (13, 160)

        movl      %ecx, %r13d

        movq      %r14, 152(%rsp)

        cfi_offset_rel_rsp (14, 152)

        movl      %eax, %r14d

        movq      %r15, 144(%rsp)

        cfi_offset_rel_rsp (15, 144)

        cfi_remember_state

.LBL_1_6:

        btl       %r14d, %r13d

        jc        .LBL_1_12

.LBL_1_7:

        lea       1(%r14), %esi

        btl       %esi, %r13d

        jc        .LBL_1_10

.LBL_1_8:

        incb      %r12b

        addl      $2, %r14d

        cmpb      $16, %r12b

        jb        .LBL_1_6

        movups    112(%rsp), %xmm8

        movups    96(%rsp), %xmm9

        movups    80(%rsp), %xmm10

        movups    64(%rsp), %xmm11

        movups    48(%rsp), %xmm12

        movups    32(%rsp), %xmm13

        movups    16(%rsp), %xmm14

        movups    (%rsp), %xmm15

        movq      136(%rsp), %rsi

        movq      128(%rsp), %rdi

        movq      168(%rsp), %r12

        cfi_restore (%r12)

        movq      160(%rsp), %r13

        cfi_restore (%r13)

        movq      152(%rsp), %r14

        cfi_restore (%r14)

        movq      144(%rsp), %r15

        cfi_restore (%r15)

        movups    256(%rsp), %xmm0

        jmp       .LBL_1_2

.LBL_1_10:

        cfi_restore_state

        movzbl    %r12b, %r15d

        movss     196(%rsp,%r15,8), %xmm0

        call      JUMPTARGET(sinf)

        movss     %xmm0, 260(%rsp,%r15,8)

        jmp       .LBL_1_8

.LBL_1_12:

        movzbl    %r12b, %r15d

        movss     192(%rsp,%r15,8), %xmm0

        call      JUMPTARGET(sinf)

        movss     %xmm0, 256(%rsp,%r15,8)

        jmp       .LBL_1_7

END(_ZGVbN4v_sinf_sse4)