Original Version

1 files changed, 22 insertions, 0 deletions

diff --git a/src/c/matrixOperations/jmat/djmata.c b/src/c/matrixOperations/jmat/djmata.c
new file mode 100644
index 00000000..6ccff61d
--- /dev/null
+++ b/src/c/matrixOperations/jmat/djmata.c
@@ -0,0 +1,22 @@
+#include <jmat.h>
+
+void djmata(double n, double m, double *out){
+	int i,nn,mm,j;
+	mm=(int) m;
+	nn=(int) n;
+
+	
+	for (i=0;i<n;i++){/*We work by blocs*/
+		for(j=0;j<m;j++){/*then by column*/
+			/*
+			 * i*mm*nn : advance to the good column
+			 * mm*(nn-(i+1)) : take the first element of the column to be 1, when j=0
+			 * j*(mm*nn+1) : advance to the next which must be 1
+			 * */
+			out[i*mm*mm*nn+(mm*(nn-(i+1))+j*(mm*nn+1))]=1;
+		}
+	}
+	
+}
+
+