Fortran Coder

查看: 177|回复: 8

[求助] fortran 中 矩阵与向量计算问题

[复制链接]

6

帖子

2

主题

0

精华

入门

F 币
33 元
贡献
19 点
发表于 2019-6-23 14:21:49 | 显示全部楼层 |阅读模式
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)但不可以,请问该如何做呢?
回复

使用道具 举报

1429

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
784 元
贡献
577 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2019-6-23 16:06:13 | 显示全部楼层
[Fortran] 纯文本查看 复制代码
integer::a(1,3)
a(1,:)=[1,2,3]
write(*,*) matmul(transpose(a),a)

6

帖子

2

主题

0

精华

入门

F 币
33 元
贡献
19 点
 楼主| 发表于 2019-6-23 17:04:00 | 显示全部楼层
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]       

120

帖子

8

主题

0

精华

宗师

F 币
3110 元
贡献
520 点
发表于 2019-6-23 17:37:04 | 显示全部楼层
妙木山养娃人 发表于 2019-6-23 17:04
您好!您这个是用列*行得到一个3*3的数组,好像计算没问题。但如果调换顺序,计算行*列的时候,结果是14 ...

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

120

帖子

8

主题

0

精华

宗师

F 币
3110 元
贡献
520 点
发表于 2019-6-23 17:38:40 | 显示全部楼层
liudy02 发表于 2019-6-23 17:37
这就是Fortran这类强类型语言的一个大坑在此
1*3的矩阵和3*1的矩阵得到的东西,程序认为是一个1*1的矩阵 ...

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

1429

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
784 元
贡献
577 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2019-6-23 17:42:42 | 显示全部楼层
可以用
k=sum(matmul(a,transpose(a))) 呀

6

帖子

2

主题

0

精华

入门

F 币
33 元
贡献
19 点
 楼主| 发表于 2019-6-23 18:02:46 | 显示全部楼层
fcode 发表于 2019-6-23 17:42
可以用
k=sum(matmul(a,transpose(a))) 呀

好的,十分感谢

6

帖子

2

主题

0

精华

入门

F 币
33 元
贡献
19 点
 楼主| 发表于 2019-6-23 18:05:12 | 显示全部楼层
liudy02 发表于 2019-6-23 17:38
这个坑其实还不算常见,更常见的坑是n*1的或者1*n的矩阵和n个元素的向量不是一回事
初学者在这个地方被坑 ...

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

120

帖子

8

主题

0

精华

宗师

F 币
3110 元
贡献
520 点
发表于 2019-6-23 18:19:39 | 显示全部楼层
fcode 发表于 2019-6-23 17:42
可以用
k=sum(matmul(a,transpose(a))) 呀

唔,还是你这个办法简洁,我咋就没想到呢……
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2019-10-17 04:02

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表