Умножение матриц

typedef int mat44[4][4];
 
void mult_mat4(mat44 factor1, mat44 factor2, mat44 result)
{
    printf("\n");
    for(int i = 0; i < 4; i++) {
        printf("i=%i\n", i);
        for(int j = 0; j < 4; j++) {
        printf("j=%i\n", j);
            for(int k = 0; k < 4; k++) {
                printf("k=%i\n", k);
                result[i][j] += factor1[i][k] * factor2[k][i];
                printf("result=%i\n", result[i][j]);
            }
            
        }
    }
}

#include <iostream>
using namespace std;
 
void function(int arr_1[4][4], int arr_2[4][4]){
    int arr_3[4][4];
    for(int i=0;i<4;i++){
        for(int j=0;j<4;j++){
            if(arr_1[i][j]=arr_2[i][j]){
                arr_3[i][j]=arr_1[i][j]*arr_2[i][j];
        
            }
        }
    }
    for(int i=0;i<4;i++){
        for(int j=0;j<4;j++){
            cout<<arr_3[i][j]<<" ";
        }
        cout<<endl;
    }
 
}
 
int main(){
    int arr_1[4][4];
    int arr_2[4][4];
    for(int i=0;i<4;i++){
        for(int j=0;j<4;j++){
            cin>>arr_1[i][j];
        }
    }
    for(int i=0;i<4;i++){
        for(int j=0;j<4;j++){
            cin>>arr_2[i][j];
        }
    }
    
    function(arr_1,arr_2);
    
    system("pause");
    return 0;
}

void mult(float a[4][4], float b[4][4], float c[4][4])
{
int i,j,k;
for (i = 0; i < 4; i++)
   for (j = 0; j < 4; j++)
   { 
      c[i][j] = 0; // чтобы было к чему прибавлять, если массив "с" исходно забит мусором
      for (k = 0; k < 4; k++)
         c[i][j] += a[i][k]*b[k][j];
   }   
}