From f481f7f62a5ea6ae9d197d79c9043a6bcbb2673c Mon Sep 17 00:00:00 2001
From: kevin <kevin@administrators-MacBook.local>
Date: Fri, 3 Jul 2015 18:22:08 -0400
Subject: [PATCH] cmComplexTypes.h : Added complex arithmetic vector
 operations.

---
 cmComplexTypes.c | 42 ++++++++++++++++++++++++++++++++++++++++++
 cmComplexTypes.h |  6 ++++++
 2 files changed, 48 insertions(+)
 create mode 100644 cmComplexTypes.c

diff --git a/cmComplexTypes.c b/cmComplexTypes.c
new file mode 100644
index 0000000..2929400
--- /dev/null
+++ b/cmComplexTypes.c
@@ -0,0 +1,42 @@
+
+
+void cmVOCR_MultVVV( cmComplexR_t* y, const cmComplex_t* x0, const cmComplexR_t* x1, unsigned n )
+{
+  unsigned i;
+  for(i=0; i<n; ++i)
+  {
+    cmReal_t ab = x0[i].r * x1[i].r;
+    cmReal_t bd = x0[i].i * x1[i].i;
+    cmReal_t bc = x0[i].i * x1[i].r;
+    cmReal_t ad = x0[i].r * x1[i].i;
+    y[i].r = ab - bd;
+    y[i].i = bc + ad;    
+  }
+
+}
+
+void cmVOCR_Abs(     cmReal_t* y, const cmComplexR_t* x, unsigned n )
+{
+  for(i=0; i<n; ++i)
+    y[i] = cmAbsR(x[i]);
+}
+
+void cmVOCR_MultR_VV(  cmComplexR_t* y, const cmReal_t* x, unsigned n )
+{
+  unsigned i;
+  for(i=0; i<n; ++i)
+  {
+    y[i].r *= x[i];
+    y[i].i *= x[i];
+  }
+}
+
+void cmVOCR_DivR_VV(   cmComplexR_t* y, const cmReal_t* x, unsigned n )
+{
+  unsigned i;
+  for(i=0; i<n; ++i)
+  {
+    y[i].r /= x[i];
+    y[i].i /= x[i];
+  }  
+}
diff --git a/cmComplexTypes.h b/cmComplexTypes.h
index e35b0db..c6711df 100644
--- a/cmComplexTypes.h
+++ b/cmComplexTypes.h
@@ -84,4 +84,10 @@
 
 #endif
 
+void cmVOCR_MultVVV( cmComplexR_t* y, const cmComplexR_t* x0, const cmComplexR_t* x1, unsigned n );
+void cmVOCR_Abs(     cmReal_t*     y, const cmComplexR_t* x, unsigned n );
+void cmVOCR_MultVS(  cmComplexR_t* y, cmReal_t v, unsigned n );
+void cmVOCR_DivVS(   cmComplexR_t* y, cmReal_t v, unsigned n );
+
+
 #endif