Calculation of spring constraint energy
Click here for full thread
Clicked A Few Times
10:31:11 AM PDT - Tue, Apr 15th 2014
I don't see any particular errors here. The term k*(r-r0)**2 is added to overall energy and forces are calculated from it consistently. It is true, however, that there is a misprint in documentation stating that 0.5*k*(r-r0)**2 term is added.