Numerical determination of anomalies in multifrequency electrical impedance tomography

The multifrequency electrical impedance tomography consists in retrieving the conductivity distribution of a sample by injecting a finite number of currents with multiple frequencies. In this paper we consider the case where the conductivity distribution is piecewise constant, takes a constant value outside a single smooth anomaly, and a frequency dependent function inside the anomaly itself. Using an original spectral decomposition of the solution of the forward conductivity problem in terms of Poincar\'e variational eigenelements, we retrieve the Cauchy data corresponding to the extreme case of a perfect conductor, and the conductivity profile. We then reconstruct the anomaly from the Cauchy data. The numerical experiments are conducted using gradient descent optimization algorithms.