summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/matrixOperations/jmat/sjmata.c
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/c/matrixOperations/jmat/sjmata.c')
-rw-r--r--2.3-1/src/c/matrixOperations/jmat/sjmata.c22
1 files changed, 22 insertions, 0 deletions
diff --git a/2.3-1/src/c/matrixOperations/jmat/sjmata.c b/2.3-1/src/c/matrixOperations/jmat/sjmata.c
new file mode 100644
index 00000000..b969fd7a
--- /dev/null
+++ b/2.3-1/src/c/matrixOperations/jmat/sjmata.c
@@ -0,0 +1,22 @@
+#include <jmat.h>
+
+void sjmata(float n, float m, float *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 to be 1
+ * */
+ out[i*mm*mm*nn+(mm*(nn-(i+1))+j*(mm*nn+1))]=1;
+ }
+ }
+
+}
+
+