Some memory optimizations

Author: sadko4u

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

@@ -119,12 +119,32 @@ LSP_DSP_LIB_SYMBOL(void, pcomplex_c2r, float *dst, const float *src, size_t coun
 
 /** Get module for complex numbers: mod = sqrt(re*re + im*im)
  *
- * @param dst_mod array to sore module
+ * @param dst_mod array to store module
  * @param src packed complex number data
  * @param count count number of elements to process
  */
 LSP_DSP_LIB_SYMBOL(void, pcomplex_mod, float *dst_mod, const float *src, size_t count);
 
+/**
+ * Get module for complex numbers and add to destination:
+ *   dst[i] = dst[i] + sqrt(re[i]*re[i] + im[i]*im[i])
+ *
+ * @param dst array to add module
+ * @param src packed complex number data
+ * @param count count number of elements to process
+ */
+LSP_DSP_LIB_SYMBOL(void, pcomplex_mod_add2, float *dst, const float *src, size_t count);
+
+/**
+ * Get module for complex numbers and add to destination:
+ *   dst[i] = src1[i] + sqrt(re[i]*re[i] + im[i]*im[i])
+ *
+ * @param dst array to add module
+ * @param src packed complex number data
+ * @param count count number of elements to process
+ */
+LSP_DSP_LIB_SYMBOL(void, pcomplex_mod_add3, float *dst, const float *src, const float *reim, size_t count);
+
 /** Convert packed complex number to polar form
  *
  * @param mod module of the complex number

include/private/dsp/arch/generic/pcomplex.h

@@ -205,6 +205,33 @@ namespace lsp
             }
         }
 
+        void pcomplex_mod_add2(float *dst, const float *src, size_t count)
+        {
+            while (count--)
+            {
+                float re        = src[0];
+                float im        = src[1];
+                *dst           += sqrtf(re*re + im*im);
+
+                src            += 2;
+                ++dst;
+            }
+        }
+
+        void pcomplex_mod_add3(float *dst, const float * src, const float *reim, size_t count)
+        {
+            while (count--)
+            {
+                float re        = reim[0];
+                float im        = reim[1];
+                *dst            = *src + sqrtf(re*re + im*im);
+
+                reim           += 2;
+                ++dst;
+                ++src;
+            }
+        }
+
         void pcomplex_arg(float *dst, const float *src, size_t count)
         {
             for (; count > 0; --count, src += 2)

src/main/generic/generic.cpp

@@ -354,6 +354,8 @@ namespace lsp
             EXPORT1(pcomplex_r2c);
             EXPORT1(pcomplex_c2r);
             EXPORT1(pcomplex_mod);
+            EXPORT1(pcomplex_mod_add2);
+            EXPORT1(pcomplex_mod_add3);
             EXPORT1(pcomplex_arg);
             EXPORT1(pcomplex_modarg);
             EXPORT1(pcomplex_corr);