在下面的代码中:
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
变量。