Electrical impedance tomography is a noninvasive imaging technique for recovering the admittivity distribution inside a body from boundary measurements of current and voltage. In practice, impedance tomography suffers from inaccurate modelling of the measurement setting: The exact electrode locations and the shape of the imaged object are not necessarily known precisely. In this work, we tackle the problem with imperfect electrode information by introducing the Fréchet derivative of the boundary measurement map of impedance tomography with respect to the electrode shapes and locations. This enables us to include the fine-tuning of the information on the electrode positions as a part of a Newton-type output least squares reconstruction algorithm; we demonstrate that this approach is feasible via a two-dimensional numerical example based on simulated data. The impedance tomography measurements are modelled by the complete electrode model, which is in good agreement with real-life electrode measurements.
We consider and analyze applying a spectral inverse iteration algorithm and its subspace iteration variant for computing eigenpairs of an elliptic operator with random coefficients. With these iterative algorithms the solution is sought from a finite dimensional space formed as the tensor product of the approximation space for the underlying stochastic function space, and the approximation space for the underlying spatial function space. Sparse polynomial approximation is employed to obtain the first one, while classical finite elements are employed to obtain the latter. An error analysis is presented for the asymptotic convergence of the spectral inverse iteration to the smallest eigenvalue and the associated eigenvector of the problem. A series of detailed numerical experiments supports the conclusions of this analysis.