7楼: Originally posted by lsloneil at 2016-02-24 03:22:14
这很简单，只要你知道了你键能还有受力（也就是键能对距离的导数取负值），直接把ebond和fbond的函数形式改掉就行了。以bond_harmonic.cpp为例，在compute()这个函数里，你可以看到以下代码
    rsq = delx*delx  ...

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++，还是补一点基础知识再去看源代码。