Question

嗨，我已经遇到过一段代码，我很难理解它为什么会起作用，如果有人能解释我，我会很感激这是代码这实际上是真实代码的一个非常简化的版本，可以在这里找到http://acts.nersc.gov/hypre/example10/ex10.cxx.html

#include <iostream>
int main(int argc, char *argv[]){
double ***elemStiff = new double**[10];
int i, k;

for (i = 0; i < 10; i++){
    elemStiff[i] = new double*[4];
    for (k = 0; k < 4; k++){
        elemStiff[i][k] = new double[4];
    }

    double **A = elemStiff[i];
    for (k = 0; k < 4; k++)
        A[k][k] = 2/3.;

    A[0][1] = A[1][0] = -1/6.;
    A[0][2] = A[2][0] = -1/3.;
    A[0][3] = A[3][0] = -1/6.;
    A[1][2] = A[2][1] = -1/6.;
    A[1][3] = A[3][1] = -1/3.;
        A[2][3] = A[3][2] = -1/6.;
    }
    std::cout << elemStiff[1][0][0] << std::endl;
}

简而言之，这里我们初始化并填充elemStiff数组的值，该数组有10行，每行是4x4矩阵。

对我来说非常奇怪的是，代码在双倍的帮助下填充了elemStiff数组的条目** A = elemStiff [i];

当我第一次看到我虽然它会给出错误但是它有效！传递给A的值，它们在elemStiff数组上同时传递，但我不明白它为什么会起作用

谢谢

Answer 1

让我们从更简单的版本开始：

int* array = new int[10];
int* ptr = array;            // equivalent to ptr = &array[0]

for (int i = 0; i < 10; ++i)
    ptr[i] = i;              // writing into the memory that ptr points to

for (int i = 0; i < 10; ++i)
    std::cout << array[i];   // outputs: 0123456789

在这个例子中，我们创建一个指针ptr，指向array的第一个元素。然后，当我们写入ptr指向的内存时，我们实际上正在写入array。因此，当我们检查array元素的值时，我们可以看到它们已正确初始化。

您的代码几乎相同，只有10个二维数组的数组：

double ***elemStiff = new double**[10];

然后就像我在我的简单示例中使用ptr一样，您使用临时指针A使其指向这10个二维数组中的每一个：

for (i = 0; i < 10; i++)
    double **A = elemStiff[i];

了解指针的使用

1 个答案: