From db464f35f5a10b58d9ed1085e0b462689adee583 Mon Sep 17 00:00:00 2001
From: Siddhesh Wani
Date: Mon, 25 May 2015 14:46:31 +0530
Subject: Original Version

---
 src/c/signalProcessing/conv2d/dconv2da.c | 35 ++++++++++++++++++++++++++++++++
 1 file changed, 35 insertions(+)
 create mode 100644 src/c/signalProcessing/conv2d/dconv2da.c

(limited to 'src/c/signalProcessing/conv2d/dconv2da.c')

diff --git a/src/c/signalProcessing/conv2d/dconv2da.c b/src/c/signalProcessing/conv2d/dconv2da.c
new file mode 100644
index 00000000..6606e81d
--- /dev/null
+++ b/src/c/signalProcessing/conv2d/dconv2da.c
@@ -0,0 +1,35 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008 - INRIA - Arnaud TORSET
+ *
+ *  This file must be used under the terms of the CeCILL.
+ *  This source file is licensed as described in the file COPYING, which
+ *  you should have received as part of this distribution.  The terms
+ *  are also available at
+ *  http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+ *
+ */
+
+#include "conv2d.h"
+
+
+void dconv2da(double *in1, int lines1, int columns1, double *in2, int lines2, int columns2, double *out){
+	int i,j,k1,k2;
+	double accu;
+	
+	for (i=0;i<lines1+lines2-1;i++){
+		for (j=0;j<columns1+columns2-1;j++){
+			accu=0;
+			for (k1=0;k1<i+1;k1++){
+				if ((k1<lines1)&&((i-k1)<lines2)){
+					for(k2=0;k2<j+1;k2++){
+						if ((k2<columns1)&&((j-k2)<columns2))
+							accu += in1[k1+k2*lines1]*in2[i-k1+(j-k2)*lines2];
+					}
+				}
+			}
+			out[i+j*(lines1+lines2-1)]=accu;
+		}
+	}		
+}
+
-- 
cgit