Optimizations of clamp function for 32-bit ARM

Author: sadko4u

include/private/dsp/arch/arm/neon-d32/pmath.h

@@ -27,6 +27,7 @@
 #endif /* PRIVATE_DSP_ARCH_ARM_NEON_D32_IMPL */
 
 #include <private/dsp/arch/arm/neon-d32/pmath/abs_vv.h>
+#include <private/dsp/arch/arm/neon-d32/pmath/clamp.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/cos.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/exp.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/fmop_kx.h>
@@ -37,6 +38,7 @@
 #include <private/dsp/arch/arm/neon-d32/pmath/minmax.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/op_kx.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/op_vv.h>
+#include <private/dsp/arch/arm/neon-d32/pmath/pmix.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/pow.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/sin.h>
 #include <private/dsp/arch/arm/neon-d32/pmath/sqr.h>

include/private/dsp/arch/arm/neon-d32/pmath/clamp.h

@@ -0,0 +1,40 @@
+/*
+ * Copyright (C) 2025 Linux Studio Plugins Project <https://lsp-plug.in/>
+ *           (C) 2025 Vladimir Sadovnikov <sadko4u@gmail.com>
+ *
+ * This file is part of lsp-dsp-lib
+ * Created on: 27 нояб. 2025 г.
+ *
+ * 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 PRIVATE_DSP_ARCH_ARM_NEON_D32_PMATH_PMIX_H_
+#define PRIVATE_DSP_ARCH_ARM_NEON_D32_PMATH_PMIX_H_
+
+#ifndef PRIVATE_DSP_ARCH_ARM_NEON_D32_IMPL
+    #error "This header should not be included directly"
+#endif /* PRIVATE_DSP_ARCH_ARM_NEON_D32_IMPL */
+
+namespace lsp
+{
+    namespace neon_d32
+    {
+
+    } /* namespace neon_d32 */
+} /* namespace lsp */
+
+
+
+
+#endif /* PRIVATE_DSP_ARCH_ARM_NEON_D32_PMATH_PMIX_H_ */

include/private/dsp/arch/arm/neon-d32/pmath/pmix.h

@@ -436,6 +436,11 @@
                 EXPORT1(uexpander_x1_curve);
                 EXPORT1(dexpander_x1_gain);
                 EXPORT1(dexpander_x1_curve);
+
+                EXPORT1(clamp_vv1);
+                EXPORT1(clamp_vv2);
+                EXPORT1(clamp_kk1);
+                EXPORT1(clamp_kk2);
             }
         } /* namespace neon_d32 */
     } /* namespace lsp */

src/main/arm/neon-d32.cpp

@@ -109,7 +109,7 @@ PTEST_BEGIN("dsp.pmath", clamp_kk1, 5, 1000)
             IF_ARCH_X86(CALL(sse::clamp_kk1));
             IF_ARCH_X86(CALL(avx::clamp_kk1));
             IF_ARCH_X86(CALL(avx512::clamp_kk1));
-//            IF_ARCH_ARM(CALL(neon_d32::clamp_kk1));
+            IF_ARCH_ARM(CALL(neon_d32::clamp_kk1));
 //            IF_ARCH_AARCH64(CALL(asimd::clamp_kk1));
             PTEST_SEPARATOR;
         }

src/test/ptest/pmath/clamp_kk1.cpp

@@ -109,7 +109,7 @@ PTEST_BEGIN("dsp.pmath", clamp_kk2, 5, 1000)
             IF_ARCH_X86(CALL(sse::clamp_kk2));
             IF_ARCH_X86(CALL(avx::clamp_kk2));
             IF_ARCH_X86(CALL(avx512::clamp_kk2));
-//            IF_ARCH_ARM(CALL(neon_d32::clamp_kk2));
+            IF_ARCH_ARM(CALL(neon_d32::clamp_kk2));
 //            IF_ARCH_AARCH64(CALL(asimd::clamp_kk2));
             PTEST_SEPARATOR;
         }

src/test/ptest/pmath/clamp_kk2.cpp

@@ -113,7 +113,7 @@ PTEST_BEGIN("dsp.pmath", clamp_vv2, 5, 1000)
             IF_ARCH_X86(CALL(sse::clamp_vv2));
             IF_ARCH_X86(CALL(avx::clamp_vv2));
             IF_ARCH_X86(CALL(avx512::clamp_vv2));
-//            IF_ARCH_ARM(CALL(neon_d32::clamp_vv2));
+            IF_ARCH_ARM(CALL(neon_d32::clamp_vv2));
 //            IF_ARCH_AARCH64(CALL(asimd::clamp_vv2));
             PTEST_SEPARATOR;
         }

src/test/ptest/pmath/clamp_vv1.cpp

@@ -113,7 +113,7 @@ PTEST_BEGIN("dsp.pmath", clamp_vv1, 5, 1000)
             IF_ARCH_X86(CALL(sse::clamp_vv1));
             IF_ARCH_X86(CALL(avx::clamp_vv1));
             IF_ARCH_X86(CALL(avx512::clamp_vv1));
-//            IF_ARCH_ARM(CALL(neon_d32::clamp_vv1));
+            IF_ARCH_ARM(CALL(neon_d32::clamp_vv1));
 //            IF_ARCH_AARCH64(CALL(asimd::clamp_vv1));
             PTEST_SEPARATOR;
         }

src/test/ptest/pmath/clamp_vv2.cpp

@@ -118,7 +118,7 @@ UTEST_BEGIN("dsp.pmath", clamp_kk1)
         IF_ARCH_X86(CALL(generic::clamp_kk1, sse::clamp_kk1, 16));
         IF_ARCH_X86(CALL(generic::clamp_kk1, avx::clamp_kk1, 32));
         IF_ARCH_X86(CALL(generic::clamp_kk1, avx512::clamp_kk1, 32));
-//        IF_ARCH_ARM(CALL(generic::clamp_kk1, neon_d32::clamp_kk1, 16));
+        IF_ARCH_ARM(CALL(generic::clamp_kk1, neon_d32::clamp_kk1, 16));
 //        IF_ARCH_AARCH64(CALL(generic::clamp_kk1, asimd::clamp_kk1, 16));
     }
 UTEST_END

src/test/utest/pmath/clamp_kk1.cpp

@@ -118,7 +118,7 @@ UTEST_BEGIN("dsp.pmath", clamp_kk2)
         IF_ARCH_X86(CALL(generic::clamp_kk2, sse::clamp_kk2, 16));
         IF_ARCH_X86(CALL(generic::clamp_kk2, avx::clamp_kk2, 32));
         IF_ARCH_X86(CALL(generic::clamp_kk2, avx512::clamp_kk2, 32));
-//        IF_ARCH_ARM(CALL(generic::clamp_kk2, neon_d32::clamp_kk2, 16));
+        IF_ARCH_ARM(CALL(generic::clamp_kk2, neon_d32::clamp_kk2, 16));
 //        IF_ARCH_AARCH64(CALL(generic::clamp_kk2, asimd::clamp_kk2, 16));
     }
 UTEST_END