Changes between Version 6 and Version 7 of Ex02課題2015

Timestamp:

Apr 20, 2015 5:38:57 PM (11 years ago)

Author:

nakasato

Comment:

--

Ex02課題2015

@@ -85 +85 @@
 
 = 課題5：行列積 =
+以下のC言語で書かれた行列乗算のプログラムを参考にして、行列積C = A x Bを計算するプログラムを作ってください。
 
+== 行列A ==
+{{{
+0  1  0  0
+2  0  0  0 
+0  0  0  3
+0  0  4  0
+}}}
+
+== 行列B ==
+{{{
+1  2  3  4
+5  6  7  8 
+9 10 11 12
+13 14 15 16
+}}}
+
+== C言語バージョン ==
+{{{
+
+#include <stdio.h>
+main()
+{
+  static int mat1[4][4] = {
+    { 1, 0, 0, 0 },
+    { 0, 1, 0, 0 },
+    { 0, 0, 1, 0 },
+    { 0, 0, 0, 1 },
+  };
+  static int mat2[4][4] = {
+    {  1,  2,  3,  4 },
+    {  5,  6,  7,  8 },
+    {  9, 10, 11, 12 },
+    { 13, 14, 15, 16 },
+  };
+  static int result[4][4];
+
+  int  i, j, k;
+  int  s;
+  int  cnt, mask, m1, m2;
+
+  /* 行列の乗算 */
+  for( i = 0; i < 4; i++ ) {
+    for( j = 0; j < 4; j++ ) {
+      s = 0;
+      for( k = 0; k < 4; k++ ) {
+        mask = 1;
+        m1 = mat1[i][k];  /* 被乗数 */
+        m2 = mat2[k][j];  /* 乗数 */
+        /* 本来は下のループで32回繰り返すべきだが、
+           乗数が32未満なので、計算時間短縮のため5回だけとした */
+        for( cnt = 0; cnt < 5; cnt++ ) {
+          if( m2 & mask ) {
+            s += m1;    
+          }
+          mask <<= 1;
+          m1 <<= 1;
+        }
+      }
+      result[i][j] = s;
+    }
+  }
+
+  /* 結果の表示 */
+  for( i = 0; i < 4; i++ ) {
+    for( j = 0; j < 4; j++ ) {
+      printf("%3d", result[i][j]); 
+    }
+    printf("\n");
+  }
+}
+}}}
+
+== テンプレート ==
 作成ファイル名は「ex02_p5.s」としてください。
+
+{{{
+        .data 
+A:      .word 0 
+        .word 1 
+        .word 0 
+        .word 0 
+    ... 配列Aの要素を書く .... 
+
+B:      .word 1 
+        .word 2 
+        .word 3 
+        .word 4 
+    ... 配列Bの要素を書く .... 
+
+C:      .space 64 ## 結果行列の保存用
+
+       ... 必要な変数の領域を追加する ....
+        .text 
+main:       ... プログラムを書く .... 
+
+
+exit:   j exit
+
+}}}