fortran 中 矩阵与向量计算问题

liudy02

妙木山养娃人 发表于 2019-6-23 17:04
您好！您这个是用列*行得到一个3*3的数组，好像计算没问题。但如果调换顺序，计算行*列的时候，结果是14 ...

这就是Fortran这类强类型语言的一个大坑在此
1*3的矩阵和3*1的矩阵得到的东西，程序认为是一个1*1的矩阵而不是一个数
程序并不认为1*1的矩阵和一个数是一回事
对此，你可以设一个1*1的矩阵变量接收矩阵乘法的结果，然后再取这个矩阵第一行第一列的元素即可

liudy02

liudy02 发表于 2019-6-23 17:37
这就是Fortran这类强类型语言的一个大坑在此
1*3的矩阵和3*1的矩阵得到的东西，程序认为是一个1*1的矩阵 ...

这个坑其实还不算常见，更常见的坑是n*1的或者1*n的矩阵和n个元素的向量不是一回事
初学者在这个地方被坑的更多些……

liudy02

fcode 发表于 2019-6-23 17:42
可以用
k=sum(matmul(a,transpose(a))) 呀

唔，还是你这个办法简洁，我咋就没想到呢……