x86-64: Add vector log10/log10f implementation to libmvec

Implement vectorized log10/log10f containing SSE, AVX, AVX2 and
AVX512 versions for libmvec as per vector ABI.  It also contains
accuracy and ABI tests for vector log10/log10f with regenerated ulps.

Reviewed-by: H.J. Lu <hjl.tools@gmail.com>

50 files changed, 3758 insertions, 1 deletions

diff --git a/bits/libm-simd-decl-stubs.h b/bits/libm-simd-decl-stubs.h
index 31878bf..4ad584c 100644
--- a/bits/libm-simd-decl-stubs.h
+++ b/bits/libm-simd-decl-stubs.h
@@ -219,4 +219,15 @@
 #define __DECL_SIMD_atan2f32x
 #define __DECL_SIMD_atan2f64x
 #define __DECL_SIMD_atan2f128x
+
+#define __DECL_SIMD_log10
+#define __DECL_SIMD_log10f
+#define __DECL_SIMD_log10l
+#define __DECL_SIMD_log10f16
+#define __DECL_SIMD_log10f32
+#define __DECL_SIMD_log10f64
+#define __DECL_SIMD_log10f128
+#define __DECL_SIMD_log10f32x
+#define __DECL_SIMD_log10f64x
+#define __DECL_SIMD_log10f128x
 #endif
diff --git a/math/bits/mathcalls.h b/math/bits/mathcalls.h
index 1bd4911..f213847 100644
--- a/math/bits/mathcalls.h
+++ b/math/bits/mathcalls.h
@@ -104,7 +104,7 @@ __MATHCALL (ldexp,, (_Mdouble_ __x, int __exponent));
 __MATHCALL_VEC (log,, (_Mdouble_ __x));
 
 /* Base-ten logarithm of X.  */
-__MATHCALL (log10,, (_Mdouble_ __x));
+__MATHCALL_VEC (log10,, (_Mdouble_ __x));
 
 /* Break VALUE into integral and fractional parts.  */
 __MATHCALL (modf,, (_Mdouble_ __x, _Mdouble_ *__iptr)) __nonnull ((2));
diff --git a/sysdeps/unix/sysv/linux/x86_64/libmvec.abilist b/sysdeps/unix/sysv/linux/x86_64/libmvec.abilist
index 2b3b8d3..8108a2a 100644
--- a/sysdeps/unix/sysv/linux/x86_64/libmvec.abilist
+++ b/sysdeps/unix/sysv/linux/x86_64/libmvec.abilist
@@ -54,6 +54,7 @@ GLIBC_2.35 _ZGVbN2v_cosh F
 GLIBC_2.35 _ZGVbN2v_exp10 F
 GLIBC_2.35 _ZGVbN2v_exp2 F
 GLIBC_2.35 _ZGVbN2v_expm1 F
+GLIBC_2.35 _ZGVbN2v_log10 F
 GLIBC_2.35 _ZGVbN2v_sinh F
 GLIBC_2.35 _ZGVbN2vv_atan2 F
 GLIBC_2.35 _ZGVbN2vv_hypot F
@@ -65,6 +66,7 @@ GLIBC_2.35 _ZGVbN4v_coshf F
 GLIBC_2.35 _ZGVbN4v_exp10f F
 GLIBC_2.35 _ZGVbN4v_exp2f F
 GLIBC_2.35 _ZGVbN4v_expm1f F
+GLIBC_2.35 _ZGVbN4v_log10f F
 GLIBC_2.35 _ZGVbN4v_sinhf F
 GLIBC_2.35 _ZGVbN4vv_atan2f F
 GLIBC_2.35 _ZGVbN4vv_hypotf F
@@ -76,6 +78,7 @@ GLIBC_2.35 _ZGVcN4v_cosh F
 GLIBC_2.35 _ZGVcN4v_exp10 F
 GLIBC_2.35 _ZGVcN4v_exp2 F
 GLIBC_2.35 _ZGVcN4v_expm1 F
+GLIBC_2.35 _ZGVcN4v_log10 F
 GLIBC_2.35 _ZGVcN4v_sinh F
 GLIBC_2.35 _ZGVcN4vv_atan2 F
 GLIBC_2.35 _ZGVcN4vv_hypot F
@@ -87,6 +90,7 @@ GLIBC_2.35 _ZGVcN8v_coshf F
 GLIBC_2.35 _ZGVcN8v_exp10f F
 GLIBC_2.35 _ZGVcN8v_exp2f F
 GLIBC_2.35 _ZGVcN8v_expm1f F
+GLIBC_2.35 _ZGVcN8v_log10f F
 GLIBC_2.35 _ZGVcN8v_sinhf F
 GLIBC_2.35 _ZGVcN8vv_atan2f F
 GLIBC_2.35 _ZGVcN8vv_hypotf F
@@ -98,6 +102,7 @@ GLIBC_2.35 _ZGVdN4v_cosh F
 GLIBC_2.35 _ZGVdN4v_exp10 F
 GLIBC_2.35 _ZGVdN4v_exp2 F
 GLIBC_2.35 _ZGVdN4v_expm1 F
+GLIBC_2.35 _ZGVdN4v_log10 F
 GLIBC_2.35 _ZGVdN4v_sinh F
 GLIBC_2.35 _ZGVdN4vv_atan2 F
 GLIBC_2.35 _ZGVdN4vv_hypot F
@@ -109,6 +114,7 @@ GLIBC_2.35 _ZGVdN8v_coshf F
 GLIBC_2.35 _ZGVdN8v_exp10f F
 GLIBC_2.35 _ZGVdN8v_exp2f F
 GLIBC_2.35 _ZGVdN8v_expm1f F
+GLIBC_2.35 _ZGVdN8v_log10f F
 GLIBC_2.35 _ZGVdN8v_sinhf F
 GLIBC_2.35 _ZGVdN8vv_atan2f F
 GLIBC_2.35 _ZGVdN8vv_hypotf F
@@ -120,6 +126,7 @@ GLIBC_2.35 _ZGVeN16v_coshf F
 GLIBC_2.35 _ZGVeN16v_exp10f F
 GLIBC_2.35 _ZGVeN16v_exp2f F
 GLIBC_2.35 _ZGVeN16v_expm1f F
+GLIBC_2.35 _ZGVeN16v_log10f F
 GLIBC_2.35 _ZGVeN16v_sinhf F
 GLIBC_2.35 _ZGVeN16vv_atan2f F
 GLIBC_2.35 _ZGVeN16vv_hypotf F
@@ -131,6 +138,7 @@ GLIBC_2.35 _ZGVeN8v_cosh F
 GLIBC_2.35 _ZGVeN8v_exp10 F
 GLIBC_2.35 _ZGVeN8v_exp2 F
 GLIBC_2.35 _ZGVeN8v_expm1 F
+GLIBC_2.35 _ZGVeN8v_log10 F
 GLIBC_2.35 _ZGVeN8v_sinh F
 GLIBC_2.35 _ZGVeN8vv_atan2 F
 GLIBC_2.35 _ZGVeN8vv_hypot F
diff --git a/sysdeps/x86/fpu/bits/math-vector.h b/sysdeps/x86/fpu/bits/math-vector.h
index 62f2890..64e80ad 100644
--- a/sysdeps/x86/fpu/bits/math-vector.h
+++ b/sysdeps/x86/fpu/bits/math-vector.h
@@ -102,6 +102,10 @@
 #  define __DECL_SIMD_atan2 __DECL_SIMD_x86_64
 #  undef __DECL_SIMD_atan2f
 #  define __DECL_SIMD_atan2f __DECL_SIMD_x86_64
+#  undef __DECL_SIMD_log10
+#  define __DECL_SIMD_log10 __DECL_SIMD_x86_64
+#  undef __DECL_SIMD_log10f
+#  define __DECL_SIMD_log10f __DECL_SIMD_x86_64
 
 # endif
 #endif
diff --git a/sysdeps/x86/fpu/finclude/math-vector-fortran.h b/sysdeps/x86/fpu/finclude/math-vector-fortran.h
index 2269b74..f5050c6 100644
--- a/sysdeps/x86/fpu/finclude/math-vector-fortran.h
+++ b/sysdeps/x86/fpu/finclude/math-vector-fortran.h
@@ -50,6 +50,8 @@
 !GCC$ builtin (cbrtf) attributes simd (notinbranch) if('x86_64')
 !GCC$ builtin (atan2) attributes simd (notinbranch) if('x86_64')
 !GCC$ builtin (atan2f) attributes simd (notinbranch) if('x86_64')
+!GCC$ builtin (log10) attributes simd (notinbranch) if('x86_64')
+!GCC$ builtin (log10f) attributes simd (notinbranch) if('x86_64')
 
 !GCC$ builtin (cos) attributes simd (notinbranch) if('x32')
 !GCC$ builtin (cosf) attributes simd (notinbranch) if('x32')
@@ -85,3 +87,5 @@
 !GCC$ builtin (cbrtf) attributes simd (notinbranch) if('x32')
 !GCC$ builtin (atan2) attributes simd (notinbranch) if('x32')
 !GCC$ builtin (atan2f) attributes simd (notinbranch) if('x32')
+!GCC$ builtin (log10) attributes simd (notinbranch) if('x32')
+!GCC$ builtin (log10f) attributes simd (notinbranch) if('x32')
diff --git a/sysdeps/x86_64/fpu/Makeconfig b/sysdeps/x86_64/fpu/Makeconfig
index 96a4085..ba37044 100644
--- a/sysdeps/x86_64/fpu/Makeconfig
+++ b/sysdeps/x86_64/fpu/Makeconfig
@@ -35,6 +35,7 @@ libmvec-funcs = \
   expm1 \
   hypot \
   log \
+  log10 \
   pow \
   sin \
   sincos \
diff --git a/sysdeps/x86_64/fpu/Versions b/sysdeps/x86_64/fpu/Versions
index f58c98e..8beaf07 100644
--- a/sysdeps/x86_64/fpu/Versions
+++ b/sysdeps/x86_64/fpu/Versions
@@ -22,6 +22,7 @@ libmvec {
     _ZGVbN2v_exp10; _ZGVcN4v_exp10; _ZGVdN4v_exp10; _ZGVeN8v_exp10;
     _ZGVbN2v_exp2; _ZGVcN4v_exp2; _ZGVdN4v_exp2; _ZGVeN8v_exp2;
     _ZGVbN2v_expm1; _ZGVcN4v_expm1; _ZGVdN4v_expm1; _ZGVeN8v_expm1;
+    _ZGVbN2v_log10; _ZGVcN4v_log10; _ZGVdN4v_log10; _ZGVeN8v_log10;
     _ZGVbN2v_sinh; _ZGVcN4v_sinh; _ZGVdN4v_sinh; _ZGVeN8v_sinh;
     _ZGVbN2vv_atan2; _ZGVcN4vv_atan2; _ZGVdN4vv_atan2; _ZGVeN8vv_atan2;
     _ZGVbN2vv_hypot; _ZGVcN4vv_hypot; _ZGVdN4vv_hypot; _ZGVeN8vv_hypot;
@@ -33,6 +34,7 @@ libmvec {
     _ZGVbN4v_exp10f; _ZGVcN8v_exp10f; _ZGVdN8v_exp10f; _ZGVeN16v_exp10f;
     _ZGVbN4v_exp2f; _ZGVcN8v_exp2f; _ZGVdN8v_exp2f; _ZGVeN16v_exp2f;
     _ZGVbN4v_expm1f; _ZGVcN8v_expm1f; _ZGVdN8v_expm1f; _ZGVeN16v_expm1f;
+    _ZGVbN4v_log10f; _ZGVcN8v_log10f; _ZGVdN8v_log10f; _ZGVeN16v_log10f;
     _ZGVbN4v_sinhf; _ZGVcN8v_sinhf; _ZGVdN8v_sinhf; _ZGVeN16v_sinhf;
     _ZGVbN4vv_atan2f; _ZGVcN8vv_atan2f; _ZGVdN8vv_atan2f; _ZGVeN16vv_atan2f;
     _ZGVbN4vv_hypotf; _ZGVcN8vv_hypotf; _ZGVdN8vv_hypotf; _ZGVeN16vv_hypotf;
diff --git a/sysdeps/x86_64/fpu/libm-test-ulps b/sysdeps/x86_64/fpu/libm-test-ulps
index 6f59c61..b0cd9d6 100644
--- a/sysdeps/x86_64/fpu/libm-test-ulps
+++ b/sysdeps/x86_64/fpu/libm-test-ulps
@@ -1641,6 +1641,26 @@ float: 2
 float128: 1
 ldouble: 1
 
+Function: "log10_vlen16":
+float: 1
+
+Function: "log10_vlen2":
+double: 1
+
+Function: "log10_vlen4":
+double: 1
+float: 1
+
+Function: "log10_vlen4_avx2":
+double: 1
+
+Function: "log10_vlen8":
+double: 1
+float: 1
+
+Function: "log10_vlen8_avx2":
+float: 1
+
 Function: "log1p":
 double: 1
 float: 1
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core-sse2.S b/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core-sse2.S
new file mode 100644
index 0000000..e654db6
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core-sse2.S
@@ -0,0 +1,20 @@
+/* SSE2 version of vectorized log10, vector length is 2.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define _ZGVbN2v_log10 _ZGVbN2v_log10_sse2
+#include "../svml_d_log102_core.S"
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core.c b/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core.c
new file mode 100644
index 0000000..1c775f3
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core.c
@@ -0,0 +1,27 @@
+/* Multiple versions of vectorized log10, vector length is 2.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define SYMBOL_NAME _ZGVbN2v_log10
+#include "ifunc-mathvec-sse4_1.h"
+
+libc_ifunc_redirected (REDIRECT_NAME, SYMBOL_NAME, IFUNC_SELECTOR ());
+
+#ifdef SHARED
+__hidden_ver1 (_ZGVbN2v_log10, __GI__ZGVbN2v_log10, __redirect__ZGVbN2v_log10)
+  __attribute__ ((visibility ("hidden")));
+#endif
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core_sse4.S b/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core_sse4.S
new file mode 100644
index 0000000..f82364c
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log102_core_sse4.S
@@ -0,0 +1,1089 @@
+/* Function log10 vectorized with SSE4.
+   Copyright (C) 2021 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
+   https://www.gnu.org/licenses/.  */
+
+/*
+ * ALGORITHM DESCRIPTION:
+ *
+ *    Get short reciprocal approximation Rcp ~ 1/mantissa(x)
+ *    R = Rcp*x - 1.0
+ *    log10(x) = k*log10(2.0) - log10(Rcp) + poly_approximation(R)
+ *       log10(Rcp) is tabulated
+ *
+ *
+ */
+
+/* Offsets for data table __svml_dlog10_data_internal
+ */
+#define Log_HA_table                  	0
+#define Log_LA_table                  	4112
+#define poly_coeff                    	8224
+#define ExpMask                       	8304
+#define Two10                         	8320
+#define MinNorm                       	8336
+#define MaxNorm                       	8352
+#define HalfMask                      	8368
+#define One                           	8384
+#define Threshold                     	8400
+#define Bias                          	8416
+#define Bias1                         	8432
+#define L2                            	8448
+
+/* Lookup bias for data table __svml_dlog10_data_internal.  */
+#define Table_Lookup_Bias               -0x406ff0
+
+#include <sysdep.h>
+
+        .text
+	.section .text.sse4,"ax",@progbits
+ENTRY(_ZGVbN2v_log10_sse4)
+        pushq     %rbp
+        cfi_def_cfa_offset(16)
+        movq      %rsp, %rbp
+        cfi_def_cfa(6, 16)
+        cfi_offset(6, -16)
+        andq      $-32, %rsp
+        subq      $64, %rsp
+
+/* exponent bits */
+        movaps    %xmm0, %xmm5
+
+/* preserve mantissa, set input exponent to 2^(-10) */
+        movups    ExpMask+__svml_dlog10_data_internal(%rip), %xmm1
+        psrlq     $20, %xmm5
+        andps     %xmm0, %xmm1
+        lea       Table_Lookup_Bias+__svml_dlog10_data_internal(%rip), %rsi
+        orps      Two10+__svml_dlog10_data_internal(%rip), %xmm1
+
+/* check range */
+        movaps    %xmm0, %xmm8
+
+/* reciprocal approximation good to at least 11 bits */
+        cvtpd2ps  %xmm1, %xmm2
+        cmpltpd   MinNorm+__svml_dlog10_data_internal(%rip), %xmm8
+        movlhps   %xmm2, %xmm2
+        movaps    %xmm0, %xmm7
+        rcpps     %xmm2, %xmm3
+        cmpnlepd  MaxNorm+__svml_dlog10_data_internal(%rip), %xmm7
+        cvtps2pd  %xmm3, %xmm12
+
+/* round reciprocal to nearest integer, will have 1+9 mantissa bits */
+        movups    .FLT_12(%rip), %xmm4
+        orps      %xmm7, %xmm8
+        addpd     %xmm4, %xmm12
+
+/* combine and get argument value range mask */
+        movmskpd  %xmm8, %edx
+
+/* argument reduction */
+        movups    HalfMask+__svml_dlog10_data_internal(%rip), %xmm9
+        subpd     %xmm4, %xmm12
+        andps     %xmm1, %xmm9
+
+/*
+ * prepare table index
+ * table lookup
+ */
+        movaps    %xmm12, %xmm10
+        subpd     %xmm9, %xmm1
+        mulpd     %xmm12, %xmm9
+        mulpd     %xmm12, %xmm1
+        subpd     One+__svml_dlog10_data_internal(%rip), %xmm9
+        addpd     %xmm9, %xmm1
+
+/* polynomial */
+        movups    poly_coeff+__svml_dlog10_data_internal(%rip), %xmm14
+        psrlq     $40, %xmm10
+        mulpd     %xmm1, %xmm14
+        movd      %xmm10, %eax
+        pshufd    $2, %xmm10, %xmm11
+        movaps    %xmm1, %xmm10
+        movups    poly_coeff+32+__svml_dlog10_data_internal(%rip), %xmm15
+        mulpd     %xmm1, %xmm10
+        addpd     poly_coeff+16+__svml_dlog10_data_internal(%rip), %xmm14
+        mulpd     %xmm1, %xmm15
+        mulpd     %xmm10, %xmm14
+        addpd     poly_coeff+48+__svml_dlog10_data_internal(%rip), %xmm15
+        movd      %xmm11, %ecx
+
+/* exponent*log(2.0) */
+        movups    Threshold+__svml_dlog10_data_internal(%rip), %xmm13
+        addpd     %xmm14, %xmm15
+        cmpltpd   %xmm12, %xmm13
+        mulpd     %xmm15, %xmm10
+        pshufd    $221, %xmm5, %xmm6
+        movups    poly_coeff+64+__svml_dlog10_data_internal(%rip), %xmm11
+
+/* biased exponent in DP format */
+        cvtdq2pd  %xmm6, %xmm3
+        mulpd     %xmm1, %xmm11
+        andps     Bias+__svml_dlog10_data_internal(%rip), %xmm13
+        orps      Bias1+__svml_dlog10_data_internal(%rip), %xmm13
+        subpd     %xmm13, %xmm3
+        addpd     %xmm10, %xmm11
+        mulpd     L2+__svml_dlog10_data_internal(%rip), %xmm3
+        movslq    %eax, %rax
+        movslq    %ecx, %rcx
+        movsd     (%rsi,%rax), %xmm2
+        movhpd    (%rsi,%rcx), %xmm2
+
+/* reconstruction */
+        addpd     %xmm11, %xmm2
+        addpd     %xmm2, %xmm3
+        testl     %edx, %edx
+
+/* Go to special inputs processing branch */
+        jne       L(SPECIAL_VALUES_BRANCH)
+                                # LOE rbx r12 r13 r14 r15 edx xmm0 xmm3
+
+/* Restore registers
+ * and exit the function
+ */
+
+L(EXIT):
+        movaps    %xmm3, %xmm0
+        movq      %rbp, %rsp
+        popq      %rbp
+        cfi_def_cfa(7, 8)
+        cfi_restore(6)
+        ret
+        cfi_def_cfa(6, 16)
+        cfi_offset(6, -16)
+
+/* Branch to process
+ * special inputs
+ */
+
+L(SPECIAL_VALUES_BRANCH):
+        movups    %xmm0, 32(%rsp)
+        movups    %xmm3, 48(%rsp)
+                                # LOE rbx r12 r13 r14 r15 edx
+
+        xorl      %eax, %eax
+        movq      %r12, 16(%rsp)
+        /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -48; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xd0, 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: -32; DW_OP_and; DW_OP_const4s: -56; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xc8, 0xff, 0xff, 0xff, 0x22
+        movl      %edx, %r13d
+        movq      %r14, (%rsp)
+        /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -64; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x22
+                                # LOE rbx r15 r12d r13d
+
+/* Range mask
+ * bits check
+ */
+
+L(RANGEMASK_CHECK):
+        btl       %r12d, %r13d
+
+/* Call scalar math function */
+        jc        L(SCALAR_MATH_CALL)
+                                # LOE rbx r15 r12d r13d
+
+/* Special inputs
+ * processing loop
+ */
+
+L(SPECIAL_VALUES_LOOP):
+        incl      %r12d
+        cmpl      $2, %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)
+        movups    48(%rsp), %xmm3
+
+/* Go to exit */
+        jmp       L(EXIT)
+        /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -48; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xd0, 0xff, 0xff, 0xff, 0x22
+        /*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -56; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xc8, 0xff, 0xff, 0xff, 0x22
+        /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -64; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xc0, 0xff, 0xff, 0xff, 0x22
+                                # LOE rbx r12 r13 r14 r15 xmm3
+
+/* Scalar math fucntion call
+ * to process special input
+ */
+
+L(SCALAR_MATH_CALL):
+        movl      %r12d, %r14d
+        movsd     32(%rsp,%r14,8), %xmm0
+        call      log10@PLT
+                                # LOE rbx r14 r15 r12d r13d xmm0
+
+        movsd     %xmm0, 48(%rsp,%r14,8)
+
+/* Process special inputs in loop */
+        jmp       L(SPECIAL_VALUES_LOOP)
+                                # LOE rbx r15 r12d r13d
+END(_ZGVbN2v_log10_sse4)
+
+        .section .rodata, "a"
+        .align 16
+
+#ifdef __svml_dlog10_data_internal_typedef
+typedef unsigned int VUINT32;
+typedef struct {
+        __declspec(align(16)) VUINT32 Log_HA_table[(1<<9)+2][2];
+        __declspec(align(16)) VUINT32 Log_LA_table[(1<<9)+1][2];
+        __declspec(align(16)) VUINT32 poly_coeff[5][2][2];
+        __declspec(align(16)) VUINT32 ExpMask[2][2];
+        __declspec(align(16)) VUINT32 Two10[2][2];
+        __declspec(align(16)) VUINT32 MinNorm[2][2];
+        __declspec(align(16)) VUINT32 MaxNorm[2][2];
+        __declspec(align(16)) VUINT32 HalfMask[2][2];
+        __declspec(align(16)) VUINT32 One[2][2];
+        __declspec(align(16)) VUINT32 Threshold[2][2];
+        __declspec(align(16)) VUINT32 Bias[2][2];
+        __declspec(align(16)) VUINT32 Bias1[2][2];
+        __declspec(align(16)) VUINT32 L2[2][2];
+} __svml_dlog10_data_internal;
+#endif
+__svml_dlog10_data_internal:
+        /* Log_HA_table */
+        .quad 0xc0733a7146f6b080, 0xbe1e707ce619c200
+        .quad 0xc0733a7547771970, 0xbe1e79c6c06d6f51
+        .quad 0xc0733a7945aacb70, 0xbe1e78e225fad29c
+        .quad 0xc0733a7d41946970, 0xbe1e76d607f9693b
+        .quad 0xc0733a813b3691f0, 0xbe1e7704b3e0685b
+        .quad 0xc0733a853293df00, 0xbe1e79c1216a27fa
+        .quad 0xc0733a8927aee660, 0xbe1e76dce5734a81
+        .quad 0xc0733a8d1a8a3920, 0xbe1e782ee2ca4dba
+        .quad 0xc0733a910b286430, 0xbe1e7812d1a0a61f
+        .quad 0xc0733a94f98bf010, 0xbe1e77e1b5ecbc61
+        .quad 0xc0733a98e5b76100, 0xbe1e76635cac1586
+        .quad 0xc0733a9ccfad36f0, 0xbe1e7638f7968f32
+        .quad 0xc0733aa0b76feda0, 0xbe1e7840ee76e365
+        .quad 0xc0733aa49d01fcb0, 0xbe1e79f3fd01907e
+        .quad 0xc0733aa88065d7a0, 0xbe1e77bbb3a9c38a
+        .quad 0xc0733aac619dedb0, 0xbe1e7742719bf41d
+        .quad 0xc0733ab040acaa20, 0xbe1e79bcedaf79cb
+        .quad 0xc0733ab41d947450, 0xbe1e762d63cb7ca0
+        .quad 0xc0733ab7f857af50, 0xbe1e77a07be83403
+        .quad 0xc0733abbd0f8ba80, 0xbe1e7763ff836ad0
+        .quad 0xc0733abfa779f130, 0xbe1e7737720ead39
+        .quad 0xc0733ac37bddaad0, 0xbe1e7776a08e55e7
+        .quad 0xc0733ac74e263af0, 0xbe1e793e3c52dd36
+        .quad 0xc0733acb1e55f160, 0xbe1e788a94695051
+        .quad 0xc0733aceec6f1a10, 0xbe1e76508114a813
+        .quad 0xc0733ad2b873fd20, 0xbe1e76909457d23e
+        .quad 0xc0733ad68266df10, 0xbe1e7664a24f9ca4
+        .quad 0xc0733ada4a4a0090, 0xbe1e7a07b3d44b18
+        .quad 0xc0733ade101f9ee0, 0xbe1e76d87594704d
+        .quad 0xc0733ae1d3e9f340, 0xbe1e79563595a182
+        .quad 0xc0733ae595ab33b0, 0xbe1e771880c3c6ab
+        .quad 0xc0733ae955659250, 0xbe1e78c171f517d4
+        .quad 0xc0733aed131b3df0, 0xbe1e77eac3874666
+        .quad 0xc0733af0cece61b0, 0xbe1e790db479d8f6
+        .quad 0xc0733af488812550, 0xbe1e7965d1aa5c90
+        .quad 0xc0733af84035ad10, 0xbe1e78ceb398ba47
+        .quad 0xc0733afbf5ee19c0, 0xbe1e779cc0dcb5aa
+        .quad 0xc0733affa9ac88c0, 0xbe1e7871053953ed
+        .quad 0xc0733b035b731420, 0xbe1e7a082cffa71a
+        .quad 0xc0733b070b43d2a0, 0xbe1e7904b4382fad
+        .quad 0xc0733b0ab920d790, 0xbe1e79b458d0b4f3
+        .quad 0xc0733b0e650c3310, 0xbe1e79d0ded414c6
+        .quad 0xc0733b120f07f200, 0xbe1e763c357a1943
+        .quad 0xc0733b15b7161dd0, 0xbe1e78b80ba6daaa
+        .quad 0xc0733b195d38bd00, 0xbe1e7998e23b8ffd
+        .quad 0xc0733b1d0171d2c0, 0xbe1e7974aa65ee8c
+        .quad 0xc0733b20a3c35f20, 0xbe1e76ccfde752ab
+        .quad 0xc0733b24442f5ef0, 0xbe1e77b4ff19debb
+        .quad 0xc0733b27e2b7cc10, 0xbe1e7772ee478542
+        .quad 0xc0733b2b7f5e9d30, 0xbe1e781d81b58b44
+        .quad 0xc0733b2f1a25c600, 0xbe1e78350d967565
+        .quad 0xc0733b32b30f3720, 0xbe1e783888e48152
+        .quad 0xc0733b364a1cde30, 0xbe1e78367bf7c111
+        .quad 0xc0733b39df50a5d0, 0xbe1e7959e57ca47d
+        .quad 0xc0733b3d72ac75c0, 0xbe1e777322423222
+        .quad 0xc0733b41043232b0, 0xbe1e767ce42a60aa
+        .quad 0xc0733b4493e3be70, 0xbe1e781d445aea19
+        .quad 0xc0733b4821c2f800, 0xbe1e7922fca18e18
+        .quad 0xc0733b4badd1bb80, 0xbe1e76fed3d40647
+        .quad 0xc0733b4f3811e210, 0xbe1e793948c9eabc
+        .quad 0xc0733b52c0854240, 0xbe1e76e487656b8c
+        .quad 0xc0733b56472daf90, 0xbe1e780ab2f71223
+        .quad 0xc0733b59cc0cfaf0, 0xbe1e77189120b09c
+        .quad 0xc0733b5d4f24f270, 0xbe1e7644a0343a12
+        .quad 0xc0733b60d0776160, 0xbe1e78f2a3e4733d
+        .quad 0xc0733b6450061080, 0xbe1e7913b2f73ae5
+        .quad 0xc0733b67cdd2c5c0, 0xbe1e7882d08393b5
+        .quad 0xc0733b6b49df4470, 0xbe1e765e1b209979
+        .quad 0xc0733b6ec42d4d20, 0xbe1e785c9c4620d4
+        .quad 0xc0733b75b394f240, 0xbe1e78878cd0e956
+        .quad 0xc0733b7c9c178630, 0xbe1e789a4112d90b
+        .quad 0xc0733b837dc2b0f0, 0xbe1e79050b8a1766
+        .quad 0xc0733b8a58a3f220, 0xbe1e7790dffc47aa
+        .quad 0xc0733b912cc8a180, 0xbe1e77174593b06a
+        .quad 0xc0733b97fa3defb0, 0xbe1e7677de2d2ecc
+        .quad 0xc0733b9ec110e6b0, 0xbe1e76cff477ca18
+        .quad 0xc0733ba5814e6a80, 0xbe1e78f8644dec7b
+        .quad 0xc0733bac3b0339d0, 0xbe1e764e1361788d
+        .quad 0xc0733bb2ee3bee30, 0xbe1e78c913e738de
+        .quad 0xc0733bb99b04fd30, 0xbe1e76666f5bddaa
+        .quad 0xc0733bc0416ab850, 0xbe1e77e87cbd8ab6
+        .quad 0xc0733bc6e1794e10, 0xbe1e76f18ba1c966
+        .quad 0xc0733bcd7b3cca10, 0xbe1e777c9461b8db
+        .quad 0xc0733bd40ec115d0, 0xbe1e78b78526ffac
+        .quad 0xc0733bda9c11f920, 0xbe1e7942abecfede
+        .quad 0xc0733be1233b1aa0, 0xbe1e76d8a684fd8c
+        .quad 0xc0733be7a4480010, 0xbe1e79622b539ac9
+        .quad 0xc0733bee1f440f30, 0xbe1e7978e7cc20ea
+        .quad 0xc0733bf4943a8de0, 0xbe1e765c9c9de825
+        .quad 0xc0733bfb0336a290, 0xbe1e775d8b138ee2
+        .quad 0xc0733c016c435500, 0xbe1e78bf33465c2f
+        .quad 0xc0733c07cf6b8e80, 0xbe1e78164f7cc441
+        .quad 0xc0733c0e2cba1a50, 0xbe1e7824e64d0b23
+        .quad 0xc0733c148439a630, 0xbe1e78373ae7dd81
+        .quad 0xc0733c1ad5f4c2c0, 0xbe1e7704513e0afe
+        .quad 0xc0733c2121f5e3d0, 0xbe1e7914aa84200f
+        .quad 0xc0733c2768476110, 0xbe1e76b1cde25cf6
+        .quad 0xc0733c2da8f37600, 0xbe1e796120e3862d
+        .quad 0xc0733c33e40442e0, 0xbe1e78ec836d7e7b
+        .quad 0xc0733c3a1983cca0, 0xbe1e77fb13b7dabb
+        .quad 0xc0733c40497bfd70, 0xbe1e783c6fcb2404
+        .quad 0xc0733c4673f6a530, 0xbe1e7628bb93dce8
+        .quad 0xc0733c4c98fd7990, 0xbe1e7857a47b5001
+        .quad 0xc0733c52b89a16d0, 0xbe1e76708dc2831f
+        .quad 0xc0733c58d2d5ffa0, 0xbe1e77b6038651f1
+        .quad 0xc0733c5ee7ba9de0, 0xbe1e792e855bb5b2
+        .quad 0xc0733c64f75142d0, 0xbe1e776cacd5c105
+        .quad 0xc0733c6b01a32740, 0xbe1e77f8a8011315
+        .quad 0xc0733c7106b96c30, 0xbe1e765cf3efcfde
+        .quad 0xc0733c77069d1ad0, 0xbe1e78d837d2efac
+        .quad 0xc0733c7d01572530, 0xbe1e78b615cf772c
+        .quad 0xc0733c82f6f06640, 0xbe1e7650bbbd7a25
+        .quad 0xc0733c88e771a220, 0xbe1e78bcf3495872
+        .quad 0xc0733c8ed2e386c0, 0xbe1e792266832e84
+        .quad 0xc0733c94b94eabd0, 0xbe1e79c1c3c2ca52
+        .quad 0xc0733c9a9abb9340, 0xbe1e78aa61e5807d
+        .quad 0xc0733ca07732a970, 0xbe1e7620fc4cf156
+        .quad 0xc0733ca64ebc4570, 0xbe1e76b914a832c5
+        .quad 0xc0733cac2160a970, 0xbe1e79227f72020e
+        .quad 0xc0733cb1ef280300, 0xbe1e77ac972cc008
+        .quad 0xc0733cb7b81a6b10, 0xbe1e798089be41f4
+        .quad 0xc0733cbd7c3fe6a0, 0xbe1e77942ae037fe
+        .quad 0xc0733cc33ba06690, 0xbe1e7956ae6463d9
+        .quad 0xc0733cc8f643c850, 0xbe1e7918a50c7942
+        .quad 0xc0733cceac31d5d0, 0xbe1e78308eeab604
+        .quad 0xc0733cd45d7245e0, 0xbe1e76dd4ea88445
+        .quad 0xc0733cda0a0cbc60, 0xbe1e77e7c1aa5909
+        .quad 0xc0733cdfb208caa0, 0xbe1e7804b9d20e54
+        .quad 0xc0733ce5556def70, 0xbe1e78f88e99d49c
+        .quad 0xc0733ceaf4439780, 0xbe1e787d74682d68
+        .quad 0xc0733cf08e911d80, 0xbe1e76edc24fe6e7
+        .quad 0xc0733cf6245dca50, 0xbe1e79b347ec86d2
+        .quad 0xc0733cfbb5b0d580, 0xbe1e797cceb2c39b
+        .quad 0xc0733d0142916530, 0xbe1e783adbdc6aa1
+        .quad 0xc0733d06cb068e70, 0xbe1e76e4c20e3d9e
+        .quad 0xc0733d0c4f175570, 0xbe1e77070bf3cf61
+        .quad 0xc0733d11cecaadc0, 0xbe1e781c43502734
+        .quad 0xc0733d174a277a80, 0xbe1e78b11268ea72
+        .quad 0xc0733d1cc1348e90, 0xbe1e7754b83bfc7d
+        .quad 0xc0733d2233f8acb0, 0xbe1e7756c29bf5e9
+        .quad 0xc0733d27a27a87d0, 0xbe1e7952fc1d9333
+        .quad 0xc0733d2d0cc0c350, 0xbe1e778c76ae6077
+        .quad 0xc0733d3272d1f2e0, 0xbe1e7a1896ba8f43
+        .quad 0xc0733d37d4b49b30, 0xbe1e76dafdf432d8
+        .quad 0xc0733d3d326f3180, 0xbe1e795330184013
+        .quad 0xc0733d428c081c80, 0xbe1e763cc774d30f
+        .quad 0xc0733d47e185b3d0, 0xbe1e77030a779c0a
+        .quad 0xc0733d4d32ee40b0, 0xbe1e7908af2a2d7e
+        .quad 0xc0733d528047fe00, 0xbe1e78c4953b797d
+        .quad 0xc0733d57c9991850, 0xbe1e78b43b096579
+        .quad 0xc0733d5d0ee7ae30, 0xbe1e7824ae0a4804
+        .quad 0xc0733d625039d040, 0xbe1e79d2b2fbb740
+        .quad 0xc0733d678d958190, 0xbe1e7662de59a1a6
+        .quad 0xc0733d6cc700b760, 0xbe1e76b251d59aaa
+        .quad 0xc0733d71fc8159b0, 0xbe1e7a00cfd1f487
+        .quad 0xc0733d772e1d4360, 0xbe1e77f4d246167e
+        .quad 0xc0733d7c5bda4200, 0xbe1e767a4ee8e6fc
+        .quad 0xc0733d8185be1640, 0xbe1e777ccf0a8aed
+        .quad 0xc0733d86abce7420, 0xbe1e767d7e279ada
+        .quad 0xc0733d8bce1102d0, 0xbe1e7a05cef4bb90
+        .quad 0xc0733d90ec8b5d40, 0xbe1e78f75369be5b
+        .quad 0xc0733d96074311d0, 0xbe1e77b9612e8c8a
+        .quad 0xc0733d9b1e3da2b0, 0xbe1e794518b9adeb
+        .quad 0xc0733da031808620, 0xbe1e7810626fb934
+        .quad 0xc0733da541112650, 0xbe1e76d87223fa6d
+        .quad 0xc0733daa4cf4e1a0, 0xbe1e794c5e7ca3b5
+        .quad 0xc0733daf55310af0, 0xbe1e789856ef816f
+        .quad 0xc0733db459cae970, 0xbe1e77d2004effbd
+        .quad 0xc0733db95ac7b8f0, 0xbe1e78467d31eb9c
+        .quad 0xc0733dbe582caa00, 0xbe1e79aaa4e25787
+        .quad 0xc0733dc351fee220, 0xbe1e762de8f107bf
+        .quad 0xc0733dc848437b90, 0xbe1e7670670a63fe
+        .quad 0xc0733dcd3aff85d0, 0xbe1e795ca237c6cc
+        .quad 0xc0733dd22a3805b0, 0xbe1e77e55c53c1d9
+        .quad 0xc0733dd715f1f520, 0xbe1e78a806213ac4
+        .quad 0xc0733ddbfe3243b0, 0xbe1e77743a2bc615
+        .quad 0xc0733de0e2fdd660, 0xbe1e78b8b45b0b7d
+        .quad 0xc0733de5c4598800, 0xbe1e78d635f2f4b9
+        .quad 0xc0733deaa24a2920, 0xbe1e7758c396a11e
+        .quad 0xc0733def7cd48020, 0xbe1e7a17a8cc454c
+        .quad 0xc0733df453fd49a0, 0xbe1e783caa73f616
+        .quad 0xc0733df927c93820, 0xbe1e7932cfa29664
+        .quad 0xc0733dfdf83cf490, 0xbe1e777d265c72a6
+        .quad 0xc0733e02c55d1e10, 0xbe1e7775e7c03c60
+        .quad 0xc0733e078f2e4a40, 0xbe1e79f65d52d232
+        .quad 0xc0733e0c55b50570, 0xbe1e76e7e7464b4e
+        .quad 0xc0733e1118f5d250, 0xbe1e77be81cad877
+        .quad 0xc0733e15d8f52a80, 0xbe1e79dd25b5fb3a
+        .quad 0xc0733e1a95b77e80, 0xbe1e78e45f1418ef
+        .quad 0xc0733e1f4f4135a0, 0xbe1e78eb7289505b
+        .quad 0xc0733e240596ae50, 0xbe1e78a468c07cad
+        .quad 0xc0733e28b8bc3e20, 0xbe1e776b558a4009
+        .quad 0xc0733e2d68b631d0, 0xbe1e77412eb9941e
+        .quad 0xc0733e321588cd80, 0xbe1e76b2853f845e
+        .quad 0xc0733e36bf384cb0, 0xbe1e76aa7184273c
+        .quad 0xc0733e3b65c8e260, 0xbe1e7832027f78fa
+        .quad 0xc0733e40093eb930, 0xbe1e7a1c7da131f5
+        .quad 0xc0733e44a99df380, 0xbe1e76a0bc2ae4bc
+        .quad 0xc0733e4946eaab30, 0xbe1e78dff13b6f5d
+        .quad 0xc0733e4de128f250, 0xbe1e765a226dea2c
+        .quad 0xc0733e52785cd290, 0xbe1e78509b989111
+        .quad 0xc0733e570c8a4de0, 0xbe1e7916a4e9803d
+        .quad 0xc0733e5b9db55e30, 0xbe1e7950c15758cc
+        .quad 0xc0733e602be1f5a0, 0xbe1e7922ba1ad420
+        .quad 0xc0733e64b713fe90, 0xbe1e794cbaabcef6
+        .quad 0xc0733e693f4f5bc0, 0xbe1e7837bf883fed
+        .quad 0xc0733e6dc497e850, 0xbe1e76f198ddbbdf
+        .quad 0xc0733e7246f177d0, 0xbe1e7a18c1067764
+        .quad 0xc0733e76c65fd6a0, 0xbe1e76b845a8fd9d
+        .quad 0xc0733e7b42e6c970, 0xbe1e7714012df506
+        .quad 0xc0733e7fbc8a0de0, 0xbe1e7765612922cd
+        .quad 0xc0733e84334d5a50, 0xbe1e7688f5424a00
+        .quad 0xc0733e88a7345df0, 0xbe1e769d011f6663
+        .quad 0xc0733e8d1842c0e0, 0xbe1e79914acbfaf7
+        .quad 0xc0733e91867c2460, 0xbe1e79a85e189bd7
+        .quad 0xc0733e95f1e422a0, 0xbe1e79ea7c726432
+        .quad 0xc0733e9a5a7e4f10, 0xbe1e768a6fbb8e6e
+        .quad 0xc0733e9ec04e3620, 0xbe1e793c75bcc9fc
+        .quad 0xc0733ea323575dd0, 0xbe1e797f78da13d4
+        .quad 0xc0733ea7839d4550, 0xbe1e78d8c9cda978
+        .quad 0xc0733eabe1236540, 0xbe1e77028d480fff
+        .quad 0xc0733eb03bed2fa0, 0xbe1e7a0d0f74ff7c
+        .quad 0xc0733eb493fe1040, 0xbe1e76732e8a35fb
+        .quad 0xc0733eb8e9596c30, 0xbe1e77220caeabeb
+        .quad 0xc0733ebd3c02a260, 0xbe1e797438b645ef
+        .quad 0xc0733ec18bfd0b80, 0xbe1e79207c5fd6e8
+        .quad 0xc0733ec5d94bf9f0, 0xbe1e781c7df8f946
+        .quad 0xc0733eca23f2b9f0, 0xbe1e76736284e2db
+        .quad 0xc0733ece6bf49190, 0xbe1e7a109cc0c3f5
+        .quad 0xc0733ed2b154c120, 0xbe1e767f14a16d50
+        .quad 0xc0733ed6f4168290, 0xbe1e789cd22acaf0
+        .quad 0xc0733edb343d0a40, 0xbe1e764355ca28ad
+        .quad 0xc0733edf71cb8660, 0xbe1e79e4c7a81c45
+        .quad 0xc0733ee3acc51fb0, 0xbe1e761e26b644c2
+        .quad 0xc0733ee7e52cf8c0, 0xbe1e793e9f8fbdd3
+        .quad 0xc0733eec1b062ed0, 0xbe1e78c432991c20
+        .quad 0xc0733ef04e53d940, 0xbe1e78cdd025f4d8
+        .quad 0xc0733ef47f1909f0, 0xbe1e778310c6446e
+        .quad 0xc0733ef8ad58cd20, 0xbe1e7871af3d6e17
+        .quad 0xc0733efcd91629b0, 0xbe1e77e0e906f697
+        .quad 0xc0733f01025420f0, 0xbe1e7a1ae9b27892
+        .quad 0xc0733f052915af00, 0xbe1e76ac64c88f9d
+        .quad 0xc0733f094d5dca60, 0xbe1e779a815589c4
+        .quad 0xc0733f0d6f2f6480, 0xbe1e788f39a4864c
+        .quad 0xc0733f118e8d6980, 0xbe1e79fc51263525
+        .quad 0xc0733f15ab7ac060, 0xbe1e783501f19e90
+        .quad 0xc0733f19c5fa4ae0, 0xbe1e767e82c327ab
+        .quad 0xc0733f1dde0ee5a0, 0xbe1e7a1785d66123
+        .quad 0xc0733f21f3bb6870, 0xbe1e7936d07203da
+        .quad 0xc0733f260702a5e0, 0xbe1e7a010a7ac699
+        .quad 0xc0733f2a17e76bb0, 0xbe1e7975e4e16312
+        .quad 0xc0733f2e266c82b0, 0xbe1e7654b5422330
+        .quad 0xc0733f323294aeb0, 0xbe1e77f8a4909d35
+        .quad 0xc0733f363c62aee0, 0xbe1e792c8e30d226
+        .quad 0xc0733f3a43d93da0, 0xbe1e76f6ac67a1ff
+        .quad 0xc0733f3e48fb1070, 0xbe1e775c2e97715a
+        .quad 0xc0733f424bcad840, 0xbe1e781cd54ae100
+        /*== Log_LA_table ==*/
+        .align 16
+        .quad 0x0000000000000000
+        .quad 0xbf4bc48a867884b7
+        .quad 0xbf5bbd9e9482af09
+        .quad 0xbf64c9096b94befd
+        .quad 0xbf6bafd47221ed26
+        .quad 0xbf714999e2ad8ea6
+        .quad 0xbf74b99563d2a1bd
+        .quad 0xbf7827de6b310350
+        .quad 0xbf7b9476a4fcd10f
+        .quad 0xbf7eff5fbaf25781
+        .quad 0xbf81344daa2d7553
+        .quad 0xbf82e8158b08d957
+        .quad 0xbf849b0851443684
+        .quad 0xbf864d26cce610dd
+        .quad 0xbf87fe71ccc4e6b0
+        .quad 0xbf89aeea1e897fdf
+        .quad 0xbf8b5e908eb13790
+        .quad 0xbf8d0d65e890405a
+        .quad 0xbf8ebb6af653e2ee
+        .quad 0xbf90345040825bad
+        .quad 0xbf910a83a8446c78
+        .quad 0xbf91e05015d30a71
+        .quad 0xbf92b5b5ec0209d3
+        .quad 0xbf938ab58d173e91
+        .quad 0xbf945f4f5acb8be0
+        .quad 0xbf953383b64bf13f
+        .quad 0xbf960753003a94ef
+        .quad 0xbf96dabd98afcc05
+        .quad 0xbf97adc3df3b1ff8
+        .quad 0xbf98806632e451d0
+        .quad 0xbf9952a4f22c5ae9
+        .quad 0xbf9a24807b0e6b5c
+        .quad 0xbf9af5f92b00e610
+        .quad 0xbf9bc70f5ef65a77
+        .quad 0xbf9c97c3735e7c0a
+        .quad 0xbf9d6815c4271775
+        .quad 0xbf9e3806acbd058f
+        .quad 0xbf9f0796880d1c19
+        .quad 0xbf9fd6c5b0851c4c
+        .quad 0xbfa052ca400a4f9b
+        .quad 0xbfa0ba01a8170000
+        .quad 0xbfa121093ce3a205
+        .quad 0xbfa187e12aad8077
+        .quad 0xbfa1ee899d74a03e
+        .quad 0xbfa25502c0fc314c
+        .quad 0xbfa2bb4cc0cafe8d
+        .quad 0xbfa32167c82bdcda
+        .quad 0xbfa38754022e18e2
+        .quad 0xbfa3ed1199a5e425
+        .quad 0xbfa452a0b92cc0ec
+        .quad 0xbfa4b8018b21ed4f
+        .quad 0xbfa51d3439aacd4a
+        .quad 0xbfa58238eeb353da
+        .quad 0xbfa5e70fd3ee6b34
+        .quad 0xbfa64bb912d65c07
+        .quad 0xbfa6b034d4ad33df
+        .quad 0xbfa71483427d2a99
+        .quad 0xbfa778a4851906f3
+        .quad 0xbfa7dc98c51c8242
+        .quad 0xbfa840602aecab3d
+        .quad 0xbfa8a3fadeb847f4
+        .quad 0xbfa90769087836e4
+        .quad 0xbfa96aaacfefcf3c
+        .quad 0xbfa9cdc05cad4042
+        .quad 0xbfaa30a9d609efea
+        .quad 0xbfaa9367632ad897
+        .quad 0xbfaaf5f92b00e610
+        .quad 0xbfab585f544951a4
+        .quad 0xbfabba9a058dfd84
+        .quad 0xbfac1ca96525cf56
+        .quad 0xbfac7e8d993509f9
+        .quad 0xbface046c7ada68d
+        .quad 0xbfad41d5164facb4
+        .quad 0xbfada338aaa98a0c
+        .quad 0xbfae0471aa1868f5
+        .quad 0xbfae658039c88690
+        .quad 0xbfaec6647eb58808
+        .quad 0xbfaf271e9daacf20
+        .quad 0xbfaf87aebb43ce06
+        .quad 0xbfafe814fbec5a77
+        .quad 0xbfb02428c1f08016
+        .quad 0xbfb054323b97a948
+        .quad 0xbfb08426fcdb1ee7
+        .quad 0xbfb0b40717932b96
+        .quad 0xbfb0e3d29d81165e
+        .quad 0xbfb11389a04f4a2e
+        .quad 0xbfb1432c31917d08
+        .quad 0xbfb172ba62c4d6de
+        .quad 0xbfb1a23445501816
+        .quad 0xbfb1d199ea83bfbe
+        .quad 0xbfb200eb639a3173
+        .quad 0xbfb23028c1b7daed
+        .quad 0xbfb25f5215eb594a
+        .quad 0xbfb28e67712d9dfc
+        .quad 0xbfb2bd68e4621371
+        .quad 0xbfb2ec568056c16f
+        .quad 0xbfb31b3055c47118
+        .quad 0xbfb349f6754ed0b4
+        .quad 0xbfb378a8ef84971e
+        .quad 0xbfb3a747d4dfa6f5
+        .quad 0xbfb3d5d335c53179
+        .quad 0xbfb4044b2285d925
+        .quad 0xbfb432afab5dd3ff
+        .quad 0xbfb46100e0750da1
+        .quad 0xbfb48f3ed1df48fb
+        .quad 0xbfb4bd698f9c41cf
+        .quad 0xbfb4eb812997cde4
+        .quad 0xbfb51985afa9fdfd
+        .quad 0xbfb5477731973e85
+        .quad 0xbfb57555bf1077f5
+        .quad 0xbfb5a32167b32f02
+        .quad 0xbfb5d0da3b09a47e
+        .quad 0xbfb5fe80488af4fd
+        .quad 0xbfb62c139f9b3837
+        .quad 0xbfb659944f8ba02d
+        .quad 0xbfb68702679a980a
+        .quad 0xbfb6b45df6f3e2c9
+        .quad 0xbfb6e1a70cb0b99a
+        .quad 0xbfb70eddb7d7ea07
+        .quad 0xbfb73c02075df3e5
+        .quad 0xbfb769140a2526fd
+        .quad 0xbfb79613cefdc07d
+        .quad 0xbfb7c30164a60836
+        .quad 0xbfb7efdcd9ca6d8f
+        .quad 0xbfb81ca63d05a44a
+        .quad 0xbfb8495d9ce0c10c
+        .quad 0xbfb8760307d355ab
+        .quad 0xbfb8a2968c438d41
+        .quad 0xbfb8cf183886480d
+        .quad 0xbfb8fb881adf3713
+        .quad 0xbfb927e64180f790
+        .quad 0xbfb95432ba8d2e2f
+        .quad 0xbfb9806d9414a209
+        .quad 0xbfb9ac96dc175776
+        .quad 0xbfb9d8aea084aa9c
+        .quad 0xbfba04b4ef3b69d8
+        .quad 0xbfba30a9d609efea
+        .quad 0xbfba5c8d62ae3dec
+        .quad 0xbfba885fa2d6151e
+        .quad 0xbfbab420a41f1076
+        .quad 0xbfbadfd07416be07
+        .quad 0xbfbb0b6f203ab82c
+        .quad 0xbfbb36fcb5f8be8a
+        .quad 0xbfbb627942aecedd
+        .quad 0xbfbb8de4d3ab3d98
+        .quad 0xbfbbb93f762cce4f
+        .quad 0xbfbbe4893762cbf7
+        .quad 0xbfbc0fc2246d20f5
+        .quad 0xbfbc3aea4a5c6eff
+        .quad 0xbfbc6601b63226cb
+        .quad 0xbfbc910874e09f98
+        .quad 0xbfbcbbfe934b2e81
+        .quad 0xbfbce6e41e463da5
+        .quad 0xbfbd11b92297632b
+        .quad 0xbfbd3c7dacf5780b
+        .quad 0xbfbd6731ca08aeb9
+        .quad 0xbfbd91d5866aa99c
+        .quad 0xbfbdbc68eea6915b
+        .quad 0xbfbde6ec0f392b05
+        .quad 0xbfbe115ef490ee07
+        .quad 0xbfbe3bc1ab0e19fe
+        .quad 0xbfbe66143f02cc5d
+        .quad 0xbfbe9056bcb315e8
+        .quad 0xbfbeba893055100b
+        .quad 0xbfbee4aba610f204
+        .quad 0xbfbf0ebe2a0125eb
+        .quad 0xbfbf38c0c8325d86
+        .quad 0xbfbf62b38ca3a706
+        .quad 0xbfbf8c9683468191
+        .quad 0xbfbfb669b7fef1a8
+        .quad 0xbfbfe02d36a3956d
+        .quad 0xbfc004f0857edc5c
+        .quad 0xbfc019c2a064b486
+        .quad 0xbfc02e8cf1dac4b8
+        .quad 0xbfc0434f7fb1f307
+        .quad 0xbfc0580a4fb4a3df
+        .quad 0xbfc06cbd67a6c3b6
+        .quad 0xbfc08168cd45d0a9
+        .quad 0xbfc0960c8648e406
+        .quad 0xbfc0aaa89860bbcf
+        .quad 0xbfc0bf3d0937c41c
+        .quad 0xbfc0d3c9de722078
+        .quad 0xbfc0e84f1dadb526
+        .quad 0xbfc0fccccc823059
+        .quad 0xbfc11142f0811357
+        .quad 0xbfc125b18f35bb8e
+        .quad 0xbfc13a18ae256b99
+        .quad 0xbfc14e7852cf5430
+        .quad 0xbfc162d082ac9d10
+        .quad 0xbfc1772143306dc6
+        .quad 0xbfc18b6a99c7f679
+        .quad 0xbfc19fac8bda7897
+        .quad 0xbfc1b3e71ec94f7b
+        .quad 0xbfc1c81a57eff8fd
+        .quad 0xbfc1dc463ca41df8
+        .quad 0xbfc1f06ad2359abd
+        .quad 0xbfc204881dee8777
+        .quad 0xbfc2189e25134081
+        .quad 0xbfc22cacece26ead
+        .quad 0xbfc240b47a950f79
+        .quad 0xbfc254b4d35e7d3c
+        .quad 0xbfc268adfc6c773e
+        .quad 0xbfc27c9ffae729c1
+        .quad 0xbfc2908ad3f13603
+        .quad 0xbfc2a46e8ca7ba2a
+        .quad 0xbfc2b84b2a225923
+        .quad 0xbfc2cc20b1734279
+        .quad 0xbfc2dfef27a73a18
+        .quad 0xbfc2f3b691c5a001
+        .quad 0xbfc30776f4d077f7
+        .quad 0xbfc31b3055c47118
+        .quad 0xbfc32ee2b998ed6e
+        .quad 0xbfc3428e2540096d
+        .quad 0x3fc331f403985097
+        .quad 0x3fc31e56798a910a
+        .quad 0x3fc30abfd8f333b6
+        .quad 0x3fc2f7301cf4e87b
+        .quad 0x3fc2e3a740b7800f
+        .quad 0x3fc2d0253f67e4cb
+        .quad 0x3fc2bcaa14381386
+        .quad 0x3fc2a935ba5f1479
+        .quad 0x3fc295c82d18f434
+        .quad 0x3fc2826167a6bc9c
+        .quad 0x3fc26f01654e6df6
+        .quad 0x3fc25ba8215af7fc
+        .quad 0x3fc24855971c3307
+        .quad 0x3fc23509c1e6d937
+        .quad 0x3fc221c49d147fb3
+        .quad 0x3fc20e8624038fed
+        .quad 0x3fc1fb4e521740f4
+        .quad 0x3fc1e81d22b790d4
+        .quad 0x3fc1d4f291513e01
+        .quad 0x3fc1c1ce9955c0c6
+        .quad 0x3fc1aeb1363b44c8
+        .quad 0x3fc19b9a637ca295
+        .quad 0x3fc1888a1c995931
+        .quad 0x3fc175805d1587c1
+        .quad 0x3fc1627d2079e731
+        .quad 0x3fc14f806253c3ed
+        .quad 0x3fc13c8a1e34f7a0
+        .quad 0x3fc1299a4fb3e306
+        .quad 0x3fc116b0f26b67bb
+        .quad 0x3fc103ce01fae223
+        .quad 0x3fc0f0f17a062353
+        .quad 0x3fc0de1b56356b04
+        .quad 0x3fc0cb4b9235619a
+        .quad 0x3fc0b88229b71227
+        .quad 0x3fc0a5bf186fe483
+        .quad 0x3fc093025a19976c
+        .quad 0x3fc0804bea723aa9
+        .quad 0x3fc06d9bc53c2941
+        .quad 0x3fc05af1e63e03b4
+        .quad 0x3fc0484e4942aa43
+        .quad 0x3fc035b0ea19373b
+        .quad 0x3fc02319c494f951
+        .quad 0x3fc01088d48d6e03
+        .quad 0x3fbffbfc2bbc7803
+        .quad 0x3fbfd6f308ce5b52
+        .quad 0x3fbfb1f6381856f4
+        .quad 0x3fbf8d05b16a6d47
+        .quad 0x3fbf68216c9cc727
+        .quad 0x3fbf4349618fa91a
+        .quad 0x3fbf1e7d882b689a
+        .quad 0x3fbef9bdd860616b
+        .quad 0x3fbed50a4a26eafc
+        .quad 0x3fbeb062d57f4de8
+        .quad 0x3fbe8bc77271b97a
+        .quad 0x3fbe6738190e394c
+        .quad 0x3fbe42b4c16caaf3
+        .quad 0x3fbe1e3d63acb3ba
+        .quad 0x3fbdf9d1f7f5b674
+        .quad 0x3fbdd5727676c959
+        .quad 0x3fbdb11ed766abf4
+        .quad 0x3fbd8cd71303bd26
+        .quad 0x3fbd689b2193f133
+        .quad 0x3fbd446afb64c7e5
+        .quad 0x3fbd204698cb42bd
+        .quad 0x3fbcfc2df223db2d
+        .quad 0x3fbcd820ffd278f3
+        .quad 0x3fbcb41fba42686d
+        .quad 0x3fbc902a19e65111
+        .quad 0x3fbc6c4017382bea
+        .quad 0x3fbc4861aab93a23
+        .quad 0x3fbc248eccf1fba6
+        .quad 0x3fbc00c7767225cb
+        .quad 0x3fbbdd0b9fd09a10
+        .quad 0x3fbbb95b41ab5ce6
+        .quad 0x3fbb95b654a78c87
+        .quad 0x3fbb721cd17157e3
+        .quad 0x3fbb4e8eb0bbf58f
+        .quad 0x3fbb2b0beb419ad0
+        .quad 0x3fbb079479c372ad
+        .quad 0x3fbae4285509950b
+        .quad 0x3fbac0c775e2fde6
+        .quad 0x3fba9d71d5258484
+        .quad 0x3fba7a276badd2c8
+        .quad 0x3fba56e8325f5c87
+        .quad 0x3fba33b4222456f1
+        .quad 0x3fba108b33edb005
+        .quad 0x3fb9ed6d60b30612
+        .quad 0x3fb9ca5aa1729f45
+        .quad 0x3fb9a752ef316149
+        .quad 0x3fb9845642fac8f0
+        .quad 0x3fb9616495e0e1e8
+        .quad 0x3fb93e7de0fc3e80
+        .quad 0x3fb91ba21d6bef77
+        .quad 0x3fb8f8d144557bdf
+        .quad 0x3fb8d60b4ee4d901
+        .quad 0x3fb8b350364c6257
+        .quad 0x3fb8909ff3c4d191
+        .quad 0x3fb86dfa808d36a0
+        .quad 0x3fb84b5fd5eaefd8
+        .quad 0x3fb828cfed29a215
+        .quad 0x3fb8064abf9b30f1
+        .quad 0x3fb7e3d04697b704
+        .quad 0x3fb7c1607b7d7e32
+        .quad 0x3fb79efb57b0f803
+        .quad 0x3fb77ca0d49cb608
+        .quad 0x3fb75a50ebb1624a
+        .quad 0x3fb7380b9665b7c8
+        .quad 0x3fb715d0ce367afc
+        .quad 0x3fb6f3a08ca67270
+        .quad 0x3fb6d17acb3e5f5e
+        .quad 0x3fb6af5f838cf654
+        .quad 0x3fb68d4eaf26d7ee
+        .quad 0x3fb66b4847a68997
+        .quad 0x3fb6494c46ac6e4d
+        .quad 0x3fb6275aa5debf81
+        .quad 0x3fb605735ee985f1
+        .quad 0x3fb5e3966b7e9295
+        .quad 0x3fb5c1c3c5557799
+        .quad 0x3fb59ffb662b815c
+        .quad 0x3fb57e3d47c3af7b
+        .quad 0x3fb55c8963e6adeb
+        .quad 0x3fb53adfb462ce16
+        .quad 0x3fb51940330c000b
+        .quad 0x3fb4f7aad9bbcbaf
+        .quad 0x3fb4d61fa2514a00
+        .quad 0x3fb4b49e86b11e5f
+        .quad 0x3fb4932780c56fe2
+        .quad 0x3fb471ba8a7de2b7
+        .quad 0x3fb450579dcf9186
+        .quad 0x3fb42efeb4b506e9
+        .quad 0x3fb40dafc92e36e2
+        .quad 0x3fb3ec6ad5407868
+        .quad 0x3fb3cb2fd2f67ef1
+        .quad 0x3fb3a9febc60540a
+        .quad 0x3fb388d78b9350ff
+        .quad 0x3fb367ba3aaa1883
+        .quad 0x3fb346a6c3c49066
+        .quad 0x3fb3259d2107db54
+        .quad 0x3fb3049d4c9e52a0
+        .quad 0x3fb2e3a740b7800f
+        .quad 0x3fb2c2baf78817b7
+        .quad 0x3fb2a1d86b49f1e2
+        .quad 0x3fb280ff963c04fc
+        .quad 0x3fb2603072a25f82
+        .quad 0x3fb23f6afac6220a
+        .quad 0x3fb21eaf28f57941
+        .quad 0x3fb1fdfcf7839804
+        .quad 0x3fb1dd5460c8b16f
+        .quad 0x3fb1bcb55f21f307
+        .quad 0x3fb19c1fecf17ee0
+        .quad 0x3fb17b94049e65d0
+        .quad 0x3fb15b11a094a1aa
+        .quad 0x3fb13a98bb450f81
+        .quad 0x3fb11a294f2569f6
+        .quad 0x3fb0f9c356b04389
+        .quad 0x3fb0d966cc6500fa
+        .quad 0x3fb0b913aac7d3a7
+        .quad 0x3fb098c9ec61b3ff
+        .quad 0x3fb078898bc05bf4
+        .quad 0x3fb0585283764178
+        .quad 0x3fb03824ce1a9101
+        .quad 0x3fb0180066492817
+        .quad 0x3fafefca8d451fd6
+        .quad 0x3fafafa6d397efdb
+        .quad 0x3faf6f9594de60f0
+        .quad 0x3faf2f96c6754aee
+        .quad 0x3faeefaa5dc2b239
+        .quad 0x3faeafd05035bd3b
+        .quad 0x3fae70089346a9e6
+        .quad 0x3fae30531c76c34a
+        .quad 0x3fadf0afe1505738
+        .quad 0x3fadb11ed766abf4
+        .quad 0x3fad719ff455f5f7
+        .quad 0x3fad32332dc34dbd
+        .quad 0x3facf2d8795ca5a5
+        .quad 0x3facb38fccd8bfdb
+        .quad 0x3fac74591df72456
+        .quad 0x3fac3534628016dd
+        .quad 0x3fabf62190448d22
+        .quad 0x3fabb7209d1e24e5
+        .quad 0x3fab78317eef1a29
+        .quad 0x3fab39542ba23d73
+        .quad 0x3faafa88992aea19
+        .quad 0x3faabbcebd84fca0
+        .quad 0x3faa7d268eb4c924
+        .quad 0x3faa3e9002c711d2
+        .quad 0x3faa000b0fd0fd6b
+        .quad 0x3fa9c197abf00dd7
+        .quad 0x3fa98335cd4a16c3
+        .quad 0x3fa944e56a0d3450
+        .quad 0x3fa906a6786fc1cb
+        .quad 0x3fa8c878eeb05074
+        .quad 0x3fa88a5cc3159e53
+        .quad 0x3fa84c51ebee8d15
+        .quad 0x3fa80e585f9218fc
+        .quad 0x3fa7d070145f4fd7
+        .quad 0x3fa7929900bd4809
+        .quad 0x3fa754d31b1b179c
+        .quad 0x3fa7171e59efcb5f
+        .quad 0x3fa6d97ab3ba5e10
+        .quad 0x3fa69be81f01af99
+        .quad 0x3fa65e6692547c4e
+        .quad 0x3fa620f604495440
+        .quad 0x3fa5e3966b7e9295
+        .quad 0x3fa5a647be9a54f6
+        .quad 0x3fa56909f44a72fe
+        .quad 0x3fa52bdd034475b8
+        .quad 0x3fa4eec0e2458f30
+        .quad 0x3fa4b1b588129203
+        .quad 0x3fa474baeb77e904
+        .quad 0x3fa437d103498eec
+        .quad 0x3fa3faf7c663060e
+        .quad 0x3fa3be2f2ba7501f
+        .quad 0x3fa381772a00e604
+        .quad 0x3fa344cfb861afae
+        .quad 0x3fa30838cdc2fbfd
+        .quad 0x3fa2cbb2612578b4
+        .quad 0x3fa28f3c69912a74
+        .quad 0x3fa252d6de1564c1
+        .quad 0x3fa21681b5c8c213
+        .quad 0x3fa1da3ce7c91bf8
+        .quad 0x3fa19e086b3b8333
+        .quad 0x3fa161e4374c37f4
+        .quad 0x3fa125d0432ea20e
+        .quad 0x3fa0e9cc861d4944
+        .quad 0x3fa0add8f759cd95
+        .quad 0x3fa071f58e2cdf9b
+        .quad 0x3fa0362241e638ec
+        .quad 0x3f9ff4be13b92920
+        .quad 0x3f9f7d57badb4ee8
+        .quad 0x3f9f061167fc31e8
+        .quad 0x3f9e8eeb09f2f6cb
+        .quad 0x3f9e17e48fa48962
+        .quad 0x3f9da0fde8038de9
+        .quad 0x3f9d2a3702105259
+        .quad 0x3f9cb38fccd8bfdb
+        .quad 0x3f9c3d0837784c41
+        .quad 0x3f9bc6a03117eb97
+        .quad 0x3f9b5057a8ee01ce
+        .quad 0x3f9ada2e8e3e546f
+        .quad 0x3f9a6424d059fc68
+        .quad 0x3f99ee3a5e9f57e8
+        .quad 0x3f99786f2879fc53
+        .quad 0x3f9902c31d62a843
+        .quad 0x3f988d362cdf359e
+        .quad 0x3f9817c846828bbd
+        .quad 0x3f97a27959ec91aa
+        .quad 0x3f972d4956ca2067
+        .quad 0x3f96b8382cd4f551
+        .quad 0x3f964345cbd3a491
+        .quad 0x3f95ce7223998b98
+        .quad 0x3f9559bd2406c3ba
+        .quad 0x3f94e526bd0814d1
+        .quad 0x3f9470aede96e7f2
+        .quad 0x3f93fc5578b93a38
+        .quad 0x3f93881a7b818f9e
+        .quad 0x3f9313fdd70ee5e8
+        .quad 0x3f929fff7b8ca79d
+        .quad 0x3f922c1f59329f1b
+        .quad 0x3f91b85d6044e9ae
+        .quad 0x3f9144b98113eac0
+        .quad 0x3f90d133abfc3f1b
+        .quad 0x3f905dcbd166b033
+        .quad 0x3f8fd503c3904f1d
+        .quad 0x3f8eeeab9b43445d
+        .quad 0x3f8e088f0b004827
+        .quad 0x3f8d22adf3f9579d
+        .quad 0x3f8c3d0837784c41
+        .quad 0x3f8b579db6dec358
+        .quad 0x3f8a726e53a6056e
+        .quad 0x3f898d79ef5eedf0
+        .quad 0x3f88a8c06bb1d2f4
+        .quad 0x3f87c441aa5e6d15
+        .quad 0x3f86dffd8d3bbf70
+        .quad 0x3f85fbf3f637ffc5
+        .quad 0x3f851824c7587eb0
+        .quad 0x3f84348fe2b99002
+        .quad 0x3f8351352a8e733f
+        .quad 0x3f826e1481213c2e
+        .quad 0x3f818b2dc8d2bb91
+        .quad 0x3f80a880e41a67f6
+        .quad 0x3f7f8c1b6b0c8d4e
+        .quad 0x3f7dc7a83f75a96d
+        .quad 0x3f7c03a80ae5e054
+        .quad 0x3f7a401a92ff827e
+        .quad 0x3f787cff9d9147a5
+        .quad 0x3f76ba56f09621bc
+        .quad 0x3f74f8205235102d
+        .quad 0x3f73365b88c0f347
+        .quad 0x3f7175085ab85ff0
+        .quad 0x3f6f684d1d8ae702
+        .quad 0x3f6be76bd77b4fc3
+        .quad 0x3f68676c71434fb9
+        .quad 0x3f64e84e793a474a
+        .quad 0x3f616a117e0d4b30
+        .quad 0x3f5bd96a1d7d9cbc
+        .quad 0x3f54e071754c98ba
+        .quad 0x3f4bd27045bfd025
+        .quad 0x3f3bcef518e29612
+        .quad 0x8000000000000000
+        /*== poly_coeff[5] ==*/
+        .align 16
+        .quad 0x3fb63C65231FBD16, 0x3fb63C65231FBD16 /* coeff5 */
+        .quad 0xbfbBCB7D4EFBE80B, 0xbfbBCB7D4EFBE80B /* coeff4 */
+        .quad 0x3fc287A7636F341E, 0x3fc287A7636F341E /* coeff3 */
+        .quad 0xbfcBCB7B1526DE36, 0xbfcBCB7B1526DE36 /* coeff2 */
+        .quad 0x3fdBCB7B1526E50E, 0x3fdBCB7B1526E50E /* coeff1 */
+        /*== ExpMask ==*/
+        .align 16
+        .quad 0x000fffffffffffff, 0x000fffffffffffff
+        /*== Two10 ==*/
+        .align 16
+        .quad 0x3f50000000000000, 0x3f50000000000000
+        /*== MinNorm ==*/
+        .align 16
+        .quad 0x0010000000000000, 0x0010000000000000
+        /*== MaxNorm ==*/
+        .align 16
+        .quad 0x7fefffffffffffff, 0x7fefffffffffffff
+        /*== HalfMask ==*/
+        .align 16
+        .quad 0xfffffffffc000000, 0xfffffffffc000000
+        /*== One ==*/
+        .align 16
+        .quad 0x3ff0000000000000, 0x3ff0000000000000
+        /*== Threshold ==*/
+        .align 16
+        .quad 0x4086a00000000000, 0x4086a00000000000
+        /*== Bias ==*/
+        .align 16
+        .quad 0x408ff80000000000, 0x408ff80000000000
+        /*== Bias1 ==*/
+        .align 16
+        .quad 0x408ff00000000000, 0x408ff00000000000
+        /*== L2 ==*/
+        .align 16
+        .quad 0x3fd34413509f79ff, 0x3fd34413509f79ff
+        .align 16
+        .type	__svml_dlog10_data_internal,@object
+        .size	__svml_dlog10_data_internal,.-__svml_dlog10_data_internal
+        .space 48, 0x00
+        .align 16
+
+.FLT_12:
+        .long	0x00000000,0x43380000,0x00000000,0x43380000
+        .type	.FLT_12,@object
+        .size	.FLT_12,16
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core-sse.S b/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core-sse.S
new file mode 100644
index 0000000..0a10166
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core-sse.S
@@ -0,0 +1,20 @@
+/* SSE version of vectorized log10, vector length is 4.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define _ZGVdN4v_log10 _ZGVdN4v_log10_sse_wrapper
+#include "../svml_d_log104_core.S"
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core.c b/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core.c
new file mode 100644
index 0000000..48c63cf
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core.c
@@ -0,0 +1,27 @@
+/* Multiple versions of vectorized log10, vector length is 4.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define SYMBOL_NAME _ZGVdN4v_log10
+#include "ifunc-mathvec-avx2.h"
+
+libc_ifunc_redirected (REDIRECT_NAME, SYMBOL_NAME, IFUNC_SELECTOR ());
+
+#ifdef SHARED
+__hidden_ver1 (_ZGVdN4v_log10, __GI__ZGVdN4v_log10, __redirect__ZGVdN4v_log10)
+  __attribute__ ((visibility ("hidden")));
+#endif
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core_avx2.S b/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core_avx2.S
new file mode 100644
index 0000000..df23926
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log104_core_avx2.S
@@ -0,0 +1,1074 @@
+/* Function log10 vectorized with AVX2.
+   Copyright (C) 2021 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
+   https://www.gnu.org/licenses/.  */
+
+/*
+ * ALGORITHM DESCRIPTION:
+ *
+ *    Get short reciprocal approximation Rcp ~ 1/mantissa(x)
+ *    R = Rcp*x - 1.0
+ *    log10(x) = k*log10(2.0) - log10(Rcp) + poly_approximation(R)
+ *       log10(Rcp) is tabulated
+ *
+ *
+ */
+
+/* Offsets for data table __svml_dlog10_data_internal
+ */
+#define Log_HA_table                  	0
+#define Log_LA_table                  	4128
+#define poly_coeff                    	8256
+#define ExpMask                       	8416
+#define Two10                         	8448
+#define MinNorm                       	8480
+#define MaxNorm                       	8512
+#define HalfMask                      	8544
+#define One                           	8576
+#define Threshold                     	8608
+#define Bias                          	8640
+#define Bias1                         	8672
+#define L2                            	8704
+
+/* Lookup bias for data table __svml_dlog10_data_internal.  */
+#define Table_Lookup_Bias               -0x406fe0
+
+#include <sysdep.h>
+
+        .text
+	.section .text.avx2,"ax",@progbits
+ENTRY(_ZGVdN4v_log10_avx2)
+        pushq     %rbp
+        cfi_def_cfa_offset(16)
+        movq      %rsp, %rbp
+        cfi_def_cfa(6, 16)
+        cfi_offset(6, -16)
+        andq      $-32, %rsp
+        subq      $96, %rsp
+        lea       Table_Lookup_Bias+__svml_dlog10_data_internal(%rip), %r8
+        vmovapd   %ymm0, %ymm3
+
+/* preserve mantissa, set input exponent to 2^(-10) */
+        vandpd    ExpMask+__svml_dlog10_data_internal(%rip), %ymm3, %ymm4
+        vorpd     Two10+__svml_dlog10_data_internal(%rip), %ymm4, %ymm2
+
+/* reciprocal approximation good to at least 11 bits */
+        vcvtpd2ps %ymm2, %xmm5
+
+/* exponent bits */
+        vpsrlq    $20, %ymm3, %ymm7
+        vmovupd   One+__svml_dlog10_data_internal(%rip), %ymm14
+        vrcpps    %xmm5, %xmm6
+
+/* check range */
+        vcmplt_oqpd MinNorm+__svml_dlog10_data_internal(%rip), %ymm3, %ymm11
+        vcmpnle_uqpd MaxNorm+__svml_dlog10_data_internal(%rip), %ymm3, %ymm12
+        vcvtps2pd %xmm6, %ymm9
+
+/* round reciprocal to nearest integer, will have 1+9 mantissa bits */
+        vroundpd  $0, %ymm9, %ymm1
+
+/* exponent*log(2.0) */
+        vmovupd   Threshold+__svml_dlog10_data_internal(%rip), %ymm9
+
+/*
+ * prepare table index
+ * table lookup
+ */
+        vpsrlq    $40, %ymm1, %ymm15
+
+/* argument reduction */
+        vfmsub213pd %ymm14, %ymm1, %ymm2
+        vcmplt_oqpd %ymm1, %ymm9, %ymm1
+        vorpd     %ymm12, %ymm11, %ymm13
+        vmovupd   poly_coeff+64+__svml_dlog10_data_internal(%rip), %ymm12
+        vfmadd213pd poly_coeff+96+__svml_dlog10_data_internal(%rip), %ymm2, %ymm12
+
+/* combine and get argument value range mask */
+        vmovmskpd %ymm13, %eax
+        vmulpd    %ymm2, %ymm2, %ymm13
+        vextractf128 $1, %ymm7, %xmm8
+        vshufps   $221, %xmm8, %xmm7, %xmm10
+
+/* biased exponent in DP format */
+        vcvtdq2pd %xmm10, %ymm0
+        vandpd    Bias+__svml_dlog10_data_internal(%rip), %ymm1, %ymm10
+        vorpd     Bias1+__svml_dlog10_data_internal(%rip), %ymm10, %ymm11
+        vsubpd    %ymm11, %ymm0, %ymm0
+        vmulpd    L2+__svml_dlog10_data_internal(%rip), %ymm0, %ymm1
+
+/* polynomial */
+        vmovupd   poly_coeff+__svml_dlog10_data_internal(%rip), %ymm0
+        vfmadd213pd poly_coeff+32+__svml_dlog10_data_internal(%rip), %ymm2, %ymm0
+        vmulpd    poly_coeff+128+__svml_dlog10_data_internal(%rip), %ymm2, %ymm2
+        vfmadd213pd %ymm12, %ymm13, %ymm0
+        vfmadd213pd %ymm2, %ymm13, %ymm0
+        vextractf128 $1, %ymm15, %xmm6
+        vmovd     %xmm15, %edx
+        vmovd     %xmm6, %esi
+        movslq    %edx, %rdx
+        vpextrd   $2, %xmm15, %ecx
+        movslq    %esi, %rsi
+        vpextrd   $2, %xmm6, %edi
+        movslq    %ecx, %rcx
+        movslq    %edi, %rdi
+        vmovsd    (%r8,%rdx), %xmm4
+        vmovsd    (%r8,%rsi), %xmm7
+        vmovhpd   (%r8,%rcx), %xmm4, %xmm5
+        vmovhpd   (%r8,%rdi), %xmm7, %xmm8
+        vinsertf128 $1, %xmm8, %ymm5, %ymm14
+
+/* reconstruction */
+        vaddpd    %ymm0, %ymm14, %ymm2
+        vaddpd    %ymm2, %ymm1, %ymm0
+        testl     %eax, %eax
+
+/* Go to special inputs processing branch */
+        jne       L(SPECIAL_VALUES_BRANCH)
+                                # LOE rbx r12 r13 r14 r15 eax ymm0 ymm3
+
+/* Restore registers
+ * and exit the function
+ */
+
+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)
+
+/* Branch to process
+ * special inputs
+ */
+
+L(SPECIAL_VALUES_BRANCH):
+        vmovupd   %ymm3, 32(%rsp)
+        vmovupd   %ymm0, 64(%rsp)
+                                # LOE rbx r12 r13 r14 r15 eax ymm0
+
+        xorl      %edx, %edx
+                                # LOE rbx r12 r13 r14 r15 eax edx
+
+        vzeroupper
+        movq      %r12, 16(%rsp)
+        /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -80; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xb0, 0xff, 0xff, 0xff, 0x22
+        movl      %edx, %r12d
+        movq      %r13, 8(%rsp)
+        /*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -88; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa8, 0xff, 0xff, 0xff, 0x22
+        movl      %eax, %r13d
+        movq      %r14, (%rsp)
+        /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -96; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa0, 0xff, 0xff, 0xff, 0x22
+                                # LOE rbx r15 r12d r13d
+
+/* Range mask
+ * bits check
+ */
+
+L(RANGEMASK_CHECK):
+        btl       %r12d, %r13d
+
+/* Call scalar math function */
+        jc        L(SCALAR_MATH_CALL)
+                                # LOE rbx r15 r12d r13d
+
+/* Special inputs
+ * processing loop
+ */
+
+L(SPECIAL_VALUES_LOOP):
+        incl      %r12d
+        cmpl      $4, %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)
+        vmovupd   64(%rsp), %ymm0
+
+/* Go to exit */
+        jmp       L(EXIT)
+        /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -80; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xb0, 0xff, 0xff, 0xff, 0x22
+        /*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -88; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa8, 0xff, 0xff, 0xff, 0x22
+        /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -96; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa0, 0xff, 0xff, 0xff, 0x22
+                                # LOE rbx r12 r13 r14 r15 ymm0
+
+/* Scalar math fucntion call
+ * to process special input
+ */
+
+L(SCALAR_MATH_CALL):
+        movl      %r12d, %r14d
+        movsd     32(%rsp,%r14,8), %xmm0
+        call      log10@PLT
+                                # LOE rbx r14 r15 r12d r13d xmm0
+
+        movsd     %xmm0, 64(%rsp,%r14,8)
+
+/* Process special inputs in loop */
+        jmp       L(SPECIAL_VALUES_LOOP)
+                                # LOE rbx r15 r12d r13d
+END(_ZGVdN4v_log10_avx2)
+
+        .section .rodata, "a"
+        .align 32
+
+#ifdef __svml_dlog10_data_internal_typedef
+typedef unsigned int VUINT32;
+typedef struct {
+        __declspec(align(32)) VUINT32 Log_HA_table[(1<<9)+2][2];
+        __declspec(align(32)) VUINT32 Log_LA_table[(1<<9)+1][2];
+        __declspec(align(32)) VUINT32 poly_coeff[5][4][2];
+        __declspec(align(32)) VUINT32 ExpMask[4][2];
+        __declspec(align(32)) VUINT32 Two10[4][2];
+        __declspec(align(32)) VUINT32 MinNorm[4][2];
+        __declspec(align(32)) VUINT32 MaxNorm[4][2];
+        __declspec(align(32)) VUINT32 HalfMask[4][2];
+        __declspec(align(32)) VUINT32 One[4][2];
+        __declspec(align(32)) VUINT32 Threshold[4][2];
+        __declspec(align(32)) VUINT32 Bias[4][2];
+        __declspec(align(32)) VUINT32 Bias1[4][2];
+        __declspec(align(32)) VUINT32 L2[4][2];
+} __svml_dlog10_data_internal;
+#endif
+__svml_dlog10_data_internal:
+        /* Log_HA_table */
+        .quad 0xc0733a7146f6b080, 0xbe1e707ce619c200
+        .quad 0xc0733a7547771970, 0xbe1e79c6c06d6f51
+        .quad 0xc0733a7945aacb70, 0xbe1e78e225fad29c
+        .quad 0xc0733a7d41946970, 0xbe1e76d607f9693b
+        .quad 0xc0733a813b3691f0, 0xbe1e7704b3e0685b
+        .quad 0xc0733a853293df00, 0xbe1e79c1216a27fa
+        .quad 0xc0733a8927aee660, 0xbe1e76dce5734a81
+        .quad 0xc0733a8d1a8a3920, 0xbe1e782ee2ca4dba
+        .quad 0xc0733a910b286430, 0xbe1e7812d1a0a61f
+        .quad 0xc0733a94f98bf010, 0xbe1e77e1b5ecbc61
+        .quad 0xc0733a98e5b76100, 0xbe1e76635cac1586
+        .quad 0xc0733a9ccfad36f0, 0xbe1e7638f7968f32
+        .quad 0xc0733aa0b76feda0, 0xbe1e7840ee76e365
+        .quad 0xc0733aa49d01fcb0, 0xbe1e79f3fd01907e
+        .quad 0xc0733aa88065d7a0, 0xbe1e77bbb3a9c38a
+        .quad 0xc0733aac619dedb0, 0xbe1e7742719bf41d
+        .quad 0xc0733ab040acaa20, 0xbe1e79bcedaf79cb
+        .quad 0xc0733ab41d947450, 0xbe1e762d63cb7ca0
+        .quad 0xc0733ab7f857af50, 0xbe1e77a07be83403
+        .quad 0xc0733abbd0f8ba80, 0xbe1e7763ff836ad0
+        .quad 0xc0733abfa779f130, 0xbe1e7737720ead39
+        .quad 0xc0733ac37bddaad0, 0xbe1e7776a08e55e7
+        .quad 0xc0733ac74e263af0, 0xbe1e793e3c52dd36
+        .quad 0xc0733acb1e55f160, 0xbe1e788a94695051
+        .quad 0xc0733aceec6f1a10, 0xbe1e76508114a813
+        .quad 0xc0733ad2b873fd20, 0xbe1e76909457d23e
+        .quad 0xc0733ad68266df10, 0xbe1e7664a24f9ca4
+        .quad 0xc0733ada4a4a0090, 0xbe1e7a07b3d44b18
+        .quad 0xc0733ade101f9ee0, 0xbe1e76d87594704d
+        .quad 0xc0733ae1d3e9f340, 0xbe1e79563595a182
+        .quad 0xc0733ae595ab33b0, 0xbe1e771880c3c6ab
+        .quad 0xc0733ae955659250, 0xbe1e78c171f517d4
+        .quad 0xc0733aed131b3df0, 0xbe1e77eac3874666
+        .quad 0xc0733af0cece61b0, 0xbe1e790db479d8f6
+        .quad 0xc0733af488812550, 0xbe1e7965d1aa5c90
+        .quad 0xc0733af84035ad10, 0xbe1e78ceb398ba47
+        .quad 0xc0733afbf5ee19c0, 0xbe1e779cc0dcb5aa
+        .quad 0xc0733affa9ac88c0, 0xbe1e7871053953ed
+        .quad 0xc0733b035b731420, 0xbe1e7a082cffa71a
+        .quad 0xc0733b070b43d2a0, 0xbe1e7904b4382fad
+        .quad 0xc0733b0ab920d790, 0xbe1e79b458d0b4f3
+        .quad 0xc0733b0e650c3310, 0xbe1e79d0ded414c6
+        .quad 0xc0733b120f07f200, 0xbe1e763c357a1943
+        .quad 0xc0733b15b7161dd0, 0xbe1e78b80ba6daaa
+        .quad 0xc0733b195d38bd00, 0xbe1e7998e23b8ffd
+        .quad 0xc0733b1d0171d2c0, 0xbe1e7974aa65ee8c
+        .quad 0xc0733b20a3c35f20, 0xbe1e76ccfde752ab
+        .quad 0xc0733b24442f5ef0, 0xbe1e77b4ff19debb
+        .quad 0xc0733b27e2b7cc10, 0xbe1e7772ee478542
+        .quad 0xc0733b2b7f5e9d30, 0xbe1e781d81b58b44
+        .quad 0xc0733b2f1a25c600, 0xbe1e78350d967565
+        .quad 0xc0733b32b30f3720, 0xbe1e783888e48152
+        .quad 0xc0733b364a1cde30, 0xbe1e78367bf7c111
+        .quad 0xc0733b39df50a5d0, 0xbe1e7959e57ca47d
+        .quad 0xc0733b3d72ac75c0, 0xbe1e777322423222
+        .quad 0xc0733b41043232b0, 0xbe1e767ce42a60aa
+        .quad 0xc0733b4493e3be70, 0xbe1e781d445aea19
+        .quad 0xc0733b4821c2f800, 0xbe1e7922fca18e18
+        .quad 0xc0733b4badd1bb80, 0xbe1e76fed3d40647
+        .quad 0xc0733b4f3811e210, 0xbe1e793948c9eabc
+        .quad 0xc0733b52c0854240, 0xbe1e76e487656b8c
+        .quad 0xc0733b56472daf90, 0xbe1e780ab2f71223
+        .quad 0xc0733b59cc0cfaf0, 0xbe1e77189120b09c
+        .quad 0xc0733b5d4f24f270, 0xbe1e7644a0343a12
+        .quad 0xc0733b60d0776160, 0xbe1e78f2a3e4733d
+        .quad 0xc0733b6450061080, 0xbe1e7913b2f73ae5
+        .quad 0xc0733b67cdd2c5c0, 0xbe1e7882d08393b5
+        .quad 0xc0733b6b49df4470, 0xbe1e765e1b209979
+        .quad 0xc0733b6ec42d4d20, 0xbe1e785c9c4620d4
+        .quad 0xc0733b75b394f240, 0xbe1e78878cd0e956
+        .quad 0xc0733b7c9c178630, 0xbe1e789a4112d90b
+        .quad 0xc0733b837dc2b0f0, 0xbe1e79050b8a1766
+        .quad 0xc0733b8a58a3f220, 0xbe1e7790dffc47aa
+        .quad 0xc0733b912cc8a180, 0xbe1e77174593b06a
+        .quad 0xc0733b97fa3defb0, 0xbe1e7677de2d2ecc
+        .quad 0xc0733b9ec110e6b0, 0xbe1e76cff477ca18
+        .quad 0xc0733ba5814e6a80, 0xbe1e78f8644dec7b
+        .quad 0xc0733bac3b0339d0, 0xbe1e764e1361788d
+        .quad 0xc0733bb2ee3bee30, 0xbe1e78c913e738de
+        .quad 0xc0733bb99b04fd30, 0xbe1e76666f5bddaa
+        .quad 0xc0733bc0416ab850, 0xbe1e77e87cbd8ab6
+        .quad 0xc0733bc6e1794e10, 0xbe1e76f18ba1c966
+        .quad 0xc0733bcd7b3cca10, 0xbe1e777c9461b8db
+        .quad 0xc0733bd40ec115d0, 0xbe1e78b78526ffac
+        .quad 0xc0733bda9c11f920, 0xbe1e7942abecfede
+        .quad 0xc0733be1233b1aa0, 0xbe1e76d8a684fd8c
+        .quad 0xc0733be7a4480010, 0xbe1e79622b539ac9
+        .quad 0xc0733bee1f440f30, 0xbe1e7978e7cc20ea
+        .quad 0xc0733bf4943a8de0, 0xbe1e765c9c9de825
+        .quad 0xc0733bfb0336a290, 0xbe1e775d8b138ee2
+        .quad 0xc0733c016c435500, 0xbe1e78bf33465c2f
+        .quad 0xc0733c07cf6b8e80, 0xbe1e78164f7cc441
+        .quad 0xc0733c0e2cba1a50, 0xbe1e7824e64d0b23
+        .quad 0xc0733c148439a630, 0xbe1e78373ae7dd81
+        .quad 0xc0733c1ad5f4c2c0, 0xbe1e7704513e0afe
+        .quad 0xc0733c2121f5e3d0, 0xbe1e7914aa84200f
+        .quad 0xc0733c2768476110, 0xbe1e76b1cde25cf6
+        .quad 0xc0733c2da8f37600, 0xbe1e796120e3862d
+        .quad 0xc0733c33e40442e0, 0xbe1e78ec836d7e7b
+        .quad 0xc0733c3a1983cca0, 0xbe1e77fb13b7dabb
+        .quad 0xc0733c40497bfd70, 0xbe1e783c6fcb2404
+        .quad 0xc0733c4673f6a530, 0xbe1e7628bb93dce8
+        .quad 0xc0733c4c98fd7990, 0xbe1e7857a47b5001
+        .quad 0xc0733c52b89a16d0, 0xbe1e76708dc2831f
+        .quad 0xc0733c58d2d5ffa0, 0xbe1e77b6038651f1
+        .quad 0xc0733c5ee7ba9de0, 0xbe1e792e855bb5b2
+        .quad 0xc0733c64f75142d0, 0xbe1e776cacd5c105
+        .quad 0xc0733c6b01a32740, 0xbe1e77f8a8011315
+        .quad 0xc0733c7106b96c30, 0xbe1e765cf3efcfde
+        .quad 0xc0733c77069d1ad0, 0xbe1e78d837d2efac
+        .quad 0xc0733c7d01572530, 0xbe1e78b615cf772c
+        .quad 0xc0733c82f6f06640, 0xbe1e7650bbbd7a25
+        .quad 0xc0733c88e771a220, 0xbe1e78bcf3495872
+        .quad 0xc0733c8ed2e386c0, 0xbe1e792266832e84
+        .quad 0xc0733c94b94eabd0, 0xbe1e79c1c3c2ca52
+        .quad 0xc0733c9a9abb9340, 0xbe1e78aa61e5807d
+        .quad 0xc0733ca07732a970, 0xbe1e7620fc4cf156
+        .quad 0xc0733ca64ebc4570, 0xbe1e76b914a832c5
+        .quad 0xc0733cac2160a970, 0xbe1e79227f72020e
+        .quad 0xc0733cb1ef280300, 0xbe1e77ac972cc008
+        .quad 0xc0733cb7b81a6b10, 0xbe1e798089be41f4
+        .quad 0xc0733cbd7c3fe6a0, 0xbe1e77942ae037fe
+        .quad 0xc0733cc33ba06690, 0xbe1e7956ae6463d9
+        .quad 0xc0733cc8f643c850, 0xbe1e7918a50c7942
+        .quad 0xc0733cceac31d5d0, 0xbe1e78308eeab604
+        .quad 0xc0733cd45d7245e0, 0xbe1e76dd4ea88445
+        .quad 0xc0733cda0a0cbc60, 0xbe1e77e7c1aa5909
+        .quad 0xc0733cdfb208caa0, 0xbe1e7804b9d20e54
+        .quad 0xc0733ce5556def70, 0xbe1e78f88e99d49c
+        .quad 0xc0733ceaf4439780, 0xbe1e787d74682d68
+        .quad 0xc0733cf08e911d80, 0xbe1e76edc24fe6e7
+        .quad 0xc0733cf6245dca50, 0xbe1e79b347ec86d2
+        .quad 0xc0733cfbb5b0d580, 0xbe1e797cceb2c39b
+        .quad 0xc0733d0142916530, 0xbe1e783adbdc6aa1
+        .quad 0xc0733d06cb068e70, 0xbe1e76e4c20e3d9e
+        .quad 0xc0733d0c4f175570, 0xbe1e77070bf3cf61
+        .quad 0xc0733d11cecaadc0, 0xbe1e781c43502734
+        .quad 0xc0733d174a277a80, 0xbe1e78b11268ea72
+        .quad 0xc0733d1cc1348e90, 0xbe1e7754b83bfc7d
+        .quad 0xc0733d2233f8acb0, 0xbe1e7756c29bf5e9
+        .quad 0xc0733d27a27a87d0, 0xbe1e7952fc1d9333
+        .quad 0xc0733d2d0cc0c350, 0xbe1e778c76ae6077
+        .quad 0xc0733d3272d1f2e0, 0xbe1e7a1896ba8f43
+        .quad 0xc0733d37d4b49b30, 0xbe1e76dafdf432d8
+        .quad 0xc0733d3d326f3180, 0xbe1e795330184013
+        .quad 0xc0733d428c081c80, 0xbe1e763cc774d30f
+        .quad 0xc0733d47e185b3d0, 0xbe1e77030a779c0a
+        .quad 0xc0733d4d32ee40b0, 0xbe1e7908af2a2d7e
+        .quad 0xc0733d528047fe00, 0xbe1e78c4953b797d
+        .quad 0xc0733d57c9991850, 0xbe1e78b43b096579
+        .quad 0xc0733d5d0ee7ae30, 0xbe1e7824ae0a4804
+        .quad 0xc0733d625039d040, 0xbe1e79d2b2fbb740
+        .quad 0xc0733d678d958190, 0xbe1e7662de59a1a6
+        .quad 0xc0733d6cc700b760, 0xbe1e76b251d59aaa
+        .quad 0xc0733d71fc8159b0, 0xbe1e7a00cfd1f487
+        .quad 0xc0733d772e1d4360, 0xbe1e77f4d246167e
+        .quad 0xc0733d7c5bda4200, 0xbe1e767a4ee8e6fc
+        .quad 0xc0733d8185be1640, 0xbe1e777ccf0a8aed
+        .quad 0xc0733d86abce7420, 0xbe1e767d7e279ada
+        .quad 0xc0733d8bce1102d0, 0xbe1e7a05cef4bb90
+        .quad 0xc0733d90ec8b5d40, 0xbe1e78f75369be5b
+        .quad 0xc0733d96074311d0, 0xbe1e77b9612e8c8a
+        .quad 0xc0733d9b1e3da2b0, 0xbe1e794518b9adeb
+        .quad 0xc0733da031808620, 0xbe1e7810626fb934
+        .quad 0xc0733da541112650, 0xbe1e76d87223fa6d
+        .quad 0xc0733daa4cf4e1a0, 0xbe1e794c5e7ca3b5
+        .quad 0xc0733daf55310af0, 0xbe1e789856ef816f
+        .quad 0xc0733db459cae970, 0xbe1e77d2004effbd
+        .quad 0xc0733db95ac7b8f0, 0xbe1e78467d31eb9c
+        .quad 0xc0733dbe582caa00, 0xbe1e79aaa4e25787
+        .quad 0xc0733dc351fee220, 0xbe1e762de8f107bf
+        .quad 0xc0733dc848437b90, 0xbe1e7670670a63fe
+        .quad 0xc0733dcd3aff85d0, 0xbe1e795ca237c6cc
+        .quad 0xc0733dd22a3805b0, 0xbe1e77e55c53c1d9
+        .quad 0xc0733dd715f1f520, 0xbe1e78a806213ac4
+        .quad 0xc0733ddbfe3243b0, 0xbe1e77743a2bc615
+        .quad 0xc0733de0e2fdd660, 0xbe1e78b8b45b0b7d
+        .quad 0xc0733de5c4598800, 0xbe1e78d635f2f4b9
+        .quad 0xc0733deaa24a2920, 0xbe1e7758c396a11e
+        .quad 0xc0733def7cd48020, 0xbe1e7a17a8cc454c
+        .quad 0xc0733df453fd49a0, 0xbe1e783caa73f616
+        .quad 0xc0733df927c93820, 0xbe1e7932cfa29664
+        .quad 0xc0733dfdf83cf490, 0xbe1e777d265c72a6
+        .quad 0xc0733e02c55d1e10, 0xbe1e7775e7c03c60
+        .quad 0xc0733e078f2e4a40, 0xbe1e79f65d52d232
+        .quad 0xc0733e0c55b50570, 0xbe1e76e7e7464b4e
+        .quad 0xc0733e1118f5d250, 0xbe1e77be81cad877
+        .quad 0xc0733e15d8f52a80, 0xbe1e79dd25b5fb3a
+        .quad 0xc0733e1a95b77e80, 0xbe1e78e45f1418ef
+        .quad 0xc0733e1f4f4135a0, 0xbe1e78eb7289505b
+        .quad 0xc0733e240596ae50, 0xbe1e78a468c07cad
+        .quad 0xc0733e28b8bc3e20, 0xbe1e776b558a4009
+        .quad 0xc0733e2d68b631d0, 0xbe1e77412eb9941e
+        .quad 0xc0733e321588cd80, 0xbe1e76b2853f845e
+        .quad 0xc0733e36bf384cb0, 0xbe1e76aa7184273c
+        .quad 0xc0733e3b65c8e260, 0xbe1e7832027f78fa
+        .quad 0xc0733e40093eb930, 0xbe1e7a1c7da131f5
+        .quad 0xc0733e44a99df380, 0xbe1e76a0bc2ae4bc
+        .quad 0xc0733e4946eaab30, 0xbe1e78dff13b6f5d
+        .quad 0xc0733e4de128f250, 0xbe1e765a226dea2c
+        .quad 0xc0733e52785cd290, 0xbe1e78509b989111
+        .quad 0xc0733e570c8a4de0, 0xbe1e7916a4e9803d
+        .quad 0xc0733e5b9db55e30, 0xbe1e7950c15758cc
+        .quad 0xc0733e602be1f5a0, 0xbe1e7922ba1ad420
+        .quad 0xc0733e64b713fe90, 0xbe1e794cbaabcef6
+        .quad 0xc0733e693f4f5bc0, 0xbe1e7837bf883fed
+        .quad 0xc0733e6dc497e850, 0xbe1e76f198ddbbdf
+        .quad 0xc0733e7246f177d0, 0xbe1e7a18c1067764
+        .quad 0xc0733e76c65fd6a0, 0xbe1e76b845a8fd9d
+        .quad 0xc0733e7b42e6c970, 0xbe1e7714012df506
+        .quad 0xc0733e7fbc8a0de0, 0xbe1e7765612922cd
+        .quad 0xc0733e84334d5a50, 0xbe1e7688f5424a00
+        .quad 0xc0733e88a7345df0, 0xbe1e769d011f6663
+        .quad 0xc0733e8d1842c0e0, 0xbe1e79914acbfaf7
+        .quad 0xc0733e91867c2460, 0xbe1e79a85e189bd7
+        .quad 0xc0733e95f1e422a0, 0xbe1e79ea7c726432
+        .quad 0xc0733e9a5a7e4f10, 0xbe1e768a6fbb8e6e
+        .quad 0xc0733e9ec04e3620, 0xbe1e793c75bcc9fc
+        .quad 0xc0733ea323575dd0, 0xbe1e797f78da13d4
+        .quad 0xc0733ea7839d4550, 0xbe1e78d8c9cda978
+        .quad 0xc0733eabe1236540, 0xbe1e77028d480fff
+        .quad 0xc0733eb03bed2fa0, 0xbe1e7a0d0f74ff7c
+        .quad 0xc0733eb493fe1040, 0xbe1e76732e8a35fb
+        .quad 0xc0733eb8e9596c30, 0xbe1e77220caeabeb
+        .quad 0xc0733ebd3c02a260, 0xbe1e797438b645ef
+        .quad 0xc0733ec18bfd0b80, 0xbe1e79207c5fd6e8
+        .quad 0xc0733ec5d94bf9f0, 0xbe1e781c7df8f946
+        .quad 0xc0733eca23f2b9f0, 0xbe1e76736284e2db
+        .quad 0xc0733ece6bf49190, 0xbe1e7a109cc0c3f5
+        .quad 0xc0733ed2b154c120, 0xbe1e767f14a16d50
+        .quad 0xc0733ed6f4168290, 0xbe1e789cd22acaf0
+        .quad 0xc0733edb343d0a40, 0xbe1e764355ca28ad
+        .quad 0xc0733edf71cb8660, 0xbe1e79e4c7a81c45
+        .quad 0xc0733ee3acc51fb0, 0xbe1e761e26b644c2
+        .quad 0xc0733ee7e52cf8c0, 0xbe1e793e9f8fbdd3
+        .quad 0xc0733eec1b062ed0, 0xbe1e78c432991c20
+        .quad 0xc0733ef04e53d940, 0xbe1e78cdd025f4d8
+        .quad 0xc0733ef47f1909f0, 0xbe1e778310c6446e
+        .quad 0xc0733ef8ad58cd20, 0xbe1e7871af3d6e17
+        .quad 0xc0733efcd91629b0, 0xbe1e77e0e906f697
+        .quad 0xc0733f01025420f0, 0xbe1e7a1ae9b27892
+        .quad 0xc0733f052915af00, 0xbe1e76ac64c88f9d
+        .quad 0xc0733f094d5dca60, 0xbe1e779a815589c4
+        .quad 0xc0733f0d6f2f6480, 0xbe1e788f39a4864c
+        .quad 0xc0733f118e8d6980, 0xbe1e79fc51263525
+        .quad 0xc0733f15ab7ac060, 0xbe1e783501f19e90
+        .quad 0xc0733f19c5fa4ae0, 0xbe1e767e82c327ab
+        .quad 0xc0733f1dde0ee5a0, 0xbe1e7a1785d66123
+        .quad 0xc0733f21f3bb6870, 0xbe1e7936d07203da
+        .quad 0xc0733f260702a5e0, 0xbe1e7a010a7ac699
+        .quad 0xc0733f2a17e76bb0, 0xbe1e7975e4e16312
+        .quad 0xc0733f2e266c82b0, 0xbe1e7654b5422330
+        .quad 0xc0733f323294aeb0, 0xbe1e77f8a4909d35
+        .quad 0xc0733f363c62aee0, 0xbe1e792c8e30d226
+        .quad 0xc0733f3a43d93da0, 0xbe1e76f6ac67a1ff
+        .quad 0xc0733f3e48fb1070, 0xbe1e775c2e97715a
+        .quad 0xc0733f424bcad840, 0xbe1e781cd54ae100
+        /*== Log_LA_table ==*/
+        .align 32
+        .quad 0x0000000000000000
+        .quad 0xbf4bc48a867884b7
+        .quad 0xbf5bbd9e9482af09
+        .quad 0xbf64c9096b94befd
+        .quad 0xbf6bafd47221ed26
+        .quad 0xbf714999e2ad8ea6
+        .quad 0xbf74b99563d2a1bd
+        .quad 0xbf7827de6b310350
+        .quad 0xbf7b9476a4fcd10f
+        .quad 0xbf7eff5fbaf25781
+        .quad 0xbf81344daa2d7553
+        .quad 0xbf82e8158b08d957
+        .quad 0xbf849b0851443684
+        .quad 0xbf864d26cce610dd
+        .quad 0xbf87fe71ccc4e6b0
+        .quad 0xbf89aeea1e897fdf
+        .quad 0xbf8b5e908eb13790
+        .quad 0xbf8d0d65e890405a
+        .quad 0xbf8ebb6af653e2ee
+        .quad 0xbf90345040825bad
+        .quad 0xbf910a83a8446c78
+        .quad 0xbf91e05015d30a71
+        .quad 0xbf92b5b5ec0209d3
+        .quad 0xbf938ab58d173e91
+        .quad 0xbf945f4f5acb8be0
+        .quad 0xbf953383b64bf13f
+        .quad 0xbf960753003a94ef
+        .quad 0xbf96dabd98afcc05
+        .quad 0xbf97adc3df3b1ff8
+        .quad 0xbf98806632e451d0
+        .quad 0xbf9952a4f22c5ae9
+        .quad 0xbf9a24807b0e6b5c
+        .quad 0xbf9af5f92b00e610
+        .quad 0xbf9bc70f5ef65a77
+        .quad 0xbf9c97c3735e7c0a
+        .quad 0xbf9d6815c4271775
+        .quad 0xbf9e3806acbd058f
+        .quad 0xbf9f0796880d1c19
+        .quad 0xbf9fd6c5b0851c4c
+        .quad 0xbfa052ca400a4f9b
+        .quad 0xbfa0ba01a8170000
+        .quad 0xbfa121093ce3a205
+        .quad 0xbfa187e12aad8077
+        .quad 0xbfa1ee899d74a03e
+        .quad 0xbfa25502c0fc314c
+        .quad 0xbfa2bb4cc0cafe8d
+        .quad 0xbfa32167c82bdcda
+        .quad 0xbfa38754022e18e2
+        .quad 0xbfa3ed1199a5e425
+        .quad 0xbfa452a0b92cc0ec
+        .quad 0xbfa4b8018b21ed4f
+        .quad 0xbfa51d3439aacd4a
+        .quad 0xbfa58238eeb353da
+        .quad 0xbfa5e70fd3ee6b34
+        .quad 0xbfa64bb912d65c07
+        .quad 0xbfa6b034d4ad33df
+        .quad 0xbfa71483427d2a99
+        .quad 0xbfa778a4851906f3
+        .quad 0xbfa7dc98c51c8242
+        .quad 0xbfa840602aecab3d
+        .quad 0xbfa8a3fadeb847f4
+        .quad 0xbfa90769087836e4
+        .quad 0xbfa96aaacfefcf3c
+        .quad 0xbfa9cdc05cad4042
+        .quad 0xbfaa30a9d609efea
+        .quad 0xbfaa9367632ad897
+        .quad 0xbfaaf5f92b00e610
+        .quad 0xbfab585f544951a4
+        .quad 0xbfabba9a058dfd84
+        .quad 0xbfac1ca96525cf56
+        .quad 0xbfac7e8d993509f9
+        .quad 0xbface046c7ada68d
+        .quad 0xbfad41d5164facb4
+        .quad 0xbfada338aaa98a0c
+        .quad 0xbfae0471aa1868f5
+        .quad 0xbfae658039c88690
+        .quad 0xbfaec6647eb58808
+        .quad 0xbfaf271e9daacf20
+        .quad 0xbfaf87aebb43ce06
+        .quad 0xbfafe814fbec5a77
+        .quad 0xbfb02428c1f08016
+        .quad 0xbfb054323b97a948
+        .quad 0xbfb08426fcdb1ee7
+        .quad 0xbfb0b40717932b96
+        .quad 0xbfb0e3d29d81165e
+        .quad 0xbfb11389a04f4a2e
+        .quad 0xbfb1432c31917d08
+        .quad 0xbfb172ba62c4d6de
+        .quad 0xbfb1a23445501816
+        .quad 0xbfb1d199ea83bfbe
+        .quad 0xbfb200eb639a3173
+        .quad 0xbfb23028c1b7daed
+        .quad 0xbfb25f5215eb594a
+        .quad 0xbfb28e67712d9dfc
+        .quad 0xbfb2bd68e4621371
+        .quad 0xbfb2ec568056c16f
+        .quad 0xbfb31b3055c47118
+        .quad 0xbfb349f6754ed0b4
+        .quad 0xbfb378a8ef84971e
+        .quad 0xbfb3a747d4dfa6f5
+        .quad 0xbfb3d5d335c53179
+        .quad 0xbfb4044b2285d925
+        .quad 0xbfb432afab5dd3ff
+        .quad 0xbfb46100e0750da1
+        .quad 0xbfb48f3ed1df48fb
+        .quad 0xbfb4bd698f9c41cf
+        .quad 0xbfb4eb812997cde4
+        .quad 0xbfb51985afa9fdfd
+        .quad 0xbfb5477731973e85
+        .quad 0xbfb57555bf1077f5
+        .quad 0xbfb5a32167b32f02
+        .quad 0xbfb5d0da3b09a47e
+        .quad 0xbfb5fe80488af4fd
+        .quad 0xbfb62c139f9b3837
+        .quad 0xbfb659944f8ba02d
+        .quad 0xbfb68702679a980a
+        .quad 0xbfb6b45df6f3e2c9
+        .quad 0xbfb6e1a70cb0b99a
+        .quad 0xbfb70eddb7d7ea07
+        .quad 0xbfb73c02075df3e5
+        .quad 0xbfb769140a2526fd
+        .quad 0xbfb79613cefdc07d
+        .quad 0xbfb7c30164a60836
+        .quad 0xbfb7efdcd9ca6d8f
+        .quad 0xbfb81ca63d05a44a
+        .quad 0xbfb8495d9ce0c10c
+        .quad 0xbfb8760307d355ab
+        .quad 0xbfb8a2968c438d41
+        .quad 0xbfb8cf183886480d
+        .quad 0xbfb8fb881adf3713
+        .quad 0xbfb927e64180f790
+        .quad 0xbfb95432ba8d2e2f
+        .quad 0xbfb9806d9414a209
+        .quad 0xbfb9ac96dc175776
+        .quad 0xbfb9d8aea084aa9c
+        .quad 0xbfba04b4ef3b69d8
+        .quad 0xbfba30a9d609efea
+        .quad 0xbfba5c8d62ae3dec
+        .quad 0xbfba885fa2d6151e
+        .quad 0xbfbab420a41f1076
+        .quad 0xbfbadfd07416be07
+        .quad 0xbfbb0b6f203ab82c
+        .quad 0xbfbb36fcb5f8be8a
+        .quad 0xbfbb627942aecedd
+        .quad 0xbfbb8de4d3ab3d98
+        .quad 0xbfbbb93f762cce4f
+        .quad 0xbfbbe4893762cbf7
+        .quad 0xbfbc0fc2246d20f5
+        .quad 0xbfbc3aea4a5c6eff
+        .quad 0xbfbc6601b63226cb
+        .quad 0xbfbc910874e09f98
+        .quad 0xbfbcbbfe934b2e81
+        .quad 0xbfbce6e41e463da5
+        .quad 0xbfbd11b92297632b
+        .quad 0xbfbd3c7dacf5780b
+        .quad 0xbfbd6731ca08aeb9
+        .quad 0xbfbd91d5866aa99c
+        .quad 0xbfbdbc68eea6915b
+        .quad 0xbfbde6ec0f392b05
+        .quad 0xbfbe115ef490ee07
+        .quad 0xbfbe3bc1ab0e19fe
+        .quad 0xbfbe66143f02cc5d
+        .quad 0xbfbe9056bcb315e8
+        .quad 0xbfbeba893055100b
+        .quad 0xbfbee4aba610f204
+        .quad 0xbfbf0ebe2a0125eb
+        .quad 0xbfbf38c0c8325d86
+        .quad 0xbfbf62b38ca3a706
+        .quad 0xbfbf8c9683468191
+        .quad 0xbfbfb669b7fef1a8
+        .quad 0xbfbfe02d36a3956d
+        .quad 0xbfc004f0857edc5c
+        .quad 0xbfc019c2a064b486
+        .quad 0xbfc02e8cf1dac4b8
+        .quad 0xbfc0434f7fb1f307
+        .quad 0xbfc0580a4fb4a3df
+        .quad 0xbfc06cbd67a6c3b6
+        .quad 0xbfc08168cd45d0a9
+        .quad 0xbfc0960c8648e406
+        .quad 0xbfc0aaa89860bbcf
+        .quad 0xbfc0bf3d0937c41c
+        .quad 0xbfc0d3c9de722078
+        .quad 0xbfc0e84f1dadb526
+        .quad 0xbfc0fccccc823059
+        .quad 0xbfc11142f0811357
+        .quad 0xbfc125b18f35bb8e
+        .quad 0xbfc13a18ae256b99
+        .quad 0xbfc14e7852cf5430
+        .quad 0xbfc162d082ac9d10
+        .quad 0xbfc1772143306dc6
+        .quad 0xbfc18b6a99c7f679
+        .quad 0xbfc19fac8bda7897
+        .quad 0xbfc1b3e71ec94f7b
+        .quad 0xbfc1c81a57eff8fd
+        .quad 0xbfc1dc463ca41df8
+        .quad 0xbfc1f06ad2359abd
+        .quad 0xbfc204881dee8777
+        .quad 0xbfc2189e25134081
+        .quad 0xbfc22cacece26ead
+        .quad 0xbfc240b47a950f79
+        .quad 0xbfc254b4d35e7d3c
+        .quad 0xbfc268adfc6c773e
+        .quad 0xbfc27c9ffae729c1
+        .quad 0xbfc2908ad3f13603
+        .quad 0xbfc2a46e8ca7ba2a
+        .quad 0xbfc2b84b2a225923
+        .quad 0xbfc2cc20b1734279
+        .quad 0xbfc2dfef27a73a18
+        .quad 0xbfc2f3b691c5a001
+        .quad 0xbfc30776f4d077f7
+        .quad 0xbfc31b3055c47118
+        .quad 0xbfc32ee2b998ed6e
+        .quad 0xbfc3428e2540096d
+        .quad 0x3fc331f403985097
+        .quad 0x3fc31e56798a910a
+        .quad 0x3fc30abfd8f333b6
+        .quad 0x3fc2f7301cf4e87b
+        .quad 0x3fc2e3a740b7800f
+        .quad 0x3fc2d0253f67e4cb
+        .quad 0x3fc2bcaa14381386
+        .quad 0x3fc2a935ba5f1479
+        .quad 0x3fc295c82d18f434
+        .quad 0x3fc2826167a6bc9c
+        .quad 0x3fc26f01654e6df6
+        .quad 0x3fc25ba8215af7fc
+        .quad 0x3fc24855971c3307
+        .quad 0x3fc23509c1e6d937
+        .quad 0x3fc221c49d147fb3
+        .quad 0x3fc20e8624038fed
+        .quad 0x3fc1fb4e521740f4
+        .quad 0x3fc1e81d22b790d4
+        .quad 0x3fc1d4f291513e01
+        .quad 0x3fc1c1ce9955c0c6
+        .quad 0x3fc1aeb1363b44c8
+        .quad 0x3fc19b9a637ca295
+        .quad 0x3fc1888a1c995931
+        .quad 0x3fc175805d1587c1
+        .quad 0x3fc1627d2079e731
+        .quad 0x3fc14f806253c3ed
+        .quad 0x3fc13c8a1e34f7a0
+        .quad 0x3fc1299a4fb3e306
+        .quad 0x3fc116b0f26b67bb
+        .quad 0x3fc103ce01fae223
+        .quad 0x3fc0f0f17a062353
+        .quad 0x3fc0de1b56356b04
+        .quad 0x3fc0cb4b9235619a
+        .quad 0x3fc0b88229b71227
+        .quad 0x3fc0a5bf186fe483
+        .quad 0x3fc093025a19976c
+        .quad 0x3fc0804bea723aa9
+        .quad 0x3fc06d9bc53c2941
+        .quad 0x3fc05af1e63e03b4
+        .quad 0x3fc0484e4942aa43
+        .quad 0x3fc035b0ea19373b
+        .quad 0x3fc02319c494f951
+        .quad 0x3fc01088d48d6e03
+        .quad 0x3fbffbfc2bbc7803
+        .quad 0x3fbfd6f308ce5b52
+        .quad 0x3fbfb1f6381856f4
+        .quad 0x3fbf8d05b16a6d47
+        .quad 0x3fbf68216c9cc727
+        .quad 0x3fbf4349618fa91a
+        .quad 0x3fbf1e7d882b689a
+        .quad 0x3fbef9bdd860616b
+        .quad 0x3fbed50a4a26eafc
+        .quad 0x3fbeb062d57f4de8
+        .quad 0x3fbe8bc77271b97a
+        .quad 0x3fbe6738190e394c
+        .quad 0x3fbe42b4c16caaf3
+        .quad 0x3fbe1e3d63acb3ba
+        .quad 0x3fbdf9d1f7f5b674
+        .quad 0x3fbdd5727676c959
+        .quad 0x3fbdb11ed766abf4
+        .quad 0x3fbd8cd71303bd26
+        .quad 0x3fbd689b2193f133
+        .quad 0x3fbd446afb64c7e5
+        .quad 0x3fbd204698cb42bd
+        .quad 0x3fbcfc2df223db2d
+        .quad 0x3fbcd820ffd278f3
+        .quad 0x3fbcb41fba42686d
+        .quad 0x3fbc902a19e65111
+        .quad 0x3fbc6c4017382bea
+        .quad 0x3fbc4861aab93a23
+        .quad 0x3fbc248eccf1fba6
+        .quad 0x3fbc00c7767225cb
+        .quad 0x3fbbdd0b9fd09a10
+        .quad 0x3fbbb95b41ab5ce6
+        .quad 0x3fbb95b654a78c87
+        .quad 0x3fbb721cd17157e3
+        .quad 0x3fbb4e8eb0bbf58f
+        .quad 0x3fbb2b0beb419ad0
+        .quad 0x3fbb079479c372ad
+        .quad 0x3fbae4285509950b
+        .quad 0x3fbac0c775e2fde6
+        .quad 0x3fba9d71d5258484
+        .quad 0x3fba7a276badd2c8
+        .quad 0x3fba56e8325f5c87
+        .quad 0x3fba33b4222456f1
+        .quad 0x3fba108b33edb005
+        .quad 0x3fb9ed6d60b30612
+        .quad 0x3fb9ca5aa1729f45
+        .quad 0x3fb9a752ef316149
+        .quad 0x3fb9845642fac8f0
+        .quad 0x3fb9616495e0e1e8
+        .quad 0x3fb93e7de0fc3e80
+        .quad 0x3fb91ba21d6bef77
+        .quad 0x3fb8f8d144557bdf
+        .quad 0x3fb8d60b4ee4d901
+        .quad 0x3fb8b350364c6257
+        .quad 0x3fb8909ff3c4d191
+        .quad 0x3fb86dfa808d36a0
+        .quad 0x3fb84b5fd5eaefd8
+        .quad 0x3fb828cfed29a215
+        .quad 0x3fb8064abf9b30f1
+        .quad 0x3fb7e3d04697b704
+        .quad 0x3fb7c1607b7d7e32
+        .quad 0x3fb79efb57b0f803
+        .quad 0x3fb77ca0d49cb608
+        .quad 0x3fb75a50ebb1624a
+        .quad 0x3fb7380b9665b7c8
+        .quad 0x3fb715d0ce367afc
+        .quad 0x3fb6f3a08ca67270
+        .quad 0x3fb6d17acb3e5f5e
+        .quad 0x3fb6af5f838cf654
+        .quad 0x3fb68d4eaf26d7ee
+        .quad 0x3fb66b4847a68997
+        .quad 0x3fb6494c46ac6e4d
+        .quad 0x3fb6275aa5debf81
+        .quad 0x3fb605735ee985f1
+        .quad 0x3fb5e3966b7e9295
+        .quad 0x3fb5c1c3c5557799
+        .quad 0x3fb59ffb662b815c
+        .quad 0x3fb57e3d47c3af7b
+        .quad 0x3fb55c8963e6adeb
+        .quad 0x3fb53adfb462ce16
+        .quad 0x3fb51940330c000b
+        .quad 0x3fb4f7aad9bbcbaf
+        .quad 0x3fb4d61fa2514a00
+        .quad 0x3fb4b49e86b11e5f
+        .quad 0x3fb4932780c56fe2
+        .quad 0x3fb471ba8a7de2b7
+        .quad 0x3fb450579dcf9186
+        .quad 0x3fb42efeb4b506e9
+        .quad 0x3fb40dafc92e36e2
+        .quad 0x3fb3ec6ad5407868
+        .quad 0x3fb3cb2fd2f67ef1
+        .quad 0x3fb3a9febc60540a
+        .quad 0x3fb388d78b9350ff
+        .quad 0x3fb367ba3aaa1883
+        .quad 0x3fb346a6c3c49066
+        .quad 0x3fb3259d2107db54
+        .quad 0x3fb3049d4c9e52a0
+        .quad 0x3fb2e3a740b7800f
+        .quad 0x3fb2c2baf78817b7
+        .quad 0x3fb2a1d86b49f1e2
+        .quad 0x3fb280ff963c04fc
+        .quad 0x3fb2603072a25f82
+        .quad 0x3fb23f6afac6220a
+        .quad 0x3fb21eaf28f57941
+        .quad 0x3fb1fdfcf7839804
+        .quad 0x3fb1dd5460c8b16f
+        .quad 0x3fb1bcb55f21f307
+        .quad 0x3fb19c1fecf17ee0
+        .quad 0x3fb17b94049e65d0
+        .quad 0x3fb15b11a094a1aa
+        .quad 0x3fb13a98bb450f81
+        .quad 0x3fb11a294f2569f6
+        .quad 0x3fb0f9c356b04389
+        .quad 0x3fb0d966cc6500fa
+        .quad 0x3fb0b913aac7d3a7
+        .quad 0x3fb098c9ec61b3ff
+        .quad 0x3fb078898bc05bf4
+        .quad 0x3fb0585283764178
+        .quad 0x3fb03824ce1a9101
+        .quad 0x3fb0180066492817
+        .quad 0x3fafefca8d451fd6
+        .quad 0x3fafafa6d397efdb
+        .quad 0x3faf6f9594de60f0
+        .quad 0x3faf2f96c6754aee
+        .quad 0x3faeefaa5dc2b239
+        .quad 0x3faeafd05035bd3b
+        .quad 0x3fae70089346a9e6
+        .quad 0x3fae30531c76c34a
+        .quad 0x3fadf0afe1505738
+        .quad 0x3fadb11ed766abf4
+        .quad 0x3fad719ff455f5f7
+        .quad 0x3fad32332dc34dbd
+        .quad 0x3facf2d8795ca5a5
+        .quad 0x3facb38fccd8bfdb
+        .quad 0x3fac74591df72456
+        .quad 0x3fac3534628016dd
+        .quad 0x3fabf62190448d22
+        .quad 0x3fabb7209d1e24e5
+        .quad 0x3fab78317eef1a29
+        .quad 0x3fab39542ba23d73
+        .quad 0x3faafa88992aea19
+        .quad 0x3faabbcebd84fca0
+        .quad 0x3faa7d268eb4c924
+        .quad 0x3faa3e9002c711d2
+        .quad 0x3faa000b0fd0fd6b
+        .quad 0x3fa9c197abf00dd7
+        .quad 0x3fa98335cd4a16c3
+        .quad 0x3fa944e56a0d3450
+        .quad 0x3fa906a6786fc1cb
+        .quad 0x3fa8c878eeb05074
+        .quad 0x3fa88a5cc3159e53
+        .quad 0x3fa84c51ebee8d15
+        .quad 0x3fa80e585f9218fc
+        .quad 0x3fa7d070145f4fd7
+        .quad 0x3fa7929900bd4809
+        .quad 0x3fa754d31b1b179c
+        .quad 0x3fa7171e59efcb5f
+        .quad 0x3fa6d97ab3ba5e10
+        .quad 0x3fa69be81f01af99
+        .quad 0x3fa65e6692547c4e
+        .quad 0x3fa620f604495440
+        .quad 0x3fa5e3966b7e9295
+        .quad 0x3fa5a647be9a54f6
+        .quad 0x3fa56909f44a72fe
+        .quad 0x3fa52bdd034475b8
+        .quad 0x3fa4eec0e2458f30
+        .quad 0x3fa4b1b588129203
+        .quad 0x3fa474baeb77e904
+        .quad 0x3fa437d103498eec
+        .quad 0x3fa3faf7c663060e
+        .quad 0x3fa3be2f2ba7501f
+        .quad 0x3fa381772a00e604
+        .quad 0x3fa344cfb861afae
+        .quad 0x3fa30838cdc2fbfd
+        .quad 0x3fa2cbb2612578b4
+        .quad 0x3fa28f3c69912a74
+        .quad 0x3fa252d6de1564c1
+        .quad 0x3fa21681b5c8c213
+        .quad 0x3fa1da3ce7c91bf8
+        .quad 0x3fa19e086b3b8333
+        .quad 0x3fa161e4374c37f4
+        .quad 0x3fa125d0432ea20e
+        .quad 0x3fa0e9cc861d4944
+        .quad 0x3fa0add8f759cd95
+        .quad 0x3fa071f58e2cdf9b
+        .quad 0x3fa0362241e638ec
+        .quad 0x3f9ff4be13b92920
+        .quad 0x3f9f7d57badb4ee8
+        .quad 0x3f9f061167fc31e8
+        .quad 0x3f9e8eeb09f2f6cb
+        .quad 0x3f9e17e48fa48962
+        .quad 0x3f9da0fde8038de9
+        .quad 0x3f9d2a3702105259
+        .quad 0x3f9cb38fccd8bfdb
+        .quad 0x3f9c3d0837784c41
+        .quad 0x3f9bc6a03117eb97
+        .quad 0x3f9b5057a8ee01ce
+        .quad 0x3f9ada2e8e3e546f
+        .quad 0x3f9a6424d059fc68
+        .quad 0x3f99ee3a5e9f57e8
+        .quad 0x3f99786f2879fc53
+        .quad 0x3f9902c31d62a843
+        .quad 0x3f988d362cdf359e
+        .quad 0x3f9817c846828bbd
+        .quad 0x3f97a27959ec91aa
+        .quad 0x3f972d4956ca2067
+        .quad 0x3f96b8382cd4f551
+        .quad 0x3f964345cbd3a491
+        .quad 0x3f95ce7223998b98
+        .quad 0x3f9559bd2406c3ba
+        .quad 0x3f94e526bd0814d1
+        .quad 0x3f9470aede96e7f2
+        .quad 0x3f93fc5578b93a38
+        .quad 0x3f93881a7b818f9e
+        .quad 0x3f9313fdd70ee5e8
+        .quad 0x3f929fff7b8ca79d
+        .quad 0x3f922c1f59329f1b
+        .quad 0x3f91b85d6044e9ae
+        .quad 0x3f9144b98113eac0
+        .quad 0x3f90d133abfc3f1b
+        .quad 0x3f905dcbd166b033
+        .quad 0x3f8fd503c3904f1d
+        .quad 0x3f8eeeab9b43445d
+        .quad 0x3f8e088f0b004827
+        .quad 0x3f8d22adf3f9579d
+        .quad 0x3f8c3d0837784c41
+        .quad 0x3f8b579db6dec358
+        .quad 0x3f8a726e53a6056e
+        .quad 0x3f898d79ef5eedf0
+        .quad 0x3f88a8c06bb1d2f4
+        .quad 0x3f87c441aa5e6d15
+        .quad 0x3f86dffd8d3bbf70
+        .quad 0x3f85fbf3f637ffc5
+        .quad 0x3f851824c7587eb0
+        .quad 0x3f84348fe2b99002
+        .quad 0x3f8351352a8e733f
+        .quad 0x3f826e1481213c2e
+        .quad 0x3f818b2dc8d2bb91
+        .quad 0x3f80a880e41a67f6
+        .quad 0x3f7f8c1b6b0c8d4e
+        .quad 0x3f7dc7a83f75a96d
+        .quad 0x3f7c03a80ae5e054
+        .quad 0x3f7a401a92ff827e
+        .quad 0x3f787cff9d9147a5
+        .quad 0x3f76ba56f09621bc
+        .quad 0x3f74f8205235102d
+        .quad 0x3f73365b88c0f347
+        .quad 0x3f7175085ab85ff0
+        .quad 0x3f6f684d1d8ae702
+        .quad 0x3f6be76bd77b4fc3
+        .quad 0x3f68676c71434fb9
+        .quad 0x3f64e84e793a474a
+        .quad 0x3f616a117e0d4b30
+        .quad 0x3f5bd96a1d7d9cbc
+        .quad 0x3f54e071754c98ba
+        .quad 0x3f4bd27045bfd025
+        .quad 0x3f3bcef518e29612
+        .quad 0x8000000000000000
+        /*== poly_coeff[5] ==*/
+        .align 32
+        .quad 0x3fb63C65231FBD16, 0x3fb63C65231FBD16, 0x3fb63C65231FBD16, 0x3fb63C65231FBD16 /* coeff5 */
+        .quad 0xbfbBCB7D4EFBE80B, 0xbfbBCB7D4EFBE80B, 0xbfbBCB7D4EFBE80B, 0xbfbBCB7D4EFBE80B /* coeff4 */
+        .quad 0x3fc287A7636F341E, 0x3fc287A7636F341E, 0x3fc287A7636F341E, 0x3fc287A7636F341E /* coeff3 */
+        .quad 0xbfcBCB7B1526DE36, 0xbfcBCB7B1526DE36, 0xbfcBCB7B1526DE36, 0xbfcBCB7B1526DE36 /* coeff2 */
+        .quad 0x3fdBCB7B1526E50E, 0x3fdBCB7B1526E50E, 0x3fdBCB7B1526E50E, 0x3fdBCB7B1526E50E /* coeff1 */
+        /*== ExpMask ==*/
+        .align 32
+        .quad 0x000fffffffffffff, 0x000fffffffffffff, 0x000fffffffffffff, 0x000fffffffffffff
+        /*== Two10 ==*/
+        .align 32
+        .quad 0x3f50000000000000, 0x3f50000000000000, 0x3f50000000000000, 0x3f50000000000000
+        /*== MinNorm ==*/
+        .align 32
+        .quad 0x0010000000000000, 0x0010000000000000, 0x0010000000000000, 0x0010000000000000
+        /*== MaxNorm ==*/
+        .align 32
+        .quad 0x7fefffffffffffff, 0x7fefffffffffffff, 0x7fefffffffffffff, 0x7fefffffffffffff
+        /*== HalfMask ==*/
+        .align 32
+        .quad 0xfffffffffc000000, 0xfffffffffc000000, 0xfffffffffc000000, 0xfffffffffc000000
+        /*== One ==*/
+        .align 32
+        .quad 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000
+        /*== Threshold ==*/
+        .align 32
+        .quad 0x4086a00000000000, 0x4086a00000000000, 0x4086a00000000000, 0x4086a00000000000
+        /*== Bias ==*/
+        .align 32
+        .quad 0x408ff80000000000, 0x408ff80000000000, 0x408ff80000000000, 0x408ff80000000000
+        /*== Bias1 ==*/
+        .align 32
+        .quad 0x408ff00000000000, 0x408ff00000000000, 0x408ff00000000000, 0x408ff00000000000
+        /*== L2 ==*/
+        .align 32
+        .quad 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff
+        .align 32
+        .type	__svml_dlog10_data_internal,@object
+        .size	__svml_dlog10_data_internal,.-__svml_dlog10_data_internal
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core-avx2.S b/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core-avx2.S
new file mode 100644
index 0000000..3432e7c
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core-avx2.S
@@ -0,0 +1,20 @@
+/* AVX2 version of vectorized log10, vector length is 8.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define _ZGVeN8v_log10 _ZGVeN8v_log10_avx2_wrapper
+#include "../svml_d_log108_core.S"
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core.c b/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core.c
new file mode 100644
index 0000000..273a0d4
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core.c
@@ -0,0 +1,27 @@
+/* Multiple versions of vectorized log10, vector length is 8.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define SYMBOL_NAME _ZGVeN8v_log10
+#include "ifunc-mathvec-avx512-skx.h"
+
+libc_ifunc_redirected (REDIRECT_NAME, SYMBOL_NAME, IFUNC_SELECTOR ());
+
+#ifdef SHARED
+__hidden_ver1 (_ZGVeN8v_log10, __GI__ZGVeN8v_log10, __redirect__ZGVeN8v_log10)
+  __attribute__ ((visibility ("hidden")));
+#endif
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core_avx512.S b/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core_avx512.S
new file mode 100644
index 0000000..0799f99
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_d_log108_core_avx512.S
@@ -0,0 +1,299 @@
+/* Function log10 vectorized with AVX-512.
+   Copyright (C) 2021 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
+   https://www.gnu.org/licenses/.  */
+
+/*
+ * ALGORITHM DESCRIPTION:
+ *
+ *    Get short reciprocal approximation Rcp ~ 1/mantissa(x)
+ *    R = Rcp*x - 1.0
+ *    log10(x) = k*log10(2.0) - log10(Rcp) + poly_approximation(R)
+ *       log10(Rcp) is tabulated
+ *
+ *
+ */
+
+/* Offsets for data table __svml_dlog10_data_internal_avx512
+ */
+#define Log_tbl                       	0
+#define One                           	128
+#define C075                          	192
+#define poly_coeff9                   	256
+#define poly_coeff8                   	320
+#define poly_coeff7                   	384
+#define poly_coeff6                   	448
+#define poly_coeff5                   	512
+#define poly_coeff4                   	576
+#define poly_coeff3                   	640
+#define poly_coeff2                   	704
+#define poly_coeff1                   	768
+#define L2                            	832
+
+#include <sysdep.h>
+
+        .text
+	.section .text.evex512,"ax",@progbits
+ENTRY(_ZGVeN8v_log10_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, %zmm7
+        vgetmantpd $8, {sae}, %zmm7, %zmm6
+        vmovups   One+__svml_dlog10_data_internal_avx512(%rip), %zmm3
+        vmovups   poly_coeff5+__svml_dlog10_data_internal_avx512(%rip), %zmm12
+        vmovups   poly_coeff3+__svml_dlog10_data_internal_avx512(%rip), %zmm13
+
+/* Start polynomial evaluation */
+        vmovups   poly_coeff9+__svml_dlog10_data_internal_avx512(%rip), %zmm10
+        vmovups   poly_coeff8+__svml_dlog10_data_internal_avx512(%rip), %zmm1
+        vmovups   poly_coeff7+__svml_dlog10_data_internal_avx512(%rip), %zmm11
+        vmovups   poly_coeff6+__svml_dlog10_data_internal_avx512(%rip), %zmm14
+
+/* Prepare exponent correction: DblRcp<0.75? */
+        vmovups   C075+__svml_dlog10_data_internal_avx512(%rip), %zmm2
+
+/* Table lookup */
+        vmovups   __svml_dlog10_data_internal_avx512(%rip), %zmm5
+
+/* GetExp(x) */
+        vgetexppd {sae}, %zmm7, %zmm0
+
+/* DblRcp ~ 1/Mantissa */
+        vrcp14pd  %zmm6, %zmm8
+
+/* x<=0? */
+        vfpclasspd $94, %zmm7, %k0
+
+/* round DblRcp to 4 fractional bits (RN mode, no Precision exception) */
+        vrndscalepd $88, {sae}, %zmm8, %zmm4
+        vmovups   poly_coeff4+__svml_dlog10_data_internal_avx512(%rip), %zmm8
+        kmovw     %k0, %edx
+
+/* Reduced argument: R = DblRcp*Mantissa - 1 */
+        vfmsub213pd {rn-sae}, %zmm3, %zmm4, %zmm6
+        vcmppd    $17, {sae}, %zmm2, %zmm4, %k1
+        vfmadd231pd {rn-sae}, %zmm6, %zmm12, %zmm8
+        vmovups   poly_coeff2+__svml_dlog10_data_internal_avx512(%rip), %zmm12
+        vfmadd231pd {rn-sae}, %zmm6, %zmm10, %zmm1
+        vfmadd231pd {rn-sae}, %zmm6, %zmm11, %zmm14
+        vmovups   poly_coeff1+__svml_dlog10_data_internal_avx512(%rip), %zmm2
+
+/* R^2 */
+        vmulpd    {rn-sae}, %zmm6, %zmm6, %zmm15
+        vfmadd231pd {rn-sae}, %zmm6, %zmm13, %zmm12
+
+/* Prepare table index */
+        vpsrlq    $48, %zmm4, %zmm9
+
+/* add 1 to Expon if DblRcp<0.75 */
+        vaddpd    {rn-sae}, %zmm3, %zmm0, %zmm0{%k1}
+        vmulpd    {rn-sae}, %zmm15, %zmm15, %zmm13
+        vfmadd213pd {rn-sae}, %zmm14, %zmm15, %zmm1
+        vfmadd213pd {rn-sae}, %zmm12, %zmm15, %zmm8
+        vpermt2pd Log_tbl+64+__svml_dlog10_data_internal_avx512(%rip), %zmm9, %zmm5
+
+/* polynomial */
+        vfmadd213pd {rn-sae}, %zmm8, %zmm13, %zmm1
+        vfmadd213pd {rn-sae}, %zmm2, %zmm6, %zmm1
+        vfmadd213pd {rn-sae}, %zmm5, %zmm1, %zmm6
+        vmovups   L2+__svml_dlog10_data_internal_avx512(%rip), %zmm1
+        vfmadd213pd {rn-sae}, %zmm6, %zmm1, %zmm0
+        testl     %edx, %edx
+
+/* Go to special inputs processing branch */
+        jne       L(SPECIAL_VALUES_BRANCH)
+                                # LOE rbx r12 r13 r14 r15 edx zmm0 zmm7
+
+/* Restore registers
+ * and exit the function
+ */
+
+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)
+
+/* Branch to process
+ * special inputs
+ */
+
+L(SPECIAL_VALUES_BRANCH):
+        vmovups   %zmm7, 64(%rsp)
+        vmovups   %zmm0, 128(%rsp)
+                                # LOE rbx r12 r13 r14 r15 edx zmm0
+
+        xorl      %eax, %eax
+                                # LOE rbx r12 r13 r14 r15 eax edx
+
+        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
+ */
+
+L(RANGEMASK_CHECK):
+        btl       %r12d, %r13d
+
+/* Call scalar math function */
+        jc        L(SCALAR_MATH_CALL)
+                                # LOE rbx r15 r12d r13d
+
+/* Special inputs
+ * processing loop
+ */
+
+L(SPECIAL_VALUES_LOOP):
+        incl      %r12d
+        cmpl      $8, %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
+ */
+
+L(SCALAR_MATH_CALL):
+        movl      %r12d, %r14d
+        movsd     64(%rsp,%r14,8), %xmm0
+        call      log10@PLT
+                                # LOE rbx r14 r15 r12d r13d xmm0
+
+        movsd     %xmm0, 128(%rsp,%r14,8)
+
+/* Process special inputs in loop */
+        jmp       L(SPECIAL_VALUES_LOOP)
+                                # LOE rbx r15 r12d r13d
+END(_ZGVeN8v_log10_skx)
+
+        .section .rodata, "a"
+        .align 64
+
+#ifdef __svml_dlog10_data_internal_avx512_typedef
+typedef unsigned int VUINT32;
+typedef struct {
+        __declspec(align(64)) VUINT32 Log_tbl[16][2];
+        __declspec(align(64)) VUINT32 One[8][2];
+        __declspec(align(64)) VUINT32 C075[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff9[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff8[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff7[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff6[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff5[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff4[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff3[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff2[8][2];
+        __declspec(align(64)) VUINT32 poly_coeff1[8][2];
+        __declspec(align(64)) VUINT32 L2[8][2];
+   } __svml_dlog10_data_internal_avx512;
+#endif
+__svml_dlog10_data_internal_avx512:
+        /*== Log_tbl ==*/
+        .quad 0x0000000000000000
+        .quad 0xbf9af5f92b00e610
+        .quad 0xbfaa30a9d609efea
+        .quad 0xbfb31b3055c47118
+        .quad 0xbfb8cf183886480d
+        .quad 0xbfbe3bc1ab0e19fe
+        .quad 0xbfc1b3e71ec94f7b
+        .quad 0xbfc42c7e7fe3fc02
+        .quad 0x3fbffbfc2bbc7803
+        .quad 0x3fbb721cd17157e3
+        .quad 0x3fb715d0ce367afc
+        .quad 0x3fb2e3a740b7800f
+        .quad 0x3fadb11ed766abf4
+        .quad 0x3fa5e3966b7e9295
+        .quad 0x3f9cb38fccd8bfdb
+        .quad 0x3f8c3d0837784c41
+        /*== One ==*/
+        .align 64
+        .quad 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000, 0x3ff0000000000000
+        /*== 0.75 ==*/
+        .align 64
+        .quad 0x3fe8000000000000, 0x3fe8000000000000, 0x3fe8000000000000, 0x3fe8000000000000, 0x3fe8000000000000, 0x3fe8000000000000, 0x3fe8000000000000, 0x3fe8000000000000
+        /*== poly_coeff9 ==*/
+        .align 64
+        .quad 0x3fa8c2d828480370, 0x3fa8c2d828480370, 0x3fa8c2d828480370, 0x3fa8c2d828480370, 0x3fa8c2d828480370, 0x3fa8c2d828480370, 0x3fa8c2d828480370, 0x3fa8c2d828480370
+        /*== poly_coeff8 ==*/
+        .align 64
+        .quad 0xbfabd80d96029814, 0xbfabd80d96029814, 0xbfabd80d96029814, 0xbfabd80d96029814, 0xbfabd80d96029814, 0xbfabd80d96029814, 0xbfabd80d96029814, 0xbfabd80d96029814
+        /*== poly_coeff7 ==*/
+        .align 64
+        .quad 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2, 0x3fafc3f6f38b58a2
+        /*== poly_coeff6 ==*/
+        .align 64
+        .quad 0xbfb287a63464dc80, 0xbfb287a63464dc80, 0xbfb287a63464dc80, 0xbfb287a63464dc80, 0xbfb287a63464dc80, 0xbfb287a63464dc80, 0xbfb287a63464dc80, 0xbfb287a63464dc80
+        /*== poly_coeff5 ==*/
+        .align 64
+        .quad 0x3fb63c62777f27d9, 0x3fb63c62777f27d9, 0x3fb63c62777f27d9, 0x3fb63c62777f27d9, 0x3fb63c62777f27d9, 0x3fb63c62777f27d9, 0x3fb63c62777f27d9, 0x3fb63c62777f27d9
+        /*== poly_coeff4 ==*/
+        .align 64
+        .quad 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3, 0xbfbbcb7b153c06a3
+        /*== poly_coeff3 ==*/
+        .align 64
+        .quad 0x3fc287a7636f428c, 0x3fc287a7636f428c, 0x3fc287a7636f428c, 0x3fc287a7636f428c, 0x3fc287a7636f428c, 0x3fc287a7636f428c, 0x3fc287a7636f428c, 0x3fc287a7636f428c
+        /*== poly_coeff2 ==*/
+        .align 64
+        .quad 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db, 0xbfcbcb7b1526e4db
+        /*== poly_coeff1 ==*/
+        .align 64
+        .quad 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e, 0x3fdbcb7b1526e50e
+        /*== L2 ==*/
+        .align 64
+        .quad 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff, 0x3fd34413509f79ff
+        .align 64
+        .type	__svml_dlog10_data_internal_avx512,@object
+        .size	__svml_dlog10_data_internal_avx512,.-__svml_dlog10_data_internal_avx512
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core-avx2.S b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core-avx2.S
new file mode 100644
index 0000000..e389e2e
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core-avx2.S
@@ -0,0 +1,20 @@
+/* AVX2 version of vectorized log10f.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define _ZGVeN16v_log10f _ZGVeN16v_log10f_avx2_wrapper
+#include "../svml_s_log10f16_core.S"
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core.c b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core.c
new file mode 100644
index 0000000..274fc7e
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core.c
@@ -0,0 +1,28 @@
+/* Multiple versions of vectorized log10f, vector length is 16.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define SYMBOL_NAME _ZGVeN16v_log10f
+#include "ifunc-mathvec-avx512-skx.h"
+
+libc_ifunc_redirected (REDIRECT_NAME, SYMBOL_NAME, IFUNC_SELECTOR ());
+
+#ifdef SHARED
+__hidden_ver1 (_ZGVeN16v_log10f, __GI__ZGVeN16v_log10f,
+	       __redirect__ZGVeN16v_log10f)
+  __attribute__ ((visibility ("hidden")));
+#endif
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core_avx512.S b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core_avx512.S
new file mode 100644
index 0000000..3dffd66
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f16_core_avx512.S
@@ -0,0 +1,238 @@
+/* Function log10f vectorized with AVX-512.
+   Copyright (C) 2021 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
+   https://www.gnu.org/licenses/.  */
+
+/*
+ * ALGORITHM DESCRIPTION:
+ *
+ *    Get short reciprocal approximation Rcp ~ 1/mantissa(x)
+ *    R = Rcp*x - 1.0
+ *    log10(x) = k*log10(2.0) - log10(Rcp) + poly_approximation(R)
+ *       log10(Rcp) is tabulated
+ *
+ *
+ */
+
+/* Offsets for data table __svml_slog10_data_internal_avx512
+ */
+#define One                           	0
+#define coeff4                        	64
+#define coeff3                        	128
+#define coeff2                        	192
+#define coeff1                        	256
+#define L2                            	320
+
+#include <sysdep.h>
+
+        .text
+	.section .text.exex512,"ax",@progbits
+ENTRY(_ZGVeN16v_log10f_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
+        vgetmantps $11, {sae}, %zmm0, %zmm3
+        vmovups   __svml_slog10_data_internal_avx512(%rip), %zmm1
+        vgetexpps {sae}, %zmm0, %zmm5
+        vmovups   L2+__svml_slog10_data_internal_avx512(%rip), %zmm10
+        vpsrld    $19, %zmm3, %zmm7
+        vgetexpps {sae}, %zmm3, %zmm6
+        vsubps    {rn-sae}, %zmm1, %zmm3, %zmm11
+        vpermps   coeff4+__svml_slog10_data_internal_avx512(%rip), %zmm7, %zmm1
+        vpermps   coeff3+__svml_slog10_data_internal_avx512(%rip), %zmm7, %zmm2
+        vsubps    {rn-sae}, %zmm6, %zmm5, %zmm9
+        vpermps   coeff2+__svml_slog10_data_internal_avx512(%rip), %zmm7, %zmm4
+        vpermps   coeff1+__svml_slog10_data_internal_avx512(%rip), %zmm7, %zmm8
+
+/* x<=0? */
+        vfpclassps $94, %zmm0, %k0
+        vfmadd213ps {rn-sae}, %zmm2, %zmm11, %zmm1
+        vmulps    {rn-sae}, %zmm10, %zmm9, %zmm12
+        vfmadd213ps {rn-sae}, %zmm4, %zmm11, %zmm1
+        kmovw     %k0, %edx
+        vfmadd213ps {rn-sae}, %zmm8, %zmm11, %zmm1
+        vfmadd213ps {rn-sae}, %zmm12, %zmm11, %zmm1
+        testl     %edx, %edx
+
+/* Go to special inputs processing branch */
+        jne       L(SPECIAL_VALUES_BRANCH)
+                                # LOE rbx r12 r13 r14 r15 edx zmm0 zmm1
+
+/* Restore registers
+ * and exit the function
+ */
+
+L(EXIT):
+        vmovaps   %zmm1, %zmm0
+        movq      %rbp, %rsp
+        popq      %rbp
+        cfi_def_cfa(7, 8)
+        cfi_restore(6)
+        ret
+        cfi_def_cfa(6, 16)
+        cfi_offset(6, -16)
+
+/* Branch to process
+ * special inputs
+ */
+
+L(SPECIAL_VALUES_BRANCH):
+        vmovups   %zmm0, 64(%rsp)
+        vmovups   %zmm1, 128(%rsp)
+                                # LOE rbx r12 r13 r14 r15 edx zmm1
+
+        xorl      %eax, %eax
+                                # LOE rbx r12 r13 r14 r15 eax edx
+
+        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
+ */
+
+L(RANGEMASK_CHECK):
+        btl       %r12d, %r13d
+
+/* Call scalar math function */
+        jc        L(SCALAR_MATH_CALL)
+                                # LOE rbx r15 r12d r13d
+
+/* Special inputs
+ * processing loop
+ */
+
+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), %zmm1
+
+/* 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 zmm1
+
+/* Scalar math fucntion call
+ * to process special input
+ */
+
+L(SCALAR_MATH_CALL):
+        movl      %r12d, %r14d
+        movss     64(%rsp,%r14,4), %xmm0
+        call      log10f@PLT
+                                # LOE rbx r14 r15 r12d r13d xmm0
+
+        movss     %xmm0, 128(%rsp,%r14,4)
+
+/* Process special inputs in loop */
+        jmp       L(SPECIAL_VALUES_LOOP)
+                                # LOE rbx r15 r12d r13d
+END(_ZGVeN16v_log10f_skx)
+
+        .section .rodata, "a"
+        .align 64
+
+#ifdef __svml_slog10_data_internal_avx512_typedef
+typedef unsigned int VUINT32;
+typedef struct {
+        __declspec(align(64)) VUINT32 One[16][1];
+        __declspec(align(64)) VUINT32 coeff4[16][1];
+        __declspec(align(64)) VUINT32 coeff3[16][1];
+        __declspec(align(64)) VUINT32 coeff2[16][1];
+        __declspec(align(64)) VUINT32 coeff1[16][1];
+        __declspec(align(64)) VUINT32 L2[16][1];
+    } __svml_slog10_data_internal_avx512;
+#endif
+__svml_slog10_data_internal_avx512:
+        /*== One ==*/
+        .long 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000
+        // c4
+        .align 64
+        .long 0xbdc9ae9b, 0xbda6fcf4
+        .long 0xbd8bac76, 0xbd6bca30
+        .long 0xbd48a99b, 0xbd2c0a9f
+        .long 0xbd1480db, 0xbd00faf2
+        .long 0xbe823aa9, 0xbe656348
+        .long 0xbe4afbb9, 0xbe346895
+        .long 0xbe20ffff, 0xbe103a0b
+        .long 0xbe01a91c, 0xbde9e84e
+        // c3
+        .align 64
+        .long 0x3e13d888, 0x3e10a87c
+        .long 0x3e0b95c3, 0x3e057f0b
+        .long 0x3dfde038, 0x3df080d9
+        .long 0x3de34c1e, 0x3dd68333
+        .long 0x3dac6e8e, 0x3dd54a51
+        .long 0x3df30f40, 0x3e04235d
+        .long 0x3e0b7033, 0x3e102c90
+        .long 0x3e12ebad, 0x3e141ff8
+        // c2
+        .align 64
+        .long 0xbe5e5a9b, 0xbe5e2677
+        .long 0xbe5d83f5, 0xbe5c6016
+        .long 0xbe5abd0b, 0xbe58a6fd
+        .long 0xbe562e02, 0xbe5362f8
+        .long 0xbe68e27c, 0xbe646747
+        .long 0xbe619a73, 0xbe5ff05a
+        .long 0xbe5f0570, 0xbe5e92d0
+        .long 0xbe5e662b, 0xbe5e5c08
+        // c1
+        .align 64
+        .long 0x3ede5bd8, 0x3ede5b45
+        .long 0x3ede57d8, 0x3ede4eb1
+        .long 0x3ede3d37, 0x3ede2166
+        .long 0x3eddf9d9, 0x3eddc5bb
+        .long 0x3ede08ed, 0x3ede32e7
+        .long 0x3ede4967, 0x3ede5490
+        .long 0x3ede597f, 0x3ede5b50
+        .long 0x3ede5bca, 0x3ede5bd9
+        /*== L2 ==*/
+        .align 64
+        .long 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b
+        .align 64
+        .type	__svml_slog10_data_internal_avx512,@object
+        .size	__svml_slog10_data_internal_avx512,.-__svml_slog10_data_internal_avx512
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core-sse2.S b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core-sse2.S
new file mode 100644
index 0000000..bb1cdee
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core-sse2.S
@@ -0,0 +1,20 @@
+/* SSE2 version of vectorized log10f, vector length is 4.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define _ZGVbN4v_log10f _ZGVbN4v_log10f_sse2
+#include "../svml_s_log10f4_core.S"
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core.c b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core.c
new file mode 100644
index 0000000..67e9e71
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core.c
@@ -0,0 +1,28 @@
+/* Multiple versions of vectorized log10f, vector length is 4.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define SYMBOL_NAME _ZGVbN4v_log10f
+#include "ifunc-mathvec-sse4_1.h"
+
+libc_ifunc_redirected (REDIRECT_NAME, SYMBOL_NAME, IFUNC_SELECTOR ());
+
+#ifdef SHARED
+__hidden_ver1 (_ZGVbN4v_log10f, __GI__ZGVbN4v_log10f,
+	       __redirect__ZGVbN4v_log10f)
+  __attribute__ ((visibility ("hidden")));
+#endif
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core_sse4.S b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core_sse4.S
new file mode 100644
index 0000000..88b3535
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f4_core_sse4.S
@@ -0,0 +1,243 @@
+/* Function log10f vectorized with SSE4.
+   Copyright (C) 2021 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
+   https://www.gnu.org/licenses/.  */
+
+/*
+ * ALGORITHM DESCRIPTION:
+ *
+ *    Get short reciprocal approximation Rcp ~ 1/mantissa(x)
+ *    R = Rcp*x - 1.0
+ *    log10(x) = k*log10(2.0) - log10(Rcp) + poly_approximation(R)
+ *       log10(Rcp) is tabulated
+ *
+ *
+ */
+
+/* Offsets for data table __svml_slog10_data_internal
+ */
+#define MinNorm                       	0
+#define MaxNorm                       	16
+#define L2H                           	32
+#define L2L                           	48
+#define iBrkValue                     	64
+#define iOffExpoMask                  	80
+#define One                           	96
+#define sPoly                         	112
+#define L2                            	256
+
+#include <sysdep.h>
+
+        .text
+	.section .text.sse4,"ax",@progbits
+ENTRY(_ZGVbN4v_log10f_sse4)
+        subq      $72, %rsp
+        cfi_def_cfa_offset(80)
+        movaps    %xmm0, %xmm1
+
+/* reduction: compute r,n */
+        movdqu    iBrkValue+__svml_slog10_data_internal(%rip), %xmm2
+        movaps    %xmm0, %xmm4
+        movdqu    iOffExpoMask+__svml_slog10_data_internal(%rip), %xmm10
+        psubd     %xmm2, %xmm1
+        pand      %xmm1, %xmm10
+        psrad     $23, %xmm1
+        paddd     %xmm2, %xmm10
+        movaps    %xmm0, %xmm3
+        movups    sPoly+__svml_slog10_data_internal(%rip), %xmm5
+        movups    sPoly+32+__svml_slog10_data_internal(%rip), %xmm6
+        movups    sPoly+64+__svml_slog10_data_internal(%rip), %xmm7
+        movups    sPoly+96+__svml_slog10_data_internal(%rip), %xmm9
+        cvtdq2ps  %xmm1, %xmm12
+        cmpltps   MinNorm+__svml_slog10_data_internal(%rip), %xmm4
+        cmpnleps  MaxNorm+__svml_slog10_data_internal(%rip), %xmm3
+        subps     One+__svml_slog10_data_internal(%rip), %xmm10
+        mulps     %xmm10, %xmm5
+        movaps    %xmm10, %xmm8
+        mulps     %xmm10, %xmm6
+        mulps     %xmm10, %xmm8
+        addps     sPoly+16+__svml_slog10_data_internal(%rip), %xmm5
+        mulps     %xmm10, %xmm7
+        addps     sPoly+48+__svml_slog10_data_internal(%rip), %xmm6
+        mulps     %xmm10, %xmm9
+        mulps     %xmm8, %xmm5
+        addps     sPoly+80+__svml_slog10_data_internal(%rip), %xmm7
+        addps     sPoly+112+__svml_slog10_data_internal(%rip), %xmm9
+        addps     %xmm5, %xmm6
+        mulps     %xmm8, %xmm6
+        orps      %xmm3, %xmm4
+
+/* combine and get argument value range mask */
+        movmskps  %xmm4, %edx
+        movups    L2L+__svml_slog10_data_internal(%rip), %xmm1
+        addps     %xmm6, %xmm7
+        mulps     %xmm12, %xmm1
+        mulps     %xmm7, %xmm8
+        movups    L2H+__svml_slog10_data_internal(%rip), %xmm11
+        addps     %xmm8, %xmm9
+        mulps     %xmm11, %xmm12
+        mulps     %xmm10, %xmm9
+        addps     sPoly+128+__svml_slog10_data_internal(%rip), %xmm9
+        mulps     %xmm9, %xmm10
+        addps     %xmm10, %xmm1
+        addps     %xmm12, %xmm1
+        testl     %edx, %edx
+
+/* Go to special inputs processing branch */
+        jne       L(SPECIAL_VALUES_BRANCH)
+                                # LOE rbx rbp r12 r13 r14 r15 edx xmm0 xmm1
+
+/* Restore registers
+ * and exit the function
+ */
+
+L(EXIT):
+        movaps    %xmm1, %xmm0
+        addq      $72, %rsp
+        cfi_def_cfa_offset(8)
+        ret
+        cfi_def_cfa_offset(80)
+
+/* Branch to process
+ * special inputs
+ */
+
+L(SPECIAL_VALUES_BRANCH):
+        movups    %xmm0, 32(%rsp)
+        movups    %xmm1, 48(%rsp)
+                                # LOE rbx rbp r12 r13 r14 r15 edx
+
+        xorl      %eax, %eax
+        movq      %r12, 16(%rsp)
+        cfi_offset(12, -64)
+        movl      %eax, %r12d
+        movq      %r13, 8(%rsp)
+        cfi_offset(13, -72)
+        movl      %edx, %r13d
+        movq      %r14, (%rsp)
+        cfi_offset(14, -80)
+                                # LOE rbx rbp r15 r12d r13d
+
+/* Range mask
+ * bits check
+ */
+
+L(RANGEMASK_CHECK):
+        btl       %r12d, %r13d
+
+/* Call scalar math function */
+        jc        L(SCALAR_MATH_CALL)
+                                # LOE rbx rbp r15 r12d r13d
+
+/* Special inputs
+ * processing loop
+ */
+
+L(SPECIAL_VALUES_LOOP):
+        incl      %r12d
+        cmpl      $4, %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    48(%rsp), %xmm1
+
+/* Go to exit */
+        jmp       L(EXIT)
+        cfi_offset(12, -64)
+        cfi_offset(13, -72)
+        cfi_offset(14, -80)
+                                # LOE rbx rbp r12 r13 r14 r15 xmm1
+
+/* Scalar math fucntion call
+ * to process special input
+ */
+
+L(SCALAR_MATH_CALL):
+        movl      %r12d, %r14d
+        movss     32(%rsp,%r14,4), %xmm0
+        call      log10f@PLT
+                                # LOE rbx rbp r14 r15 r12d r13d xmm0
+
+        movss     %xmm0, 48(%rsp,%r14,4)
+
+/* Process special inputs in loop */
+        jmp       L(SPECIAL_VALUES_LOOP)
+                                # LOE rbx rbp r15 r12d r13d
+END(_ZGVbN4v_log10f_sse4)
+
+        .section .rodata, "a"
+        .align 16
+
+#ifdef __svml_slog10_data_internal_typedef
+typedef unsigned int VUINT32;
+typedef struct {
+        __declspec(align(16)) VUINT32 MinNorm[4][1];
+        __declspec(align(16)) VUINT32 MaxNorm[4][1];
+        __declspec(align(16)) VUINT32 L2H[4][1];
+        __declspec(align(16)) VUINT32 L2L[4][1];
+        __declspec(align(16)) VUINT32 iBrkValue[4][1];
+        __declspec(align(16)) VUINT32 iOffExpoMask[4][1];
+        __declspec(align(16)) VUINT32 One[4][1];
+        __declspec(align(16)) VUINT32 sPoly[9][4][1];
+        __declspec(align(16)) VUINT32 L2[4][1];
+} __svml_slog10_data_internal;
+#endif
+__svml_slog10_data_internal:
+        /*== MinNorm ==*/
+        .long 0x00800000, 0x00800000, 0x00800000, 0x00800000
+        /*== MaxNorm ==*/
+        .align 16
+        .long 0x7f7fffff, 0x7f7fffff, 0x7f7fffff, 0x7f7fffff
+        /*== L2H ==*/
+        .align 16
+        .long 0x3e9a2100, 0x3e9a2100, 0x3e9a2100, 0x3e9a2100
+        /*== L2L ==*/
+        .align 16
+        .long 0xb64AF600, 0xb64AF600, 0xb64AF600, 0xb64AF600
+        /*== iBrkValue = SP 2/3 ==*/
+        .align 16
+        .long 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab
+        /*== iOffExpoMask = SP significand mask ==*/
+        .align 16
+        .long 0x007fffff, 0x007fffff, 0x007fffff, 0x007fffff
+        /*== sOne = SP 1.0 ==*/
+        .align 16
+        .long 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000
+        /*== spoly[9] ==*/
+        .align 16
+        .long 0x3d8063B4, 0x3d8063B4, 0x3d8063B4, 0x3d8063B4 /* coeff9 */
+        .long 0xbd890073, 0xbd890073, 0xbd890073, 0xbd890073 /* coeff8 */
+        .long 0x3d775317, 0x3d775317, 0x3d775317, 0x3d775317 /* coeff7 */
+        .long 0xbd91FB27, 0xbd91FB27, 0xbd91FB27, 0xbd91FB27 /* coeff6 */
+        .long 0x3dB20B96, 0x3dB20B96, 0x3dB20B96, 0x3dB20B96 /* coeff5 */
+        .long 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20 /* coeff4 */
+        .long 0x3e143CE5, 0x3e143CE5, 0x3e143CE5, 0x3e143CE5 /* coeff3 */
+        .long 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5 /* coeff2 */
+        .long 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9 /* coeff1 */
+        /*== L2 ==*/
+        .align 16
+        .long 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b
+        .align 16
+        .type	__svml_slog10_data_internal,@object
+        .size	__svml_slog10_data_internal,.-__svml_slog10_data_internal
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core-sse.S b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core-sse.S
new file mode 100644
index 0000000..e3467e5
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core-sse.S
@@ -0,0 +1,20 @@
+/* SSE version of vectorized log10f, vector length is 8.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define _ZGVdN8v_log10f _ZGVdN8v_log10f_sse_wrapper
+#include "../svml_s_log10f8_core.S"
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core.c b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core.c
new file mode 100644
index 0000000..bfd3ef6
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core.c
@@ -0,0 +1,28 @@
+/* Multiple versions of vectorized log10f, vector length is 8.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#define SYMBOL_NAME _ZGVdN8v_log10f
+#include "ifunc-mathvec-avx2.h"
+
+libc_ifunc_redirected (REDIRECT_NAME, SYMBOL_NAME, IFUNC_SELECTOR ());
+
+#ifdef SHARED
+__hidden_ver1 (_ZGVdN8v_log10f, __GI__ZGVdN8v_log10f,
+	       __redirect__ZGVdN8v_log10f)
+  __attribute__ ((visibility ("hidden")));
+#endif
diff --git a/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core_avx2.S b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core_avx2.S
new file mode 100644
index 0000000..58e2634
--- /dev/null
+++ b/sysdeps/x86_64/fpu/multiarch/svml_s_log10f8_core_avx2.S
@@ -0,0 +1,243 @@
+/* Function log10f vectorized with AVX2.
+   Copyright (C) 2021 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
+   https://www.gnu.org/licenses/.  */
+
+/*
+ * ALGORITHM DESCRIPTION:
+ *
+ *    Get short reciprocal approximation Rcp ~ 1/mantissa(x)
+ *    R = Rcp*x - 1.0
+ *    log10(x) = k*log10(2.0) - log10(Rcp) + poly_approximation(R)
+ *       log10(Rcp) is tabulated
+ *
+ *
+ */
+
+/* Offsets for data table __svml_slog10_data_internal
+ */
+#define MinNorm                       	0
+#define MaxNorm                       	32
+#define L2H                           	64
+#define L2L                           	96
+#define iBrkValue                     	128
+#define iOffExpoMask                  	160
+#define One                           	192
+#define sPoly                         	224
+#define L2                            	512
+
+#include <sysdep.h>
+
+        .text
+	.section .text.avx2,"ax",@progbits
+ENTRY(_ZGVdN8v_log10f_avx2)
+        pushq     %rbp
+        cfi_def_cfa_offset(16)
+        movq      %rsp, %rbp
+        cfi_def_cfa(6, 16)
+        cfi_offset(6, -16)
+        andq      $-32, %rsp
+        subq      $96, %rsp
+
+/* reduction: compute r,n */
+        vmovups   iBrkValue+__svml_slog10_data_internal(%rip), %ymm4
+        vmovups   sPoly+__svml_slog10_data_internal(%rip), %ymm15
+        vmovups   sPoly+64+__svml_slog10_data_internal(%rip), %ymm9
+        vmovups   sPoly+128+__svml_slog10_data_internal(%rip), %ymm10
+        vmovups   sPoly+192+__svml_slog10_data_internal(%rip), %ymm12
+        vpsubd    %ymm4, %ymm0, %ymm1
+        vcmplt_oqps MinNorm+__svml_slog10_data_internal(%rip), %ymm0, %ymm5
+        vcmpnle_uqps MaxNorm+__svml_slog10_data_internal(%rip), %ymm0, %ymm6
+        vpand     iOffExpoMask+__svml_slog10_data_internal(%rip), %ymm1, %ymm3
+        vpsrad    $23, %ymm1, %ymm2
+        vpaddd    %ymm4, %ymm3, %ymm8
+        vcvtdq2ps %ymm2, %ymm1
+        vsubps    One+__svml_slog10_data_internal(%rip), %ymm8, %ymm13
+        vmulps    L2L+__svml_slog10_data_internal(%rip), %ymm1, %ymm14
+        vfmadd213ps sPoly+32+__svml_slog10_data_internal(%rip), %ymm13, %ymm15
+        vfmadd213ps sPoly+96+__svml_slog10_data_internal(%rip), %ymm13, %ymm9
+        vmulps    %ymm13, %ymm13, %ymm11
+        vfmadd213ps sPoly+160+__svml_slog10_data_internal(%rip), %ymm13, %ymm10
+        vfmadd213ps sPoly+224+__svml_slog10_data_internal(%rip), %ymm13, %ymm12
+        vfmadd213ps %ymm9, %ymm11, %ymm15
+        vfmadd213ps %ymm10, %ymm11, %ymm15
+        vfmadd213ps %ymm12, %ymm11, %ymm15
+        vfmadd213ps sPoly+256+__svml_slog10_data_internal(%rip), %ymm13, %ymm15
+        vfmadd213ps %ymm14, %ymm13, %ymm15
+        vorps     %ymm6, %ymm5, %ymm7
+
+/* combine and get argument value range mask */
+        vmovmskps %ymm7, %edx
+        vfmadd132ps L2H+__svml_slog10_data_internal(%rip), %ymm15, %ymm1
+        testl     %edx, %edx
+
+/* Go to special inputs processing branch */
+        jne       L(SPECIAL_VALUES_BRANCH)
+                                # LOE rbx r12 r13 r14 r15 edx ymm0 ymm1
+
+/* Restore registers
+ * and exit the function
+ */
+
+L(EXIT):
+        vmovaps   %ymm1, %ymm0
+        movq      %rbp, %rsp
+        popq      %rbp
+        cfi_def_cfa(7, 8)
+        cfi_restore(6)
+        ret
+        cfi_def_cfa(6, 16)
+        cfi_offset(6, -16)
+
+/* Branch to process
+ * special inputs
+ */
+
+L(SPECIAL_VALUES_BRANCH):
+        vmovups   %ymm0, 32(%rsp)
+        vmovups   %ymm1, 64(%rsp)
+                                # LOE rbx r12 r13 r14 r15 edx ymm1
+
+        xorl      %eax, %eax
+                                # LOE rbx r12 r13 r14 r15 eax edx
+
+        vzeroupper
+        movq      %r12, 16(%rsp)
+        /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -80; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xb0, 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: -32; DW_OP_and; DW_OP_const4s: -88; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa8, 0xff, 0xff, 0xff, 0x22
+        movl      %edx, %r13d
+        movq      %r14, (%rsp)
+        /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -96; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa0, 0xff, 0xff, 0xff, 0x22
+                                # LOE rbx r15 r12d r13d
+
+/* Range mask
+ * bits check
+ */
+
+L(RANGEMASK_CHECK):
+        btl       %r12d, %r13d
+
+/* Call scalar math function */
+        jc        L(SCALAR_MATH_CALL)
+                                # LOE rbx r15 r12d r13d
+
+/* Special inputs
+ * processing loop
+ */
+
+L(SPECIAL_VALUES_LOOP):
+        incl      %r12d
+        cmpl      $8, %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   64(%rsp), %ymm1
+
+/* Go to exit */
+        jmp       L(EXIT)
+        /*  DW_CFA_expression: r12 (r12) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -80; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0c, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xb0, 0xff, 0xff, 0xff, 0x22
+        /*  DW_CFA_expression: r13 (r13) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -88; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0d, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa8, 0xff, 0xff, 0xff, 0x22
+        /*  DW_CFA_expression: r14 (r14) (DW_OP_lit8; DW_OP_minus; DW_OP_const4s: -32; DW_OP_and; DW_OP_const4s: -96; DW_OP_plus)  */
+        .cfi_escape 0x10, 0x0e, 0x0e, 0x38, 0x1c, 0x0d, 0xe0, 0xff, 0xff, 0xff, 0x1a, 0x0d, 0xa0, 0xff, 0xff, 0xff, 0x22
+                                # LOE rbx r12 r13 r14 r15 ymm1
+
+/* Scalar math fucntion call
+ * to process special input
+ */
+
+L(SCALAR_MATH_CALL):
+        movl      %r12d, %r14d
+        movss     32(%rsp,%r14,4), %xmm0
+        call      log10f@PLT
+                                # LOE rbx r14 r15 r12d r13d xmm0
+
+        movss     %xmm0, 64(%rsp,%r14,4)
+
+/* Process special inputs in loop */
+        jmp       L(SPECIAL_VALUES_LOOP)
+                                # LOE rbx r15 r12d r13d
+END(_ZGVdN8v_log10f_avx2)
+
+        .section .rodata, "a"
+        .align 32
+
+#ifdef __svml_slog10_data_internal_typedef
+typedef unsigned int VUINT32;
+typedef struct {
+        __declspec(align(32)) VUINT32 MinNorm[8][1];
+        __declspec(align(32)) VUINT32 MaxNorm[8][1];
+        __declspec(align(32)) VUINT32 L2H[8][1];
+        __declspec(align(32)) VUINT32 L2L[8][1];
+        __declspec(align(32)) VUINT32 iBrkValue[8][1];
+        __declspec(align(32)) VUINT32 iOffExpoMask[8][1];
+        __declspec(align(32)) VUINT32 One[8][1];
+        __declspec(align(32)) VUINT32 sPoly[9][8][1];
+        __declspec(align(32)) VUINT32 L2[8][1];
+} __svml_slog10_data_internal;
+#endif
+__svml_slog10_data_internal:
+        /*== MinNorm ==*/
+        .long 0x00800000, 0x00800000, 0x00800000, 0x00800000, 0x00800000, 0x00800000, 0x00800000, 0x00800000
+        /*== MaxNorm ==*/
+        .align 32
+        .long 0x7f7fffff, 0x7f7fffff, 0x7f7fffff, 0x7f7fffff, 0x7f7fffff, 0x7f7fffff, 0x7f7fffff, 0x7f7fffff
+        /*== L2H ==*/
+        .align 32
+        .long 0x3e9a2100, 0x3e9a2100, 0x3e9a2100, 0x3e9a2100, 0x3e9a2100, 0x3e9a2100, 0x3e9a2100, 0x3e9a2100
+        /*== L2L ==*/
+        .align 32
+        .long 0xb64AF600, 0xb64AF600, 0xb64AF600, 0xb64AF600, 0xb64AF600, 0xb64AF600, 0xb64AF600, 0xb64AF600
+        /*== iBrkValue = SP 2/3 ==*/
+        .align 32
+        .long 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab, 0x3f2aaaab
+        /*== iOffExpoMask = SP significand mask ==*/
+        .align 32
+        .long 0x007fffff, 0x007fffff, 0x007fffff, 0x007fffff, 0x007fffff, 0x007fffff, 0x007fffff, 0x007fffff
+        /*== sOne = SP 1.0 ==*/
+        .align 32
+        .long 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000
+        /*== spoly[9] ==*/
+        .align 32
+        .long 0x3d8063B4, 0x3d8063B4, 0x3d8063B4, 0x3d8063B4, 0x3d8063B4, 0x3d8063B4, 0x3d8063B4, 0x3d8063B4 /* coeff9 */
+        .long 0xbd890073, 0xbd890073, 0xbd890073, 0xbd890073, 0xbd890073, 0xbd890073, 0xbd890073, 0xbd890073 /* coeff8 */
+        .long 0x3d775317, 0x3d775317, 0x3d775317, 0x3d775317, 0x3d775317, 0x3d775317, 0x3d775317, 0x3d775317 /* coeff7 */
+        .long 0xbd91FB27, 0xbd91FB27, 0xbd91FB27, 0xbd91FB27, 0xbd91FB27, 0xbd91FB27, 0xbd91FB27, 0xbd91FB27 /* coeff6 */
+        .long 0x3dB20B96, 0x3dB20B96, 0x3dB20B96, 0x3dB20B96, 0x3dB20B96, 0x3dB20B96, 0x3dB20B96, 0x3dB20B96 /* coeff5 */
+        .long 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20, 0xbdDE6E20 /* coeff4 */
+        .long 0x3e143CE5, 0x3e143CE5, 0x3e143CE5, 0x3e143CE5, 0x3e143CE5, 0x3e143CE5, 0x3e143CE5, 0x3e143CE5 /* coeff3 */
+        .long 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5, 0xbe5E5BC5 /* coeff2 */
+        .long 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9, 0x3eDE5BD9 /* coeff1 */
+        /*== L2 ==*/
+        .align 32
+        .long 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b, 0x3e9a209b
+        .align 32
+        .type	__svml_slog10_data_internal,@object
+        .size	__svml_slog10_data_internal,.-__svml_slog10_data_internal
diff --git a/sysdeps/x86_64/fpu/svml_d_log102_core.S b/sysdeps/x86_64/fpu/svml_d_log102_core.S
new file mode 100644
index 0000000..3d0c058
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_d_log102_core.S
@@ -0,0 +1,29 @@
+/* Function log10 vectorized with SSE2.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_d_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVbN2v_log10)
+WRAPPER_IMPL_SSE2 log10
+END (_ZGVbN2v_log10)
+
+#ifndef USE_MULTIARCH
+ libmvec_hidden_def (_ZGVbN2v_log10)
+#endif
diff --git a/sysdeps/x86_64/fpu/svml_d_log104_core.S b/sysdeps/x86_64/fpu/svml_d_log104_core.S
new file mode 100644
index 0000000..9e32c62
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_d_log104_core.S
@@ -0,0 +1,29 @@
+/* Function log10 vectorized with AVX2, wrapper version.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_d_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVdN4v_log10)
+WRAPPER_IMPL_AVX _ZGVbN2v_log10
+END (_ZGVdN4v_log10)
+
+#ifndef USE_MULTIARCH
+ libmvec_hidden_def (_ZGVdN4v_log10)
+#endif
diff --git a/sysdeps/x86_64/fpu/svml_d_log104_core_avx.S b/sysdeps/x86_64/fpu/svml_d_log104_core_avx.S
new file mode 100644
index 0000000..2b073b1
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_d_log104_core_avx.S
@@ -0,0 +1,25 @@
+/* Function log10 vectorized in AVX ISA as wrapper to SSE4 ISA version.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_d_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVcN4v_log10)
+WRAPPER_IMPL_AVX _ZGVbN2v_log10
+END (_ZGVcN4v_log10)
diff --git a/sysdeps/x86_64/fpu/svml_d_log108_core.S b/sysdeps/x86_64/fpu/svml_d_log108_core.S
new file mode 100644
index 0000000..853d791
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_d_log108_core.S
@@ -0,0 +1,25 @@
+/* Function log10 vectorized with AVX-512, wrapper to AVX2.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_d_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVeN8v_log10)
+WRAPPER_IMPL_AVX512 _ZGVdN4v_log10
+END (_ZGVeN8v_log10)
diff --git a/sysdeps/x86_64/fpu/svml_s_log10f16_core.S b/sysdeps/x86_64/fpu/svml_s_log10f16_core.S
new file mode 100644
index 0000000..769603c
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_s_log10f16_core.S
@@ -0,0 +1,25 @@
+/* Function log10f vectorized with AVX-512. Wrapper to AVX2 version.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_s_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVeN16v_log10f)
+WRAPPER_IMPL_AVX512 _ZGVdN8v_log10f
+END (_ZGVeN16v_log10f)
diff --git a/sysdeps/x86_64/fpu/svml_s_log10f4_core.S b/sysdeps/x86_64/fpu/svml_s_log10f4_core.S
new file mode 100644
index 0000000..5235254
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_s_log10f4_core.S
@@ -0,0 +1,29 @@
+/* Function log10f vectorized with SSE2, wrapper version.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_s_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVbN4v_log10f)
+WRAPPER_IMPL_SSE2 log10f
+END (_ZGVbN4v_log10f)
+
+#ifndef USE_MULTIARCH
+ libmvec_hidden_def (_ZGVbN4v_log10f)
+#endif
diff --git a/sysdeps/x86_64/fpu/svml_s_log10f8_core.S b/sysdeps/x86_64/fpu/svml_s_log10f8_core.S
new file mode 100644
index 0000000..630ec76
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_s_log10f8_core.S
@@ -0,0 +1,29 @@
+/* Function log10f vectorized with AVX2, wrapper version.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_s_wrapper_impl.h"
+
+	.text
+ENTRY (_ZGVdN8v_log10f)
+WRAPPER_IMPL_AVX _ZGVbN4v_log10f
+END (_ZGVdN8v_log10f)
+
+#ifndef USE_MULTIARCH
+ libmvec_hidden_def (_ZGVdN8v_log10f)
+#endif
diff --git a/sysdeps/x86_64/fpu/svml_s_log10f8_core_avx.S b/sysdeps/x86_64/fpu/svml_s_log10f8_core_avx.S
new file mode 100644
index 0000000..374208c
--- /dev/null
+++ b/sysdeps/x86_64/fpu/svml_s_log10f8_core_avx.S
@@ -0,0 +1,25 @@
+/* Function log10f vectorized in AVX ISA as wrapper to SSE4 ISA version.
+   Copyright (C) 2021 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
+   <https://www.gnu.org/licenses/>.  */
+
+#include <sysdep.h>
+#include "svml_s_wrapper_impl.h"
+
+        .text
+ENTRY (_ZGVcN8v_log10f)
+WRAPPER_IMPL_AVX _ZGVbN4v_log10f
+END (_ZGVcN8v_log10f)
diff --git a/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx.c b/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx.c
new file mode 100644
index 0000000..770fd72
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx.c
@@ -0,0 +1 @@
+#include "test-double-libmvec-log10.c"
diff --git a/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx2.c b/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx2.c
new file mode 100644
index 0000000..770fd72
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx2.c
@@ -0,0 +1 @@
+#include "test-double-libmvec-log10.c"
diff --git a/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx512f.c b/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx512f.c
new file mode 100644
index 0000000..770fd72
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-double-libmvec-log10-avx512f.c
@@ -0,0 +1 @@
+#include "test-double-libmvec-log10.c"
diff --git a/sysdeps/x86_64/fpu/test-double-libmvec-log10.c b/sysdeps/x86_64/fpu/test-double-libmvec-log10.c
new file mode 100644
index 0000000..cb1ab36
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-double-libmvec-log10.c
@@ -0,0 +1,3 @@
+#define LIBMVEC_TYPE double
+#define LIBMVEC_FUNC log10
+#include "test-vector-abi-arg1.h"
diff --git a/sysdeps/x86_64/fpu/test-double-vlen2-wrappers.c b/sysdeps/x86_64/fpu/test-double-vlen2-wrappers.c
index 37a7a1c..3dce136 100644
--- a/sysdeps/x86_64/fpu/test-double-vlen2-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-double-vlen2-wrappers.c
@@ -38,6 +38,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1), _ZGVbN2v_expm1)
 VECTOR_WRAPPER (WRAPPER_NAME (sinh), _ZGVbN2v_sinh)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrt), _ZGVbN2v_cbrt)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2), _ZGVbN2vv_atan2)
+VECTOR_WRAPPER (WRAPPER_NAME (log10), _ZGVbN2v_log10)
 
 #define VEC_INT_TYPE __m128i
 
diff --git a/sysdeps/x86_64/fpu/test-double-vlen4-avx2-wrappers.c b/sysdeps/x86_64/fpu/test-double-vlen4-avx2-wrappers.c
index 4313f67..1852625 100644
--- a/sysdeps/x86_64/fpu/test-double-vlen4-avx2-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-double-vlen4-avx2-wrappers.c
@@ -41,6 +41,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1), _ZGVdN4v_expm1)
 VECTOR_WRAPPER (WRAPPER_NAME (sinh), _ZGVdN4v_sinh)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrt), _ZGVdN4v_cbrt)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2), _ZGVdN4vv_atan2)
+VECTOR_WRAPPER (WRAPPER_NAME (log10), _ZGVdN4v_log10)
 
 #ifndef __ILP32__
 # define VEC_INT_TYPE __m256i
diff --git a/sysdeps/x86_64/fpu/test-double-vlen4-wrappers.c b/sysdeps/x86_64/fpu/test-double-vlen4-wrappers.c
index 4b8b00f..cf9ea35 100644
--- a/sysdeps/x86_64/fpu/test-double-vlen4-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-double-vlen4-wrappers.c
@@ -38,6 +38,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1), _ZGVcN4v_expm1)
 VECTOR_WRAPPER (WRAPPER_NAME (sinh), _ZGVcN4v_sinh)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrt), _ZGVcN4v_cbrt)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2), _ZGVcN4vv_atan2)
+VECTOR_WRAPPER (WRAPPER_NAME (log10), _ZGVcN4v_log10)
 
 #define VEC_INT_TYPE __m128i
 
diff --git a/sysdeps/x86_64/fpu/test-double-vlen8-wrappers.c b/sysdeps/x86_64/fpu/test-double-vlen8-wrappers.c
index d06522a..b6457ea 100644
--- a/sysdeps/x86_64/fpu/test-double-vlen8-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-double-vlen8-wrappers.c
@@ -38,6 +38,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1), _ZGVeN8v_expm1)
 VECTOR_WRAPPER (WRAPPER_NAME (sinh), _ZGVeN8v_sinh)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrt), _ZGVeN8v_cbrt)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2), _ZGVeN8vv_atan2)
+VECTOR_WRAPPER (WRAPPER_NAME (log10), _ZGVeN8v_log10)
 
 #ifndef __ILP32__
 # define VEC_INT_TYPE __m512i
diff --git a/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx.c b/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx.c
new file mode 100644
index 0000000..04f017f
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx.c
@@ -0,0 +1 @@
+#include "test-float-libmvec-log10f.c"
diff --git a/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx2.c b/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx2.c
new file mode 100644
index 0000000..04f017f
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx2.c
@@ -0,0 +1 @@
+#include "test-float-libmvec-log10f.c"
diff --git a/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx512f.c b/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx512f.c
new file mode 100644
index 0000000..04f017f
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-float-libmvec-log10f-avx512f.c
@@ -0,0 +1 @@
+#include "test-float-libmvec-log10f.c"
diff --git a/sysdeps/x86_64/fpu/test-float-libmvec-log10f.c b/sysdeps/x86_64/fpu/test-float-libmvec-log10f.c
new file mode 100644
index 0000000..682ce1e
--- /dev/null
+++ b/sysdeps/x86_64/fpu/test-float-libmvec-log10f.c
@@ -0,0 +1,3 @@
+#define LIBMVEC_TYPE float
+#define LIBMVEC_FUNC log10f
+#include "test-vector-abi-arg1.h"
diff --git a/sysdeps/x86_64/fpu/test-float-vlen16-wrappers.c b/sysdeps/x86_64/fpu/test-float-vlen16-wrappers.c
index 0bd631b..272e754 100644
--- a/sysdeps/x86_64/fpu/test-float-vlen16-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-float-vlen16-wrappers.c
@@ -38,6 +38,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1f), _ZGVeN16v_expm1f)
 VECTOR_WRAPPER (WRAPPER_NAME (sinhf), _ZGVeN16v_sinhf)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrtf), _ZGVeN16v_cbrtf)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2f), _ZGVeN16vv_atan2f)
+VECTOR_WRAPPER (WRAPPER_NAME (log10f), _ZGVeN16v_log10f)
 
 #define VEC_INT_TYPE __m512i
 
diff --git a/sysdeps/x86_64/fpu/test-float-vlen4-wrappers.c b/sysdeps/x86_64/fpu/test-float-vlen4-wrappers.c
index 1018398..b892258 100644
--- a/sysdeps/x86_64/fpu/test-float-vlen4-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-float-vlen4-wrappers.c
@@ -38,6 +38,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1f), _ZGVbN4v_expm1f)
 VECTOR_WRAPPER (WRAPPER_NAME (sinhf), _ZGVbN4v_sinhf)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrtf), _ZGVbN4v_cbrtf)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2f), _ZGVbN4vv_atan2f)
+VECTOR_WRAPPER (WRAPPER_NAME (log10f), _ZGVbN4v_log10f)
 
 #define VEC_INT_TYPE __m128i
 
diff --git a/sysdeps/x86_64/fpu/test-float-vlen8-avx2-wrappers.c b/sysdeps/x86_64/fpu/test-float-vlen8-avx2-wrappers.c
index 42ea28f..1c6ead7 100644
--- a/sysdeps/x86_64/fpu/test-float-vlen8-avx2-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-float-vlen8-avx2-wrappers.c
@@ -41,6 +41,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1f), _ZGVdN8v_expm1f)
 VECTOR_WRAPPER (WRAPPER_NAME (sinhf), _ZGVdN8v_sinhf)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrtf), _ZGVdN8v_cbrtf)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2f), _ZGVdN8vv_atan2f)
+VECTOR_WRAPPER (WRAPPER_NAME (log10f), _ZGVdN8v_log10f)
 
 /* Redefinition of wrapper to be compatible with _ZGVdN8vvv_sincosf.  */
 #undef VECTOR_WRAPPER_fFF
diff --git a/sysdeps/x86_64/fpu/test-float-vlen8-wrappers.c b/sysdeps/x86_64/fpu/test-float-vlen8-wrappers.c
index 70a0216..71f5d8d 100644
--- a/sysdeps/x86_64/fpu/test-float-vlen8-wrappers.c
+++ b/sysdeps/x86_64/fpu/test-float-vlen8-wrappers.c
@@ -38,6 +38,7 @@ VECTOR_WRAPPER (WRAPPER_NAME (expm1f), _ZGVcN8v_expm1f)
 VECTOR_WRAPPER (WRAPPER_NAME (sinhf), _ZGVcN8v_sinhf)
 VECTOR_WRAPPER (WRAPPER_NAME (cbrtf), _ZGVcN8v_cbrtf)
 VECTOR_WRAPPER_ff (WRAPPER_NAME (atan2f), _ZGVcN8vv_atan2f)
+VECTOR_WRAPPER (WRAPPER_NAME (log10f), _ZGVcN8v_log10f)
 
 #define VEC_INT_TYPE __m128i