Блог

Недавно я начал работать с C для численных вычислений, где я хочу использовать a Struct Operators для хранения 3D-полей в ходе моделирования. Я создаю 3D-массивы в куче с помощью

 const unsigned int RES = 256;
auto arr3D = new double [RES][RES][RES];
 

потому что, судя по тому, что я протестировал, этот подход быстрее, чем использование Boost_multiarr, собственного тензора или вложенных векторов.
До сих пор это хорошо работало с моим минималистичным Test.cpp но когда я пытаюсь реализовать эти же 3D-массивы в качестве членов my Struct Operators , я больше не могу использовать auto команду:

 const unsigned int RES = 256;

struct Operators {
public:
    std::complex<double>*** wfc;         // edited, see 'spitconsumers' comment

    Operators(Settings amp;set) {        // just another structure used by Operators
        wfc = new std::complex<double> [RES][RES][RES];

        // ...Initializing wfc using Settings

};
 

В этом случае я не нашел способа объявить wfc , чтобы я не получал ошибок типа

ошибка: не удается преобразовать «std::сложный (*)[256][256]’ в «std::комплекс***» в задании

Поэтому мой вопрос заключается в том, как правильно объявить 3D-массив wfc и возможно ли вообще/полезно поддерживать такой подход к структуре. Было бы вообще быстрее получить доступ к wfc[i][j][k] if wfc , если бы он не был членом структуры? (Мне придется сделать это ~10^6 раз)

Заранее спасибо!

Сообщение об ошибке возвращает правильное объявление, std::complex<double>(*wfc)[RES][RES]; .

 const unsigned int RES = 256;
struct Settings {};
struct Operators {
public:
    std::complex<double>(*wfc)[RES][RES];         // edited, see 'spitconsumers' comment

    Operators(Settingsamp; set) {        // just another structure used by Operators
        wfc = new std::complex<double>[RES][RES][RES];

        // ...Initializing wfc using Settings
        // Setting the last element
        wfc[254][254[254] = 42;

    };
}
 

В вашем сообщении об ошибке компилятор сообщает вам, что new complex [RES][RES][RES] возвращает тип std::complex (*)[256][256] .

double wfc должно быть std::complex*** wfc .

Типы членов класса должны быть известны, и, поскольку вы не инициализируете переменную-член во время ее определения, компилятор не может определить тип переменной, поэтому auto она не работает как тип члена класса, когда вы не назначаете ей что-то немедленно. Но вы присваиваете значение непосредственно переменной в своей тестовой программе, позволяя компилятору вывести тип там.

Кроме того, new оператор возвращает указатель на массив неинициализированных std::complex<double>** указателей. Вам нужно вызвать new все неинициализированные указатели.

 std::complex<double>*** wfc;

wfc = new std::complex<double>** [RES];

for(int i = 0; i < RES; i  ) {
    wfc[i] = new std::complex<double>* [RES];

    for(int j = 0; j < RES; j  ) {
        wfc[i][j] = new std::complex<double> [RES];

        for(int k = 0; k < RES; k  ) {
            wfc[i][j][k] = /* insert default value here */;
        }
    }
}