Fixed susceptibility calculation to use magnetisation instead of spin and to work with multiple materials
Bugfix: Prevented nan in susceptibility calculation
Merged changes from upstream