Fortran转VBA：P(I,1)、P(1,J)含义及代码解析求助

核心背景先明确 #

P是一个NF×NF的二维数组，Fortran中二维数组采用列优先存储——即先存完第一列的所有元素，再存第二列，以此类推，同一行的相邻列元素在内存中会间隔NF个位置（因为每列有NF个元素）。

关于P(I,1) #

它是二维数组P的第I行、第1列的元素，但在调用INPR子程序时，它的作用是传递一个行向量的起始内存地址：
结合INPR的步长参数ia=NF，子程序内部会从这个地址开始，每次跳过NF个元素，循环NF次，最终遍历的是P数组的第I行全部元素（从P(I,1)到P(I,NF)），也就是一个长度为NF的行向量，作为点积运算的第一个输入向量。

关于P(1,J) #

它是二维数组P的第1行、第J列的元素：

• 如果是代码中实际出现的P(J,1)（推测你可能笔误），逻辑和P(I,1)完全一致——代表P数组的第J行全部元素，作为点积的第二个输入向量，运算结果存入C(I,J)，再通过C(J,I)=C(I,J)生成对称矩阵C。
• 如果你确实要让P(1,J)代表列向量（比如P的第J列），那INPR的步长参数ib需要改为1，此时子程序会遍历P(1,J)、P(2,J)……P(NF,J)，也就是第J列的全部元素，形成列向量参与点积。

转VBA的关键注意点 #

1. 存储顺序差异：VBA的二维数组是行优先存储，和Fortran的列优先完全相反，直接按索引转写会导致元素对应错误。比如Fortran中P(I,J)的内存位置对应VBA中P(J,I)（当数组维度相同时），需要调整索引逻辑。
2. INPR的转写：不需要模拟Fortran的步长寻址逻辑，直接在VBA中提取目标行/列的元素，循环计算乘积和即可。比如计算P第I行和第J行的点积，VBA代码可以写成：

Function CalcDotProduct(P As Variant, row1 As Integer, row2 As Integer, NF As Integer) As Double
    Dim dot As Double
    dot = 0
    For k = 1 To NF
        dot = dot + P(row1, k) * P(row2, k)
    Next k
    CalcDotProduct = dot
End Function

内容的提问来源于stack exchange，提问作者user3387046