Fortran OpenMPはどこに配置されますか

Question

Fortran-OpenMPと割り当て可能な配列に関する質問があります。それは簡単です：スペースはどこに割り当てられますか？私のようなものがある場合

!$omp parallel default(shared) private(arr) 
!$omp critical 
    allocate(arr(BIGNUMBER)) 
!$omp end critical 

!do calculations with many arr accesses 

!$omp critical 
    deallocate(arr) 
!$omp end critical 
!$omp end parallel 

スペースはスタックまたはヒープに割り当てられますか？それは、ヒープ上にある場合は、上記のコードと最初のコードでは、この

allocate(arr(BIGNUMBER, nThread)) 
!$omp parallel default(shared) private(localArr) 
    iThread = omp_get_thread_num() 

    localArr => arr(:, iThread) 

    !do calculations with many localArr accesses 
!$omp end parallel 

deallocate(arr) 

• のようなものとの間に差がある、2つの重要な領域が存在します。彼らは実行を遅くし、うまくスケールしないと私は仮定します。 （私は実際には、割り振りがスレッドセーブであるため、それらを残しておくことができないのかどうかは実際は分かりません）。しかし、配列がスタックに割り当てられていれば、高速アクセスのために高速になるはずです。
• 2番目のコードでは、配列がヒープ上にあることを確認していますが、これはアクセスが遅くなります。しかし、最初のコードの配列がヒープ上にも割り当てられていれば、クリティカルな試薬を保存します+それはただ1つの割り当て/割り当て解除です。速いはずですか？
• 配列のサイズはこれに任意のロールですか？
• ヒープ上に割り当てられる場合は、スタックに強制的に割り当てる方法がありますか？

短い質問は基本的に問題の最適解と思われるでしょうか？

Answer 1

OpenMPは、スタック上に自動変数（配列を含む）を割り当てる傾向があります。あなたが明示的な割り当てを行うとき、私はそれらがヒープに割り当てられると仮定しますが、Fortran標準はスタックやヒープについて全く話さないことに注意してください。例： number 1プライベート変数を割り当てているので、クリティカルセクションをそのまま残しておきます。サイズに関しては、自動配列が大きすぎるためにスタックオーバーフローが発生することがありますが、これはおそらくあなたのケースではありません。私が知らない最速のアプローチは何ですか？

編集：このプログラムは、スタック上のヒープ

integer,parameter :: BIGNUMBER = 100000000 
real,dimension(:),allocatable :: arr 

!$omp parallel default(shared) private(Arr) 
allocate(arr(BIGNUMBER)) 
    iThread = omp_get_thread_num() 

    arr = 5 

    print *, arr 

deallocate(arr) 
!$omp end parallel 


end 

と、このいずれかに配列を割り当てるように思わ

integer,parameter :: BIGNUMBER = 100000000 
real arr(BIGNUMBER) 

!$omp parallel default(shared) private(Arr) 
    iThread = omp_get_thread_num() 

    arr = 5 

    print *, arr 

!$omp end parallel 


end 

Answer 2

OKは、ウラジミールは、ほとんどの私はだろうと言う（それが失敗しました） （標準では言及されていませんが、それは実装に完全に対応していますが、どうしてあなたの秘密を守るためにクリティカルを使用していますか？）

しかし、あなたは、スタックに割り当てられたメモリへのアクセスがヒープ上のメモリへのアクセスよりもいくらか速いと考える印象を与えます。典型的な実装では、そうではありません。アクセス時間は同じです。 スタック上のメモリの割り当ては通常、ヒープ上よりも高速ですが、割り当てられた後はアクセスする時間は同じです。だから、私はクリティカルをカットしてルート1に行きます。プライベートなものを保つのは簡単ですし、ポインタが悪いですし、メモリ割り当て時間が限られている場合は、パラレル領域を並列化するには十分な作業がほとんどありませんそれは価値がある。