应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » fortran如何实现张量相乘（并乘）
101/1返回列表
查看: 4753  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

btx97

金虫 (小有名气)

[交流] fortran如何实现张量相乘（并乘）已有3人参与


如上的计算公式，Fortran是否有不用DO循环就能实现的办法。
不想用DO的原因是我的计算中会出现高阶张量，如果用DO，嵌套（用的对否？）次数最多会达到十二次。顺便也问下FORTRAN里DO嵌套的效率怎么样？有人用过12次，或7次的吗，会严重影响效率吗？

[ Last edited by btx97 on 2012-7-13 at 10:40 ]
回复此楼

» 收录本帖的淘帖专辑推荐

资源收集

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-07-13 17:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
什么是张量相乘，向量的点乘？
赞一下 回复此楼
2楼2012-07-13 19:50:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)
引用回帖:
2楼: Originally posted by snoopyzhao at 2012-07-13 12:50:41
什么是张量相乘，向量的点乘？

不是点乘，点乘的表达式是这样的：
两个二阶张量点乘结果是二阶张量，两个二阶张量并乘是一个四阶张量。如果是两个一阶的张量并乘其结果是二阶张量，可以跟列向量a通过得到矩阵进行类比.
你知道有什么办法可以避免太多的循环嵌套吗？

[ Last edited by btx97 on 2012-7-13 at 13:27 ]
赞一下 回复此楼
3楼2012-07-13 20:21:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
没有计算公式不好说。如果计算不能向量化，那么基本上也只能用循环了……
赞一下 回复此楼
4楼2012-07-13 20:47:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
看看下面这段代码与你的这个计算有没有相关性，至少不需要循环，呵呵……
CODE:
program test

   integer, dimension(5) :: a = (/1,2,3,4,5/)
   integer, dimension(5) :: b = (/5,4,3,2,1/)
   integer, dimension(5) :: c
   integer, dimension(3, 3) :: d
   integer, dimension(3, 3) :: e
   integer, dimension(3, 3) :: f
   
   d = reshape((/ 1, 2, 3, 4, 5, 6, 7, 8, 9 /), shape(d))
   e = reshape((/ 1, 2, 3, 4, 5, 6, 7, 8, 9 /), shape(e))
   
   c = a * b
   f = d * e
   write(*,*) c
   write(*,*) f

end program test
赞一下 回复此楼
5楼2012-07-13 21:09:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)
引用回帖:
4楼: Originally posted by snoopyzhao at 2012-07-13 13:47:01
没有计算公式不好说。如果计算不能向量化，那么基本上也只能用循环了……

谢谢！ 计算公式的话就是上面的，算这个除了用下面这样的，还有别的办法吗？
CODE:
do i=1,3
    do j=1,3
        do k=1,3
           do l=1,3
                 c(i, j, k, l)=a(i, j)*b(k, l)
           end do
         end do
    end do
end do
赞一下 回复此楼
6楼2012-07-13 21:38:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

supersunyan

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
forall语句好像只能赋值，也没有想到什么好的方法。
赞一下 回复此楼
默默奋斗，实现梦想
7楼2012-07-13 23:33:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)
引用回帖:
7楼: Originally posted by supersunyan at 2012-07-13 16:33:14
forall语句好像只能赋值，也没有想到什么好的方法。

看来fortran不支持并乘计算啊。
赞一下 回复此楼
8楼2012-07-15 04:45:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whu_zyg

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这个四重循环觉得应该可以减少为两重循环，因为两个向量直接相乘（product命令）就是求和了，你可以试试
赞一下 回复此楼
9楼2012-07-16 00:25:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)
引用回帖:
9楼: Originally posted by whu_zyg at 2012-07-15 17:25:16
这个四重循环觉得应该可以减少为两重循环，因为两个向量直接相乘（product命令）就是求和了，你可以试试

谢谢参与。不过这里没有求和的，是向量中的每个元素与另一向量中的每个元素相乘的算法。
PS: 才发现你前半句是对的，是可以减少为两重

[ Last edited by btx97 on 2012-7-17 at 19:53 ]
赞一下 回复此楼
10楼2012-07-16 01:31:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 btx97 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭