我试图了解用于实习的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
,该变量指向类型为real
的kind_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
的目标开始,在前两个维度的索引ibx
和iby
处,它包含第三个维度的第二个值,从索引2开始并在flds%nc * flds%nl
停止。输出数组中的前两个维被忽略了,因为这些维中的索引是简单的标量。该表达式不能在语句中单独使用,因此我假设您在问题中遗漏了该行的某些部分。