Solvers now consistent and yield same result as prototype

TODO:

• check why there is a slight difference in the iterates between de_geus(...) and muSpectre_proto
• write tests for the field_maps' mean() method