ml_util.cpp

00001 #include "ml_util.hpp"
00002 
00003 void multiply_mat_vec(float*  mat, // m * n 行列
00004                       float*  vec, // n * 1 ベクトル
00005                       float*  result, // m * 1 計算結果ベクトル
00006                       int     m,   // m
00007                       int     n)   // n
00008 {
00009     register float mat_sum;
00010     for (int mat_mul_i = 0; mat_mul_i < m; mat_mul_i++) {
00011         mat_sum = 0;
00012         /*
00013         printf("result[%d] : %p, MATRIX_AT(mat,n,mat_mul_i,mat_mul_j) : %p, vec[0] : %p", mat_mul_i, 
00014                 &(result[mat_mul_i]), &(MATRIX_AT(mat,n,mat_mul_i,0)), vec);
00015                 */
00016         for (int mat_mul_j = 0; mat_mul_j < n; mat_mul_j++) {
00017             mat_sum += MATRIX_AT(mat,n,mat_mul_i,mat_mul_j) * vec[mat_mul_j];
00018         }
00019         result[mat_mul_i] = mat_sum;
00020     }
00021 }