solver: fill nans to be sure

usadelndsoc/solver.py

@@ -388,6 +388,7 @@ class Solver:
 
         if np.isnan(Phi).any():
             _log_solve.error("Failed to solve")
+            Phi[...] = np.nan
             M = None
 
         self._M = M