In this paper, we develop a shape optimization-based algorithm for the
electrical impedance tomography (EIT) problem of determining a piecewise
constant conductivity on a polygonal partition from boundary measurements. The
key tool is to use a distributed shape derivative of a suitable cost functional
with respect to movements of the partition. Numerical simulations showing the
robustness and accuracy of the method are presented for simulated test cases in
two dimensions