OpenACCとdouble型の対応について

Question

前提・実現したいこと

c++でOpenACCを使っています．
こんな感じのN個の成分を持つdouble型の配列ｐに関してp[0],p[1],...と
並列化していってそれぞれに対して決められた回数だけ1を足し続けるという簡単な動作確認をしています．

ただここでdouble型を使っているはずがint型の境界2147483647を超すとそこでストップしてしまうようです．double型は100桁だろうと問題なく計算できるはずだと認識しており，10^10程度の計算がまさかできないというわけでもないでしょうし，また使用しているCPUやGPUは現行機の中でもスペックは高い方と思われますので処理限界というわけでもなさそうです．

OpenACCなどのGPUプログラミングは倍精度計算には対応していると聞いたはずですが,どうすれば10^10程度の数までの計算を扱えるのでしょうか？

ご指導よろしくお願いします．

該当のソースコード

#include <iostream>
#include <string> 
#include <stdint.h>
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <omp.h>
#include <numeric>
#include <cmath>
#include <iostream>
#include <math.h>
#include <iomanip>
#include <fstream>
#include <sstream>
#include <algorithm> 
#include <chrono>
#include <utility> 
#include <time.h>
#include <tuple> 
#include <cstdint> 
#include <cstdio>
#include <random>

int N, STEP;
using namespace std;

int main(int argc, char* argv[]) {
	cout << "N=?\n";
	cin >> N;
	cout << "STEP=?\n";
	cin >> STEP;
	
double p[N];
#pragma acc kernels
#pragma acc loop independent
for (long long i = 0; i < N; ++i)  {

	for (long long T = 1; T <= STEP; ++T) {
		p[i]+=1;
		
	}
	}

cout<<"p[5]="<<p[5]<<endl;
}

発生している問題・エラーメッセージ

$ ./a.out
N=?
1000
STEP=?
1000
p[5]=1000　　←OK
$ ./a.out
N=?
10000
STEP=?
20000000　　←OK
p[5]=2e+07
$ ./a.out
N=?
10000
STEP=?
2900000000
p[5]=2.14748e+09　　←？？
$ ./a.out
N=?
100
STEP=?
29000000000
p[5]=2.14748e+09　　←？？
$ ./a.out
N=?
100
STEP=?
2900000000
p[5]=2.14748e+09　　←？？
$ ./a.out
N=?
100
STEP=?
290000000
p[5]=2.9e+08　　←OK

(cout << fixed << setprecision(20);とした場合)
./a.out
N=?
5031
STEP=?
9898564565
p[5]=2147483647.00000000000000000000　　←？？

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

nvc++ 20.9-0
g++ (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0
tesla V100
Xeon

Answer 1

STEPの型がintなので、どれだけ大きい値を入力してもINT_MAXに丸められているのではないでしょうか

std::cout << N << ", " << STEP << std::endl;

などを挟んでみると分かると思います。

あと、

double p[N];
p[i]+=1;

これではpが未初期化なので、未定義動作を踏んでいると思います。とりあえずこちらが原因ではないようですが、未定義動作なら何が起きても不思議ではないです。

Comments

[Kinsho]

言われてみればまさに仰る通りで．．．
int型上限だなぁとまで気づいていながらこんなに情けないケアレスミスで長々と詰まっていたとは誠にお恥ずかしい限りです．
ありがとうございました．