нулевой выход при использовании OpenACC
Я использую PGI Community Edition 17.10 для компиляции и запуска следующего кода. почему вывод неправильный, когда я добавляю директивы OpenACC?
можете ли вы помочь мне, почему это произошло?
Заранее спасибо,
Саяд
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <assert.h>
#include <openacc.h>
#include<time.h>
#include <string.h>
#include <malloc.h>
#include <cuda_runtime_api.h>
#define NX 2
#define NY 2
#define NZ 2
int main(void)
{
static int i, j, k;
static double A[NX][NY][NZ]=2 ,B[NX][NY][NZ]=10.,C[NX]=10.,D[NY]=10.,E[NZ]=10.;
FILE *file;
file = fopen("BB-and-A.csv", "w");
#pragma acc data copy( A ,B,C,D,E,i, j, k)
{
#pragma acc kernels loop private(i, j, k)
for (i = 0; i <= NX; i++) {
for (j =0; j <= NY ; j++) {
for (k =0; k <= NZ ; k++) {
C[i]=i;
D[j]=j;
E[k]=k;
}
}
}
}
for (i = 0; i <= NX; i++) {
for (j =0; j <= NY ; j++) {
for (k =0; k <= NZ ; k++) {
fprintf(file, "%e, %e, %e \n", C[i], D[j],E[k] );
}
}
}
fclose(file);
}
Решение
У вас есть ряд проблем с этим кодом.
1) Ваши границы массива неверны. Поскольку петли идут от 1 до <= N, но массивы имеют только N членов, вы записываете конец массива.
2) Ваш цикл не распараллеливается, так как вы пишете в каждый элемент из нескольких итераций цикла. Чтобы исправить, я бы сделал эти три отдельных цикла.
3) Переменные индекса цикла не должны быть статическими. Это помещает их в глобальное хранилище и, таким образом, вызывает зависимость. Хотя вы можете исправить это, поместив их в закрытое предложение, лучше убрать статические и позволить компилятору неявно приватизировать их.
4) Нет необходимости копировать переменные индекса цикла.
Попробуйте что-то вроде следующего:
% cat test2.c
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <assert.h>
#include <openacc.h>
#include<time.h>
#include <string.h>
#include <malloc.h>
// #include <cuda_runtime_api.h>
#define NX 2
#define NY 2
#define NZ 2int main(void)
{
int i, j, k;
static double A[NX+1][NY+1][NZ+1]=2 ,B[NX+1][NY+1][NZ+1]=10.,C[NX+1]=10.,D[NY+1]=10.,E[NZ+1]=10.;
FILE *file;
file = fopen("BB-and-A.csv", "w");
#pragma acc data copy(A,B,C,D,E)
{
#pragma acc kernels
{
for (i = 0; i <= NX; i++) C[i]=i;
for (j =0; j <= NY ; j++) D[j]=j;
for (k =0; k <= NZ ; k++) E[k]=k;
} }
for (i = 0; i <= NX; i++) {
for (j =0; j <= NY ; j++) {
for (k =0; k <= NZ ; k++) {
fprintf(file,"%e, %e, %e \n", C[i], D[j],E[k] );
} } }
fclose(file);
}
% pgcc test2.c -ta=tesla:cc60 -Minfo=accel
main:
23, Generating copy(A[:][:][:],B[:][:][:],C[:],E[:],D[:])
27, Loop is parallelizable
Accelerator kernel generated
Generating Tesla code
27, #pragma acc loop gang, vector(32) /* blockIdx.x threadIdx.x */
28, Loop is parallelizable
Accelerator kernel generated
Generating Tesla code
28, #pragma acc loop gang, vector(32) /* blockIdx.x threadIdx.x */
29, Loop is parallelizable
Accelerator kernel generated
Generating Tesla code
29, #pragma acc loop gang, vector(32) /* blockIdx.x threadIdx.x */
% a.out
% cat BB-and-A.csv
0.000000e+00, 0.000000e+00, 0.000000e+00
0.000000e+00, 0.000000e+00, 1.000000e+00
0.000000e+00, 0.000000e+00, 2.000000e+00
0.000000e+00, 1.000000e+00, 0.000000e+00
0.000000e+00, 1.000000e+00, 1.000000e+00
0.000000e+00, 1.000000e+00, 2.000000e+00
0.000000e+00, 2.000000e+00, 0.000000e+00
0.000000e+00, 2.000000e+00, 1.000000e+00
0.000000e+00, 2.000000e+00, 2.000000e+00
1.000000e+00, 0.000000e+00, 0.000000e+00
1.000000e+00, 0.000000e+00, 1.000000e+00
1.000000e+00, 0.000000e+00, 2.000000e+00
1.000000e+00, 1.000000e+00, 0.000000e+00
1.000000e+00, 1.000000e+00, 1.000000e+00
1.000000e+00, 1.000000e+00, 2.000000e+00
1.000000e+00, 2.000000e+00, 0.000000e+00
1.000000e+00, 2.000000e+00, 1.000000e+00
1.000000e+00, 2.000000e+00, 2.000000e+00
2.000000e+00, 0.000000e+00, 0.000000e+00
2.000000e+00, 0.000000e+00, 1.000000e+00
2.000000e+00, 0.000000e+00, 2.000000e+00
2.000000e+00, 1.000000e+00, 0.000000e+00
2.000000e+00, 1.000000e+00, 1.000000e+00
2.000000e+00, 1.000000e+00, 2.000000e+00
2.000000e+00, 2.000000e+00, 0.000000e+00
2.000000e+00, 2.000000e+00, 1.000000e+00
2.000000e+00, 2.000000e+00, 2.000000e+00
Другие решения
Других решений пока нет …