値を追加するときにFortranアレイが自動的に増加する

Question

Fortranで拡張アレイをエミュレートする既存の方法はありますか？ C++のベクトルと似ています。私はインターネット上でこのテーマについて何も見つけていないとき、非常に驚​​いていました。

動機付けの例として、いくつかの反復関係を計算し、私が得るすべての中間の数値を保存したいとします。私の停止基準は隣接する結果の違いですので、私はこれに割り当てるメモリ量を事前に知ることができません。

Answer 1

私はこれまでこのサイトのどこかに表示されていると確信していますが、私はそれを見つけることができません。 francescalusによってコメントとして

まず、Fortran 2003の中で、あなたは簡単な

a = [a, item] 

によって一つの要素を追加することができます。これは配列を頻繁に再割り当てする可能性が高く、遅くなります。

配列数は、要素数nよりやや大きめになるようにすることができます。要素数nが配列size(a)のサイズを超えると、新しい要素をいくつかの要素（ここでは2x）で割り当ててそこの古い要素をコピーできます。残念ながら、Fortranにはrealloc()はありません。

module growing_array 
    implicit none 

    real, allocatable :: a(:) 

    integer :: n 

contains 

    subroutine add_item(item) 
    real, allocatable :: tmp(:) 
    real, intent(in) :: item 

    if (n == size(a)) then 
     !this statement is F2003, it can be avoided, but I don't see why in 2016 
     call move_alloc(a, tmp) 

     allocate(a(n*2)) 
     a(1:n) = tmp 
    end if 

    n = n + 1 

    a(n) = item 
    end subroutine 
end module 

私は最初の割り当てを省いただけで十分です。

タイプバインドされたプロシージャを使用して派生型に入れることができ、データ構造として使用できますが、これは純粋なFortran 2003であり、90を望みました。したがって、Fortran 90には欠陥がある割り当て可能な配列のための多くの方法であり、必然的に時代遅れであり、本質的に死んでいる。