在下面的代码中:
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
,但没有将 float 的值正确存储到数组 xyz
中。我觉得这是一个引用/取消引用问题,但我不清楚我做错了什么。
如果我尝试将值直接读取到 xyzs[iat][0]
而不是 &xyzs[iat][0]
,我会遇到段错误。
最佳答案
%f
用于float
变量。
%lf
用于double
变量。
关于c++ - 无法将字符串读入矩阵,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/21125032/