矩陣乘法運算

 

#include "stdafx.h"

void matrix_mutiply(int *m1, int *m2, int *r, int x, int y, int z);

int main(int argc, char* argv[])
{
    int m1[]={2, -6 , 3, 5, 1, -1};
    int m2[]={4, -2, -4, -5, -7, -3, 6, 7};
    int r[100];

    matrix_mutiply(m1, m2, r, 3, 2, 4);

    for (int t=0; t<3; t++)
    {
        for (int h=0; h<4; h++)
        {
            printf("%3d ", r[4*t+h]);
        }
        printf("\n");
    }

    return 0;
}

void matrix_mutiply(int *m1, int *m2, int *r, int x, int y, int z)
{
    for (int i=0; i<x; i++)
    {
        for (int j=0; j<z; j++)
        {
            r[i * z + j]=0;
            for (int k=0; k<y; k++)
            {
                r[i * z + j] = r[i * z + j] + m1[i * y + k] * m2[k * z +j];
            }
        }
    }
}