So far, the long form of the quadratic bending energy has been
implemented and tested, and it appears to do the right thing.

This is a snapshot commit.  The code remains out-commented until
it is finished.

[[Imported from SVN: r8403]]