Fortranのディープアレイコピー

Question

Fortran（90）で（実数）配列のディープコピーが必要ですが、参照の仕組みを完全に理解していないため、取得方法がわかりません。直感的に、私はこれは私が欲しいものを私を得ることを期待する：

do i=1,n 
    b(i) = a(i) 
end do 

しかし、最近b(1:n) = a(1:n)は上記のコードと同等であることを私に指摘されました。直観的には、b(1:n) = a(1:n)は、b(1:n)の参照がメモリ内のa(1:n)の位置を指し示すに過ぎないと予想します。

b(1:n) = a(1:n)ディープコピーですか？どうして？ b = aとは対照的に、根本的な参照で何が起こっていますか？

Answer 1

三つの方法、doループ、b(1:n) = a(1:n)とb = aを参照してください、すべての等価です。アレイaの内容をアレイbにコピーします。 aとbは単純な配列であり、空想的なポインタなどではないので、代入a = bは基本的に数式と同じです。参考までに（ユーザーが知る必要のある）魔法はありません。そのため、Fortranは学ぶためのかなり簡単な言語です。あなたはFortranでポインタ配列を使うことができますが、これはまったく別の問題です。

M MetcalfとJ ReidのFortran 90/95で説明されているは、Fortranの言語機能に関する相談には常に参考になります。ページ48から：

3.11固有の割り当てによって

アレイ割り当て、配列式は、式の各要素が割り当てられたかのように解釈されるものと同じ形状の配列変数に割り当てることができます変数の対応する要素に追加します。例えば、宣言と

real, dimension(10, 20) :: a 

割り当て

a = a + 1.0 

はi=1,2..,10とj=1,2,..,20ためa(i,j) + 1.0によってa(i,j)を置き換えます。

もスカラー式はsaclar値は、すべての配列要素にブロードキャストされる場合には、アレイに割り当てられてもよいことに留意されたいです。

これはすべて実際にどのように実装されているかという点では、私はあなたがあなたの疑問を抱いていると思っていますが、これはFortran標準によって完全には特定されていません。この種のことは、コンパイラ作成者がどのような最適化を行っても同じように行うことができるように、未定義のままです。たとえば、a = bという割り当てでは、bの要素がaにコピーされる順序は標準では指定されていないため、さまざまな方法で異なるコンパイラがこれを行うことができます。 aとbがポインタでない場合、aとbは別個の配列であり、一方の要素を変更しても他方の対応する要素は変更されないということです。ですから、a=bは "ディープコピー"で、これはbのすべてのアイテムをaのメモリ位置にコピーしていると考えることができます。

Answer 2

a = bは、配列b全体をaにコピーします。 nより大きい次元の配列の一部だけが必要な場合は、添字表記a（1：n）= b（1：n）を使用できます。つまり、Fortran 90よりも高いレベルの言語です。 "a = b"はコピーであり、ポインタ "a"は演算子 "="を使用しているため、ターゲットbと関連付けません。ポインタの関連付けは=>を使用します。

編集：コピーすると、おそらくdoループと同じマシンコードで複製されます。ポインタの関連付けは、すべての配列要素を複製するdoループを実行せずに参照を行います。

を使用すると、配列のコピーに言及http://en.wikipedia.org/wiki/Fortran_95_language_features