在下面的代码中:

  atnames = new char[natoms];
  xyzs = new double*[natoms];
  for(int iat=0;iat<natoms;iat++)
  {
    int idum1;
    xyzs[iat] = new double[3];
    for(int ii = 0; ii<3; ii++)
    {
      xyzs[iat][ii] = 0.0;
    }
    getline(finp, tline);
    sscanf(tline.c_str(), "%c(%i) %f %f %f ATOM", &atnames[iat], &idum1, &xyzs[iat][0], &xyzs[iat][1], &xyzs[iat][2]);
  }

  for(int iat=0; iat<natoms; iat++)
  {
    cout << atnames[iat] << iat+1 << ":" << xyzs[iat][0] <<  " " << xyzs[iat][1] << " " << xyzs[iat][2] << endl;
  }


我正在尝试阅读以下内容:

O(1)       1.23799   0.00000   0.00000 ATOM
N(2)      -0.75911   1.08623   0.00000 ATOM
C(3)       0.00000  -0.00000   0.00000 ATOM
H(4)      -1.75297   0.97208   0.03384 ATOM
H(5)      -0.30916   1.98585  -0.01442 ATOM
H(6)      -0.52646  -0.95018   0.03999 ATOM


但是,当我打印出这些值时,我得到:

O1:5.27341e-315 0 0
N2:1.58536e-314 5.26712e-315 0
C3:0 0 0
H4:1.59047e-314 5.26123e-315 5.05975e-315
H5:1.58004e-314 5.3044e-315 1.56185e-314
H6:1.58343e-314 1.58694e-314 5.06791e-315


因此,显然我正确地读取了第一个char,但是没有将浮点数的值正确地存储到数组xyz中。我感觉这是一个引用/引用问题,但是我不清楚我在做什么错。

如果尝试将值直接读取到xyzs[iat][0]而不是&xyzs[iat][0],则会出现段错误。

最佳答案

%f用于float变量。

%lf用于double变量。

08-28 04:31