#include /* * Note: * On multiplication of matrices, * A column with B row is must equal. * * New update: * Function concept, new algorithm, and now * support generic types for matrix */ template void multiplication(T (&R)[maxR][maxOtherC], T (&A)[maxR][maxC], T (&B)[maxOtherR][maxOtherC]); template void printArray(T (&arr)[maxR][maxC]); int main() { /* * Example matrices (int and float) * result1 = A * B * result2 = C * D */ int A[4][3] = { {1, 2, 4}, {1, 3, 3}, {3, 4, 2}, {3, 2, 1} }; int B[3][5] = { {1, 2, 2, 4, 1}, {4, 3, 1, 3, 3}, {1, 1, 3, 5, 2}, }; float C[4][3] = { {2.5f, 3.0f, 1.5f}, {1.5f, 2.0f, 2.5f}, {2.1f, 1.5f, 3.5f}, {2.5f, 2.5f, 4.5f} }; float D[3][2] = { {1.5f, 1.7f}, {2.0f, 2.0f}, {2.9f, 3.1f} }; //(result row = A row) and (result column = B column) int result1[4][5] = {}; multiplication(result1, A, B); float result2[4][2] = {}; multiplication(result2, C, D); std::cout << "Result1 :\n"; printArray(result1); std::cout << "\nResult2 :\n"; printArray(result2); return 0; } template void multiplication(T (&R)[maxR][maxOtherC], T (&A)[maxR][maxC], T (&B)[maxOtherR][maxOtherC]) { int r, c, otherR, otherC; for(r = 0; r < maxR; r++) { for(otherC = 0; otherC < maxOtherC; otherC++) { for(c = 0, otherR = 0; otherR < maxOtherR; otherR++, c++) { R[r][otherC] += A[r][c] * B[otherR][otherC]; } c = 0; otherR = 0; } } } template void printArray(T (&arr)[maxR][maxC]) { for(int r = 0; r < maxR; r++) { for(int c = 0; c < maxC; c++) { std::cout << arr[r][c] << " "; } std::cout << std::endl; } }

Multiplication of Mаtrices [Update]