Debugging :-(

Author: kbarros

examples/03_LSWT_SU3_FeI2.jl

@@ -197,8 +197,13 @@ suggest_magnetic_supercell([[0, -1/4, 1/4]])
 # [`energy_per_site`](@ref) is unchanged by the reshaping.
 
 sys_min = reshape_supercell(sys, [1 0 0; 0 1 -2; 0 1 2])
+@show energy_per_site(sys)
+@show energy_per_site(sys_min)
 randomize_spins!(sys_min)
 minimize_energy!(sys_min)
+@show energy_per_site(sys_min)
+@assert energy_per_site(sys_min) ≈ -2.9189311821844806
+@assert energy_per_site(sys) ≈ -2.9189311821844806
 @assert energy_per_site(sys_min) ≈ energy_per_site(sys)
 
 # Plot the system again, now including "ghost" spins out to 12Å. The correct