sPuReMD: fix issue with torsion & conjugation forces being incorrect due to rel_box arithmetic being wrong (incorrect directions for position vectors). Redesign atomic distance and position update functions to support both periodic boundary conditions and non-periodic boundary conditions.  Behavior change: let random_velocity parameter in control file take precedent over specifying a non-zero initial temperature value. Other code clean-up.