0

我在 MPI C 中有一个用于图像处理（pgm 文件）的 MPI 程序，我对 2D 数组使用动态分配，如下所示。

float **masterbuf;
masterbuf = arralloc(sizeof(float), 2, M, N);

当我使用

float masterbuf[M][N];

该程序提供的图像看起来不错。

问题是当我使用动态分配时，图像在其左侧丢失了一些像素。所以这些缺失的像素会形成一条黑线。就像图像向右移动了 2 个像素。我没有对图像做任何其他操作，只是读取它并再次打印它。

我用来写图像的函数是：

void pgmwrite(char *filename, void *vx, int nx, int ny)
{
  FILE *fp;

  int i, j, k, grey;

  float xmin, xmax, tmp, fval;
  float thresh = 255.0;

  float *x = (float *) vx;

  if (NULL == (fp = fopen(filename,"w")))
  {
    fprintf(stderr, "pgmwrite: cannot create <%s>\n", filename);
    exit(-1);
  }

  printf("Writing %d x %d picture into file: %s\n", nx, ny, filename);

  /*
   *  Find the max and min absolute values of the array
   */

  xmin = fabs(x[0]);
  xmax = fabs(x[0]);

  for (i=0; i < nx*ny; i++)
  {
    if (fabs(x[i]) < xmin) xmin = fabs(x[i]);
    if (fabs(x[i]) > xmax) xmax = fabs(x[i]);
  }

  if (xmin == xmax) xmin = xmax-1.0;

  fprintf(fp, "P2\n");
  fprintf(fp, "# Written by pgmwrite\n");
  fprintf(fp, "%d %d\n", nx, ny);
  fprintf(fp, "%d\n", (int) thresh);

  k = 0;

  for (j=ny-1; j >=0 ; j--)
  {
    for (i=0; i < nx; i++)
    {
      /*
       *  Access the value of x[i][j]
       */

      tmp = x[j+ny*i];

      /*
       *  Scale the value appropriately so it lies between 0 and thresh
       */

      fval = thresh*((fabs(tmp)-xmin)/(xmax-xmin))+0.5;
      grey = (int) fval;

      fprintf(fp, "%3d ", grey);

      if (0 == (k+1)%16) fprintf(fp, "\n");

      k++;
    }
  }

  if (0 != k%16) fprintf(fp, "\n");
  fclose(fp);
}
4

1 回答 1

2

您对 masterbuf 的两个定义可能都创建了 2D 数组，但它们的创建方式不同。函数arralloc()为数据和指针创建空间——不仅仅是像简单的静态数组定义那样的数据。这意味着在 pgmwrite() 中，无论使用何种方法，x[i][j] 都会返回相同的结果，但由于指针的参与，x[i] 将意味着两个不同的东西。

值得注意的是，如果您将void *vx原型更改为float *vx. 由于您立即无条件地将这个 void * 转换为 float *，因此无论如何这样做会更好。

(2nd edit:) Also, if interested, check out this response. It shows how to index using two dimensions into a single malloc'd block, without arralloc().

于 2012-11-27T10:28:51.363 回答