我有一个整数二维矩阵NUMI和一个称为PROB
的三维双矩阵。
这里的2分配:
int main ( int argc, char* argv[]){
double ***PROB;
int **NUMI;
NUMI = (int **)malloc((dim)*sizeof(int *));
PROB = (double ***)malloc((dim)*sizeof(double**));
...
for( i = 0; i < n ; i++){
PROB[ACT][ACTSTART][i] = value;
NUMI[i][i]= value2;
}
}
有多少行和列具有
NUMI
?是dim x dim matrix
吗???PROB
是一个三维矩阵。。。这里的分配是dim x dim x dim
?
您应该为每个维度分配内存:
NUMI = (int **)malloc((dim)*sizeof(int *));
for (i = 0; i < dim; i++)
NUMI[i] = (int*)malloc(dim * sizeof(int));
PROB = (double ***)malloc((dim)*sizeof(double**));
for (i = 0; i < dim; i++)
{
int j;
PROB[i] = (double**)malloc(dim* sizeof(double*));
for (j = 0; j < dim; j++)
{
PROB[i][j] = (double*)malloc(dim * sizeof(double));
}
}
NUMI没有行和列,它是一个指向int的指针,并且恰好指向一个分配的内存,该内存有空间容纳指向int的dim
指针,而不是指向dim * dim
int的指针。这相当于将其视为已声明int* NUMI[dim]
呼叫
int* NUMI;
NUMI= malloc( dim*dim*sizeof(int) );
将分配一个CCD_ 9乘CCD_。
但是,请注意,对于多维数组,例如int example[a][b]
,分配区域的大小等于int* example_ = malloc(a*b*sizeof(int))
,编译器计算出从多维索引到一维索引的转换,即example[c][d]
->example_[c*a+d]
所以当你做时
int* NUMI;
NUMI= malloc( dim*dim*sizeof(int) );
//...
NUMI[i][i]= value2;
编译器没有从多维转换为等效的一维所需的信息。
PROB是类似的,它指向一个内存区域,该内存区域为dim
指针提供了空间,指向指向双精度指针,而不是dim * dim * dim
双精度指针。要得到dim
立方矩阵,需要
double *PROB;
PROB = (double *)malloc( dim*dim*dim*sizeof(double) );
并在多维索引中遇到同样的问题。
如果dim是一个编译时常数,那么您可以直接声明多维数组,而不需要malloc
double PROB[dim][dim][dim];
int NUMI[dim][dim];
main((末尾的循环现在应该可以按预期工作了。
如果您必须使用malloc,请使用如上所述的malloc:
NUMI= (int *) malloc( dim*dim*sizeof(int) );
PROB = (double *)malloc( dim*dim*dim*sizeof(double) );
并将环路主体修改为
PROB[(ACT * dim * dim) + (ACTSTART * dim ) + i] = value;
NUMI[i + dim * i]= value2;
或者,如Alexey和Paolo所述,在多个循环中调用malloc。在我的解决方案中,每个变量有一个malloc()
调用,因此每个变量都引用一个连续的内存区域。循环中有多个malloc()
调用时,您有多个分配的内存区域,这些区域不太可能是连续的。
NUMI
"是"指向int
的指针的l+1
数组。PROB
"是"指向double
的指针的l+1
数组。
这更接近于您正在寻找的:
#include <cstdlib>
int main( int argc, char* argv[] )
{
size_t dim = 100;
int **NUMI = (int **)malloc( dim * sizeof(int*) );
for( size_t i = 0; i < dim; ++i )
NUMI[i] = (int*)malloc( dim * sizeof(int) );
double*** PROB = (double ***)malloc( dim * sizeof(double**) );
for( size_t i = 0; i < dim; ++i )
{
PROB[i] = (double**)malloc( dim * sizeof(double*) );
for( size_t j = 0; j < dim; ++j )
PROB[i][j] = (double*)malloc( dim * sizeof(double) );
}
/// ...
size_t ACT = 0, ACTSTART = 0;
for( size_t i = 0 ; i < dim; i++ )
{
PROB[ACT][ACTSTART][i] = 1;
NUMI[i][i] = 2;
}
}
int main() {
int **p;
int m=4,n=4,i;
/* Allocate memory */
p = (int *) malloc(sizeof(int *) * m); /* Row pointers */
for(i = 0; i < m; i++)
{
p[i] = (int) malloc(sizeof(int) * n); /* Rows */
}
这将动态分配大小为[4][4]的数组。类似地,也可以对三维阵列进行分配。