вложенные классы с указателями с openACC

У меня довольно большой код на C ++, мне пришлось интегрировать какой-то новый класс в базовый класс, как показано ниже.

class A
    {
    int N;
    B  b;
    double *__restrict__ w;
    construct();
    }



  A::construct()
    {
    w=new double[N];

    #pragma acc data enter create(this)
    #pragma acc update device(this)
    #pragma acc data enter create(w)
    // allocate  class A
    b.construct()
    }



  class B
    {
    double *__restrict__ u;
    double *__restrict__ v;
    B(){};
    construct();
    }
    B::construct()
    {
    u=new double[N];
    v=new double[N];
    #pragma acc data enter create(this)
    #pragma acc update device(this)
    #pragma acc data enter create(u)
    #pragma acc data enter create(v)

    }

Я думаю, что столкнулся с проблемой глубокого копирования, поскольку указатели класса B недействительны и, следовательно, поведение кода на GPU i undefined. Я был бы признателен за отзывы о том, как выполнить включение класса в другой класс, не вдаваясь в проблему глубокого копирования. Я подозреваю, что устройство обновления (это) каким-то образом вызывает это.


gpu openacc pgi pgi-accelerator
person JimBamFeng    schedule 27.12.2018    source источник

Ответы (1)


arrow_upward
1
arrow_downward

У вас есть полный пример, который воссоздает ошибку, которую вы видите? Я написал небольшой тестовый пример, используя ваш фрагмент кода, и он отлично сработал. (Увидеть ниже)

Если бы вы обновляли указатель «this» после создания массивов, это было бы проблемой, поскольку вы бы перезаписывали указатели устройств указателями узлов. Но, как вы показали выше, это не должно быть проблемой.

% cat test.cpp
#include <iostream>

class B
{
        public:
                int N;
                double *__restrict__ u;
                double *__restrict__ v;
                void construct(int);
};

void B::construct(int _N)
{
        N=_N;
        u=new double[N];
        v=new double[N];
#pragma acc enter data create(this)
#pragma acc update device(this)
#pragma acc enter data create(u[:N])
#pragma acc enter data create(v[:N])
}

class A
{
        public:
                int N;
                B  b;
                double *__restrict__ w;
                void construct(int);
};

void A::construct(int _N)
{
        N=_N;
        w=new double[N];
#pragma acc enter data create(this)
#pragma acc update device(this)
#pragma acc enter data create(w[:N])

        // allocate  class A
        b.construct(N);
}


int main() {

        A myA;
        int N=32;
        myA.construct(N);

#pragma acc parallel loop present(myA)
        for (int i=0; i<N; ++i) {
                myA.w[i] = i;
                myA.b.u[i] = i;
                myA.b.v[i] = i;
        }
#pragma acc update host( myA.w[:N], myA.b.u[:N], myA.b.v[:N])
        for (int i=0; i<N; ++i) {
                std::cout << myA.w[i] << ":" << myA.b.u[i] << ":" << myA.b.v[i] << std::endl;
        }
        return 0;
}
% pgc++ test.cpp -Minfo=accel -V18.10 -ta=tesla; a.out
main:
     49, Generating present(myA)
         Accelerator kernel generated
         Generating Tesla code
         52, #pragma acc loop gang, vector(32) /* blockIdx.x threadIdx.x */
     56, Generating update self(myA.b.u[:N],myA.w[:N],myA.b.v[:N])
B::construct(int):
     21, Generating update device(this[:1])
         Generating enter data create(this[:1],v[:N],u[:N])
A::construct(int):
     41, Generating update device(this[:1])
         Generating enter data create(w[:N],this[:1])
0:0:0
1:1:1
2:2:2
3:3:3
4:4:4
5:5:5
6:6:6
7:7:7
8:8:8
9:9:9
10:10:10
11:11:11
12:12:12
13:13:13
14:14:14
15:15:15
16:16:16
17:17:17
18:18:18
19:19:19
20:20:20
21:21:21
22:22:22
23:23:23
24:24:24
25:25:25
26:26:26
27:27:27
28:28:28
29:29:29
30:30:30
31:31:31
person Mat Colgrove    schedule 27.12.2018
comment
Конечно. Отправьте его в службу поддержки PGI (trs @ pgroup, com) и попросите Алекса переслать его мне. - person Mat Colgrove; 28.12.2018
comment
Большое спасибо, поскольку я боролся с тем, почему возникает эта ситуация, поскольку оба класса разрабатываются отдельно и работают нормально. когда я их соединяю, возникает проблема. Я только что отправил это. Еще раз спасибо за ваше время и помощь - person JimBamFeng; 28.12.2018