用c分配整数矩阵



我有一个整数二维矩阵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]的数组。类似地,也可以对三维阵列进行分配。

相关内容

  • 没有找到相关文章

最新更新