如何用gfortran实现包含参数化派生类型的自定义类型？

gfortran在处理嵌套延迟长度参数化派生类型时存在部分内部编译错误（ICE），以下是符合Fortran标准的实现方案，同时避开这些问题：

核心实现代码 #

module param_types
    implicit none

    ! 参数化方阵类型：N为方阵阶数
    type :: sqmat_t(N)
        integer, len :: N
        real :: mat(N,N)
    end type sqmat_t

    ! 包含两个同阶方阵的嵌套参数化类型
    type :: two_sqmat_t(M)
        integer, len :: M
        type(sqmat_t(M)) :: mat1, mat2
    end type two_sqmat_t

contains

    ! 初始化延迟参数化类型变量的子程序
    subroutine init_two_sqmat(two_mat, order, val1, val2)
        type(two_sqmat_t(:)), allocatable, intent(out) :: two_mat
        integer, intent(in) :: order
        real, intent(in) :: val1, val2

        allocate(two_sqmat_t(order) :: two_mat)
        two_mat%mat1%mat = val1
        two_mat%mat2%mat = val2
    end subroutine init_two_sqmat

end module param_types

program test_param_types
    use param_types
    implicit none

    type(two_sqmat_t(:)), allocatable :: dyn_twosqmat
    type(two_sqmat_t(2)) :: fixed_twosqmat
    integer :: i, j

    ! 初始化动态参数变量
    call init_two_sqmat(dyn_twosqmat, 3, 1.0, 2.0)

    ! 初始化固定参数变量（用循环替代直接数组赋值，规避ICE）
    do i = 1, 2
        do j = 1, 2
            fixed_twosqmat%mat1%mat(i,j) = 1.0
            fixed_twosqmat%mat2%mat(i,j) = 3.0
        end do
    end do

    ! 验证输出
    print *, "动态参数mat1第一行：", dyn_twosqmat%mat1%mat(1,:)
    print *, "固定参数mat1第一行：", fixed_twosqmat%mat1%mat(1,:)

    if(allocated(dyn_twosqmat)) deallocate(dyn_twosqmat)

end program test_param_types

关键注意事项 #

• 延迟参数变量必须先分配再使用：直接声明type(two_sqmat_t(:)), allocatable :: var后，必须通过allocate(two_sqmat_t(阶数) :: var)指定参数值，否则会触发gfortran的ICE；
• 规避直接数组赋值：对固定参数类型的内部数组直接做var%mat1%mat(:,:) = 1会触发ICE，改用显式循环或者子程序初始化；
• 封装初始化逻辑：将初始化逻辑放到模块子程序中，既符合Fortran模块化编程规范，也能绕过编译器的bug。

关于gfortran的编译错误 #

你遇到的f951: internal compiler error: Segmentation fault属于gfortran对嵌套参数化类型处理的已知bug，尤其是在直接赋值嵌套参数化类型的数组成员时。上述代码通过间接初始化的方式完全规避了这些问题，同时兼容Fortran 2008及以上标准。

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