diff mbox series

[v3,4/7] x86: Optimize svml_s_tanhf16_core_avx512.S

Message ID 20220609165835.2527760-4-goldstein.w.n@gmail.com
State New
Headers show
Series [v3,1/7] x86: Improve svml_s_atanhf16_core_avx512.S | expand

Commit Message

Noah Goldstein June 9, 2022, 4:58 p.m. UTC
Optimizations are:
    1. Reduce code size (-67 bytes).
    2. Remove redundant move instructions.
    3. Slightly improve instruction selection/scheduling where
       possible.
    4. Reduce rodata usage (-448 bytes).

Result is roughly a 14% speedup:

       Function, New Time, Old Time, New / Old
_ZGVeN16v_tanhf,    0.649,    0.752,     0.863
---
 .../multiarch/svml_s_tanhf16_core_avx512.S    | 527 ++++++++++--------
 1 file changed, 287 insertions(+), 240 deletions(-)

Comments

H.J. Lu June 9, 2022, 5:07 p.m. UTC | #1
On Thu, Jun 9, 2022 at 9:58 AM Noah Goldstein <goldstein.w.n@gmail.com> wrote:
>
> Optimizations are:
>     1. Reduce code size (-67 bytes).
>     2. Remove redundant move instructions.
>     3. Slightly improve instruction selection/scheduling where
>        possible.
>     4. Reduce rodata usage (-448 bytes).
>
> Result is roughly a 14% speedup:
>
>        Function, New Time, Old Time, New / Old
> _ZGVeN16v_tanhf,    0.649,    0.752,     0.863
> ---
>  .../multiarch/svml_s_tanhf16_core_avx512.S    | 527 ++++++++++--------
>  1 file changed, 287 insertions(+), 240 deletions(-)
>
> diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S b/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S
> index 5b1f9f151c..7edc74a116 100644
> --- a/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S
> +++ b/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S
> @@ -70,310 +70,357 @@
>   *
>   */
>
> -/* Offsets for data table __svml_stanh_data_internal
> +/* Offsets for data table __svml_stanh_data_internal_avx512. Ordered
> +   by use in the function. On cold-starts this might help the
> +   prefetcher. Possibly a better idea is to interleave start/end so
> +   that the prefetcher is less likely to detect a stream and pull
> +   irrelivant lines into cache.  */
> +
> +/* Offsets for data table __svml_stanh_data_internal. 4 bytes each.
>   */
> -#define _sC                            0
> -#define _sP0                           128
> -#define _sP2                           256
> -#define _sP3                           384
> -#define _sP4                           512
> -#define _sP5                           640
> -#define _sP6                           768
> -#define _sP7                           896
> -#define _iExpMantMask_UISA             1024
> -#define _iMinIdxOfsMask_UISA           1088
> -#define _iMaxIdxMask_UISA              1152
> -#define _sSignMask                     1216
> -#define _sAbsMask                      1280
> -#define _iExpMantMask                  1344
> -#define _iExpMask                      1408
> -#define _iMinIdxOfsMask                        1472
> -#define _iMaxIdxMask                   1536
> +#define _iExpMantMask_UISA             0
> +#define _iMinIdxOfsMask_UISA           4
> +#define _iMaxIdxMask_UISA              8
> +#define _iExpMask                      12
> +
> +/* Offsets for data table __svml_stanh_data_internal_al64. 64 bytes
> +   each.  */
> +#define _sC_lo                         0
> +#define _sC_hi                         64
> +#define _sP7_lo                                128
> +#define _sP7_hi                                192
> +#define _sSignMask                     256
> +#define _sP6_lo                                320
> +#define _sP6_hi                                384
> +#define _sP5_lo                                448
> +#define _sP5_hi                                512
> +#define _sP4_lo                                576
> +#define _sP4_hi                                640
> +#define _sP3_lo                                704
> +#define _sP3_hi                                768
> +#define _sP2_lo                                832
> +#define _sP2_hi                                896
> +#define _sP0_lo                                960
> +#define _sP0_hi                                1024
>
>  #include <sysdep.h>
> +#define TANHF_DATA(x)                  ((x)+__svml_stanh_data_internal_al64)
> +#define TANHF_DATA_UNALIGNED(x)                ((x)+__svml_stanh_data_internal)
>
>         .section .text.exex512, "ax", @progbits
>  ENTRY(_ZGVeN16v_tanhf_skx)
> -       pushq   %rbp
> -       cfi_def_cfa_offset(16)
> -       movq    %rsp, %rbp
> -       cfi_def_cfa(6, 16)
> -       cfi_offset(6, -16)
> -       andq    $-64, %rsp
> -       subq    $192, %rsp
> -       vmovaps %zmm0, %zmm1
> -       vmovups __svml_stanh_data_internal(%rip), %zmm9
> -       vmovups _sP6+__svml_stanh_data_internal(%rip), %zmm11
> -       vmovups _sP5+__svml_stanh_data_internal(%rip), %zmm12
> -       vmovups _sP4+__svml_stanh_data_internal(%rip), %zmm13
> -       vmovups _sP3+__svml_stanh_data_internal(%rip), %zmm14
> -       vmovups _sP2+__svml_stanh_data_internal(%rip), %zmm15
> -       vpternlogd $255, %zmm2, %zmm2, %zmm2
> -       vandps  _sAbsMask+__svml_stanh_data_internal(%rip), %zmm1, %zmm8
> -       vandps  _sSignMask+__svml_stanh_data_internal(%rip), %zmm1, %zmm0
> -
>         /* Here huge arguments, INF and NaNs are filtered out to callout. */
> -       vpandd  _iExpMantMask_UISA+__svml_stanh_data_internal(%rip), %zmm1, %zmm3
> -       vpsubd  _iMinIdxOfsMask_UISA+__svml_stanh_data_internal(%rip), %zmm3, %zmm4
> -       vpcmpd  $2, _iExpMask+__svml_stanh_data_internal(%rip), %zmm3, %k1
> +       vpandd  TANHF_DATA_UNALIGNED(_iExpMantMask_UISA)(%rip){1to16}, %zmm0, %zmm1
> +       vpsubd  TANHF_DATA_UNALIGNED(_iMinIdxOfsMask_UISA)(%rip){1to16}, %zmm1, %zmm2
>
> -       /*
> -        *  small table specific variables *
> -        *  Constant loading
> -        */
> -       vpxord  %zmm5, %zmm5, %zmm5
> -
> -       /* if VMIN, VMAX is defined for I type */
> -       vpmaxsd %zmm5, %zmm4, %zmm6
> -       vpminsd _iMaxIdxMask_UISA+__svml_stanh_data_internal(%rip), %zmm6, %zmm7
> -       vpsrld  $21, %zmm7, %zmm10
> -       vmovups _sP7+__svml_stanh_data_internal(%rip), %zmm4
> -       vpermt2ps _sC+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm9
> -       vpermt2ps _sP6+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm11
> -       vpermt2ps _sP7+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm4
> -       vpermt2ps _sP5+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm12
> -       vpermt2ps _sP4+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm13
> -       vpermt2ps _sP3+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm14
> -       vpermt2ps _sP2+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm15
> -       vpandnd %zmm3, %zmm3, %zmm2{%k1}
> -       vptestmd %zmm2, %zmm2, %k0
> -       vmovups _sP0+__svml_stanh_data_internal(%rip), %zmm3
> -       vsubps  {rn-sae}, %zmm9, %zmm8, %zmm2
> -       kmovw   %k0, %edx
> -       vfmadd213ps {rn-sae}, %zmm11, %zmm2, %zmm4
> -       vpermt2ps _sP0+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm3
> -       vfmadd213ps {rn-sae}, %zmm12, %zmm2, %zmm4
> -       vfmadd213ps {rn-sae}, %zmm13, %zmm2, %zmm4
> -       vfmadd213ps {rn-sae}, %zmm14, %zmm2, %zmm4
> -       vfmadd213ps {rn-sae}, %zmm15, %zmm2, %zmm4
> -       vfmadd213ps {rn-sae}, %zmm3, %zmm2, %zmm4
> -       vorps   %zmm0, %zmm4, %zmm0
> -       testl   %edx, %edx
> +       /* Selection arguments between [0, 0x03e00000] into zmm3.  */
> +       vpxord  %zmm3, %zmm3, %zmm3
> +       vpmaxsd %zmm3, %zmm2, %zmm3
> +       vpminsd TANHF_DATA_UNALIGNED(_iMaxIdxMask_UISA)(%rip){1to16}, %zmm3, %zmm3
>
> -       /* Go to special inputs processing branch */
> -       jne     L(SPECIAL_VALUES_BRANCH)
> -       # LOE rbx r12 r13 r14 r15 edx zmm0 zmm1
> +       /* Setup permute indices in zmm3.  */
> +       vpsrld  $21, %zmm3, %zmm3
>
> -       /* Restore registers
> -        * and exit the function
> -        */
> +       /* Store if there are any special cases in k1.  */
> +       vpcmpd  $6, TANHF_DATA_UNALIGNED(_iExpMask)(%rip){1to16}, %zmm1, %k1
>
> -L(EXIT):
> -       movq    %rbp, %rsp
> -       popq    %rbp
> -       cfi_def_cfa(7, 8)
> -       cfi_restore(6)
> -       ret
> -       cfi_def_cfa(6, 16)
> -       cfi_offset(6, -16)
> +       vmovaps TANHF_DATA(_sC_lo)(%rip), %zmm5
> +       vpermt2ps TANHF_DATA(_sC_hi)(%rip), %zmm3, %zmm5
>
> -       /* Branch to process
> -        * special inputs
> -        */
> +       vmovaps TANHF_DATA(_sP7_lo)(%rip), %zmm2
> +       vpermt2ps TANHF_DATA(_sP7_hi)(%rip), %zmm3, %zmm2
>
> -L(SPECIAL_VALUES_BRANCH):
> -       vmovups %zmm1, 64(%rsp)
> -       vmovups %zmm0, 128(%rsp)
> -       # LOE rbx r12 r13 r14 r15 edx zmm0
> +       /* Store absolute values of inputs in zmm1.  */
> +       vmovaps TANHF_DATA(_sSignMask)(%rip), %zmm4
> +       vandnps %zmm0, %zmm4, %zmm1
> +       vsubps  {rn-sae}, %zmm5, %zmm1, %zmm1
>
> -       xorl    %eax, %eax
> -       # LOE rbx r12 r13 r14 r15 eax edx
> +       vmovaps TANHF_DATA(_sP6_lo)(%rip), %zmm5
> +       vpermt2ps TANHF_DATA(_sP6_hi)(%rip), %zmm3, %zmm5
>
> -       vzeroupper
> -       movq    %r12, 16(%rsp)
> -       /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -176; DW_OP_plus)  */
> -       .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x50, 0xff, 0xff, 0xff, 0x22
> -       movl    %eax, %r12d
> -       movq    %r13, 8(%rsp)
> -       /*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -184; DW_OP_plus)  */
> -       .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x48, 0xff, 0xff, 0xff, 0x22
> -       movl    %edx, %r13d
> -       movq    %r14, (%rsp)
> -       /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -192; DW_OP_plus)  */
> -       .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x40, 0xff, 0xff, 0xff, 0x22
> -       # LOE rbx r15 r12d r13d
> -
> -       /* Range mask
> -        * bits check
> -        */
> +       vmovaps TANHF_DATA(_sP5_lo)(%rip), %zmm6
> +       vpermt2ps TANHF_DATA(_sP5_hi)(%rip), %zmm3, %zmm6
> +
> +       vfmadd213ps {rn-sae}, %zmm5, %zmm1, %zmm2
> +       vfmadd213ps {rn-sae}, %zmm6, %zmm1, %zmm2
>
> -L(RANGEMASK_CHECK):
> -       btl     %r12d, %r13d
> +       vmovaps TANHF_DATA(_sP4_lo)(%rip), %zmm7
> +       vpermt2ps TANHF_DATA(_sP4_hi)(%rip), %zmm3, %zmm7
>
> -       /* Call scalar math function */
> -       jc      L(SCALAR_MATH_CALL)
> -       # LOE rbx r15 r12d r13d
> +       vmovaps TANHF_DATA(_sP3_lo)(%rip), %zmm8
> +       vpermt2ps TANHF_DATA(_sP3_hi)(%rip), %zmm3, %zmm8
>
> -       /* Special inputs
> -        * processing loop
> +       vfmadd213ps {rn-sae}, %zmm7, %zmm1, %zmm2
> +       vfmadd213ps {rn-sae}, %zmm8, %zmm1, %zmm2
> +
> +       vmovaps TANHF_DATA(_sP2_lo)(%rip), %zmm9
> +       vpermt2ps TANHF_DATA(_sP2_hi)(%rip), %zmm3, %zmm9
> +
> +       vmovaps TANHF_DATA(_sP0_lo)(%rip), %zmm10
> +       vpermt2ps TANHF_DATA(_sP0_hi)(%rip), %zmm3, %zmm10
> +
> +       vfmadd213ps {rn-sae}, %zmm9, %zmm1, %zmm2
> +       vfmadd213ps {rn-sae}, %zmm10, %zmm1, %zmm2
> +
> +       kmovw   %k1, %edx
> +       testl   %edx, %edx
> +
> +       /* Go to special inputs processing branch.  */
> +       jne     L(SPECIAL_VALUES_BRANCH)
> +       # LOE rbx r12 r13 r14 r15 zmm0 zmm2 zmm4
> +       /* Wait until after branch of write over zmm0.  */
> +       vpternlogd $0xec, %zmm4, %zmm2, %zmm0
> +
> +       /* No stack restoration on the fastpath.  */
> +       ret
> +
> +       /* Cold case. edx has 1s where there was a special value that
> +          needs to be handled by a tanhf call. Optimize for code size
> +          more so than speed here. */
> +L(SPECIAL_VALUES_BRANCH):
> +       # LOE rbx rdx r12 r13 r14 r15 zmm0 zmm2 zmm4
> +    /* Use r13 to save/restore the stack. This allows us to use rbp as
> +       callee save register saving code size. */
> +       pushq   %r13
> +       cfi_adjust_cfa_offset(8)
> +       cfi_offset(r13, -16)
> +       /* Need to callee save registers to preserve state across tanhf calls.
>          */
> +       pushq   %rbx
> +       cfi_adjust_cfa_offset(8)
> +       cfi_offset(rbx, -24)
> +       pushq   %rbp
> +       cfi_adjust_cfa_offset(8)
> +       cfi_offset(rbp, -32)
> +       movq    %rsp, %r13
> +       cfi_def_cfa_register(r13)
> +
> +       /* Align stack and make room for 2x zmm vectors.  */
> +       andq    $-64, %rsp
> +       addq    $-128, %rsp
> +
> +       /* Save original input (zmm0 unchanged up to this point).  */
> +       vmovaps %zmm0, 64(%rsp)
> +       /* Save all already computed inputs.  */
> +       vpternlogd $0xec, %zmm4, %zmm2, %zmm0
> +       vmovaps %zmm0, (%rsp)
>
> +       vzeroupper
> +
> +       /* edx has 1s where there was a special value that needs to be handled
> +          by a tanhf call.  */
> +       movl    %edx, %ebx
>  L(SPECIAL_VALUES_LOOP):
> -       incl    %r12d
> -       cmpl    $16, %r12d
> -
> -       /* Check bits in range mask */
> -       jl      L(RANGEMASK_CHECK)
> -       # LOE rbx r15 r12d r13d
> -
> -       movq    16(%rsp), %r12
> -       cfi_restore(12)
> -       movq    8(%rsp), %r13
> -       cfi_restore(13)
> -       movq    (%rsp), %r14
> -       cfi_restore(14)
> -       vmovups 128(%rsp), %zmm0
> -
> -       /* Go to exit */
> -       jmp     L(EXIT)
> -       /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -176; DW_OP_plus)  */
> -       .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x50, 0xff, 0xff, 0xff, 0x22
> -       /*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -184; DW_OP_plus)  */
> -       .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x48, 0xff, 0xff, 0xff, 0x22
> -       /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -192; DW_OP_plus)  */
> -       .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x40, 0xff, 0xff, 0xff, 0x22
> -       # LOE rbx r12 r13 r14 r15 zmm0
> -
> -       /* Scalar math fucntion call
> -        * to process special input
> -        */
> +       # LOE rbx rbp r12 r13 r14 r15
> +       /* use rbp as index for special value that is saved across calls to
> +          tanhf. We technically don't need a callee save register here as offset
> +          to rsp is always [0, 56] so we can restore rsp by realigning to 64.
> +          Essentially the tradeoff is 1 extra save/restore vs 2 extra instructions
> +          in the loop. Realigning also costs more code size.  */
> +       xorl    %ebp, %ebp
> +       tzcntl  %ebx, %ebp
>
> -L(SCALAR_MATH_CALL):
> -       movl    %r12d, %r14d
> -       movss   64(%rsp, %r14, 4), %xmm0
> +       /* Scalar math fucntion call to process special input.  */
> +       movss   64(%rsp, %rbp, 4), %xmm0
>         call    tanhf@PLT
> -       # LOE rbx r14 r15 r12d r13d xmm0
>
> -       movss   %xmm0, 128(%rsp, %r14, 4)
> +       /* No good way to avoid the store-forwarding fault this will cause on
> +          return. `lfence` avoids the SF fault but at greater cost as it
> +          serialized stack/callee save restoration.  */
> +       movss   %xmm0, (%rsp, %rbp, 4)
>
> -       /* Process special inputs in loop */
> -       jmp     L(SPECIAL_VALUES_LOOP)
> -       # LOE rbx r15 r12d r13d
> +       blsrl   %ebx, %ebx
> +       jnz     L(SPECIAL_VALUES_LOOP)
> +       # LOE r12 r13 r14 r15
> +
> +       /* All results have been written to (%rsp).  */
> +       vmovaps (%rsp), %zmm0
> +       /* Restore rsp.  */
> +       movq    %r13, %rsp
> +       cfi_def_cfa_register(rsp)
> +       /* Restore callee save registers.  */
> +       popq    %rbp
> +       cfi_adjust_cfa_offset(-8)
> +       cfi_restore(rbp)
> +       popq    %rbx
> +       cfi_adjust_cfa_offset(-8)
> +       cfi_restore(rbp)
> +       popq    %r13
> +       cfi_adjust_cfa_offset(-8)
> +       cfi_restore(r13)
> +       ret
>  END(_ZGVeN16v_tanhf_skx)
>
>         .section .rodata, "a"
> -       .align  64
> -
> +       .align  16
>  #ifdef __svml_stanh_data_internal_typedef
>  typedef unsigned int VUINT32;
> -typedef struct {
> -       __declspec(align(64)) VUINT32 _sC[32][1];
> -       __declspec(align(64)) VUINT32 _sP0[32][1];
> -       __declspec(align(64)) VUINT32 _sP2[32][1];
> -       __declspec(align(64)) VUINT32 _sP3[32][1];
> -       __declspec(align(64)) VUINT32 _sP4[32][1];
> -       __declspec(align(64)) VUINT32 _sP5[32][1];
> -       __declspec(align(64)) VUINT32 _sP6[32][1];
> -       __declspec(align(64)) VUINT32 _sP7[32][1];
> -       __declspec(align(64)) VUINT32 _iExpMantMask_UISA[16][1];
> -       __declspec(align(64)) VUINT32 _iMinIdxOfsMask_UISA[16][1];
> -       __declspec(align(64)) VUINT32 _iMaxIdxMask_UISA[16][1];
> +typedef struct
> +       {
> +       __declspec(align(4)) VUINT32 _iExpMantMask_UISA[1][1];
> +       __declspec(align(4)) VUINT32 _iMinIdxOfsMask_UISA[1][1];
> +       __declspec(align(4)) VUINT32 _iMaxIdxMask_UISA[1][1];
> +       __declspec(align(4)) VUINT32 _iExpMask[1][1];
> +       __declspec(align(64)) VUINT32 _sC_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sC_hi[16][1];
> +       __declspec(align(64)) VUINT32 _sP7_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP7_hi[16][1];
>         __declspec(align(64)) VUINT32 _sSignMask[16][1];
> -       __declspec(align(64)) VUINT32 _sAbsMask[16][1];
> -       __declspec(align(64)) VUINT32 _iExpMantMask[16][1];
> -       __declspec(align(64)) VUINT32 _iExpMask[16][1];
> -       __declspec(align(64)) VUINT32 _iMinIdxOfsMask[16][1];
> -       __declspec(align(64)) VUINT32 _iMaxIdxMask[16][1];
> +       __declspec(align(64)) VUINT32 _sP6_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP6_hi[16][1];
> +       __declspec(align(64)) VUINT32 _sP5_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP5_hi[16][1];
> +       __declspec(align(64)) VUINT32 _sP4_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP4_hi[16][1];
> +       __declspec(align(64)) VUINT32 _sP3_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP3_hi[16][1];
> +       __declspec(align(64)) VUINT32 _sP2_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP2_hi[16][1];
> +       __declspec(align(64)) VUINT32 _sP0_lo[16][1];
> +       __declspec(align(64)) VUINT32 _sP0_hi[16][1];
>  } __svml_stanh_data_internal;
>  #endif
> +
>  __svml_stanh_data_internal:
> -       /* _sC */
> +       .align  4
> +       /* _iExpMantMask_UISA */
> +       .long   0x7fe00000
> +
> +       .align  4
> +       /* _iMinIdxOfsMask_UISA */
> +       .long   0x3d400000
> +
> +       .align  4
> +       /* _iMaxIdxMask_UISA */
> +       .long   0x03e00000
> +
> +       .align  4
> +       /* _iExpMask */
> +       .long   0x7f000000
> +
> +       .align  64
> +__svml_stanh_data_internal_al64:
> +       .align  64
> +       /* _sC_lo */
>         .long   0x00000000, 0x3d700000, 0x3d900000, 0x3db00000
>         .long   0x3dd00000, 0x3df00000, 0x3e100000, 0x3e300000
>         .long   0x3e500000, 0x3e700000, 0x3e900000, 0x3eb00000
>         .long   0x3ed00000, 0x3ef00000, 0x3f100000, 0x3f300000
> +
> +       .align  64
> +       /* _sC_hi */
>         .long   0x3f500000, 0x3f700000, 0x3f900000, 0x3fb00000
>         .long   0x3fd00000, 0x3ff00000, 0x40100000, 0x40300000
>         .long   0x40500000, 0x40700000, 0x40900000, 0x40b00000
>         .long   0x40d00000, 0x40f00000, 0x41100000, 0x00000000
> -       /* p0 */
> -       .align  64
> -       .long   0x00000000, 0x3d6fb9c9, 0x3d8fc35f, 0x3daf9169
> -       .long   0x3dcf49ab, 0x3deee849, 0x3e0f0ee8, 0x3e2e4984
> -       .long   0x3e4d2f8e, 0x3e6bb32e, 0x3e8c51cd, 0x3ea96163
> -       .long   0x3ec543f1, 0x3edfd735, 0x3f028438, 0x3f18abf0
> -       .long   0x3f2bc480, 0x3f3bec1c, 0x3f4f2e5b, 0x3f613c53
> -       .long   0x3f6ce37d, 0x3f743c4f, 0x3f7a5feb, 0x3f7dea85
> -       .long   0x3f7f3b3d, 0x3f7fb78c, 0x3f7fefd4, 0x3f7ffdd0
> -       .long   0x3f7fffb4, 0x3f7ffff6, 0x3f7fffff, 0x3f800000
> -       /* p2 */
> -       .align  64
> -       .long   0x3f800000, 0x3f7f1f84, 0x3f7ebd11, 0x3f7e1e5f
> -       .long   0x3f7d609f, 0x3f7c842d, 0x3f7b00e5, 0x3f789580
> -       .long   0x3f75b8ad, 0x3f726fd9, 0x3f6cc59b, 0x3f63fb92
> -       .long   0x3f59ff97, 0x3f4f11d7, 0x3f3d7573, 0x3f24f360
> -       .long   0x3f0cbfe7, 0x3eec1a69, 0x3eb0a801, 0x3e6753a2
> -       .long   0x3e132f1a, 0x3db7e7d3, 0x3d320845, 0x3c84d3d4
> -       .long   0x3bc477b7, 0x3b10d3da, 0x3a01601e, 0x388c1a3b
> -       .long   0x3717b0da, 0x35a43bce, 0x338306c6, 0x00000000
> -       /* p3 */
> +
>         .align  64
> -       .long   0xb0343c7b, 0xbd6ee69d, 0xbd8f0da7, 0xbdae477d
> -       .long   0xbdcd2a1f, 0xbdeba80d, 0xbe0c443b, 0xbe293cf3
> -       .long   0xbe44f282, 0xbe5f3651, 0xbe81c7c0, 0xbe96d7ca
> -       .long   0xbea7fb8e, 0xbeb50e9e, 0xbec12efe, 0xbec4be92
> -       .long   0xbebce070, 0xbead510e, 0xbe8ef7d6, 0xbe4b8704
> -       .long   0xbe083237, 0xbdaf7449, 0xbd2e1ec4, 0xbc83bf06
> -       .long   0xbbc3e0b5, 0xbb10aadc, 0xba0157db, 0xb88c18f2
> -       .long   0xb717b096, 0xb5a43bae, 0xb383012c, 0x00000000
> -       /* p4 */
> +       /* _sP7_lo */
> +       .long   0xbc0e2f66, 0x460bda12, 0x43d638ef, 0xc3e11c3e
> +       .long   0xc2baa4e9, 0xc249da2d, 0xc1859b82, 0x40dd5b57
> +       .long   0x40494640, 0x40c730a8, 0xbf0f160e, 0x3e30e76f
> +       .long   0xbea81387, 0xbdb26a1c, 0xbd351e57, 0xbb4c01a0
> +
>         .align  64
> -       .long   0xbeaaaaa5, 0xbeab0612, 0xbea7f01f, 0xbea4e120
> -       .long   0xbea387b7, 0xbea15962, 0xbe9d57f7, 0xbe976b5a
> -       .long   0xbe90230d, 0xbe880dff, 0xbe7479b3, 0xbe4c3d88
> -       .long   0xbe212482, 0xbdeb8cba, 0xbd5e78ad, 0x3c6b5e6e
> -       .long   0x3d839143, 0x3dc21ee1, 0x3de347af, 0x3dcbec96
> -       .long   0x3d99ef2d, 0x3d542ea1, 0x3cdde701, 0x3c2cca67
> -       .long   0x3b81cb27, 0x3ac073a1, 0x39ac3032, 0x383a94d9
> -       .long   0x36ca081d, 0x355abd4c, 0x332b3cb6, 0x00000000
> -       /* p5 */
> +       /* _sP7_hi */
> +       .long   0x3c1d7bfb, 0x3c722cd1, 0x3c973f1c, 0x3c33a31b
> +       .long   0x3b862ef4, 0x3a27b3d0, 0xba3b5907, 0xba0efc22
> +       .long   0xb97f9f0f, 0xb8c8af50, 0xb7bdddfb, 0xb64f2950
> +       .long   0xb4e085b1, 0xb3731dfa, 0xb15a1f04, 0x00000000
> +
>         .align  64
> -       .long   0xb76dd6b9, 0xbe1c276d, 0x3c1dcf2f, 0x3dc1a78d
> -       .long   0x3d96f985, 0x3da2b61b, 0x3dc13397, 0x3dd2f670
> -       .long   0x3df48a0a, 0x3e06c5a8, 0x3e1a3aba, 0x3e27c405
> -       .long   0x3e2e78d0, 0x3e2c3e44, 0x3e1d3097, 0x3df4a8f4
> -       .long   0x3da38508, 0x3d31416a, 0x3b562657, 0xbcaeeac9
> -       .long   0xbcce9419, 0xbcaaeac4, 0xbc49e7d0, 0xbba71ddd
> -       .long   0xbb003b0e, 0xba3f9a05, 0xb92c08a7, 0xb7ba9232
> -       .long   0xb64a0b0f, 0xb4dac169, 0xb2ab78ac, 0x00000000
> -       /* p6 */
> +       /* _sSignMask */
> +       .long   0x80000000, 0x80000000, 0x80000000, 0x80000000
> +       .long   0x80000000, 0x80000000, 0x80000000, 0x80000000
> +       .long   0x80000000, 0x80000000, 0x80000000, 0x80000000
> +       .long   0x80000000, 0x80000000, 0x80000000, 0x80000000
> +
>         .align  64
> +       /* _sP6_lo */
>         .long   0x3e0910e9, 0x43761143, 0x4165ecdc, 0xc190f756
>         .long   0xc08c097d, 0xc02ba813, 0xbf7f6bda, 0x3f2b1dc0
>         .long   0x3ece105d, 0x3f426a94, 0xbadb0dc4, 0x3da43b17
>         .long   0xbd51ab88, 0xbcaea23d, 0xbd3b6d8d, 0xbd6caaad
> +
> +       .align  64
> +       /* _sP6_hi */
>         .long   0xbd795bed, 0xbd5fddda, 0xbd038f3b, 0xbc1cad63
>         .long   0x3abb4766, 0x3b95f10b, 0x3b825873, 0x3afaea66
>         .long   0x3a49f878, 0x39996bf3, 0x388f3e6c, 0x371bb0e3
>         .long   0x35a8a5e6, 0x34369b17, 0x322487b0, 0x00000000
> -       /* p7 */
> +
>         .align  64
> -       .long   0xbc0e2f66, 0x460bda12, 0x43d638ef, 0xc3e11c3e
> -       .long   0xc2baa4e9, 0xc249da2d, 0xc1859b82, 0x40dd5b57
> -       .long   0x40494640, 0x40c730a8, 0xbf0f160e, 0x3e30e76f
> -       .long   0xbea81387, 0xbdb26a1c, 0xbd351e57, 0xbb4c01a0
> -       .long   0x3c1d7bfb, 0x3c722cd1, 0x3c973f1c, 0x3c33a31b
> -       .long   0x3b862ef4, 0x3a27b3d0, 0xba3b5907, 0xba0efc22
> -       .long   0xb97f9f0f, 0xb8c8af50, 0xb7bdddfb, 0xb64f2950
> -       .long   0xb4e085b1, 0xb3731dfa, 0xb15a1f04, 0x00000000
> +       /* _sP5_lo */
> +       .long   0xb76dd6b9, 0xbe1c276d, 0x3c1dcf2f, 0x3dc1a78d
> +       .long   0x3d96f985, 0x3da2b61b, 0x3dc13397, 0x3dd2f670
> +       .long   0x3df48a0a, 0x3e06c5a8, 0x3e1a3aba, 0x3e27c405
> +       .long   0x3e2e78d0, 0x3e2c3e44, 0x3e1d3097, 0x3df4a8f4
> +
>         .align  64
> -       .long   0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000 /* _iExpMantMask_UISA */
> +       /* _sP5_hi */
> +       .long   0x3da38508, 0x3d31416a, 0x3b562657, 0xbcaeeac9
> +       .long   0xbcce9419, 0xbcaaeac4, 0xbc49e7d0, 0xbba71ddd
> +       .long   0xbb003b0e, 0xba3f9a05, 0xb92c08a7, 0xb7ba9232
> +       .long   0xb64a0b0f, 0xb4dac169, 0xb2ab78ac, 0x00000000
> +
>         .align  64
> -       .long   0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000 /* _iMinIdxOfsMask_UISA */
> +       /* _sP4_lo */
> +       .long   0xbeaaaaa5, 0xbeab0612, 0xbea7f01f, 0xbea4e120
> +       .long   0xbea387b7, 0xbea15962, 0xbe9d57f7, 0xbe976b5a
> +       .long   0xbe90230d, 0xbe880dff, 0xbe7479b3, 0xbe4c3d88
> +       .long   0xbe212482, 0xbdeb8cba, 0xbd5e78ad, 0x3c6b5e6e
> +
>         .align  64
> -       .long   0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000 /* _iMaxIdxMask_UISA */
> +       /* _sP4_hi */
> +       .long   0x3d839143, 0x3dc21ee1, 0x3de347af, 0x3dcbec96
> +       .long   0x3d99ef2d, 0x3d542ea1, 0x3cdde701, 0x3c2cca67
> +       .long   0x3b81cb27, 0x3ac073a1, 0x39ac3032, 0x383a94d9
> +       .long   0x36ca081d, 0x355abd4c, 0x332b3cb6, 0x00000000
> +
>         .align  64
> -       .long   0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000 /* _sSignMask */
> +       /* _sP3_lo */
> +       .long   0xb0343c7b, 0xbd6ee69d, 0xbd8f0da7, 0xbdae477d
> +       .long   0xbdcd2a1f, 0xbdeba80d, 0xbe0c443b, 0xbe293cf3
> +       .long   0xbe44f282, 0xbe5f3651, 0xbe81c7c0, 0xbe96d7ca
> +       .long   0xbea7fb8e, 0xbeb50e9e, 0xbec12efe, 0xbec4be92
> +
>         .align  64
> -       .long   0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff /* _sAbsMask */
> +       /* _sP3_hi */
> +       .long   0xbebce070, 0xbead510e, 0xbe8ef7d6, 0xbe4b8704
> +       .long   0xbe083237, 0xbdaf7449, 0xbd2e1ec4, 0xbc83bf06
> +       .long   0xbbc3e0b5, 0xbb10aadc, 0xba0157db, 0xb88c18f2
> +       .long   0xb717b096, 0xb5a43bae, 0xb383012c, 0x00000000
> +
>         .align  64
> -       .long   0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000 /* _iExpMantMask */
> +       /* _sP2_lo */
> +       .long   0x3f800000, 0x3f7f1f84, 0x3f7ebd11, 0x3f7e1e5f
> +       .long   0x3f7d609f, 0x3f7c842d, 0x3f7b00e5, 0x3f789580
> +       .long   0x3f75b8ad, 0x3f726fd9, 0x3f6cc59b, 0x3f63fb92
> +       .long   0x3f59ff97, 0x3f4f11d7, 0x3f3d7573, 0x3f24f360
> +
>         .align  64
> -       .long   0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000 /* _iExpMask */
> +       /* _sP2_hi */
> +       .long   0x3f0cbfe7, 0x3eec1a69, 0x3eb0a801, 0x3e6753a2
> +       .long   0x3e132f1a, 0x3db7e7d3, 0x3d320845, 0x3c84d3d4
> +       .long   0x3bc477b7, 0x3b10d3da, 0x3a01601e, 0x388c1a3b
> +       .long   0x3717b0da, 0x35a43bce, 0x338306c6, 0x00000000
> +
>         .align  64
> -       .long   0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000 /* _iMinIdxOfsMask */
> +       /* _sP0_lo */
> +       .long   0x00000000, 0x3d6fb9c9, 0x3d8fc35f, 0x3daf9169
> +       .long   0x3dcf49ab, 0x3deee849, 0x3e0f0ee8, 0x3e2e4984
> +       .long   0x3e4d2f8e, 0x3e6bb32e, 0x3e8c51cd, 0x3ea96163
> +       .long   0x3ec543f1, 0x3edfd735, 0x3f028438, 0x3f18abf0
> +
>         .align  64
> -       .long   0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000 /* _iMaxIdxMask */
> +       /* _sP0_hi */
> +       .long   0x3f2bc480, 0x3f3bec1c, 0x3f4f2e5b, 0x3f613c53
> +       .long   0x3f6ce37d, 0x3f743c4f, 0x3f7a5feb, 0x3f7dea85
> +       .long   0x3f7f3b3d, 0x3f7fb78c, 0x3f7fefd4, 0x3f7ffdd0
> +       .long   0x3f7fffb4, 0x3f7ffff6, 0x3f7fffff, 0x3f800000
> +
>         .align  64
> +       .type   __svml_stanh_data_internal_al64, @object
> +       .size   __svml_stanh_data_internal_al64, .-__svml_stanh_data_internal_al64
>         .type   __svml_stanh_data_internal, @object
>         .size   __svml_stanh_data_internal, .-__svml_stanh_data_internal
> --
> 2.34.1
>

LGTM.

Thanks.
diff mbox series

Patch

diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S b/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S
index 5b1f9f151c..7edc74a116 100644
--- a/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_tanhf16_core_avx512.S
@@ -70,310 +70,357 @@ 
  *
  */
 
-/* Offsets for data table __svml_stanh_data_internal
+/* Offsets for data table __svml_stanh_data_internal_avx512. Ordered
+   by use in the function. On cold-starts this might help the
+   prefetcher. Possibly a better idea is to interleave start/end so
+   that the prefetcher is less likely to detect a stream and pull
+   irrelivant lines into cache.  */
+
+/* Offsets for data table __svml_stanh_data_internal. 4 bytes each.
  */
-#define _sC				0
-#define _sP0				128
-#define _sP2				256
-#define _sP3				384
-#define _sP4				512
-#define _sP5				640
-#define _sP6				768
-#define _sP7				896
-#define _iExpMantMask_UISA		1024
-#define _iMinIdxOfsMask_UISA		1088
-#define _iMaxIdxMask_UISA		1152
-#define _sSignMask			1216
-#define _sAbsMask			1280
-#define _iExpMantMask			1344
-#define _iExpMask			1408
-#define _iMinIdxOfsMask			1472
-#define _iMaxIdxMask			1536
+#define _iExpMantMask_UISA		0
+#define _iMinIdxOfsMask_UISA		4
+#define _iMaxIdxMask_UISA		8
+#define _iExpMask			12
+
+/* Offsets for data table __svml_stanh_data_internal_al64. 64 bytes
+   each.  */
+#define _sC_lo				0
+#define _sC_hi				64
+#define _sP7_lo				128
+#define _sP7_hi				192
+#define _sSignMask			256
+#define _sP6_lo				320
+#define _sP6_hi				384
+#define _sP5_lo				448
+#define _sP5_hi				512
+#define _sP4_lo				576
+#define _sP4_hi				640
+#define _sP3_lo				704
+#define _sP3_hi				768
+#define _sP2_lo				832
+#define _sP2_hi				896
+#define _sP0_lo				960
+#define _sP0_hi				1024
 
 #include <sysdep.h>
+#define TANHF_DATA(x)			((x)+__svml_stanh_data_internal_al64)
+#define TANHF_DATA_UNALIGNED(x)		((x)+__svml_stanh_data_internal)
 
 	.section .text.exex512, "ax", @progbits
 ENTRY(_ZGVeN16v_tanhf_skx)
-	pushq	%rbp
-	cfi_def_cfa_offset(16)
-	movq	%rsp, %rbp
-	cfi_def_cfa(6, 16)
-	cfi_offset(6, -16)
-	andq	$-64, %rsp
-	subq	$192, %rsp
-	vmovaps	%zmm0, %zmm1
-	vmovups	__svml_stanh_data_internal(%rip), %zmm9
-	vmovups	_sP6+__svml_stanh_data_internal(%rip), %zmm11
-	vmovups	_sP5+__svml_stanh_data_internal(%rip), %zmm12
-	vmovups	_sP4+__svml_stanh_data_internal(%rip), %zmm13
-	vmovups	_sP3+__svml_stanh_data_internal(%rip), %zmm14
-	vmovups	_sP2+__svml_stanh_data_internal(%rip), %zmm15
-	vpternlogd $255, %zmm2, %zmm2, %zmm2
-	vandps	_sAbsMask+__svml_stanh_data_internal(%rip), %zmm1, %zmm8
-	vandps	_sSignMask+__svml_stanh_data_internal(%rip), %zmm1, %zmm0
-
 	/* Here huge arguments, INF and NaNs are filtered out to callout. */
-	vpandd	_iExpMantMask_UISA+__svml_stanh_data_internal(%rip), %zmm1, %zmm3
-	vpsubd	_iMinIdxOfsMask_UISA+__svml_stanh_data_internal(%rip), %zmm3, %zmm4
-	vpcmpd	$2, _iExpMask+__svml_stanh_data_internal(%rip), %zmm3, %k1
+	vpandd	TANHF_DATA_UNALIGNED(_iExpMantMask_UISA)(%rip){1to16}, %zmm0, %zmm1
+	vpsubd	TANHF_DATA_UNALIGNED(_iMinIdxOfsMask_UISA)(%rip){1to16}, %zmm1, %zmm2
 
-	/*
-	 *  small table specific variables *
-	 *  Constant loading
-	 */
-	vpxord	%zmm5, %zmm5, %zmm5
-
-	/* if VMIN, VMAX is defined for I type */
-	vpmaxsd	%zmm5, %zmm4, %zmm6
-	vpminsd	_iMaxIdxMask_UISA+__svml_stanh_data_internal(%rip), %zmm6, %zmm7
-	vpsrld	$21, %zmm7, %zmm10
-	vmovups	_sP7+__svml_stanh_data_internal(%rip), %zmm4
-	vpermt2ps _sC+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm9
-	vpermt2ps _sP6+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm11
-	vpermt2ps _sP7+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm4
-	vpermt2ps _sP5+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm12
-	vpermt2ps _sP4+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm13
-	vpermt2ps _sP3+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm14
-	vpermt2ps _sP2+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm15
-	vpandnd	%zmm3, %zmm3, %zmm2{%k1}
-	vptestmd %zmm2, %zmm2, %k0
-	vmovups	_sP0+__svml_stanh_data_internal(%rip), %zmm3
-	vsubps	{rn-sae}, %zmm9, %zmm8, %zmm2
-	kmovw	%k0, %edx
-	vfmadd213ps {rn-sae}, %zmm11, %zmm2, %zmm4
-	vpermt2ps _sP0+64+__svml_stanh_data_internal(%rip), %zmm10, %zmm3
-	vfmadd213ps {rn-sae}, %zmm12, %zmm2, %zmm4
-	vfmadd213ps {rn-sae}, %zmm13, %zmm2, %zmm4
-	vfmadd213ps {rn-sae}, %zmm14, %zmm2, %zmm4
-	vfmadd213ps {rn-sae}, %zmm15, %zmm2, %zmm4
-	vfmadd213ps {rn-sae}, %zmm3, %zmm2, %zmm4
-	vorps	%zmm0, %zmm4, %zmm0
-	testl	%edx, %edx
+	/* Selection arguments between [0, 0x03e00000] into zmm3.  */
+	vpxord	%zmm3, %zmm3, %zmm3
+	vpmaxsd	%zmm3, %zmm2, %zmm3
+	vpminsd	TANHF_DATA_UNALIGNED(_iMaxIdxMask_UISA)(%rip){1to16}, %zmm3, %zmm3
 
-	/* Go to special inputs processing branch */
-	jne	L(SPECIAL_VALUES_BRANCH)
-	# LOE rbx r12 r13 r14 r15 edx zmm0 zmm1
+	/* Setup permute indices in zmm3.  */
+	vpsrld	$21, %zmm3, %zmm3
 
-	/* Restore registers
-	 * and exit the function
-	 */
+	/* Store if there are any special cases in k1.  */
+	vpcmpd	$6, TANHF_DATA_UNALIGNED(_iExpMask)(%rip){1to16}, %zmm1, %k1
 
-L(EXIT):
-	movq	%rbp, %rsp
-	popq	%rbp
-	cfi_def_cfa(7, 8)
-	cfi_restore(6)
-	ret
-	cfi_def_cfa(6, 16)
-	cfi_offset(6, -16)
+	vmovaps	TANHF_DATA(_sC_lo)(%rip), %zmm5
+	vpermt2ps TANHF_DATA(_sC_hi)(%rip), %zmm3, %zmm5
 
-	/* Branch to process
-	 * special inputs
-	 */
+	vmovaps	TANHF_DATA(_sP7_lo)(%rip), %zmm2
+	vpermt2ps TANHF_DATA(_sP7_hi)(%rip), %zmm3, %zmm2
 
-L(SPECIAL_VALUES_BRANCH):
-	vmovups	%zmm1, 64(%rsp)
-	vmovups	%zmm0, 128(%rsp)
-	# LOE rbx r12 r13 r14 r15 edx zmm0
+	/* Store absolute values of inputs in zmm1.  */
+	vmovaps	TANHF_DATA(_sSignMask)(%rip), %zmm4
+	vandnps	%zmm0, %zmm4, %zmm1
+	vsubps	{rn-sae}, %zmm5, %zmm1, %zmm1
 
-	xorl	%eax, %eax
-	# LOE rbx r12 r13 r14 r15 eax edx
+	vmovaps	TANHF_DATA(_sP6_lo)(%rip), %zmm5
+	vpermt2ps TANHF_DATA(_sP6_hi)(%rip), %zmm3, %zmm5
 
-	vzeroupper
-	movq	%r12, 16(%rsp)
-	/*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -176; DW_OP_plus)  */
-	.cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x50, 0xff, 0xff, 0xff, 0x22
-	movl	%eax, %r12d
-	movq	%r13, 8(%rsp)
-	/*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -184; DW_OP_plus)  */
-	.cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x48, 0xff, 0xff, 0xff, 0x22
-	movl	%edx, %r13d
-	movq	%r14, (%rsp)
-	/*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -192; DW_OP_plus)  */
-	.cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x40, 0xff, 0xff, 0xff, 0x22
-	# LOE rbx r15 r12d r13d
-
-	/* Range mask
-	 * bits check
-	 */
+	vmovaps	TANHF_DATA(_sP5_lo)(%rip), %zmm6
+	vpermt2ps TANHF_DATA(_sP5_hi)(%rip), %zmm3, %zmm6
+
+	vfmadd213ps {rn-sae}, %zmm5, %zmm1, %zmm2
+	vfmadd213ps {rn-sae}, %zmm6, %zmm1, %zmm2
 
-L(RANGEMASK_CHECK):
-	btl	%r12d, %r13d
+	vmovaps	TANHF_DATA(_sP4_lo)(%rip), %zmm7
+	vpermt2ps TANHF_DATA(_sP4_hi)(%rip), %zmm3, %zmm7
 
-	/* Call scalar math function */
-	jc	L(SCALAR_MATH_CALL)
-	# LOE rbx r15 r12d r13d
+	vmovaps	TANHF_DATA(_sP3_lo)(%rip), %zmm8
+	vpermt2ps TANHF_DATA(_sP3_hi)(%rip), %zmm3, %zmm8
 
-	/* Special inputs
-	 * processing loop
+	vfmadd213ps {rn-sae}, %zmm7, %zmm1, %zmm2
+	vfmadd213ps {rn-sae}, %zmm8, %zmm1, %zmm2
+
+	vmovaps	TANHF_DATA(_sP2_lo)(%rip), %zmm9
+	vpermt2ps TANHF_DATA(_sP2_hi)(%rip), %zmm3, %zmm9
+
+	vmovaps	TANHF_DATA(_sP0_lo)(%rip), %zmm10
+	vpermt2ps TANHF_DATA(_sP0_hi)(%rip), %zmm3, %zmm10
+
+	vfmadd213ps {rn-sae}, %zmm9, %zmm1, %zmm2
+	vfmadd213ps {rn-sae}, %zmm10, %zmm1, %zmm2
+
+	kmovw	%k1, %edx
+	testl	%edx, %edx
+
+	/* Go to special inputs processing branch.  */
+	jne	L(SPECIAL_VALUES_BRANCH)
+	# LOE rbx r12 r13 r14 r15 zmm0 zmm2 zmm4
+	/* Wait until after branch of write over zmm0.  */
+	vpternlogd $0xec, %zmm4, %zmm2, %zmm0
+
+	/* No stack restoration on the fastpath.  */
+	ret
+
+	/* Cold case. edx has 1s where there was a special value that
+	   needs to be handled by a tanhf call. Optimize for code size
+	   more so than speed here. */
+L(SPECIAL_VALUES_BRANCH):
+	# LOE rbx rdx r12 r13 r14 r15 zmm0 zmm2 zmm4
+    /* Use r13 to save/restore the stack. This allows us to use rbp as
+       callee save register saving code size. */
+	pushq	%r13
+	cfi_adjust_cfa_offset(8)
+	cfi_offset(r13, -16)
+	/* Need to callee save registers to preserve state across tanhf calls.
 	 */
+	pushq	%rbx
+	cfi_adjust_cfa_offset(8)
+	cfi_offset(rbx, -24)
+	pushq	%rbp
+	cfi_adjust_cfa_offset(8)
+	cfi_offset(rbp, -32)
+	movq	%rsp, %r13
+	cfi_def_cfa_register(r13)
+
+	/* Align stack and make room for 2x zmm vectors.  */
+	andq	$-64, %rsp
+	addq	$-128, %rsp
+
+	/* Save original input (zmm0 unchanged up to this point).  */
+	vmovaps	%zmm0, 64(%rsp)
+	/* Save all already computed inputs.  */
+	vpternlogd $0xec, %zmm4, %zmm2, %zmm0
+	vmovaps	%zmm0, (%rsp)
 
+	vzeroupper
+
+	/* edx has 1s where there was a special value that needs to be handled
+	   by a tanhf call.  */
+	movl	%edx, %ebx
 L(SPECIAL_VALUES_LOOP):
-	incl	%r12d
-	cmpl	$16, %r12d
-
-	/* Check bits in range mask */
-	jl	L(RANGEMASK_CHECK)
-	# LOE rbx r15 r12d r13d
-
-	movq	16(%rsp), %r12
-	cfi_restore(12)
-	movq	8(%rsp), %r13
-	cfi_restore(13)
-	movq	(%rsp), %r14
-	cfi_restore(14)
-	vmovups	128(%rsp), %zmm0
-
-	/* Go to exit */
-	jmp	L(EXIT)
-	/*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -176; DW_OP_plus)  */
-	.cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x50, 0xff, 0xff, 0xff, 0x22
-	/*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -184; DW_OP_plus)  */
-	.cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x48, 0xff, 0xff, 0xff, 0x22
-	/*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -64; DW_OP_and; DW_OP_const4s: -192; DW_OP_plus)  */
-	.cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0x40, 0xff, 0xff, 0xff, 0x22
-	# LOE rbx r12 r13 r14 r15 zmm0
-
-	/* Scalar math fucntion call
-	 * to process special input
-	 */
+	# LOE rbx rbp r12 r13 r14 r15
+	/* use rbp as index for special value that is saved across calls to
+	   tanhf. We technically don't need a callee save register here as offset
+	   to rsp is always [0, 56] so we can restore rsp by realigning to 64.
+	   Essentially the tradeoff is 1 extra save/restore vs 2 extra instructions
+	   in the loop. Realigning also costs more code size.  */
+	xorl	%ebp, %ebp
+	tzcntl	%ebx, %ebp
 
-L(SCALAR_MATH_CALL):
-	movl	%r12d, %r14d
-	movss	64(%rsp, %r14, 4), %xmm0
+	/* Scalar math fucntion call to process special input.  */
+	movss	64(%rsp, %rbp, 4), %xmm0
 	call	tanhf@PLT
-	# LOE rbx r14 r15 r12d r13d xmm0
 
-	movss	%xmm0, 128(%rsp, %r14, 4)
+	/* No good way to avoid the store-forwarding fault this will cause on
+	   return. `lfence` avoids the SF fault but at greater cost as it
+	   serialized stack/callee save restoration.  */
+	movss	%xmm0, (%rsp, %rbp, 4)
 
-	/* Process special inputs in loop */
-	jmp	L(SPECIAL_VALUES_LOOP)
-	# LOE rbx r15 r12d r13d
+	blsrl   %ebx, %ebx
+	jnz	L(SPECIAL_VALUES_LOOP)
+	# LOE r12 r13 r14 r15
+
+	/* All results have been written to (%rsp).  */
+	vmovaps	(%rsp), %zmm0
+	/* Restore rsp.  */
+	movq	%r13, %rsp
+	cfi_def_cfa_register(rsp)
+	/* Restore callee save registers.  */
+	popq	%rbp
+	cfi_adjust_cfa_offset(-8)
+	cfi_restore(rbp)
+	popq	%rbx
+	cfi_adjust_cfa_offset(-8)
+	cfi_restore(rbp)
+	popq	%r13
+	cfi_adjust_cfa_offset(-8)
+	cfi_restore(r13)
+	ret
 END(_ZGVeN16v_tanhf_skx)
 
 	.section .rodata, "a"
-	.align	64
-
+	.align	16
 #ifdef __svml_stanh_data_internal_typedef
 typedef unsigned int VUINT32;
-typedef struct {
-	__declspec(align(64)) VUINT32 _sC[32][1];
-	__declspec(align(64)) VUINT32 _sP0[32][1];
-	__declspec(align(64)) VUINT32 _sP2[32][1];
-	__declspec(align(64)) VUINT32 _sP3[32][1];
-	__declspec(align(64)) VUINT32 _sP4[32][1];
-	__declspec(align(64)) VUINT32 _sP5[32][1];
-	__declspec(align(64)) VUINT32 _sP6[32][1];
-	__declspec(align(64)) VUINT32 _sP7[32][1];
-	__declspec(align(64)) VUINT32 _iExpMantMask_UISA[16][1];
-	__declspec(align(64)) VUINT32 _iMinIdxOfsMask_UISA[16][1];
-	__declspec(align(64)) VUINT32 _iMaxIdxMask_UISA[16][1];
+typedef struct
+	{
+	__declspec(align(4)) VUINT32 _iExpMantMask_UISA[1][1];
+	__declspec(align(4)) VUINT32 _iMinIdxOfsMask_UISA[1][1];
+	__declspec(align(4)) VUINT32 _iMaxIdxMask_UISA[1][1];
+	__declspec(align(4)) VUINT32 _iExpMask[1][1];
+	__declspec(align(64)) VUINT32 _sC_lo[16][1];
+	__declspec(align(64)) VUINT32 _sC_hi[16][1];
+	__declspec(align(64)) VUINT32 _sP7_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP7_hi[16][1];
 	__declspec(align(64)) VUINT32 _sSignMask[16][1];
-	__declspec(align(64)) VUINT32 _sAbsMask[16][1];
-	__declspec(align(64)) VUINT32 _iExpMantMask[16][1];
-	__declspec(align(64)) VUINT32 _iExpMask[16][1];
-	__declspec(align(64)) VUINT32 _iMinIdxOfsMask[16][1];
-	__declspec(align(64)) VUINT32 _iMaxIdxMask[16][1];
+	__declspec(align(64)) VUINT32 _sP6_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP6_hi[16][1];
+	__declspec(align(64)) VUINT32 _sP5_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP5_hi[16][1];
+	__declspec(align(64)) VUINT32 _sP4_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP4_hi[16][1];
+	__declspec(align(64)) VUINT32 _sP3_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP3_hi[16][1];
+	__declspec(align(64)) VUINT32 _sP2_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP2_hi[16][1];
+	__declspec(align(64)) VUINT32 _sP0_lo[16][1];
+	__declspec(align(64)) VUINT32 _sP0_hi[16][1];
 } __svml_stanh_data_internal;
 #endif
+
 __svml_stanh_data_internal:
-	/* _sC */
+	.align	4
+	/* _iExpMantMask_UISA */
+	.long	0x7fe00000
+
+	.align	4
+	/* _iMinIdxOfsMask_UISA */
+	.long	0x3d400000
+
+	.align	4
+	/* _iMaxIdxMask_UISA */
+	.long	0x03e00000
+
+	.align	4
+	/* _iExpMask */
+	.long	0x7f000000
+
+	.align	64
+__svml_stanh_data_internal_al64:
+	.align	64
+	/* _sC_lo */
 	.long	0x00000000, 0x3d700000, 0x3d900000, 0x3db00000
 	.long	0x3dd00000, 0x3df00000, 0x3e100000, 0x3e300000
 	.long	0x3e500000, 0x3e700000, 0x3e900000, 0x3eb00000
 	.long	0x3ed00000, 0x3ef00000, 0x3f100000, 0x3f300000
+
+	.align	64
+	/* _sC_hi */
 	.long	0x3f500000, 0x3f700000, 0x3f900000, 0x3fb00000
 	.long	0x3fd00000, 0x3ff00000, 0x40100000, 0x40300000
 	.long	0x40500000, 0x40700000, 0x40900000, 0x40b00000
 	.long	0x40d00000, 0x40f00000, 0x41100000, 0x00000000
-	/* p0 */
-	.align	64
-	.long	0x00000000, 0x3d6fb9c9, 0x3d8fc35f, 0x3daf9169
-	.long	0x3dcf49ab, 0x3deee849, 0x3e0f0ee8, 0x3e2e4984
-	.long	0x3e4d2f8e, 0x3e6bb32e, 0x3e8c51cd, 0x3ea96163
-	.long	0x3ec543f1, 0x3edfd735, 0x3f028438, 0x3f18abf0
-	.long	0x3f2bc480, 0x3f3bec1c, 0x3f4f2e5b, 0x3f613c53
-	.long	0x3f6ce37d, 0x3f743c4f, 0x3f7a5feb, 0x3f7dea85
-	.long	0x3f7f3b3d, 0x3f7fb78c, 0x3f7fefd4, 0x3f7ffdd0
-	.long	0x3f7fffb4, 0x3f7ffff6, 0x3f7fffff, 0x3f800000
-	/* p2 */
-	.align	64
-	.long	0x3f800000, 0x3f7f1f84, 0x3f7ebd11, 0x3f7e1e5f
-	.long	0x3f7d609f, 0x3f7c842d, 0x3f7b00e5, 0x3f789580
-	.long	0x3f75b8ad, 0x3f726fd9, 0x3f6cc59b, 0x3f63fb92
-	.long	0x3f59ff97, 0x3f4f11d7, 0x3f3d7573, 0x3f24f360
-	.long	0x3f0cbfe7, 0x3eec1a69, 0x3eb0a801, 0x3e6753a2
-	.long	0x3e132f1a, 0x3db7e7d3, 0x3d320845, 0x3c84d3d4
-	.long	0x3bc477b7, 0x3b10d3da, 0x3a01601e, 0x388c1a3b
-	.long	0x3717b0da, 0x35a43bce, 0x338306c6, 0x00000000
-	/* p3 */
+
 	.align	64
-	.long	0xb0343c7b, 0xbd6ee69d, 0xbd8f0da7, 0xbdae477d
-	.long	0xbdcd2a1f, 0xbdeba80d, 0xbe0c443b, 0xbe293cf3
-	.long	0xbe44f282, 0xbe5f3651, 0xbe81c7c0, 0xbe96d7ca
-	.long	0xbea7fb8e, 0xbeb50e9e, 0xbec12efe, 0xbec4be92
-	.long	0xbebce070, 0xbead510e, 0xbe8ef7d6, 0xbe4b8704
-	.long	0xbe083237, 0xbdaf7449, 0xbd2e1ec4, 0xbc83bf06
-	.long	0xbbc3e0b5, 0xbb10aadc, 0xba0157db, 0xb88c18f2
-	.long	0xb717b096, 0xb5a43bae, 0xb383012c, 0x00000000
-	/* p4 */
+	/* _sP7_lo */
+	.long	0xbc0e2f66, 0x460bda12, 0x43d638ef, 0xc3e11c3e
+	.long	0xc2baa4e9, 0xc249da2d, 0xc1859b82, 0x40dd5b57
+	.long	0x40494640, 0x40c730a8, 0xbf0f160e, 0x3e30e76f
+	.long	0xbea81387, 0xbdb26a1c, 0xbd351e57, 0xbb4c01a0
+
 	.align	64
-	.long	0xbeaaaaa5, 0xbeab0612, 0xbea7f01f, 0xbea4e120
-	.long	0xbea387b7, 0xbea15962, 0xbe9d57f7, 0xbe976b5a
-	.long	0xbe90230d, 0xbe880dff, 0xbe7479b3, 0xbe4c3d88
-	.long	0xbe212482, 0xbdeb8cba, 0xbd5e78ad, 0x3c6b5e6e
-	.long	0x3d839143, 0x3dc21ee1, 0x3de347af, 0x3dcbec96
-	.long	0x3d99ef2d, 0x3d542ea1, 0x3cdde701, 0x3c2cca67
-	.long	0x3b81cb27, 0x3ac073a1, 0x39ac3032, 0x383a94d9
-	.long	0x36ca081d, 0x355abd4c, 0x332b3cb6, 0x00000000
-	/* p5 */
+	/* _sP7_hi */
+	.long	0x3c1d7bfb, 0x3c722cd1, 0x3c973f1c, 0x3c33a31b
+	.long	0x3b862ef4, 0x3a27b3d0, 0xba3b5907, 0xba0efc22
+	.long	0xb97f9f0f, 0xb8c8af50, 0xb7bdddfb, 0xb64f2950
+	.long	0xb4e085b1, 0xb3731dfa, 0xb15a1f04, 0x00000000
+
 	.align	64
-	.long	0xb76dd6b9, 0xbe1c276d, 0x3c1dcf2f, 0x3dc1a78d
-	.long	0x3d96f985, 0x3da2b61b, 0x3dc13397, 0x3dd2f670
-	.long	0x3df48a0a, 0x3e06c5a8, 0x3e1a3aba, 0x3e27c405
-	.long	0x3e2e78d0, 0x3e2c3e44, 0x3e1d3097, 0x3df4a8f4
-	.long	0x3da38508, 0x3d31416a, 0x3b562657, 0xbcaeeac9
-	.long	0xbcce9419, 0xbcaaeac4, 0xbc49e7d0, 0xbba71ddd
-	.long	0xbb003b0e, 0xba3f9a05, 0xb92c08a7, 0xb7ba9232
-	.long	0xb64a0b0f, 0xb4dac169, 0xb2ab78ac, 0x00000000
-	/* p6 */
+	/* _sSignMask */
+	.long	0x80000000, 0x80000000, 0x80000000, 0x80000000
+	.long	0x80000000, 0x80000000, 0x80000000, 0x80000000
+	.long	0x80000000, 0x80000000, 0x80000000, 0x80000000
+	.long	0x80000000, 0x80000000, 0x80000000, 0x80000000
+
 	.align	64
+	/* _sP6_lo */
 	.long	0x3e0910e9, 0x43761143, 0x4165ecdc, 0xc190f756
 	.long	0xc08c097d, 0xc02ba813, 0xbf7f6bda, 0x3f2b1dc0
 	.long	0x3ece105d, 0x3f426a94, 0xbadb0dc4, 0x3da43b17
 	.long	0xbd51ab88, 0xbcaea23d, 0xbd3b6d8d, 0xbd6caaad
+
+	.align	64
+	/* _sP6_hi */
 	.long	0xbd795bed, 0xbd5fddda, 0xbd038f3b, 0xbc1cad63
 	.long	0x3abb4766, 0x3b95f10b, 0x3b825873, 0x3afaea66
 	.long	0x3a49f878, 0x39996bf3, 0x388f3e6c, 0x371bb0e3
 	.long	0x35a8a5e6, 0x34369b17, 0x322487b0, 0x00000000
-	/* p7 */
+
 	.align	64
-	.long	0xbc0e2f66, 0x460bda12, 0x43d638ef, 0xc3e11c3e
-	.long	0xc2baa4e9, 0xc249da2d, 0xc1859b82, 0x40dd5b57
-	.long	0x40494640, 0x40c730a8, 0xbf0f160e, 0x3e30e76f
-	.long	0xbea81387, 0xbdb26a1c, 0xbd351e57, 0xbb4c01a0
-	.long	0x3c1d7bfb, 0x3c722cd1, 0x3c973f1c, 0x3c33a31b
-	.long	0x3b862ef4, 0x3a27b3d0, 0xba3b5907, 0xba0efc22
-	.long	0xb97f9f0f, 0xb8c8af50, 0xb7bdddfb, 0xb64f2950
-	.long	0xb4e085b1, 0xb3731dfa, 0xb15a1f04, 0x00000000
+	/* _sP5_lo */
+	.long	0xb76dd6b9, 0xbe1c276d, 0x3c1dcf2f, 0x3dc1a78d
+	.long	0x3d96f985, 0x3da2b61b, 0x3dc13397, 0x3dd2f670
+	.long	0x3df48a0a, 0x3e06c5a8, 0x3e1a3aba, 0x3e27c405
+	.long	0x3e2e78d0, 0x3e2c3e44, 0x3e1d3097, 0x3df4a8f4
+
 	.align	64
-	.long	0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000, 0x7fe00000 /* _iExpMantMask_UISA */
+	/* _sP5_hi */
+	.long	0x3da38508, 0x3d31416a, 0x3b562657, 0xbcaeeac9
+	.long	0xbcce9419, 0xbcaaeac4, 0xbc49e7d0, 0xbba71ddd
+	.long	0xbb003b0e, 0xba3f9a05, 0xb92c08a7, 0xb7ba9232
+	.long	0xb64a0b0f, 0xb4dac169, 0xb2ab78ac, 0x00000000
+
 	.align	64
-	.long	0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000, 0x3d400000 /* _iMinIdxOfsMask_UISA */
+	/* _sP4_lo */
+	.long	0xbeaaaaa5, 0xbeab0612, 0xbea7f01f, 0xbea4e120
+	.long	0xbea387b7, 0xbea15962, 0xbe9d57f7, 0xbe976b5a
+	.long	0xbe90230d, 0xbe880dff, 0xbe7479b3, 0xbe4c3d88
+	.long	0xbe212482, 0xbdeb8cba, 0xbd5e78ad, 0x3c6b5e6e
+
 	.align	64
-	.long	0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000, 0x03e00000 /* _iMaxIdxMask_UISA */
+	/* _sP4_hi */
+	.long	0x3d839143, 0x3dc21ee1, 0x3de347af, 0x3dcbec96
+	.long	0x3d99ef2d, 0x3d542ea1, 0x3cdde701, 0x3c2cca67
+	.long	0x3b81cb27, 0x3ac073a1, 0x39ac3032, 0x383a94d9
+	.long	0x36ca081d, 0x355abd4c, 0x332b3cb6, 0x00000000
+
 	.align	64
-	.long	0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000, 0x80000000 /* _sSignMask */
+	/* _sP3_lo */
+	.long	0xb0343c7b, 0xbd6ee69d, 0xbd8f0da7, 0xbdae477d
+	.long	0xbdcd2a1f, 0xbdeba80d, 0xbe0c443b, 0xbe293cf3
+	.long	0xbe44f282, 0xbe5f3651, 0xbe81c7c0, 0xbe96d7ca
+	.long	0xbea7fb8e, 0xbeb50e9e, 0xbec12efe, 0xbec4be92
+
 	.align	64
-	.long	0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff, 0x7fffffff /* _sAbsMask */
+	/* _sP3_hi */
+	.long	0xbebce070, 0xbead510e, 0xbe8ef7d6, 0xbe4b8704
+	.long	0xbe083237, 0xbdaf7449, 0xbd2e1ec4, 0xbc83bf06
+	.long	0xbbc3e0b5, 0xbb10aadc, 0xba0157db, 0xb88c18f2
+	.long	0xb717b096, 0xb5a43bae, 0xb383012c, 0x00000000
+
 	.align	64
-	.long	0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000, 0x7ff80000 /* _iExpMantMask */
+	/* _sP2_lo */
+	.long	0x3f800000, 0x3f7f1f84, 0x3f7ebd11, 0x3f7e1e5f
+	.long	0x3f7d609f, 0x3f7c842d, 0x3f7b00e5, 0x3f789580
+	.long	0x3f75b8ad, 0x3f726fd9, 0x3f6cc59b, 0x3f63fb92
+	.long	0x3f59ff97, 0x3f4f11d7, 0x3f3d7573, 0x3f24f360
+
 	.align	64
-	.long	0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000, 0x7f000000 /* _iExpMask */
+	/* _sP2_hi */
+	.long	0x3f0cbfe7, 0x3eec1a69, 0x3eb0a801, 0x3e6753a2
+	.long	0x3e132f1a, 0x3db7e7d3, 0x3d320845, 0x3c84d3d4
+	.long	0x3bc477b7, 0x3b10d3da, 0x3a01601e, 0x388c1a3b
+	.long	0x3717b0da, 0x35a43bce, 0x338306c6, 0x00000000
+
 	.align	64
-	.long	0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000, 0x3cf80000 /* _iMinIdxOfsMask */
+	/* _sP0_lo */
+	.long	0x00000000, 0x3d6fb9c9, 0x3d8fc35f, 0x3daf9169
+	.long	0x3dcf49ab, 0x3deee849, 0x3e0f0ee8, 0x3e2e4984
+	.long	0x3e4d2f8e, 0x3e6bb32e, 0x3e8c51cd, 0x3ea96163
+	.long	0x3ec543f1, 0x3edfd735, 0x3f028438, 0x3f18abf0
+
 	.align	64
-	.long	0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000, 0x04280000 /* _iMaxIdxMask */
+	/* _sP0_hi */
+	.long	0x3f2bc480, 0x3f3bec1c, 0x3f4f2e5b, 0x3f613c53
+	.long	0x3f6ce37d, 0x3f743c4f, 0x3f7a5feb, 0x3f7dea85
+	.long	0x3f7f3b3d, 0x3f7fb78c, 0x3f7fefd4, 0x3f7ffdd0
+	.long	0x3f7fffb4, 0x3f7ffff6, 0x3f7fffff, 0x3f800000
+
 	.align	64
+	.type	__svml_stanh_data_internal_al64, @object
+	.size	__svml_stanh_data_internal_al64, .-__svml_stanh_data_internal_al64
 	.type	__svml_stanh_data_internal, @object
 	.size	__svml_stanh_data_internal, .-__svml_stanh_data_internal