Topology optimization method utilizing iterative solvers with subspace recycling applied to high-resolution electrical impedance tomography