我试图了解用于实习的Fortran代码,但我是一名C ++程序员,但我被困在此代码位的作用。有人可以启发我吗?

real(kind=kind_real), pointer :: c(:,:,:)  (I think this created a pointer to a matrix )

self%c => self%gfld3d(:,:,ioff+1:) (this is in the constructor of the class)

flds%c(ibx,iby,2:flds%nc*flds%nl:2) (this is a function that calls the c variable of the class but I don't understand how the parameters work and what does the colon operator do as a a parameter)

最佳答案

real(kind=kind_real), pointer :: c(:,:,:)


正如您已经写过的,此行声明了一个指针变量c,该变量指向类型为realkind_real数字的3维数组。这里的冒号表示每个尺寸的大小未知。假设它在type块内,则c是该类型的成员。请注意,在Fortran中,()既用于数组下标又用于函数调用。


self%c => self%gfld3d(:,:,ioff+1:)


在构造函数中,为指针分配了一个目标,该目标是3维数组self%gfld3d的一部分,该数组也是相同类型的成员。 Fortran中的指针包含有关它们指向的数组切片的尺寸的信息。在这种情况下,前两个冒号表示数组切片跨越了目标数组的整个前两个维度。 ioff+1:表示在第三维中,指针目标从给定值self%gfld3d开始包括ioff+1中的所有值。这是常规切片符号的简化情况,其中a(start:end:step)表示数组a的切片,从索引start开始,以end结束,以step的步长。


flds%c(ibx,iby,2:flds%nc*flds%nl:2)


基于以上所述,这不是函数调用。它是一个表达式,计算结果为一维非连续数组切片。从flds%c的目标开始,在前两个维度的索引ibxiby处,它包含第三个维度的第二个值,从索引2开始并在flds%nc * flds%nl停止。输出数组中的前两个维被忽略了,因为这些维中的索引是简单的标量。该表达式不能在语句中单独使用,因此我假设您在问题中遗漏了该行的某些部分。

07-27 15:12