在堆上存储 3D 数组作为结构成员

Question

我最近开始使用 C++ 进行数值计算，我想在模拟过程中使用 Struct Operators 来存储 3D 场。 我用

在堆上创建3D数组

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

因为根据我的测试，这种方法比使用 Boost_multiarr、Eigen Tensor 或嵌套向量更快。 到目前为止，这适用于我的简约 Test.cpp，但是当我尝试将这些相同的 3D 数组实现为 Struct Operators 的成员时，我无法再使用 auto 命令：

const unsigned int RES = 256;

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

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

        // ...Initializing wfc using Settings

};

在这种情况下，我找不到声明 wfc 的方法，这样我就不会收到该类型的错误

错误：无法在赋值中将 'std::complex (*)[256][256]' 转换为 'std::complex***'

所以我的问题是如何正确声明 3D 数组wfc 以及维护这种结构方法是否完全可能/有用。如果wfc 不是结构的成员，访问wfc[i][j][k] 通常会更快吗？ （我将不得不这样做〜10 ^ 6次）

提前致谢！

Answer 1

错误消息返回正确的声明，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(Settings& 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;

    };
}

Answer 2

在您的错误消息中，编译器告诉您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 */;
        }
    }
}