指针的 Fortran 副本

Question

我有一个类型 var，它包含一个指针 p。 我需要将 var 复制到与 var 相同类型的另一个变量 var1 上（通过做 var1 "=" var，在引号中，因为我不知道这是否正确，见下文）。

在我的实现中 var 和 var1 被传递给一个需要将 var1 初始化为 var 然后修改 var1 的子例程。 如果我改变了 var1%p 的值，那么 var%p 也会被修改。 所以，我需要改变var1%p指向的内存区域的值，而不修改var%p指向的内存区域。 在 Fortran (2003, Intel) 中这样做的正确方法是什么？

代码示例：

type variable 
     real, dimension(:), pointer:: p
end type variable


subroutine init(var,var1) !This of course will not work
  type (variable):: var, var1
  var1=var
end subroutine

现在，我需要能够在不影响 var%p 的情况下执行与更改 var1%p 等价的操作。我无法更改变量（我正在修改现有代码）。

---

以下解决方案（如果我是正确的）： 我将 var1 声明为：

Type(variable), allocatable:: Var1 

然后分配和初始化：

allocate(var1, source=var) 
var1 = var  

然后

call somesub(var1)

修改 var1%p 的结果是修改 var%p 使得 var%p 等于 var1%p。

编辑 3： 正在做：

 subroutine init(var,var1)
 type(variable), intent(in):: var
 type(variable), allocatable, intent(inout):: var1
 allocate(var1%p, source=var%p)
 var1%p = 2 
 end subroutine init

 type variable:: var
 type variable, allocatable:: var1
 call init(var, var1)

以上编译但崩溃（无限挂起）。 改为使用：

allocate(var1, source=var) 

运行但将 var%p 更新为 1。

Answer 1

当您对派生类型变量进行内在赋值时，任何指针都会使用指针赋值=> 复制（它只是复制现有内存块的地址和其他属性）。您必须确保在内存中分配了一个新目标，并将数组的值复制到那里。

您可以为给定类型创建自己的重载赋值，然后使用语法var1 = var1，或者您可以将正确的赋值编码到您的子例程中，如果足够的话：

subroutine init(var,var1) !This of course will not work
 type(variable) :: var, var1  !DO NOT FORGET THIS WHEN POSTING EXAMPLES NEXT TIME!

 allocate( var1%p(lbound(var%p,1):ubound(var%p,1)) )

 var1%p = var%p
end subroutine

正如 HighPerformance Mark 评论的那样，如果您没有其他理由使用指针，那么可分配组件会好得多。

type variable 
     real, dimension(:), allocatable :: p
end type variable


subroutine init(var,var1) !This WILL work
  type (variable):: var, var1
  var1 = var
end subroutine

---

在您的编辑中，您遇到了同样的问题：

var1 = var 

这和之前的问题一模一样，你什么都没改！

还有

allocate(var1, source=var) 

应该和原代码一样的效果，相信你需要

allocate(var1%p, source=var%p) 

但我已经评论说，最常见的 gfortran 版本不会接受这一点，并且需要明确指定形状。

变量var 和var1 不必是可分配的。但如果出于其他原因需要，它们也可以。

---

可以编写一次用户定义的赋值，然后重载内部赋值 (=)。

   type variable
     real :: pointer :: p(:)
   contains
     procedure :: assign
     generic :: assignment(=) => assign
   end type

   subroutine assign(out, in)
      class(variable), intent(out) :: out
      class(variable), intent(in) :: in

      allocate( out%p(lbound(in%p,1):ubound(in%p,1)) )

      out%p = in%p
   end subroutine

然后你可以随时写var1 = var，它会自动调用子程序assign。