fortran 中 矩阵与向量计算问题

妙木山养娃人

implicit none
integer::a(3,1)

a(1,1)=1
a(2,1)=2
a(3,1)=3

假如有一列向量a，应该如何计算a向量的转置乘a。我想用matmul(transpose(a),a)但不可以，请问该如何做呢？

妙木山养娃人

fcode 发表于 2019-6-23 16:06
[mw_shl_code=fortran,true]integer::a(1,3)
a(1,:)=[1,2,3]
write(*,*) matmul(transpose(a),a)[/mw_shl_c ...

您好！您这个是用列*行得到一个3*3的数组，好像计算没问题。但如果调换顺序，计算行*列的时候，结果是14，直接write是可以的，但用变量表示就不行。请问是什么原因呢？
integer::k
k=matmul(a,transpose(a))
write(*,*) k

提示错误The shapes of the array expressions do not conform.   [K]       

妙木山养娃人

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

好的，十分感谢

妙木山养娃人

liudy02 发表于 2019-6-23 17:38
这个坑其实还不算常见，更常见的坑是n*1的或者1*n的矩阵和n个元素的向量不是一回事
初学者在这个地方被坑 ...

谢谢 很惭愧，以前学程序很马虎，也不常用，现在需要自己编程的时候才发觉捉襟见肘