数组下标不是整数答案

【问题标题】：Array subscript is not an integer数组下标不是整数
【发布时间】：2011-06-05 19:53:52
【问题描述】：

按照上一个问题Malloc Memory Corruption in C，现在我有另一个问题。我有相同的代码。现在我正在尝试将数组 A * vc 中包含的值相乘并存储在 res 中。然后将 A 设置为零，我与 res 和 vc 进行第二次乘法，并将值存储在 A 中。（A 和 Q 是方阵，mc 和 vc 是 N 行两列矩阵或数组）。这是我的代码：

    int jacobi_gpu(double A[], double Q[], 
           double tol, long int dim){
  int nrot, p, q, k, tid;
  double c, s;
  double *mc, *vc, *res;
  int i,kc;
  double vc1, vc2;

  mc   = (double *)malloc(2 * dim * sizeof(double));
  vc   = (double *)malloc(2 * dim * sizeof(double));
  vc   = (double *)malloc(dim * dim * sizeof(double));

  if( mc == NULL || vc == NULL){
    fprintf(stderr, "pb allocation matricre\n");
    exit(1);
  }

  nrot = 0;

  for(k = 0; k < dim - 1; k++){

    eye(mc, dim);
    eye(vc, dim);

    for(tid = 0; tid < floor(dim /2); tid++){
      p = (tid + k)%(dim - 1);
      if(tid != 0)
    q = (dim - tid + k - 1)%(dim - 1);
      else
    q = dim - 1;

      printf("p = %d | q = %d\n", p, q);
      if(fabs(A[p + q*dim]) > tol){

    nrot++;
    symschur2(A, dim, p, q, &c, &s);


    mc[2*tid] = p;                                               vc[2 * tid] = c;
    mc[2*tid + 1] = q;                                           vc[2*tid + 1] = -s;
    mc[2*tid + 2*(dim - 2*tid) - 2] = p;                         vc[2*tid + 2*(dim - 2*tid)   - 2 ] = s;
    mc[2*tid + 2*(dim - 2*tid) - 1] = q;                         vc[2 * tid + 2*(dim - 2*tid) - 1 ] = c;



    }
    }

    for( i = 0; i< dim; i++){
      for(kc=0; kc < dim; kc++){
    if( kc < floor(dim/2)) {
      vc1 = vc[2*kc + i*dim];
      vc2 = vc[2*kc + 2*(dim - 2*kc) - 2];
    }else {
      vc1 = vc[2*kc+1 + i*dim];
      vc2 = vc[2*kc - 2*(dim - 2*kc) - 1];
    }
    res[kc + i*dim] = A[mc[2*kc] + i*dim]*vc1 + A[mc[2*kc + 1] + i*dim]*vc2;
      }
    }

    zero(A, dim);

    for( i = 0; i< dim; i++){
      for(kc=0; kc < dim; k++){
    if( k < floor(dim/2)){
      vc1 = vc[2*kc + i*dim];
      vc2 = vc[2*kc + 2*(dim - 2*kc) - 2];
    }else {
      vc1 = vc[2*kc+1 + i*dim];
      vc2 = vc[2*kc - 2*(dim - 2*kc) - 1];
    }
    A[kc + i*dim] = res[mc[2*kc] + i*dim]*vc1 + res[mc[2*kc + 1] + i*dim]*vc2;
      }
    }


    affiche(mc,dim,2,"Matrice creuse");
    affiche(vc,dim,2,"Valeur creuse");

  }

  free(mc);
  free(vc);
  free(res);
  return nrot;
}

当我尝试编译时，我有这个错误：

    jacobi_gpu.c: In function ‘jacobi_gpu’:
jacobi_gpu.c:103: error: array subscript is not an integer
jacobi_gpu.c:103: error: array subscript is not an integer
jacobi_gpu.c:118: error: array subscript is not an integer
jacobi_gpu.c:118: error: array subscript is not an integer
make: *** [jacobi_gpu.o] Erreur 1

相应的行是我将结果存储在 res 和 A 中的位置：

res[kc + i*dim] = A[mc[2*kc] + i*dim]*vc1 + A[mc[2*kc + 1] + i*dim]*vc2;

和

A[kc + i*dim] = res[mc[2*kc] + i*dim]*vc1 + res[mc[2*kc + 1] + i*dim]*vc2;

谁能解释一下这个错误是什么，我该如何纠正它？谢谢你的帮助。 ;)

【问题讨论】：

标签： c arrays cuda gpu

【解决方案1】：

mc 是双精度类型。必须是整型

【讨论】：

【解决方案2】：

编译器正在抱怨，因为您用作数组索引的表达式的计算结果为类型 double。

换句话说，表达式：

mc[2*kc] + i*dim

...会给你一个double类型的结果。如果您不明白为什么此表达式的计算结果为 double，您可能需要查看 C 中 usual arithmetic type conversions 的规则。

问题是数组索引必须是整数类型，例如int 或long。这是因为 C 中的数组下标运算符基本上是指针运算的简写。换句话说，说array[N] 与说*(array + N) 是一样的。但是你不能用像float或double这样的非整数类型来做指针运算，所以当然数组下标运算符也不会那样工作。

要解决此问题，您需要将数组索引表达式的结果转换为整数类型。

【讨论】：

【解决方案3】：

mc 是指向双精度的指针。

A[mc[2*kc + 1]

在上面，您使用 mc（双数组）中的值对 A 进行索引。而且，还有其他类似的案例。如果您确定这些值，请转换为 int

【讨论】：

【解决方案4】：

你的mc声明：

mc   = (double *)malloc(2 * dim * sizeof(double));

然后您在数组访问中多次使用mc。例如：

A[mc[2*kc + 1] ...]

您能否将mc 更改为int 数组而不是double？

【讨论】：

【解决方案5】：

看起来您正在使用 mc 中的条目，它们是双精度值，作为数组下标的一部分，从而使整个下标成为双精度值。

如果您打算这样做，请尝试转换回整数。我不知道这个问题的背景是什么，但我会好好看看你在做什么，以确保你真的想使用 mc 的内容作为下标。

【讨论】：

【解决方案6】：

mc 是一个双精度数组，浮点值不能用于索引数组。我注意到，在您的代码中，除了整数之外，您没有将任何其他内容分配给 mc。您应该考虑将 mc 的类型更改为整数数组。

【讨论】：