fortran 中 矩阵乘法

 

001、

[root@rocky01 fortran]# ls
test01.f90
[root@rocky01 fortran]# cat test01.f90       ## 测试程序
program test01_fortran
        implicit none
        integer :: ay1(2,2)
        integer :: ay2(2,2)
        integer :: ay3(2,2)
        ay1 = reshape((/2,3,4,8/), (/2,2/))    ## 二维数组赋值
        ay2 = reshape((/1,0,2,4/), (/2,2/))    ## 二维数组赋值
        ay3 = matmul(ay1,ay2)                  ## 矩阵相乘
        print *, ay3(1,:)
        print *, ay3(2,:)
end program test01_fortran
[root@rocky01 fortran]# gfortran test01.f90 -o kkk && ./kkk
           2          20
           3          38

 

。