solver: fix solve_many with terminals

usadelndsoc/solver.py

@@ -561,7 +561,6 @@ class Solver:
         self._t_solve_min = np.inf
         self._M = None
         self._Phi[...] = 0
-        self._core.fix_terminals(Phi)
 
         for j in js:
             self.solve(omega=omega[j], **kw)