指针的 Fortran 副本

我有一个包含指针 p 的 var 类型。 我需要将 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

现在，我需要能够做一些相当于更改 var1%p 而不影响 var%p 的事情。我无法更改变量（我正在修改现有代码）。

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

Type(variable), allocatable:: Var1 

然后分配并初始化：

allocate(var1, source=var) 
var1 = var  

Then

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。

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

您可以为给定类型创建自己的重载赋值，然后使用语法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

正如 High Performance 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 

这和以前的问题一模一样，你没有改变任何东西！

Also

allocate(var1, source=var) 

应该与原始代码具有相同的效果，我相信你需要

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

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

变量var and 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自动地。