2017-03-31 183 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); 
}