6楼: Originally posted by frcon at 2016-02-23 04:30:42
好的，谢谢你，我已经能把代码和势函数对应上了，还有个疑问麻烦了，就是力的公司是怎么写进去的，因为我不是学习分子动力学的，所以不是很清楚，了解的还不够

    rsq = delx*delx + dely*dely + delz*delz;
    r = sqrt(rsq);
    dr = r - r0[type];
    rk = k[type] * dr;

    // force & energy

    if (r > 0.0) fbond = -2.0*rk/r;
    else fbond = 0.0;

    if (eflag) ebond = rk*dr;

2楼: Originally posted by lsloneil at 2016-02-23 15:27:18
这很好理解啊，模拟中可以定义多种类型的键，每一种类型的键都有相应的力常数k和平衡距离r0， type就表示健的种类。

p.s. 如果不懂c++，还是补一点基础知识再去看源代码。

2楼: Originally posted by lsloneil at 2016-02-23 15:27:18
这很好理解啊，模拟中可以定义多种类型的键，每一种类型的键都有相应的力常数k和平衡距离r0， type就表示健的种类。

p.s. 如果不懂c++，还是补一点基础知识再去看源代码。