On the non-convexity issue in the radial Calderón problem

A classical approach to the Calderón problem is to estimate the unknown conductivity by solving a nonlinear least-squares problem. It is generally believed that it leads to a nonconvex optimization problem which is riddled with bad local minimums. This has motivated the development of reconstruction methods based on convex optimization, one recent contribution being the nonlinear convex semidefinite programming approach of Harrach (2023). In this work, we investigate the computational viability of this convex approach in a simple setting where the conductivities are piecewise constant and radial. We implement this convex reconstruction method and compare it extensively to the least squares approach. Our experiments suggest that this convex programming approach only allows to accurately estimate the unknown for problems with a very small size. Moreover, surprisingly, it is consistently outperformed by Newton-type least squares solvers, which are also faster and require less measurements. We revisit the issue of nonconvexity in this piecewise constant radial setting and prove that, contrary to previous claims, there are no local minimums in the case of two scalar unknowns with no measurement noise. We also provide a partial proof of this result in the general setting which holds under a numerically verifiable assumption.