Implementing periodic boundaries for rigid bodies

Lots of the changes here are entirely cosmetic or else were intended entirely for my use - e.g. most of the new comments and the changes to pymolwrapper.py which were designed to allow visualisation of small rigid-body periodic systems.

The important changes are the ones where I've defined new or overloaded functions and classes. In particular the new _otp_bulk module, the new classes in minpermdist_stochastic.py and periodic_exact_match.py, and the changes to the cpp files/their wrappers. I've also put in a fairly extensive test routine.
This only works for equally weighted rigid bodies.  It is treating all rigid bodies as if the have the same mass.  We'll have to think about what we actually want to do here.  Do we want the center of mass or the center of geometry?