OpenACCで2次元配列と構造体配列の受け渡しができません

Question

前提・実現したいこと

運動シミュレーションをOpenACCで並列化を試みています．

単純な並列化

double p=0;
for(lon long i=0;i<10000000000;i++){
p+=1;
}
cout<<p;

では処理時間短縮が確認され動作環境には問題ないようですが，csvから作った2次元ベクトルの成分や構造体配列参照しようとするとエラーを返されてしまいます．
例として以下のようなコードを考えます．
構造体配列のa[0]の成分p.vxに並列化してひたすらBr[4][4]を足し続けるだけの至ってシンプルな内容です．

データ移動ができていないのかなと思い,こちらの記事を参考にcopyinやpresentなどいろいろ試しましたが結果は変わりませんでしたし，そもそもGenerating implicit copy(Br,a)と書かれているのでそれも関係なさそうです．

ご教示よろしくお願いします．

補足情報（FW/ツールのバージョンなど）

GPUはnvidiaのtesla，nvc++は最新版にしてあります．
コンプラ絡みの諸事情で解決後に質問文は大幅編集しました。ご容赦ください。

Comments

[thkana]

https://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q13233822608

どんなに慌てているのかはわかりませんが、ルールとして https://teratail.com/help#posted-otherservice ということにはなっているので然るべく処置してください。

Answer 1

色々試した結果vector型が対応していないようです．
a[]とかa[][]みたいのは通って計算結果確認もできたので，うまい方法が分かるまではそちらで代用します

Answer 2

警告は知りませんけど、Segmentation faultの方は

int N; //静的寿命の変数なので初期値0

vector<agent> a(N);//Nには特に値が与えられて無いので0
であれば
a[0].p.vxにアクセスしようとしたら破綻するような気がします。

#include <vector>
int N;
struct coordinate {
    double x;
    double y;
    double z;
    double vx;
    double vy;
    double vz;
};

struct agent {
    int category;
    struct coordinate p;
};//categoryとxyz方向の6成分を表す構造体

int main() {
  std::vector<agent> a(N);
  a[0].p.vx=100;
  return 0;
}

これで再現。

---

質問が変わっちゃって、全然知らない分野に突入しそうだけど...

そもそも

   a[0].p.vx=0;
   for(long long i=0;i<100;i++){
    a[0].p.vx+=Br[4][4];
  }

は
a[0].p.vx+=Br[4][4]*100;
と等価ではないのですか?
等価なら並列で計算する意味はなさそうだし、等価でない状況なら並列計算が成り立たないのでは。

Comments

[Kinsho]

質問前に慌てていてつい誤った文をコンパイルしていました．
エラー文とともに訂正いたしましたのでよろしくお願いします．

[Kinsho]

そうですね．ただこれはサンプルとして作っただけで，実際はいろんなデータを大量に足し合わせるイメージです．

元のコードをそのまま掲載すると1000行超えてて誰も読まないと思うので．．．

[Kinsho]

また等価でない状況がーとの懸念ですが，一番上に掲載したようなひたすら1を足し合わせるやつを10^13回でやってみても1秒以内で終わったので内部でうまく並列化してくれているようです．（並列化したやつを最後に足し合わせているとかだと思います）

[thkana]

> サンプルとして作っただけで
サンプルとして役に立っていないのでは...

並列計算に詳しくないのですが
p=p+1
みたいに、これまでの結果が今回の演算に影響するようなものはそもそも並列計算に出来ないと思うのだけれど。n回めの計算とn+1回めの計算が独立だから並列化できるんじゃないのですか? (これぐらい簡単に論破出来ないようだったらちょっと苦しいんじゃないかと思う。急いでいるなら相手しなくていいですけれど)