fortran - 类型绑定赋值将成员重新初始化为零

Question

我是 fortran 的新手，并尝试使用动态内存分配和指针在其中实现一个矩阵类。当我试图重载赋值运算符时，我观察到对象成员在赋值子例程的入口处被重新初始化为零。以下是重现问题的简单程序。

module matModule

    implicit none
    ! Define the data kind for integer and real numbers
    integer, private, parameter :: iKind = selected_int_kind(8), rKind = selected_real_kind(15,307);

    type mat
        integer(kind=iKind), public :: nrows,ncols,npage;
        real(kind=rKind), pointer, public :: vals(:,:,:) => null();
    contains
        procedure, private, pass :: assign_scal_r
        generic :: assignment(=) => assign_scal_r
    end type mat

    contains
        subroutine assign_scal_r(this,value)
            implicit none
            real, intent(in) :: value
            class(mat), intent(out) :: this
            integer(kind=iKind) :: err
            print *, loc(this%nrows), loc(this%ncols), loc(this%npage)
            print *, associated(this%vals), this%nrows, this%ncols, this%npage
            this%nrows = 1; this%ncols = 1; this%npage = 1;
            if ( associated(this%vals) ) then
                print *, "The array on the left is already allocated."
                deallocate(this%vals, stat=err)
                if (err /= 0) print *, "Error deallocating the variable."
            end if
            allocate(this%vals(this%nrows,this%ncols,this%npage), stat=err)
            if (err /= 0) print *, "Error allocating the variable."
            this%vals = value;
        end subroutine
end module matModule

program main
    use matModule
    implicit none
    type(mat) :: k
    k = 5.0;
    k = 1.0;

end program main

这段代码的输出是

140700383959168       140700383959172       140700383959176
F           0           0           0
140700383959168       140700383959172       140700383959176
F           0           0           0

请注意，输出是从赋值子例程中生成的。主程序中的变量“k”经过两次赋值操作。可以理解的是，成员在第一次分配期间默认初始化为零。但是，当分配子例程第二次执行时，成员包含 0。我还打印了成员的内存地址，据我了解，这表明成员的值在调用分配子例程期间被覆盖。

我不确定这是否是预期的行为。如果这是预期的行为，那么我的问题是指针是否会在赋值子例程的条目上自动释放？我在 Intel Fortran 编译器和 gfortran 上都试过了，都产生了相同的结果。

如果这整个问题是由于我在某处的代码中犯了一些愚蠢的错误引起的，请原谅我。

标签： fortrangfortranintel-fortran