Release 1.0.19

* AVX2 optimization of search functions for maximum and minimum.
* Implemented SIMD-optimized gate functions.
* Fixed bun in AVX-512 implementation of mid/side conversion functions.
* AVX512 optimization of packed complex functions.
* Additional optimizations of parallel mathematics functions with AVX-512 instruction set.

Author: sadko4u

CHANGELOG

@@ -2,6 +2,13 @@
 * RECENT CHANGES
 *******************************************************************************
 
+=== 1.0.19 ===
+* AVX2 optimization of search functions for maximum and minimum.
+* Implemented SIMD-optimized gate functions.
+* Fixed bun in AVX-512 implementation of mid/side conversion functions.
+* AVX512 optimization of packed complex functions.
+* Additional optimizations of parallel mathematics functions with AVX-512 instruction set.
+
 === 1.0.18 ===
 * Fixed compilation regression for 32-bit Clang compiler.
 

include/lsp-plug.in/dsp/common/dynamics.h

@@ -26,6 +26,7 @@
 
 #include <lsp-plug.in/dsp/common/dynamics/types.h>
 #include <lsp-plug.in/dsp/common/dynamics/compressor.h>
+#include <lsp-plug.in/dsp/common/dynamics/gate.h>
 
 
 #endif /* LSP_PLUG_IN_DSP_COMMON_DYNAMICS_H_ */

include/lsp-plug.in/dsp/common/dynamics/gate.h

@@ -0,0 +1,33 @@
+/*
+ * Copyright (C) 2023 Linux Studio Plugins Project <https://lsp-plug.in/>
+ *           (C) 2023 Vladimir Sadovnikov <sadko4u@gmail.com>
+ *
+ * This file is part of lsp-dsp-lib
+ * Created on: 19 окт. 2023 г.
+ *
+ * lsp-dsp-lib 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 3 of the License, or
+ * any later version.
+ *
+ * lsp-dsp-lib 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 lsp-dsp-lib. If not, see <https://www.gnu.org/licenses/>.
+ */
+
+#ifndef LSP_PLUG_IN_DSP_COMMON_DYNAMICS_GATE_H_
+#define LSP_PLUG_IN_DSP_COMMON_DYNAMICS_GATE_H_
+
+#include <lsp-plug.in/dsp/common/types.h>
+#include <lsp-plug.in/dsp/common/dynamics/types.h>
+
+LSP_DSP_LIB_SYMBOL(void, gate_x1_gain, float *dst, const float *src, const LSP_DSP_LIB_TYPE(gate_knee_t) *c, size_t count);
+
+LSP_DSP_LIB_SYMBOL(void, gate_x1_curve, float *dst, const float *src, const LSP_DSP_LIB_TYPE(gate_knee_t) *c, size_t count);
+
+
+#endif /* LSP_PLUG_IN_DSP_COMMON_DYNAMICS_GATE_H_ */

include/lsp-plug.in/dsp/common/dynamics/types.h

@@ -29,7 +29,7 @@ LSP_DSP_LIB_BEGIN_NAMESPACE
 #pragma pack(push, 1)
 
 /**
- * Compressor knee is a curve that constists of three parts:
+ * Compressor knee is a curve that consists of three parts:
  *   1. Part with constant gain amplification in the range [-inf .. start] dB
  *   2. Soft compression knee in the range (start .. end) dB present by the quadratic function (2nd-order polynom)
  *   3. Gain reduction part in the range [end .. +inf] dB present by the linear function (1st-order polynom)
@@ -40,7 +40,7 @@ LSP_DSP_LIB_BEGIN_NAMESPACE
  *   3. Compute the natural logarithm of the x: lx = logf(x).
  *   4. If x < end then compute the gain using the 2nd-order polynom: gain = (herm[0]*lx + herm[1])*lx + herm[2]
  *   5. Otherwise compute the gain using the 1st-order polynom: gain = tilt[0]*lx + tilt[1]
- *   6. return expf(gain)
+ *   6. return expf(gain) * x
  */
 typedef struct LSP_DSP_LIB_TYPE(compressor_knee_t)
 {
@@ -61,6 +61,30 @@ typedef struct LSP_DSP_LIB_TYPE(compressor_x2_t)
     LSP_DSP_LIB_TYPE(compressor_knee_t)   k[2];
 } LSP_DSP_LIB_TYPE(compressor_x2_t);
 
+
+/**
+ * Gate knee is a curve that consists of three parts:
+ *   1. Part with constant gain amplification in the range [-inf .. start] dB
+ *   2. Transition zone in the range (start .. end) dB present by the quadratic function (2nd-order polynom)
+ *   3. Part with constant gain amplification in the range [end .. +inf] dB
+ *
+ * The typical algorithm of computing the gate's curve:
+ *   1. Take absolute value of the sample: x = fabfs(in)
+ *   2. If x <= start then return gain_start * x
+ *   3. If x <= end then return gain_end * x
+ *   4. Compute the natural logarithm of the x: lx = logf(x).
+ *   5. Compute the gain using the 3rd-order polynom: gain = ((herm[0]*lx + herm[1])*lx + herm[2]*lx) + herm[3]
+ *   6. return expf(gain) * x
+ */
+typedef struct LSP_DSP_LIB_TYPE(gate_knee_t)
+{
+    float       start;          // The start of the knee, in gain units
+    float       end;            // The end of the knee, in gain units
+    float       gain_start;     // Gain below the start threshold
+    float       gain_end;       // Gain above the end threshold
+    float       herm[4];        // Hermite interpolation of the knee with the 3rd-order polynom
+} LSP_DSP_LIB_TYPE(gate_knee_t);
+
 #pragma pack(pop)
 
 LSP_DSP_LIB_END_NAMESPACE

include/lsp-plug.in/dsp/version.h

@@ -25,7 +25,7 @@
 // Define version of headers
 #define LSP_DSP_LIB_MAJOR           1
 #define LSP_DSP_LIB_MINOR           0
-#define LSP_DSP_LIB_MICRO           18
+#define LSP_DSP_LIB_MICRO           19
 
 #if defined(__WINDOWS__) || defined(__WIN32__) || defined(__WIN64__) || defined(_WIN64) || defined(_WIN32) || defined(__WINNT) || defined(__WINNT__)
     #define LSP_DSP_LIB_EXPORT_MODIFIER     __declspec(dllexport)

include/private/dsp/arch/aarch64/asimd/dynamics.h

@@ -27,6 +27,6 @@
 #endif /* PRIVATE_DSP_ARCH_AARCH64_ASIMD_IMPL */
 
 #include <private/dsp/arch/aarch64/asimd/dynamics/compressor.h>
-
+#include <private/dsp/arch/aarch64/asimd/dynamics/gate.h>
 
 #endif /* PRIVATE_DSP_ARCH_AARCH64_ASIMD_DYNAMICS_H_ */

include/private/dsp/arch/aarch64/asimd/dynamics/compressor.h

@@ -294,6 +294,11 @@ namespace lsp
             );
         }
 
+    #undef KNEE_LOAD
+    #undef PROCESS_KNEE_SINGLE_X4
+    #undef PROCESS_KNEE_SINGLE_X8
+    #undef PROCESS_COMP_FULL_X8
+    #undef PROCESS_COMP_FULL_X4
 
     } /* namespace asimd */
 } /* namespace lsp */