Bugfix: Exchange energy calculation in C++ missed a factor of Ms. This
did not affect any calculations because E_exch is calculated using M.dotSum(H) in exchange_field.py...