From 124ef41ebab756797803b30c4c07ce37213a6477 Mon Sep 17 00:00:00 2001
From: Brijeshcr
Date: Fri, 7 Jul 2017 16:43:27 +0530
Subject: Toeplitz Added

---
 src/c/matrixOperations/toeplitz/ctoeplitza.c | 36 ++++++++++++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100644 src/c/matrixOperations/toeplitz/ctoeplitza.c

(limited to 'src/c/matrixOperations/toeplitz/ctoeplitza.c')

diff --git a/src/c/matrixOperations/toeplitz/ctoeplitza.c b/src/c/matrixOperations/toeplitz/ctoeplitza.c
new file mode 100644
index 0000000..048a367
--- /dev/null
+++ b/src/c/matrixOperations/toeplitz/ctoeplitza.c
@@ -0,0 +1,36 @@
+#include <stdio.h>
+#include "toeplitz.h"
+#include "floatComplex.h"
+#include "stdlib.h"
+#include "string.h"
+#include "cat.h"
+
+/*Function to build a Toeplitz Matrix for inputs of SingleComplex datatype*/
+
+void ctoeplitza(floatComplex* inp1,int size1,floatComplex* inp2,int size2,floatComplex* oup)
+{
+    if (inp1[0]!=inp2[0])
+    {
+        printf("Error!The first elements of the Vectors are not equal."); // First element of both input vectors must be equal for Toeplitz.
+        return;
+    }
+    int i, j;
+
+    for(i=0;i<size1*size2;i++) oup[i] = 0;       // Initializing the output matrix with zeros.
+
+    for (i = 0; i<size1; i++)
+    {
+        for (j = 0; j<size2; j++)
+        {
+            oup[j*size1] = inp2[j];              // Elements of the second input vector are copied to the first row of the Toeplitx Matrix.
+        }
+        oup[i] = inp1[i];                        // Elements of the first input vector are copied to the first column of the Toeplitx Matrix.
+    }
+    for (i = size2+1; i<size1*size2; i++)        // Loop to build the rest of the Toeplitz matrix.
+    {
+        if (oup[i] == 0)
+        oup[i] = oup[i-size2-1];
+    }		
+}
+
+
-- 
cgit 


From 4ebc92d414c6e87cea84770f4137434b922a76ca Mon Sep 17 00:00:00 2001
From: Brijeshcr
Date: Fri, 7 Jul 2017 22:29:47 +0530
Subject: Changes to Toeplitz

---
 src/c/matrixOperations/toeplitz/ctoeplitza.c | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

(limited to 'src/c/matrixOperations/toeplitz/ctoeplitza.c')

diff --git a/src/c/matrixOperations/toeplitz/ctoeplitza.c b/src/c/matrixOperations/toeplitz/ctoeplitza.c
index 048a367..56c3ef8 100644
--- a/src/c/matrixOperations/toeplitz/ctoeplitza.c
+++ b/src/c/matrixOperations/toeplitz/ctoeplitza.c
@@ -9,14 +9,14 @@
 
 void ctoeplitza(floatComplex* inp1,int size1,floatComplex* inp2,int size2,floatComplex* oup)
 {
-    if (inp1[0]!=inp2[0])
+    if ((creals(inp1[0])!=creals(inp2[0]))&&(cimags(inp1[0])!=cimags(inp2[0])))
     {
         printf("Error!The first elements of the Vectors are not equal."); // First element of both input vectors must be equal for Toeplitz.
         return;
     }
     int i, j;
 
-    for(i=0;i<size1*size2;i++) oup[i] = 0;       // Initializing the output matrix with zeros.
+    for(i=0;i<size1*size2;i++) oup[i] = FloatComplex(0,0);       // Initializing the output matrix with zeros.
 
     for (i = 0; i<size1; i++)
     {
@@ -28,7 +28,7 @@ void ctoeplitza(floatComplex* inp1,int size1,floatComplex* inp2,int size2,floatC
     }
     for (i = size2+1; i<size1*size2; i++)        // Loop to build the rest of the Toeplitz matrix.
     {
-        if (oup[i] == 0)
+        if (creals(oup[i]) == 0 && cimags(oup[i]) == 0)
         oup[i] = oup[i-size2-1];
     }		
 }
-- 
cgit