aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorSunil K Pandey <skpgkp2@gmail.com>2022-03-07 10:47:10 -0800
committerSunil K Pandey <skpgkp2@gmail.com>2022-03-07 21:14:10 -0800
commit1447e84caf034d6547ca52ad7ef195a1681f3730 (patch)
tree02cf8788ef942eb1797b68921bcdf40781e439ea
parentf55b59764a360162eb5fb9cfad8a00864df14ed5 (diff)
downloadglibc-1447e84caf034d6547ca52ad7ef195a1681f3730.zip
glibc-1447e84caf034d6547ca52ad7ef195a1681f3730.tar.gz
glibc-1447e84caf034d6547ca52ad7ef195a1681f3730.tar.bz2
x86_64: Fix svml_d_atan22_core_sse4.S code formatting
This commit contains following formatting changes 1. Instructions proceeded by a tab. 2. Instruction less than 8 characters in length have a tab between it and the first operand. 3. Instruction greater than 7 characters in length have a space between it and the first operand. 4. Tabs after `#define`d names and their value. 5. 8 space at the beginning of line replaced by tab. 6. Indent comments with code. 7. Remove redundent .text section. 8. 1 space between line content and line comment. 9. Space after all commas. Reviewed-by: Noah Goldstein <goldstein.w.n@gmail.com>
-rw-r--r--sysdeps/x86_64/fpu/multiarch/svml_d_atan22_core_sse4.S861
1 files changed, 430 insertions, 431 deletions
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_atan22_core_sse4.S b/sysdeps/x86_64/fpu/multiarch/svml_d_atan22_core_sse4.S
index 138ff2f..e1298b5 100644
--- a/sysdeps/x86_64/fpu/multiarch/svml_d_atan22_core_sse4.S
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_atan22_core_sse4.S
@@ -30,467 +30,466 @@
/* Offsets for data table __svml_datan2_data_internal
*/
-#define dPI 0
-#define dPIO2 16
-#define dA19 32
-#define dA18 48
-#define dA17 64
-#define dA16 80
-#define dA15 96
-#define dA14 112
-#define dA13 128
-#define dA12 144
-#define dA11 160
-#define dA10 176
-#define dA09 192
-#define dA08 208
-#define dA07 224
-#define dA06 240
-#define dA05 256
-#define dA04 272
-#define dA03 288
-#define dA02 304
-#define dA01 320
-#define dA00 336
-#define dSIGN_MASK 352
-#define iCHK_WORK_SUB 368
-#define iCHK_WORK_CMP 384
-#define dABS_MASK 400
-#define dZERO 416
+#define dPI 0
+#define dPIO2 16
+#define dA19 32
+#define dA18 48
+#define dA17 64
+#define dA16 80
+#define dA15 96
+#define dA14 112
+#define dA13 128
+#define dA12 144
+#define dA11 160
+#define dA10 176
+#define dA09 192
+#define dA08 208
+#define dA07 224
+#define dA06 240
+#define dA05 256
+#define dA04 272
+#define dA03 288
+#define dA02 304
+#define dA01 320
+#define dA00 336
+#define dSIGN_MASK 352
+#define iCHK_WORK_SUB 368
+#define iCHK_WORK_CMP 384
+#define dABS_MASK 400
+#define dZERO 416
#include <sysdep.h>
- .text
- .section .text.sse4,"ax",@progbits
+ .section .text.sse4, "ax", @progbits
ENTRY(_ZGVbN2vv_atan2_sse4)
- subq $88, %rsp
- cfi_def_cfa_offset(96)
- movaps %xmm1, %xmm11
-
-/*
- * #define NO_VECTOR_ZERO_ATAN2_ARGS
- * Declarations
- * Variables
- * Constants
- * The end of declarations
- * Implementation
- * Get r0~=1/B
- * Cannot be replaced by VQRCP(D, dR0, dB);
- * Argument Absolute values
- */
- movups dABS_MASK+__svml_datan2_data_internal(%rip), %xmm1
- movaps %xmm0, %xmm10
- movaps %xmm1, %xmm9
- andps %xmm10, %xmm1
- andps %xmm11, %xmm9
- movaps %xmm1, %xmm4
- cmpnltpd %xmm9, %xmm4
-
-/* Argument signs */
- movups dSIGN_MASK+__svml_datan2_data_internal(%rip), %xmm5
- movaps %xmm4, %xmm0
- movaps %xmm5, %xmm8
- movaps %xmm5, %xmm7
-
-/*
- * 1) If y<x then a= y, b=x, PIO2=0
- * 2) If y>x then a=-x, b=y, PIO2=Pi/2
- */
- orps %xmm9, %xmm5
- andnps %xmm1, %xmm0
- andps %xmm4, %xmm5
- andps %xmm11, %xmm8
- movups dPIO2+__svml_datan2_data_internal(%rip), %xmm6
- orps %xmm5, %xmm0
- movaps %xmm4, %xmm5
- andps %xmm4, %xmm6
- andnps %xmm9, %xmm5
- andps %xmm1, %xmm4
- orps %xmm4, %xmm5
- andps %xmm10, %xmm7
- divpd %xmm5, %xmm0
- movq iCHK_WORK_SUB+__svml_datan2_data_internal(%rip), %xmm2
- xorl %edx, %edx
-
-/* Check if y and x are on main path. */
- pshufd $221, %xmm9, %xmm3
- xorl %eax, %eax
- pshufd $221, %xmm1, %xmm13
- psubd %xmm2, %xmm3
- psubd %xmm2, %xmm13
- movdqa %xmm3, %xmm4
- movq iCHK_WORK_CMP+__svml_datan2_data_internal(%rip), %xmm12
- movdqa %xmm13, %xmm14
- pcmpgtd %xmm12, %xmm4
- pcmpeqd %xmm12, %xmm3
- pcmpgtd %xmm12, %xmm14
- pcmpeqd %xmm12, %xmm13
-
-/* Polynomial. */
- movaps %xmm0, %xmm12
- por %xmm3, %xmm4
- mulpd %xmm0, %xmm12
-
-/* P = A19*R2 + A18 */
- movups dA19+__svml_datan2_data_internal(%rip), %xmm15
- movaps %xmm11, %xmm2
- mulpd %xmm12, %xmm15
- addpd dA18+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A17 */
- mulpd %xmm12, %xmm15
- addpd dA17+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A16 */
- mulpd %xmm12, %xmm15
- addpd dA16+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A15 */
- mulpd %xmm12, %xmm15
- addpd dA15+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A14 */
- mulpd %xmm12, %xmm15
- addpd dA14+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A13 */
- mulpd %xmm12, %xmm15
- addpd dA13+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A12 */
- mulpd %xmm12, %xmm15
- addpd dA12+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A11 */
- mulpd %xmm12, %xmm15
- addpd dA11+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A10 */
- mulpd %xmm12, %xmm15
- addpd dA10+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A09 */
- mulpd %xmm12, %xmm15
- addpd dA09+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A08 */
- mulpd %xmm12, %xmm15
- addpd dA08+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A07 */
- mulpd %xmm12, %xmm15
- addpd dA07+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A06 */
- mulpd %xmm12, %xmm15
- addpd dA06+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A05 */
- mulpd %xmm12, %xmm15
- addpd dA05+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A04 */
- mulpd %xmm12, %xmm15
- addpd dA04+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A03 */
- mulpd %xmm12, %xmm15
- addpd dA03+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A02 */
- mulpd %xmm12, %xmm15
- addpd dA02+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 + A01 */
- mulpd %xmm12, %xmm15
- addpd dA01+__svml_datan2_data_internal(%rip), %xmm15
-
-/* P = P*R2 */
- mulpd %xmm15, %xmm12
-
-/*
- * Reconstruction.
- * dP=(R+R*dP) + dPIO2
- */
- mulpd %xmm0, %xmm12
- addpd %xmm12, %xmm0
-
-/* if x<0, dPI = Pi, else dPI =0 */
- movups dZERO+__svml_datan2_data_internal(%rip), %xmm3
- por %xmm13, %xmm14
- cmplepd %xmm3, %xmm2
- addpd %xmm6, %xmm0
- andps __svml_datan2_data_internal(%rip), %xmm2
- orps %xmm8, %xmm0
- addpd %xmm2, %xmm0
- por %xmm14, %xmm4
- orps %xmm7, %xmm0
- movmskps %xmm4, %ecx
-
-/* Special branch for fast (vector) processing of zero arguments */
- testb $3, %cl
-
-/* Go to auxilary branch */
- jne L(AUX_BRANCH)
- # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0 xmm1 xmm2 xmm3 xmm4 xmm5 xmm6 xmm7 xmm8 xmm9 xmm10 xmm11
-
-/* Return from auxilary branch
- * for out of main path inputs
- */
+ subq $88, %rsp
+ cfi_def_cfa_offset(96)
+ movaps %xmm1, %xmm11
+
+ /*
+ * #define NO_VECTOR_ZERO_ATAN2_ARGS
+ * Declarations
+ * Variables
+ * Constants
+ * The end of declarations
+ * Implementation
+ * Get r0~=1/B
+ * Cannot be replaced by VQRCP(D, dR0, dB);
+ * Argument Absolute values
+ */
+ movups dABS_MASK+__svml_datan2_data_internal(%rip), %xmm1
+ movaps %xmm0, %xmm10
+ movaps %xmm1, %xmm9
+ andps %xmm10, %xmm1
+ andps %xmm11, %xmm9
+ movaps %xmm1, %xmm4
+ cmpnltpd %xmm9, %xmm4
+
+ /* Argument signs */
+ movups dSIGN_MASK+__svml_datan2_data_internal(%rip), %xmm5
+ movaps %xmm4, %xmm0
+ movaps %xmm5, %xmm8
+ movaps %xmm5, %xmm7
+
+ /*
+ * 1) If y<x then a= y, b=x, PIO2=0
+ * 2) If y>x then a=-x, b=y, PIO2=Pi/2
+ */
+ orps %xmm9, %xmm5
+ andnps %xmm1, %xmm0
+ andps %xmm4, %xmm5
+ andps %xmm11, %xmm8
+ movups dPIO2+__svml_datan2_data_internal(%rip), %xmm6
+ orps %xmm5, %xmm0
+ movaps %xmm4, %xmm5
+ andps %xmm4, %xmm6
+ andnps %xmm9, %xmm5
+ andps %xmm1, %xmm4
+ orps %xmm4, %xmm5
+ andps %xmm10, %xmm7
+ divpd %xmm5, %xmm0
+ movq iCHK_WORK_SUB+__svml_datan2_data_internal(%rip), %xmm2
+ xorl %edx, %edx
+
+ /* Check if y and x are on main path. */
+ pshufd $221, %xmm9, %xmm3
+ xorl %eax, %eax
+ pshufd $221, %xmm1, %xmm13
+ psubd %xmm2, %xmm3
+ psubd %xmm2, %xmm13
+ movdqa %xmm3, %xmm4
+ movq iCHK_WORK_CMP+__svml_datan2_data_internal(%rip), %xmm12
+ movdqa %xmm13, %xmm14
+ pcmpgtd %xmm12, %xmm4
+ pcmpeqd %xmm12, %xmm3
+ pcmpgtd %xmm12, %xmm14
+ pcmpeqd %xmm12, %xmm13
+
+ /* Polynomial. */
+ movaps %xmm0, %xmm12
+ por %xmm3, %xmm4
+ mulpd %xmm0, %xmm12
+
+ /* P = A19*R2 + A18 */
+ movups dA19+__svml_datan2_data_internal(%rip), %xmm15
+ movaps %xmm11, %xmm2
+ mulpd %xmm12, %xmm15
+ addpd dA18+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A17 */
+ mulpd %xmm12, %xmm15
+ addpd dA17+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A16 */
+ mulpd %xmm12, %xmm15
+ addpd dA16+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A15 */
+ mulpd %xmm12, %xmm15
+ addpd dA15+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A14 */
+ mulpd %xmm12, %xmm15
+ addpd dA14+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A13 */
+ mulpd %xmm12, %xmm15
+ addpd dA13+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A12 */
+ mulpd %xmm12, %xmm15
+ addpd dA12+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A11 */
+ mulpd %xmm12, %xmm15
+ addpd dA11+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A10 */
+ mulpd %xmm12, %xmm15
+ addpd dA10+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A09 */
+ mulpd %xmm12, %xmm15
+ addpd dA09+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A08 */
+ mulpd %xmm12, %xmm15
+ addpd dA08+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A07 */
+ mulpd %xmm12, %xmm15
+ addpd dA07+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A06 */
+ mulpd %xmm12, %xmm15
+ addpd dA06+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A05 */
+ mulpd %xmm12, %xmm15
+ addpd dA05+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A04 */
+ mulpd %xmm12, %xmm15
+ addpd dA04+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A03 */
+ mulpd %xmm12, %xmm15
+ addpd dA03+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A02 */
+ mulpd %xmm12, %xmm15
+ addpd dA02+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 + A01 */
+ mulpd %xmm12, %xmm15
+ addpd dA01+__svml_datan2_data_internal(%rip), %xmm15
+
+ /* P = P*R2 */
+ mulpd %xmm15, %xmm12
+
+ /*
+ * Reconstruction.
+ * dP=(R+R*dP) + dPIO2
+ */
+ mulpd %xmm0, %xmm12
+ addpd %xmm12, %xmm0
+
+ /* if x<0, dPI = Pi, else dPI =0 */
+ movups dZERO+__svml_datan2_data_internal(%rip), %xmm3
+ por %xmm13, %xmm14
+ cmplepd %xmm3, %xmm2
+ addpd %xmm6, %xmm0
+ andps __svml_datan2_data_internal(%rip), %xmm2
+ orps %xmm8, %xmm0
+ addpd %xmm2, %xmm0
+ por %xmm14, %xmm4
+ orps %xmm7, %xmm0
+ movmskps %xmm4, %ecx
+
+ /* Special branch for fast (vector) processing of zero arguments */
+ testb $3, %cl
+
+ /* Go to auxilary branch */
+ jne L(AUX_BRANCH)
+ # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0 xmm1 xmm2 xmm3 xmm4 xmm5 xmm6 xmm7 xmm8 xmm9 xmm10 xmm11
+
+ /* Return from auxilary branch
+ * for out of main path inputs
+ */
L(AUX_BRANCH_RETURN):
-/*
- * Special branch for fast (vector) processing of zero arguments
- * The end of implementation
- */
- testl %edx, %edx
+ /*
+ * Special branch for fast (vector) processing of zero arguments
+ * The end of implementation
+ */
+ testl %edx, %edx
-/* Go to special inputs processing branch */
- jne L(SPECIAL_VALUES_BRANCH)
- # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0 xmm10 xmm11
+ /* Go to special inputs processing branch */
+ jne L(SPECIAL_VALUES_BRANCH)
+ # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0 xmm10 xmm11
-/* Restore registers
- * and exit the function
- */
+ /* Restore registers
+ * and exit the function
+ */
L(EXIT):
- addq $88, %rsp
- cfi_def_cfa_offset(8)
- ret
- cfi_def_cfa_offset(96)
+ addq $88, %rsp
+ cfi_def_cfa_offset(8)
+ ret
+ cfi_def_cfa_offset(96)
-/* Branch to process
- * special inputs
- */
+ /* Branch to process
+ * special inputs
+ */
L(SPECIAL_VALUES_BRANCH):
- movups %xmm10, 32(%rsp)
- movups %xmm11, 48(%rsp)
- movups %xmm0, 64(%rsp)
- # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0
-
- movq %r12, 16(%rsp)
- cfi_offset(12, -80)
- movl %eax, %r12d
- movq %r13, 8(%rsp)
- cfi_offset(13, -88)
- movl %edx, %r13d
- movq %r14, (%rsp)
- cfi_offset(14, -96)
- # LOE rbx rbp r15 r12d r13d
-
-/* Range mask
- * bits check
- */
+ movups %xmm10, 32(%rsp)
+ movups %xmm11, 48(%rsp)
+ movups %xmm0, 64(%rsp)
+ # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0
+
+ movq %r12, 16(%rsp)
+ cfi_offset(12, -80)
+ movl %eax, %r12d
+ movq %r13, 8(%rsp)
+ cfi_offset(13, -88)
+ movl %edx, %r13d
+ movq %r14, (%rsp)
+ cfi_offset(14, -96)
+ # LOE rbx rbp r15 r12d r13d
+
+ /* Range mask
+ * bits check
+ */
L(RANGEMASK_CHECK):
- btl %r12d, %r13d
+ btl %r12d, %r13d
-/* Call scalar math function */
- jc L(SCALAR_MATH_CALL)
- # LOE rbx rbp r15 r12d r13d
+ /* Call scalar math function */
+ jc L(SCALAR_MATH_CALL)
+ # LOE rbx rbp r15 r12d r13d
-/* Special inputs
- * processing loop
- */
+ /* Special inputs
+ * processing loop
+ */
L(SPECIAL_VALUES_LOOP):
- incl %r12d
- cmpl $2, %r12d
-
-/* Check bits in range mask */
- jl L(RANGEMASK_CHECK)
- # LOE rbx rbp r15 r12d r13d
-
- movq 16(%rsp), %r12
- cfi_restore(12)
- movq 8(%rsp), %r13
- cfi_restore(13)
- movq (%rsp), %r14
- cfi_restore(14)
- movups 64(%rsp), %xmm0
-
-/* Go to exit */
- jmp L(EXIT)
- cfi_offset(12, -80)
- cfi_offset(13, -88)
- cfi_offset(14, -96)
- # LOE rbx rbp r12 r13 r14 r15 xmm0
-
-/* Scalar math fucntion call
- * to process special input
- */
+ incl %r12d
+ cmpl $2, %r12d
+
+ /* Check bits in range mask */
+ jl L(RANGEMASK_CHECK)
+ # LOE rbx rbp r15 r12d r13d
+
+ movq 16(%rsp), %r12
+ cfi_restore(12)
+ movq 8(%rsp), %r13
+ cfi_restore(13)
+ movq (%rsp), %r14
+ cfi_restore(14)
+ movups 64(%rsp), %xmm0
+
+ /* Go to exit */
+ jmp L(EXIT)
+ cfi_offset(12, -80)
+ cfi_offset(13, -88)
+ cfi_offset(14, -96)
+ # LOE rbx rbp r12 r13 r14 r15 xmm0
+
+ /* Scalar math fucntion call
+ * to process special input
+ */
L(SCALAR_MATH_CALL):
- movl %r12d, %r14d
- movsd 32(%rsp,%r14,8), %xmm0
- movsd 48(%rsp,%r14,8), %xmm1
- call atan2@PLT
- # LOE rbx rbp r14 r15 r12d r13d xmm0
-
- movsd %xmm0, 64(%rsp,%r14,8)
-
-/* Process special inputs in loop */
- jmp L(SPECIAL_VALUES_LOOP)
- cfi_restore(12)
- cfi_restore(13)
- cfi_restore(14)
- # LOE rbx rbp r15 r12d r13d
-
-/* Auxilary branch
- * for out of main path inputs
- */
+ movl %r12d, %r14d
+ movsd 32(%rsp, %r14, 8), %xmm0
+ movsd 48(%rsp, %r14, 8), %xmm1
+ call atan2@PLT
+ # LOE rbx rbp r14 r15 r12d r13d xmm0
-L(AUX_BRANCH):
-/* Check if both X & Y are not NaNs: iXYnotNAN */
- movaps %xmm11, %xmm13
- movaps %xmm10, %xmm12
- cmpordpd %xmm11, %xmm13
- cmpordpd %xmm10, %xmm12
+ movsd %xmm0, 64(%rsp, %r14, 8)
-/* Check if at least on of Y or Y is zero: iAXAYZERO */
- cmpeqpd %xmm3, %xmm9
- cmpeqpd %xmm3, %xmm1
+ /* Process special inputs in loop */
+ jmp L(SPECIAL_VALUES_LOOP)
+ cfi_restore(12)
+ cfi_restore(13)
+ cfi_restore(14)
+ # LOE rbx rbp r15 r12d r13d
-/*
- * Path for zero arguments (at least one of both)
- * Check if both args are zeros (den. is zero)
- */
- cmpeqpd %xmm3, %xmm5
- andps %xmm12, %xmm13
- orps %xmm1, %xmm9
- pshufd $221, %xmm9, %xmm1
- pshufd $221, %xmm13, %xmm9
-
-/* Check if at least on of Y or Y is zero and not NaN: iAXAYZEROnotNAN */
- pand %xmm9, %xmm1
-
-/* Exclude from previous callout mask zero (and not NaN) arguments */
- movdqa %xmm1, %xmm14
- pandn %xmm4, %xmm14
-
-/* Set sPIO2 to zero if den. is zero */
- movaps %xmm5, %xmm4
- andnps %xmm6, %xmm4
- andps %xmm3, %xmm5
-
-/* Res = sign(Y)*(X<0)?(PIO2+PI):PIO2 */
- pshufd $221, %xmm3, %xmm3
- orps %xmm5, %xmm4
- pshufd $221, %xmm11, %xmm5
- orps %xmm8, %xmm4
- pcmpgtd %xmm5, %xmm3
- pshufd $80, %xmm3, %xmm6
- andps %xmm2, %xmm6
- addpd %xmm6, %xmm4
-
-/* Go to callout */
- movmskps %xmm14, %edx
-
-/* Merge results from main and spec path */
- pshufd $80, %xmm1, %xmm2
- orps %xmm7, %xmm4
- movdqa %xmm2, %xmm7
- andps %xmm2, %xmm4
- andnps %xmm0, %xmm7
- andl $3, %edx
- movaps %xmm7, %xmm0
- orps %xmm4, %xmm0
-
-/* Return to main vector processing path */
- jmp L(AUX_BRANCH_RETURN)
- # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0 xmm10 xmm11
+ /* Auxilary branch
+ * for out of main path inputs
+ */
+
+L(AUX_BRANCH):
+ /* Check if both X & Y are not NaNs: iXYnotNAN */
+ movaps %xmm11, %xmm13
+ movaps %xmm10, %xmm12
+ cmpordpd %xmm11, %xmm13
+ cmpordpd %xmm10, %xmm12
+
+ /* Check if at least on of Y or Y is zero: iAXAYZERO */
+ cmpeqpd %xmm3, %xmm9
+ cmpeqpd %xmm3, %xmm1
+
+ /*
+ * Path for zero arguments (at least one of both)
+ * Check if both args are zeros (den. is zero)
+ */
+ cmpeqpd %xmm3, %xmm5
+ andps %xmm12, %xmm13
+ orps %xmm1, %xmm9
+ pshufd $221, %xmm9, %xmm1
+ pshufd $221, %xmm13, %xmm9
+
+ /* Check if at least on of Y or Y is zero and not NaN: iAXAYZEROnotNAN */
+ pand %xmm9, %xmm1
+
+ /* Exclude from previous callout mask zero (and not NaN) arguments */
+ movdqa %xmm1, %xmm14
+ pandn %xmm4, %xmm14
+
+ /* Set sPIO2 to zero if den. is zero */
+ movaps %xmm5, %xmm4
+ andnps %xmm6, %xmm4
+ andps %xmm3, %xmm5
+
+ /* Res = sign(Y)*(X<0)?(PIO2+PI):PIO2 */
+ pshufd $221, %xmm3, %xmm3
+ orps %xmm5, %xmm4
+ pshufd $221, %xmm11, %xmm5
+ orps %xmm8, %xmm4
+ pcmpgtd %xmm5, %xmm3
+ pshufd $80, %xmm3, %xmm6
+ andps %xmm2, %xmm6
+ addpd %xmm6, %xmm4
+
+ /* Go to callout */
+ movmskps %xmm14, %edx
+
+ /* Merge results from main and spec path */
+ pshufd $80, %xmm1, %xmm2
+ orps %xmm7, %xmm4
+ movdqa %xmm2, %xmm7
+ andps %xmm2, %xmm4
+ andnps %xmm0, %xmm7
+ andl $3, %edx
+ movaps %xmm7, %xmm0
+ orps %xmm4, %xmm0
+
+ /* Return to main vector processing path */
+ jmp L(AUX_BRANCH_RETURN)
+ # LOE rbx rbp r12 r13 r14 r15 eax edx xmm0 xmm10 xmm11
END(_ZGVbN2vv_atan2_sse4)
- .section .rodata, "a"
- .align 16
+ .section .rodata, "a"
+ .align 16
#ifdef __svml_datan2_data_internal_typedef
typedef unsigned int VUINT32;
typedef struct {
- __declspec(align(16)) VUINT32 dPI[2][2];
- __declspec(align(16)) VUINT32 dPIO2[2][2];
- __declspec(align(16)) VUINT32 dA19[2][2];
- __declspec(align(16)) VUINT32 dA18[2][2];
- __declspec(align(16)) VUINT32 dA17[2][2];
- __declspec(align(16)) VUINT32 dA16[2][2];
- __declspec(align(16)) VUINT32 dA15[2][2];
- __declspec(align(16)) VUINT32 dA14[2][2];
- __declspec(align(16)) VUINT32 dA13[2][2];
- __declspec(align(16)) VUINT32 dA12[2][2];
- __declspec(align(16)) VUINT32 dA11[2][2];
- __declspec(align(16)) VUINT32 dA10[2][2];
- __declspec(align(16)) VUINT32 dA09[2][2];
- __declspec(align(16)) VUINT32 dA08[2][2];
- __declspec(align(16)) VUINT32 dA07[2][2];
- __declspec(align(16)) VUINT32 dA06[2][2];
- __declspec(align(16)) VUINT32 dA05[2][2];
- __declspec(align(16)) VUINT32 dA04[2][2];
- __declspec(align(16)) VUINT32 dA03[2][2];
- __declspec(align(16)) VUINT32 dA02[2][2];
- __declspec(align(16)) VUINT32 dA01[2][2];
- __declspec(align(16)) VUINT32 dA00[2][2];
- __declspec(align(16)) VUINT32 dSIGN_MASK[2][2];
- __declspec(align(16)) VUINT32 iCHK_WORK_SUB[4][1];
- __declspec(align(16)) VUINT32 iCHK_WORK_CMP[4][1];
- __declspec(align(16)) VUINT32 dABS_MASK[2][2];
- __declspec(align(16)) VUINT32 dZERO[2][2];
+ __declspec(align(16)) VUINT32 dPI[2][2];
+ __declspec(align(16)) VUINT32 dPIO2[2][2];
+ __declspec(align(16)) VUINT32 dA19[2][2];
+ __declspec(align(16)) VUINT32 dA18[2][2];
+ __declspec(align(16)) VUINT32 dA17[2][2];
+ __declspec(align(16)) VUINT32 dA16[2][2];
+ __declspec(align(16)) VUINT32 dA15[2][2];
+ __declspec(align(16)) VUINT32 dA14[2][2];
+ __declspec(align(16)) VUINT32 dA13[2][2];
+ __declspec(align(16)) VUINT32 dA12[2][2];
+ __declspec(align(16)) VUINT32 dA11[2][2];
+ __declspec(align(16)) VUINT32 dA10[2][2];
+ __declspec(align(16)) VUINT32 dA09[2][2];
+ __declspec(align(16)) VUINT32 dA08[2][2];
+ __declspec(align(16)) VUINT32 dA07[2][2];
+ __declspec(align(16)) VUINT32 dA06[2][2];
+ __declspec(align(16)) VUINT32 dA05[2][2];
+ __declspec(align(16)) VUINT32 dA04[2][2];
+ __declspec(align(16)) VUINT32 dA03[2][2];
+ __declspec(align(16)) VUINT32 dA02[2][2];
+ __declspec(align(16)) VUINT32 dA01[2][2];
+ __declspec(align(16)) VUINT32 dA00[2][2];
+ __declspec(align(16)) VUINT32 dSIGN_MASK[2][2];
+ __declspec(align(16)) VUINT32 iCHK_WORK_SUB[4][1];
+ __declspec(align(16)) VUINT32 iCHK_WORK_CMP[4][1];
+ __declspec(align(16)) VUINT32 dABS_MASK[2][2];
+ __declspec(align(16)) VUINT32 dZERO[2][2];
} __svml_datan2_data_internal;
#endif
__svml_datan2_data_internal:
- .quad 0x400921FB54442D18, 0x400921FB54442D18 //dPI
- .align 16
- .quad 0x3FF921FB54442D18, 0x3FF921FB54442D18 //dPIO2
- .align 16
- .quad 0xBEF4FDB537ABC7A3, 0xBEF4FDB537ABC7A3 // dA19
- .align 16
- .quad 0x3F2CED0A36665209, 0x3F2CED0A36665209 // dA18
- .align 16
- .quad 0xBF52E67C93954C23, 0xBF52E67C93954C23 // dA17
- .align 16
- .quad 0x3F6F5A1DAE82AFB3, 0x3F6F5A1DAE82AFB3 // dA16
- .align 16
- .quad 0xBF82B2EC618E4BAD, 0xBF82B2EC618E4BAD // dA15
- .align 16
- .quad 0x3F914F4C661116A5, 0x3F914F4C661116A5 // dA14
- .align 16
- .quad 0xBF9A5E83B081F69C, 0xBF9A5E83B081F69C // dA13
- .align 16
- .quad 0x3FA169980CB6AD4F, 0x3FA169980CB6AD4F // dA12
- .align 16
- .quad 0xBFA4EFA2E563C1BC, 0xBFA4EFA2E563C1BC // dA11
- .align 16
- .quad 0x3FA7EC0FBC50683B, 0x3FA7EC0FBC50683B // dA10
- .align 16
- .quad 0xBFAAD261EAA09954, 0xBFAAD261EAA09954 // dA09
- .align 16
- .quad 0x3FAE1749BD612DCF, 0x3FAE1749BD612DCF // dA08
- .align 16
- .quad 0xBFB11084009435E0, 0xBFB11084009435E0 // dA07
- .align 16
- .quad 0x3FB3B12A49295651, 0x3FB3B12A49295651 // dA06
- .align 16
- .quad 0xBFB745D009BADA94, 0xBFB745D009BADA94 // dA05
- .align 16
- .quad 0x3FBC71C707F7D5B5, 0x3FBC71C707F7D5B5 // dA04
- .align 16
- .quad 0xBFC2492491EE55C7, 0xBFC2492491EE55C7 // dA03
- .align 16
- .quad 0x3FC999999997EE34, 0x3FC999999997EE34 // dA02
- .align 16
- .quad 0xBFD55555555553C5, 0xBFD55555555553C5 // dA01
- .align 16
- .quad 0x3FF0000000000000, 0x3FF0000000000000 // dA00
- .align 16
- .quad 0x8000000000000000, 0x8000000000000000 //dSIGN_MASK
- .align 16
- .long 0x80300000, 0x80300000, 0x80300000, 0x80300000 //iCHK_WORK_SUB
- .align 16
- .long 0xfdd00000, 0xfdd00000, 0xfdd00000, 0xfdd00000 //iCHK_WORK_CMP
- .align 16
- .quad 0x7fffffffffffffff, 0x7fffffffffffffff //dABS_MASK
- .align 16
- .quad 0x0000000000000000, 0x0000000000000000 //dZERO
- .align 16
- .type __svml_datan2_data_internal,@object
- .size __svml_datan2_data_internal,.-__svml_datan2_data_internal
+ .quad 0x400921FB54442D18, 0x400921FB54442D18 // dPI
+ .align 16
+ .quad 0x3FF921FB54442D18, 0x3FF921FB54442D18 // dPIO2
+ .align 16
+ .quad 0xBEF4FDB537ABC7A3, 0xBEF4FDB537ABC7A3 // dA19
+ .align 16
+ .quad 0x3F2CED0A36665209, 0x3F2CED0A36665209 // dA18
+ .align 16
+ .quad 0xBF52E67C93954C23, 0xBF52E67C93954C23 // dA17
+ .align 16
+ .quad 0x3F6F5A1DAE82AFB3, 0x3F6F5A1DAE82AFB3 // dA16
+ .align 16
+ .quad 0xBF82B2EC618E4BAD, 0xBF82B2EC618E4BAD // dA15
+ .align 16
+ .quad 0x3F914F4C661116A5, 0x3F914F4C661116A5 // dA14
+ .align 16
+ .quad 0xBF9A5E83B081F69C, 0xBF9A5E83B081F69C // dA13
+ .align 16
+ .quad 0x3FA169980CB6AD4F, 0x3FA169980CB6AD4F // dA12
+ .align 16
+ .quad 0xBFA4EFA2E563C1BC, 0xBFA4EFA2E563C1BC // dA11
+ .align 16
+ .quad 0x3FA7EC0FBC50683B, 0x3FA7EC0FBC50683B // dA10
+ .align 16
+ .quad 0xBFAAD261EAA09954, 0xBFAAD261EAA09954 // dA09
+ .align 16
+ .quad 0x3FAE1749BD612DCF, 0x3FAE1749BD612DCF // dA08
+ .align 16
+ .quad 0xBFB11084009435E0, 0xBFB11084009435E0 // dA07
+ .align 16
+ .quad 0x3FB3B12A49295651, 0x3FB3B12A49295651 // dA06
+ .align 16
+ .quad 0xBFB745D009BADA94, 0xBFB745D009BADA94 // dA05
+ .align 16
+ .quad 0x3FBC71C707F7D5B5, 0x3FBC71C707F7D5B5 // dA04
+ .align 16
+ .quad 0xBFC2492491EE55C7, 0xBFC2492491EE55C7 // dA03
+ .align 16
+ .quad 0x3FC999999997EE34, 0x3FC999999997EE34 // dA02
+ .align 16
+ .quad 0xBFD55555555553C5, 0xBFD55555555553C5 // dA01
+ .align 16
+ .quad 0x3FF0000000000000, 0x3FF0000000000000 // dA00
+ .align 16
+ .quad 0x8000000000000000, 0x8000000000000000 // dSIGN_MASK
+ .align 16
+ .long 0x80300000, 0x80300000, 0x80300000, 0x80300000 // iCHK_WORK_SUB
+ .align 16
+ .long 0xfdd00000, 0xfdd00000, 0xfdd00000, 0xfdd00000 // iCHK_WORK_CMP
+ .align 16
+ .quad 0x7fffffffffffffff, 0x7fffffffffffffff // dABS_MASK
+ .align 16
+ .quad 0x0000000000000000, 0x0000000000000000 // dZERO
+ .align 16
+ .type __svml_datan2_data_internal, @object
+ .size __svml_datan2_data_internal, .-__svml_datan2_data_internal