Added implementation of second order uniaxial anisotropy
Changed default element to Fe to avoid errors in rasmol/povray output.
Bugfix: Added global reduction of max_torque
Bugfix: Corrected applied field resolution in static hysteresis program.
Bugfix: Corrected field units in static hysteresis loop
Initial implementation of Lagrange multiplier constrained minimisation
Implemented ellipsoid particles shapes
Bugfix: fixed parallel compilation error in statistics calculation
Implemented periodic boundary conditions for parallel simulation

Modified generation of halo atoms to include replication of system generating periodic boundaries. Initialisation of total number of super cells has been necessarily moved to master create function. Neighbourlist calculation has been modified to accept negative supercell indices. Tidied up tabs.
Implemented macrocell data output including position, magnetisation and demagnetisation field data.
Implemented centre of mass demag calculation 

Improved macrocell volume approximation. Added MPI timer for accurate timing for parallel version.
Bugfix: Fixes in calculation of macrocells

Error in calculation of centre of mass corrected and shift added for fence post problem when assigning atoms to macrocells. Also added reduction operation for correct output of macrocell demagnetisation fields.
Implemented temperature dependent lattice anisotropy
Changed definition of fourth order anisotropy terms
Changed default values for lattice anisotropy to prevent NaN
Corrected sign error in energy calculations
Implemented coupling of materials to the phonon temperature
Merge branch 'develop' into vio-refactor