arrays - Fortran:从矩阵中提取列和行并将它们用作乘法的一维数组

Question

假设我有矩阵 A 和 B 大小都为 n x n 。现在我想将 A 的每一列视为列矩阵 ColA(i) [i=1,n] 并将 B 矩阵的每一行视为行矩阵 RowB(j) [j=1,n]。

现在我可以使用 matmul 函数构造一个名为 C 的新矩阵，其元素是通过将所有可能的 i 和 j 值与 ColA 和 RowB 相乘得到的吗？有点像

do i = 1, n
 do j= 1, n
    C(i,j) = matmul( ColA(i), RowB(j) ) 
 end do
end do

Answer 1

好吧，对于您的问题，我可以看到三个可能的答案，但我不确定您指的是哪一个...

• A 的 i 列和 B 的 j 行的点积，等于矩阵-矩阵乘法 B*A:

C = matmul( B,A )

• 外积:C=colA * rowB^T，但是每个组合 i,j 都会产生一个单独的矩阵 C。

C = matmul( A(:,i), B(j,:) )

• 逐元素乘法，这意味着 C_ij = A_ij * B_ij。这在 Fortran 中非常简单:

C = A*B

你选吧!