Refactor initialization to initialize_fe()
Refactor the matrix A variables, document better
Precalculate only once
Reuse the electronic density
Print the FFT and FE energies in the same way
Disable the hard checks for imaginary parts
Do 10 MD steps