Oh, I forgot to mention, but I've been making some changes to the simulation unrelated to the testing. Well, some of them have been related…
In the tension model utilities, you can now pick x_min and x_max when plotting energy landscapes E(x) (that was for plotting Kramers’ cusp potential). You can also pick F_max when plotting F(x).
I also ripped apart the tension-group idea and put it back together. The initial implementation was a bit embarassing. Now you can do things that tension groups are intended for, like having two seperate WLC groups (for, say, folded- and unfolded-domains).
Subversion is lots of fun, by the way ;).