Calculation of spring constraint energy


Click here for full thread
Clicked A Few Times
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.