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

问题描述

我是 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,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%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 上都试过这个，都产生了相同的结果。

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

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

fortran gfortran intel-fortran