无法将字符串读入矩阵



在下面的代码中:

  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],我会得到一个段错误。

float变量使用 %f

double变量使用 %lf

相关内容

  • 没有找到相关文章

最新更新