multithreading tests from 152 lab 5

[riscv-tests.git] / mt / am_matmul / matmul2.c

diff --git a/mt/am_matmul/matmul2.c b/mt/am_matmul/matmul2.c
new file mode 100644 (file)
index 0000000..30c705d
--- /dev/null
+++ b/mt/am_matmul/matmul2.c
@@ -0,0 +1,73 @@
+/*size_t i;
+  size_t j;
+  size_t max_dim = lda*lda;
+  if (coreid==0){
+    for (i=0; i<max_dim/(ncores*2); i+=8){
+      data_t elementA1 = A[i];
+      data_t elementA12 = A[i+1];
+      data_t elementA13 = A[i+2];
+      data_t elementA14 = A[i+3];
+      data_t elementA15 = A[i+4];
+      data_t elementA16 = A[i+5];
+      data_t elementA17 = A[i+6];
+      data_t elementA18 = A[i+7];
+      data_t elementA2 = A[i+32*8];
+      data_t elementA21 = A[i+32*8+1];
+      data_t elementA22 = A[i+32*8+2];
+      data_t elementA23 = A[i+32*8+3];
+      data_t elementA24 = A[i+32*8+4];
+      data_t elementA25 = A[i+32*8+5];
+      data_t elementA26 = A[i+32*8+6];
+      data_t elementA27 = A[i+32*8+7];
+      int row= (int)(i/32)*32;
+      int row2 = row+8*32;
+      int column1 = i%32*32;
+      int column12 = (i+1)%32*32;
+      int column13 = (i+2)%32*32;
+      int column14 = (i+3)%32*32;
+      int column15 = (i+4)%32*32;
+      int column16 = (i+5)%32*32;
+      int column17 = (i+6)%32*32;
+      int column18 = (i+7)%32*32;
+    
+      for (j=0; j<lda; j++){
+       C[row+j]+=elementA1*B[column1+j]+elementA12*B[column12+j]+elementA13*B[column13+j]+elementA14*B[column14+j]+elementA15*B[column15+j]+elementA16*B[column16+j]+elementA17*B[column17+j]+elementA18*B[column18+j]
+ 
+       C[row2+j]+=elementA2*B[column1+j]+elementA21*B[column12+j]+elementA22*B[column13+j]+elementA23*B[column14+j]+elementA24*B[column15+j]+elementA25*B[column16+j]+elementA26*B[column17+j]+elementA27*B[column18+j];
+      }
+  }}else{
+    for (i=max_dim/2; i<(max_dim/(ncores*2)+max_dim/2); i+=8){
+      data_t elementA1 = A[i];
+      data_t elementA12 = A[i+1];
+      data_t elementA13 = A[i+2];
+      data_t elementA14 = A[i+3];
+      data_t elementA15 = A[i+4];
+      data_t elementA16 = A[i+5];
+      data_t elementA17 = A[i+6];
+      data_t elementA18 = A[i+7];
+      data_t elementA2 = A[i+32*8];
+      data_t elementA21 = A[i+32*8+1];
+      data_t elementA22 = A[i+32*8+2];
+      data_t elementA23 = A[i+32*8+3];
+      data_t elementA24 = A[i+32*8+4];
+      data_t elementA25 = A[i+32*8+5];
+      data_t elementA26 = A[i+32*8+6];
+      data_t elementA27 = A[i+32*8+7];
+      int row= (int)(i/32)*32;
+      int row2 = row+8*32;
+      int column1 = i%32*32;
+      int column12 = (i+1)%32*32;
+      int column13 = (i+2)%32*32;
+      int column14 = (i+3)%32*32;
+      int column15 = (i+4)%32*32;
+      int column16 = (i+5)%32*32;
+      int column17 = (i+6)%32*32;
+      int column18 = (i+7)%32*32;
+    
+      for (j=0; j<lda; j++){
+       C[row+j]+=elementA1*B[column1+j]+elementA12*B[column12+j]+elementA13*B[column13+j]+elementA14*B[column14+j]+elementA15*B[column15+j]+elementA16*B[column16+j]+elementA17*B[column17+j]+elementA18*B[column18+j];
+       C[row2+j]+=elementA2*B[column1+j]+elementA21*B[column12+j]+elementA22*B[column13+j]+elementA23*B[column14+j]+elementA24*B[column15+j]+elementA25*B[column16+j]+elementA26*B[column17+j]+elementA27*B[column18+j];
+
+      }
+    }
+    }*/