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

2楼: Originally posted by snoopyzhao at 2012-07-13 12:50:41
什么是张量相乘，向量的点乘？

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

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

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

7楼: Originally posted by supersunyan at 2012-07-13 16:33:14
forall语句好像只能赋值，也没有想到什么好的方法。

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