C++计算矩阵乘法

遇到一个矩阵乘法的问题，目的是计算A'A和AA'

设A 是3行3列的矩阵

计算AA'如下

int sk[9]={0};

for (i=0;i<3;i++)

{

for (j=0;j<3;j++)

{

for (k=0;k<3;k++)

{

sk[i*3+j]+=a[i*3+k]*a[j*3+k];

}

}

}

计算A'A如下：

for (i=0;i<3;i++)

{

for (j=0;j<3;j++)

{

for (k=0;k<3;k++)

{

sk[i*3+j]+=a[k*3+i]*a[k*3+j];

}

}

}