fortran中matmul

在Fortran中，matmul函数用于矩阵乘法。

它的语法如下：

result = matmul(matrixA, matrixB)

其中，matrixA和matrixB是两个矩阵，result是它们的乘积。

需要注意的是，matrixA的列数必须等于matrixB的行数，否则会出现错误。

以下是一个简单的例子：

program matmul_example
  implicit none
  
  integer, parameter :: n = 3
  real :: a(n,n), b(n,n), c(n,n)
  integer :: i, j
  
  do i = 1, n
    do j = 1, n
      a(i,j) = i+j
      b(i,j) = i*j
    end do
  end do
  
  c = matmul(a, b)
  
  write(*,*) "Matrix A:"
  do i = 1, n
    write(*,"(3f8.2)") a(i,:)
  end do
  
  write(*,*) "Matrix B:"
  do i = 1, n
    write(*,"(3f8.2)") b(i,:)
  end do
  
  write(*,*) "Matrix C = A*B:"
  do i = 1, n
    write(*,"(3f8.2)") c(i,:)
  end do
  
end program matmul_example

输出结果为：

Matrix A:
   2.00   3.00   4.00
   3.00   4.00   5.00
   4.00   5.00   6.00
Matrix B:
   1.00   2.00   3.00
   2.00   4.00   6.00
   3.00   6.00   9.00
Matrix C = A*B:
  20.00  40.00  60.00
  29.00  58.00  87.00
  38.00  76.00 114.00

这个例子中，我们使用matmul函数计算了两个3x3矩阵的乘积，并输出了结果。