The natural spaces for boundary voltages and currents are and (see Sobolev and Trace Spaces). They can be computed discretely by spectral interpolation between and on the boundary.
Solve the generalised eigenproblem with the Boundary Mass and Stiffness Matrices
Clip tiny negative eigenvalues caused by round-off to zero. The columns of are the discrete Laplace–Beltrami eigenfunctions (on a circle: Fourier modes, with ). For define
- : , the norm.
- : , the norm.
- : the discrete norms. On a circle, .
Duality. Since , one gets and hence
So the matrix is the -dual of the matrix, as in the continuous setting.
Uses in EIT. Weighting the voltage misfit in or the current misfit in matches the mapping properties of the Neumann-to-Dirichlet Map and changes how high-frequency boundary data are weighted relative to low-frequency data (see Data Fidelity Terms).
On uniform rectangle grids the eigenvectors are cosine modes, and the same construction works for σ in the interior at cost; see Spectral Sobolev Norms on Rectangles.
References
- M. Arioli, D. Loghin (2009). Discrete Interpolation Norms with Applications. SIAM J. Numer. Anal. 47(4), 2924–2951. doi:10.1137/080729360
- W. McLean (2000). Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press. ISBN 978-0-521-66375-5