diff --git a/include/private/dsp/arch/x86/avx512/mix.h b/include/private/dsp/arch/x86/avx512/mix.h
index 68de9e3..93b3e8e 100644
--- a/include/private/dsp/arch/x86/avx512/mix.h
+++ b/include/private/dsp/arch/x86/avx512/mix.h
@@ -663,6 +663,680 @@ namespace lsp
             );
         }
 
+        void mix4(float *a, const float *b, const float *c, const float *d, float k1, float k2, float k3, float k4, size_t count)
+        {
+            IF_ARCH_X86(size_t off);
+            ARCH_X86_ASM
+            (
+                // 64x blocks
+                __ASM_EMIT("xor             %[off], %[off]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("vbroadcastss    %[k1], %%zmm6")
+                __ASM_EMIT("vbroadcastss    %[k2], %%zmm7")
+                __ASM_EMIT("vbroadcastss    %[k3], %%zmm5")
+                __ASM_EMIT("vbroadcastss    %[k4], %%zmm4")
+                __ASM_EMIT("jb              2f")
+                __ASM_EMIT("1:")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a],%[off]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vmulps          0x080(%[a],%[off]), %%zmm6, %%zmm2")
+                __ASM_EMIT("vmulps          0x0c0(%[a],%[off]), %%zmm6, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b],%[off]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[b],%[off]), %%zmm7, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[b],%[off]), %%zmm7, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c],%[off]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[c],%[off]), %%zmm5, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[c],%[off]), %%zmm5, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d],%[off]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[d],%[off]), %%zmm4, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[d],%[off]), %%zmm4, %%zmm3")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[a],%[off])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[a],%[off])")
+                __ASM_EMIT("vmovups         %%zmm2, 0x080(%[a],%[off])")
+                __ASM_EMIT("vmovups         %%zmm3, 0x0c0(%[a],%[off])")
+                __ASM_EMIT("add             $0x100, %[off]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("jae             1b")
+                // 32x block
+                __ASM_EMIT("2:")
+                __ASM_EMIT("add             $32, %[count]")
+                __ASM_EMIT("jl              4f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a],%[off]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b],%[off]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c],%[off]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d],%[off]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[a],%[off])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[a],%[off])")
+                __ASM_EMIT("sub             $32, %[count]")
+                __ASM_EMIT("add             $0x80, %[off]")
+                // 16x block
+                __ASM_EMIT("4:")
+                __ASM_EMIT("add             $16, %[count]")
+                __ASM_EMIT("jl              6f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%ymm6, %%ymm0")
+                __ASM_EMIT("vmulps          0x020(%[a],%[off]), %%ymm6, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%ymm7, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[b],%[off]), %%ymm7, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%ymm5, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[c],%[off]), %%ymm5, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%ymm4, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[d],%[off]), %%ymm4, %%ymm1")
+                __ASM_EMIT("vmovups         %%ymm0, 0x000(%[a],%[off])")
+                __ASM_EMIT("vmovups         %%ymm1, 0x020(%[a],%[off])")
+                __ASM_EMIT("sub             $16, %[count]")
+                __ASM_EMIT("add             $0x040, %[off]")
+                // 8x block
+                __ASM_EMIT("6:")
+                __ASM_EMIT("add             $8, %[count]")
+                __ASM_EMIT("jl              8f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vmulps          0x010(%[a],%[off]), %%xmm6, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[b],%[off]), %%xmm7, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[c],%[off]), %%xmm5, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[d],%[off]), %%xmm4, %%xmm1")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[a],%[off])")
+                __ASM_EMIT("vmovups         %%xmm1, 0x010(%[a],%[off])")
+                __ASM_EMIT("sub             $8, %[count]")
+                __ASM_EMIT("add             $0x020, %[off]")
+                // 4x block
+                __ASM_EMIT("8:")
+                __ASM_EMIT("add             $4, %[count]")
+                __ASM_EMIT("jl              10f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[a],%[off])")
+                __ASM_EMIT("sub             $4, %[count]")
+                __ASM_EMIT("add             $0x010, %[off]")
+                // 1x blocks
+                __ASM_EMIT("10:")
+                __ASM_EMIT("add             $3, %[count]")
+                __ASM_EMIT("jl              12f")
+                __ASM_EMIT("11:")
+                __ASM_EMIT("vmulss          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovss          %%xmm0, 0x00(%[a],%[off])")
+                __ASM_EMIT("add             $0x04, %[off]")
+                __ASM_EMIT("dec             %[count]")
+                __ASM_EMIT("jge             11b")
+                // End
+                __ASM_EMIT("12:")
+                : [count] "+r" (count),
+                  [off] "=&r" (off)
+                : [a] "r" (a), [b] "r" (b), [c] "r" (c), [d] "r" (d),
+                  [k1] "m" (k1), [k2] "m" (k2), [k3] "m" (k3), [k4] "m" (k4)
+                : "cc", "memory",
+                  "%xmm0", "%xmm1", "%xmm2", "%xmm3",
+                  "%xmm4", "%xmm5", "%xmm6", "%xmm7"
+            );
+        }
+
+    #if defined(ARCH_I386)
+        void mix_copy4(float *dst, const float *a, const float *b, const float *c, const float *d, float k1, float k2, float k3, float k4, size_t count)
+        {
+            ARCH_X86_ASM
+            (
+                // 64x blocks
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("vbroadcastss    %[k1], %%zmm6")
+                __ASM_EMIT("vbroadcastss    %[k2], %%zmm7")
+                __ASM_EMIT("vbroadcastss    %[k3], %%zmm5")
+                __ASM_EMIT("vbroadcastss    %[k4], %%zmm4")
+                __ASM_EMIT("jb              2f")
+                __ASM_EMIT("1:")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vmulps          0x080(%[a]), %%zmm6, %%zmm2")
+                __ASM_EMIT("vmulps          0x0c0(%[a]), %%zmm6, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[b]), %%zmm7, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[b]), %%zmm7, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[c]), %%zmm5, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[c]), %%zmm5, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[d]), %%zmm4, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[d]), %%zmm4, %%zmm3")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm2, 0x080(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm3, 0x0c0(%[dst])")
+                __ASM_EMIT("add             $0x100, %[a]")
+                __ASM_EMIT("add             $0x100, %[b]")
+                __ASM_EMIT("add             $0x100, %[c]")
+                __ASM_EMIT("add             $0x100, %[d]")
+                __ASM_EMIT("add             $0x100, %[dst]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("jae             1b")
+                // 32x block
+                __ASM_EMIT("2:")
+                __ASM_EMIT("add             $32, %[count]")
+                __ASM_EMIT("jl              4f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst])")
+                __ASM_EMIT("sub             $32, %[count]")
+                __ASM_EMIT("add             $0x080, %[a]")
+                __ASM_EMIT("add             $0x080, %[b]")
+                __ASM_EMIT("add             $0x080, %[c]")
+                __ASM_EMIT("add             $0x080, %[d]")
+                __ASM_EMIT("add             $0x080, %[dst]")
+                // 16x block
+                __ASM_EMIT("4:")
+                __ASM_EMIT("add             $16, %[count]")
+                __ASM_EMIT("jl              6f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%ymm6, %%ymm0")
+                __ASM_EMIT("vmulps          0x020(%[a]), %%ymm6, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%ymm7, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[b]), %%ymm7, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%ymm5, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[c]), %%ymm5, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%ymm4, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[d]), %%ymm4, %%ymm1")
+                __ASM_EMIT("vmovups         %%ymm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%ymm1, 0x020(%[dst])")
+                __ASM_EMIT("sub             $16, %[count]")
+                __ASM_EMIT("add             $0x040, %[a]")
+                __ASM_EMIT("add             $0x040, %[b]")
+                __ASM_EMIT("add             $0x040, %[c]")
+                __ASM_EMIT("add             $0x040, %[d]")
+                __ASM_EMIT("add             $0x040, %[dst]")
+                // 8x block
+                __ASM_EMIT("6:")
+                __ASM_EMIT("add             $8, %[count]")
+                __ASM_EMIT("jl              8f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vmulps          0x010(%[a]), %%xmm6, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[b]), %%xmm7, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[c]), %%xmm5, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[d]), %%xmm4, %%xmm1")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%xmm1, 0x010(%[dst])")
+                __ASM_EMIT("sub             $8, %[count]")
+                __ASM_EMIT("add             $0x020, %[a]")
+                __ASM_EMIT("add             $0x020, %[b]")
+                __ASM_EMIT("add             $0x020, %[c]")
+                __ASM_EMIT("add             $0x020, %[d]")
+                __ASM_EMIT("add             $0x020, %[dst]")
+                // 4x block
+                __ASM_EMIT("8:")
+                __ASM_EMIT("add             $4, %[count]")
+                __ASM_EMIT("jl              10f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst])")
+                __ASM_EMIT("sub             $4, %[count]")
+                __ASM_EMIT("add             $0x010, %[a]")
+                __ASM_EMIT("add             $0x010, %[b]")
+                __ASM_EMIT("add             $0x010, %[c]")
+                __ASM_EMIT("add             $0x010, %[d]")
+                __ASM_EMIT("add             $0x010, %[dst]")
+                // 1x blocks
+                __ASM_EMIT("10:")
+                __ASM_EMIT("add             $3, %[count]")
+                __ASM_EMIT("jl              12f")
+                __ASM_EMIT("11:")
+                __ASM_EMIT("vmulss          0x000(%[a]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[b]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[c]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[d]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovss          %%xmm0, 0x00(%[dst])")
+                __ASM_EMIT("add             $0x004, %[a]")
+                __ASM_EMIT("add             $0x004, %[b]")
+                __ASM_EMIT("add             $0x004, %[c]")
+                __ASM_EMIT("add             $0x004, %[d]")
+                __ASM_EMIT("add             $0x004, %[dst]")
+                __ASM_EMIT("dec             %[count]")
+                __ASM_EMIT("jge             11b")
+                // End
+                __ASM_EMIT("12:")
+                : [count] "+r" (count),
+                  [dst] "+r" (dst), [a] "+r" (a), [b] "+r" (b), [c] "+r" (c), [d] "+r" (d)
+                : [k1] "m" (k1), [k2] "m" (k2), [k3] "m" (k3), [k4] "m" (k4)
+                : "cc", "memory",
+                  "%xmm0", "%xmm1", "%xmm2", "%xmm3",
+                  "%xmm4", "%xmm5", "%xmm6", "%xmm7"
+            );
+        }
+
+        void mix_add4(float *dst, const float *a, const float *b, const float *c, const float *d, float k1, float k2, float k3, float k4, size_t count)
+        {
+            ARCH_X86_ASM
+            (
+                // 64x blocks
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("vbroadcastss    %[k1], %%zmm6")
+                __ASM_EMIT("vbroadcastss    %[k2], %%zmm7")
+                __ASM_EMIT("vbroadcastss    %[k3], %%zmm5")
+                __ASM_EMIT("vbroadcastss    %[k4], %%zmm4")
+                __ASM_EMIT("jb              2f")
+                __ASM_EMIT("1:")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vmulps          0x080(%[a]), %%zmm6, %%zmm2")
+                __ASM_EMIT("vmulps          0x0c0(%[a]), %%zmm6, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[b]), %%zmm7, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[b]), %%zmm7, %%zmm3")
+                __ASM_EMIT("vaddps          0x000(%[dst]), %%zmm0, %%zmm0")
+                __ASM_EMIT("vaddps          0x040(%[dst]), %%zmm1, %%zmm1")
+                __ASM_EMIT("vaddps          0x080(%[dst]), %%zmm2, %%zmm2")
+                __ASM_EMIT("vaddps          0x0c0(%[dst]), %%zmm3, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[c]), %%zmm5, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[c]), %%zmm5, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[d]), %%zmm4, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[d]), %%zmm4, %%zmm3")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm2, 0x080(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm3, 0x0c0(%[dst])")
+                __ASM_EMIT("add             $0x100, %[a]")
+                __ASM_EMIT("add             $0x100, %[b]")
+                __ASM_EMIT("add             $0x100, %[c]")
+                __ASM_EMIT("add             $0x100, %[d]")
+                __ASM_EMIT("add             $0x100, %[dst]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("jae             1b")
+                // 32x block
+                __ASM_EMIT("2:")
+                __ASM_EMIT("add             $32, %[count]")
+                __ASM_EMIT("jl              4f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vaddps          0x000(%[dst]), %%zmm0, %%zmm0")
+                __ASM_EMIT("vaddps          0x040(%[dst]), %%zmm1, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst])")
+                __ASM_EMIT("sub             $32, %[count]")
+                __ASM_EMIT("add             $0x080, %[a]")
+                __ASM_EMIT("add             $0x080, %[b]")
+                __ASM_EMIT("add             $0x080, %[c]")
+                __ASM_EMIT("add             $0x080, %[d]")
+                __ASM_EMIT("add             $0x080, %[dst]")
+                // 16x block
+                __ASM_EMIT("4:")
+                __ASM_EMIT("add             $16, %[count]")
+                __ASM_EMIT("jl              6f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%ymm6, %%ymm0")
+                __ASM_EMIT("vmulps          0x020(%[a]), %%ymm6, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%ymm7, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[b]), %%ymm7, %%ymm1")
+                __ASM_EMIT("vaddps          0x000(%[dst]), %%ymm0, %%ymm0")
+                __ASM_EMIT("vaddps          0x020(%[dst]), %%ymm1, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%ymm5, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[c]), %%ymm5, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%ymm4, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[d]), %%ymm4, %%ymm1")
+                __ASM_EMIT("vmovups         %%ymm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%ymm1, 0x020(%[dst])")
+                __ASM_EMIT("sub             $16, %[count]")
+                __ASM_EMIT("add             $0x040, %[a]")
+                __ASM_EMIT("add             $0x040, %[b]")
+                __ASM_EMIT("add             $0x040, %[c]")
+                __ASM_EMIT("add             $0x040, %[d]")
+                __ASM_EMIT("add             $0x040, %[dst]")
+                // 8x block
+                __ASM_EMIT("6:")
+                __ASM_EMIT("add             $8, %[count]")
+                __ASM_EMIT("jl              8f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vmulps          0x010(%[a]), %%xmm6, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[b]), %%xmm7, %%xmm1")
+                __ASM_EMIT("vaddps          0x000(%[dst]), %%xmm0, %%xmm0")
+                __ASM_EMIT("vaddps          0x010(%[dst]), %%xmm1, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[c]), %%xmm5, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[d]), %%xmm4, %%xmm1")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst])")
+                __ASM_EMIT("vmovups         %%xmm1, 0x010(%[dst])")
+                __ASM_EMIT("sub             $8, %[count]")
+                __ASM_EMIT("add             $0x020, %[a]")
+                __ASM_EMIT("add             $0x020, %[b]")
+                __ASM_EMIT("add             $0x020, %[c]")
+                __ASM_EMIT("add             $0x020, %[d]")
+                __ASM_EMIT("add             $0x020, %[dst]")
+                // 4x block
+                __ASM_EMIT("8:")
+                __ASM_EMIT("add             $4, %[count]")
+                __ASM_EMIT("jl              10f")
+                __ASM_EMIT("vmulps          0x000(%[a]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vaddps          0x000(%[dst]), %%xmm0, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst])")
+                __ASM_EMIT("sub             $4, %[count]")
+                __ASM_EMIT("add             $0x010, %[a]")
+                __ASM_EMIT("add             $0x010, %[b]")
+                __ASM_EMIT("add             $0x010, %[c]")
+                __ASM_EMIT("add             $0x010, %[d]")
+                __ASM_EMIT("add             $0x010, %[dst]")
+                // 1x blocks
+                __ASM_EMIT("10:")
+                __ASM_EMIT("add             $3, %[count]")
+                __ASM_EMIT("jl              12f")
+                __ASM_EMIT("11:")
+                __ASM_EMIT("vmulss          0x000(%[a]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[b]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vaddss          0x000(%[dst]), %%xmm0, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[c]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[d]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovss          %%xmm0, 0x00(%[dst])")
+                __ASM_EMIT("add             $0x004, %[a]")
+                __ASM_EMIT("add             $0x004, %[b]")
+                __ASM_EMIT("add             $0x004, %[c]")
+                __ASM_EMIT("add             $0x004, %[d]")
+                __ASM_EMIT("add             $0x004, %[dst]")
+                __ASM_EMIT("dec             %[count]")
+                __ASM_EMIT("jge             11b")
+                // End
+                __ASM_EMIT("12:")
+                : [count] "+r" (count),
+                  [dst] "+r" (dst), [a] "+r" (a), [b] "+r" (b), [c] "+r" (c), [d] "+r" (d)
+                : [k1] "m" (k1), [k2] "m" (k2), [k3] "m" (k3), [k4] "m" (k4)
+                : "cc", "memory",
+                  "%xmm0", "%xmm1", "%xmm2", "%xmm3",
+                  "%xmm4", "%xmm5", "%xmm6", "%xmm7"
+            );
+        }
+
+    #else
+        void mix_copy4(float *dst, const float *a, const float *b, const float *c, const float *d, float k1, float k2, float k3, float k4, size_t count)
+        {
+            IF_ARCH_X86(size_t off);
+            ARCH_X86_ASM
+            (
+                // 64x blocks
+                __ASM_EMIT("xor             %[off], %[off]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("vbroadcastss    %[k1], %%zmm6")
+                __ASM_EMIT("vbroadcastss    %[k2], %%zmm7")
+                __ASM_EMIT("vbroadcastss    %[k3], %%zmm5")
+                __ASM_EMIT("vbroadcastss    %[k4], %%zmm4")
+                __ASM_EMIT("jb              2f")
+                __ASM_EMIT("1:")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a],%[off]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vmulps          0x080(%[a],%[off]), %%zmm6, %%zmm2")
+                __ASM_EMIT("vmulps          0x0c0(%[a],%[off]), %%zmm6, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b],%[off]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[b],%[off]), %%zmm7, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[b],%[off]), %%zmm7, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c],%[off]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[c],%[off]), %%zmm5, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[c],%[off]), %%zmm5, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d],%[off]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[d],%[off]), %%zmm4, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[d],%[off]), %%zmm4, %%zmm3")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm2, 0x080(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm3, 0x0c0(%[dst],%[off])")
+                __ASM_EMIT("add             $0x100, %[off]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("jae             1b")
+                // 32x block
+                __ASM_EMIT("2:")
+                __ASM_EMIT("add             $32, %[count]")
+                __ASM_EMIT("jl              4f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a],%[off]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b],%[off]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c],%[off]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d],%[off]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst],%[off])")
+                __ASM_EMIT("sub             $32, %[count]")
+                __ASM_EMIT("add             $0x80, %[off]")
+                // 16x block
+                __ASM_EMIT("4:")
+                __ASM_EMIT("add             $16, %[count]")
+                __ASM_EMIT("jl              6f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%ymm6, %%ymm0")
+                __ASM_EMIT("vmulps          0x020(%[a],%[off]), %%ymm6, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%ymm7, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[b],%[off]), %%ymm7, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%ymm5, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[c],%[off]), %%ymm5, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%ymm4, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[d],%[off]), %%ymm4, %%ymm1")
+                __ASM_EMIT("vmovups         %%ymm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%ymm1, 0x020(%[dst],%[off])")
+                __ASM_EMIT("sub             $16, %[count]")
+                __ASM_EMIT("add             $0x040, %[off]")
+                // 8x block
+                __ASM_EMIT("6:")
+                __ASM_EMIT("add             $8, %[count]")
+                __ASM_EMIT("jl              8f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vmulps          0x010(%[a],%[off]), %%xmm6, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[b],%[off]), %%xmm7, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[c],%[off]), %%xmm5, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[d],%[off]), %%xmm4, %%xmm1")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%xmm1, 0x010(%[dst],%[off])")
+                __ASM_EMIT("sub             $8, %[count]")
+                __ASM_EMIT("add             $0x020, %[off]")
+                // 4x block
+                __ASM_EMIT("8:")
+                __ASM_EMIT("add             $4, %[count]")
+                __ASM_EMIT("jl              10f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("sub             $4, %[count]")
+                __ASM_EMIT("add             $0x010, %[off]")
+                // 1x blocks
+                __ASM_EMIT("10:")
+                __ASM_EMIT("add             $3, %[count]")
+                __ASM_EMIT("jl              12f")
+                __ASM_EMIT("11:")
+                __ASM_EMIT("vmulss          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovss          %%xmm0, 0x00(%[dst],%[off])")
+                __ASM_EMIT("add             $0x04, %[off]")
+                __ASM_EMIT("dec             %[count]")
+                __ASM_EMIT("jge             11b")
+                // End
+                __ASM_EMIT("12:")
+                : [count] "+r" (count),
+                  [off] "=&r" (off)
+                : [dst] "r" (dst), [a] "r" (a), [b] "r" (b), [c] "r" (c), [d] "r" (d),
+                  [k1] "m" (k1), [k2] "m" (k2), [k3] "m" (k3), [k4] "m" (k4)
+                : "cc", "memory",
+                  "%xmm0", "%xmm1", "%xmm2", "%xmm3",
+                  "%xmm4", "%xmm5", "%xmm6", "%xmm7"
+            );
+        }
+
+        void mix_add4(float *dst, const float *a, const float *b, const float *c, const float *d, float k1, float k2, float k3, float k4, size_t count)
+        {
+            IF_ARCH_X86(size_t off);
+            ARCH_X86_ASM
+            (
+                // 64x blocks
+                __ASM_EMIT("xor             %[off], %[off]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("vbroadcastss    %[k1], %%zmm6")
+                __ASM_EMIT("vbroadcastss    %[k2], %%zmm7")
+                __ASM_EMIT("vbroadcastss    %[k3], %%zmm5")
+                __ASM_EMIT("vbroadcastss    %[k4], %%zmm4")
+                __ASM_EMIT("jb              2f")
+                __ASM_EMIT("1:")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a],%[off]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vmulps          0x080(%[a],%[off]), %%zmm6, %%zmm2")
+                __ASM_EMIT("vmulps          0x0c0(%[a],%[off]), %%zmm6, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b],%[off]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[b],%[off]), %%zmm7, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[b],%[off]), %%zmm7, %%zmm3")
+                __ASM_EMIT("vaddps          0x000(%[dst],%[off]), %%zmm0, %%zmm0")
+                __ASM_EMIT("vaddps          0x040(%[dst],%[off]), %%zmm1, %%zmm1")
+                __ASM_EMIT("vaddps          0x080(%[dst],%[off]), %%zmm2, %%zmm2")
+                __ASM_EMIT("vaddps          0x0c0(%[dst],%[off]), %%zmm3, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c],%[off]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[c],%[off]), %%zmm5, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[c],%[off]), %%zmm5, %%zmm3")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d],%[off]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x080(%[d],%[off]), %%zmm4, %%zmm2")
+                __ASM_EMIT("vfmadd231ps     0x0c0(%[d],%[off]), %%zmm4, %%zmm3")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm2, 0x080(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm3, 0x0c0(%[dst],%[off])")
+                __ASM_EMIT("add             $0x100, %[off]")
+                __ASM_EMIT("sub             $64, %[count]")
+                __ASM_EMIT("jae             1b")
+                // 32x block
+                __ASM_EMIT("2:")
+                __ASM_EMIT("add             $32, %[count]")
+                __ASM_EMIT("jl              4f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%zmm6, %%zmm0")
+                __ASM_EMIT("vmulps          0x040(%[a],%[off]), %%zmm6, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%zmm7, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[b],%[off]), %%zmm7, %%zmm1")
+                __ASM_EMIT("vaddps          0x000(%[dst],%[off]), %%zmm0, %%zmm0")
+                __ASM_EMIT("vaddps          0x040(%[dst],%[off]), %%zmm1, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%zmm5, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[c],%[off]), %%zmm5, %%zmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%zmm4, %%zmm0")
+                __ASM_EMIT("vfmadd231ps     0x040(%[d],%[off]), %%zmm4, %%zmm1")
+                __ASM_EMIT("vmovups         %%zmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%zmm1, 0x040(%[dst],%[off])")
+                __ASM_EMIT("sub             $32, %[count]")
+                __ASM_EMIT("add             $0x80, %[off]")
+                // 16x block
+                __ASM_EMIT("4:")
+                __ASM_EMIT("add             $16, %[count]")
+                __ASM_EMIT("jl              6f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%ymm6, %%ymm0")
+                __ASM_EMIT("vmulps          0x020(%[a],%[off]), %%ymm6, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%ymm7, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[b],%[off]), %%ymm7, %%ymm1")
+                __ASM_EMIT("vaddps          0x000(%[dst],%[off]), %%ymm0, %%ymm0")
+                __ASM_EMIT("vaddps          0x020(%[dst],%[off]), %%ymm1, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%ymm5, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[c],%[off]), %%ymm5, %%ymm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%ymm4, %%ymm0")
+                __ASM_EMIT("vfmadd231ps     0x020(%[d],%[off]), %%ymm4, %%ymm1")
+                __ASM_EMIT("vmovups         %%ymm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%ymm1, 0x020(%[dst],%[off])")
+                __ASM_EMIT("sub             $16, %[count]")
+                __ASM_EMIT("add             $0x040, %[off]")
+                // 8x block
+                __ASM_EMIT("6:")
+                __ASM_EMIT("add             $8, %[count]")
+                __ASM_EMIT("jl              8f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vmulps          0x010(%[a],%[off]), %%xmm6, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[b],%[off]), %%xmm7, %%xmm1")
+                __ASM_EMIT("vaddps          0x000(%[dst],%[off]), %%xmm0, %%xmm0")
+                __ASM_EMIT("vaddps          0x010(%[dst],%[off]), %%xmm1, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[c],%[off]), %%xmm5, %%xmm1")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x010(%[d],%[off]), %%xmm4, %%xmm1")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("vmovups         %%xmm1, 0x010(%[dst],%[off])")
+                __ASM_EMIT("sub             $8, %[count]")
+                __ASM_EMIT("add             $0x020, %[off]")
+                // 4x block
+                __ASM_EMIT("8:")
+                __ASM_EMIT("add             $4, %[count]")
+                __ASM_EMIT("jl              10f")
+                __ASM_EMIT("vmulps          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vaddps          0x000(%[dst],%[off]), %%xmm0, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ps     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovups         %%xmm0, 0x000(%[dst],%[off])")
+                __ASM_EMIT("sub             $4, %[count]")
+                __ASM_EMIT("add             $0x010, %[off]")
+                // 1x blocks
+                __ASM_EMIT("10:")
+                __ASM_EMIT("add             $3, %[count]")
+                __ASM_EMIT("jl              12f")
+                __ASM_EMIT("11:")
+                __ASM_EMIT("vmulss          0x000(%[a],%[off]), %%xmm6, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[b],%[off]), %%xmm7, %%xmm0")
+                __ASM_EMIT("vaddss          0x000(%[dst],%[off]), %%xmm0, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[c],%[off]), %%xmm5, %%xmm0")
+                __ASM_EMIT("vfmadd231ss     0x000(%[d],%[off]), %%xmm4, %%xmm0")
+                __ASM_EMIT("vmovss          %%xmm0, 0x00(%[dst],%[off])")
+                __ASM_EMIT("add             $0x04, %[off]")
+                __ASM_EMIT("dec             %[count]")
+                __ASM_EMIT("jge             11b")
+                // End
+                __ASM_EMIT("12:")
+                : [count] "+r" (count),
+                  [off] "=&r" (off)
+                : [dst] "r" (dst), [a] "r" (a), [b] "r" (b), [c] "r" (c), [d] "r" (d),
+                  [k1] "m" (k1), [k2] "m" (k2), [k3] "m" (k3), [k4] "m" (k4)
+                : "cc", "memory",
+                  "%xmm0", "%xmm1", "%xmm2", "%xmm3",
+                  "%xmm4", "%xmm5", "%xmm6", "%xmm7"
+            );
+        }
+    #endif
+
     } /* namespace avx512 */
 } /* namespace lsp */
 
diff --git a/src/main/x86/avx512.cpp b/src/main/x86/avx512.cpp
index 0e7a63d..830b772 100644
--- a/src/main/x86/avx512.cpp
+++ b/src/main/x86/avx512.cpp
@@ -319,6 +319,10 @@
                 CEXPORT1(vl, mix3);
                 CEXPORT1(vl, mix_copy3);
                 CEXPORT1(vl, mix_add3);
+
+                CEXPORT1(vl, mix4);
+                CEXPORT1(vl, mix_copy4);
+                CEXPORT1(vl, mix_add4);
             }
         } /* namespace avx2 */
     } /* namespace lsp */
diff --git a/src/test/ptest/mix4.cpp b/src/test/ptest/mix4.cpp
index f6bd859..1d0652d 100644
--- a/src/test/ptest/mix4.cpp
+++ b/src/test/ptest/mix4.cpp
@@ -54,9 +54,9 @@ namespace lsp
 
         namespace avx512
         {
-//            void mix4(float *dst, const float *src1, const float *src2, const float *src3, float k1, float k2, float k3, float k4, size_t count);
-//            void mix_copy4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
-//            void mix_add4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
+            void mix4(float *dst, const float *src1, const float *src2, const float *src3, float k1, float k2, float k3, float k4, size_t count);
+            void mix_copy4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
+            void mix_add4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
         }
     )
 
@@ -141,7 +141,7 @@ PTEST_MAIN
         CALL(generic::mix4);
         IF_ARCH_X86(CALL(sse::mix4));
         IF_ARCH_X86(CALL(avx::mix4));
-//        IF_ARCH_X86(CALL(avx512::mix4));
+        IF_ARCH_X86(CALL(avx512::mix4));
         IF_ARCH_ARM(CALL(neon_d32::mix4));
         IF_ARCH_AARCH64(CALL(asimd::mix4));
         PTEST_SEPARATOR;
@@ -149,7 +149,7 @@ PTEST_MAIN
         CALL(generic::mix_copy4);
         IF_ARCH_X86(CALL(sse::mix_copy4));
         IF_ARCH_X86(CALL(avx::mix_copy4));
-//        IF_ARCH_X86(CALL(avx512::mix_copy4));
+        IF_ARCH_X86(CALL(avx512::mix_copy4));
         IF_ARCH_ARM(CALL(neon_d32::mix_copy4));
         IF_ARCH_AARCH64(CALL(asimd::mix_copy4));
         PTEST_SEPARATOR;
@@ -157,7 +157,7 @@ PTEST_MAIN
         CALL(generic::mix_add4);
         IF_ARCH_X86(CALL(sse::mix_add4));
         IF_ARCH_X86(CALL(avx::mix_add4));
-//        IF_ARCH_X86(CALL(avx512::mix_add4));
+        IF_ARCH_X86(CALL(avx512::mix_add4));
         IF_ARCH_ARM(CALL(neon_d32::mix_add4));
         IF_ARCH_AARCH64(CALL(asimd::mix_add4));
         PTEST_SEPARATOR2;
diff --git a/src/test/utest/mix4.cpp b/src/test/utest/mix4.cpp
index aab144c..2a8ed88 100644
--- a/src/test/utest/mix4.cpp
+++ b/src/test/utest/mix4.cpp
@@ -49,9 +49,9 @@ namespace lsp
 
         namespace avx512
         {
-//            void mix4(float *dst, const float *src1, const float *src2, const float *src3, float k1, float k2, float k3, float k4, size_t count);
-//            void mix_copy4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
-//            void mix_add4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
+            void mix4(float *dst, const float *src1, const float *src2, const float *src3, float k1, float k2, float k3, float k4, size_t count);
+            void mix_copy4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
+            void mix_add4(float *dst, const float *src1, const float *src2, const float *src3, const float *src4, float k1, float k2, float k3, float k4, size_t count);
         }
     )
 
@@ -181,9 +181,9 @@ UTEST_BEGIN("dsp", mix4)
         IF_ARCH_X86(CALL(generic::mix_copy4, avx::mix_copy4, 32));
         IF_ARCH_X86(CALL(generic::mix_add4, avx::mix_add4, 32));
 
-//        IF_ARCH_X86(CALL(generic::mix4, avx512::mix4, 64));
-//        IF_ARCH_X86(CALL(generic::mix_copy4, avx512::mix_copy4, 64));
-//        IF_ARCH_X86(CALL(generic::mix_add4, avx512::mix_add4, 64));
+        IF_ARCH_X86(CALL(generic::mix4, avx512::mix4, 64));
+        IF_ARCH_X86(CALL(generic::mix_copy4, avx512::mix_copy4, 64));
+        IF_ARCH_X86(CALL(generic::mix_add4, avx512::mix_add4, 64));
 
         IF_ARCH_ARM(CALL(generic::mix4, neon_d32::mix4, 16));
         IF_ARCH_ARM(CALL(generic::mix_copy4, neon_d32::mix_copy4, 16));