0

我正在从一个文件中扫描一个双精度数组。我可以打印这个数组来检索文件中每一行的正确值。

当我在此过程中尝试将另一个文件扫描到另一个阵列中时，就会出现我的问题。

我有：

//
//

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

int main (int argc, char **argv) {

  double c[13], d[13];
  char filename[20];
  double x1, x2, y1, y2, z1, z2, d1, d2, au, aux, auy;
  double dx, dy, dz, nnid, chk, nn;
  int n, b;
  char *in;
  in=argv[1];
  FILE* infile;
  FILE* inbfile;
  infile = fopen(in, "r");
  inbfile = fopen("copy.bt", "r");

  while (!feof(infile)) {

    fscanf(infile, "%f %f %e %e %e %e %e %e %e %e %e %e %e %f",c,c+1,c+2,c+3,c+4,c+5,c+6,c+7,c+8,c+9,c+10,c+11,c+12,c+13);
    printf("Selected Particle A: %f\n",c[0]);
    n=0;

    while (!feof(infile)) { 
      printf("%f\n",c[0]);
      fscanf(infile, "%f %f %e %e %e %e %e %e %e %e %e %e %e %f",d,d+1,d+2,d+3,d+4,d+5,d+6,d+7,d+8,d+9,d+10,d+11,d+12,d+13);
      printf("%f\n",c[0]);
      printf("Selected Particle B: %f\n",d[0]);

      /**/
      printf("%f = %f ?\n",c[0],d[0]);
      if (c[0]==d[0]) {
        printf("Same Particle SKIP...\n");
      }
      else {  

        dx = (d[4])-(c[4]);
        dy = (d[5])-(c[5]);
        dz = (d[6])-(c[6]);

        printf("dx dy dz %e %e %e\n",d[4],d[5],d[6]);

        /**/
        if (n == 0) {  
          au=pow(((dx*dx*dx)+(dy*dy*dy)+(dz*dz*dz)),(1.0/3.0));
          printf("%f is %e from %f\n",c[0],au,d[0]);
        }
        else { 
          aux=pow(((dx*dx)+(dy*dy)+(dz*dz)),(1.0/3.0));
          printf("%f is %e from %f\n",c[0],aux,d[0]);
          if (aux < au) {
            au = aux;
            nnid = d[0];
          }
        }
        /**/

      }
      n++;
      nn=d[1];
      /**/

    }

    printf("%f Is Particle %f At %e\n", c[1], nnid, au);

  }

  fclose(infile);
  fclose(inbfile);
  return 0;

}

那么为什么值会改变呢？我所做的只是将另一个文件的第一行扫描到一个单独的数组中。

4

2 回答 2

3

您试图将 14 个值读入具有 13 个元素的数组中，从而导致未定义的行为。将数组声明为float c[14], d[14]。

当您声明为时d，float d[13]指向d+12数组中的最后一个元素，并且d+13超过了数组的末尾，因此不允许读取它。

于 2012-12-02T16:50:29.280 回答
1

添加到 Interjay 的答案：您正在另一个地方调用未定义的行为。您正在尝试double使用格式说明符扫描 s %f，而正确的格式说明符double是%lf.

编辑：看来你给我们的代码错误，所以上面的陈述不太正确。（但是，我会保留此作为未来访问者的参考。）

于 2012-12-02T16:53:58.430 回答