2017-03-31 225 views
-1

我的任務有問題,我有三個矩陣。我應該將第一個(a)乘以第二個(b),並將結果返回第三個(c)。還有值得說的是我只能通過指針移動,我不能在循環中使用索引。在C++中將矩陣乘以另一個矩陣

主要問題是在dim1 == dim和matrix不是square的情況下。所有應該用C++編寫。我的代碼返回幾乎好的結果,但它看起來像在我傳遞垃圾值的第一行。看起來我沒有填充這些值。我會很感激任何幫助。謝謝

void multiply(int* a[], int* b[], int* c[], int dim1, int dim2, int 
     dim3){ 

    for(int i = 0; i < dim1; i++){ 
     for(int j = 0; j < dim3; j++){ 
      for(int k = 0; k < dim2; k++){ 
       *(*(c+i)+j) += *(*(a+i)+k) * *(*(b+k)+j); 
      } 
     } 
    } 

    for(int i = 0; i < dim1; i++){ 
     for(int j = 0; j < dim3; j++){ 
      cout << *(*(c+i)+j) << "\t"; 
     } 
     cout << endl; 
    } 
} 

int main(){ 
    const int dim1 = 3; 
    const int dim2 = 2; 
    const int dim3 = 3; 


    int a[dim1][dim2] = {{1,2},{3,4},{5,6}}; 
    int b[dim2][dim3] = {{1,2,3},{4,5,6}}; 
    int c[dim1][dim3] = {}; 

    int *tab1 [dim1]; 
    int *tab2 [dim2]; 
    int *tab3 [dim3]; 


    for(int i = 0; i < dim1; i++){ 
     *(tab1+i) = *(a+i); 
     *(tab2+i) = *(b+i); 
     *(tab3+i) = *(c+i); 
    } 

    multiply(tab1,tab2,tab3,dim1,dim2,dim3); 
}                   
+1

你應該在'k'循環之前設置'*(*(c + i)+ j)= 0'。 – LutzL

回答

0

除了tab初始化,一切都正確。你迭代到dim1=3,但tab2dim2=2大小。試試這個:

const int dim1 = 3; 
const int dim2 = 2; 
const int dim3 = 3; 

void multiply(int a[dim1][dim2], int b[dim2][dim3], int c[dim1][dim3], int dim1, int dim2, int 
     dim3){ 

    for(int i = 0; i < dim1; i++){ 
     for(int j = 0; j < dim3; j++){ 
      for(int k = 0; k < dim2; k++){ 
       *(*(c+i)+j) += *(*(a+i)+k) * *(*(b+k)+j); 
      } 
     } 
    } 

    for(int i = 0; i < dim1; i++){ 
     for(int j = 0; j < dim3; j++){ 
      cout << *(*(c+i)+j) << "\t"; 
     } 
     cout << endl; 
    } 
} 

int main() 
{ 
    int a[dim1][dim2] = {{1,2},{3,4},{5,6}}; 
    int b[dim2][dim3] = {{1,2,3},{4,5,6}}; 
    int c[dim1][dim3] = {}; 

    multiply(a,b,c,dim1,dim2,dim3); 
}