Source author record

Stefan Kindermann

Stefan Kindermann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

8works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Regularization of linear inverse problems by rational Krylov methods

For approximately solving linear ill-posed problems in Hilbert spaces, we investigate the regularization properties of the aggregation method and the RatCG method. These recent algorithms use previously calculated solutions of Tikhonov regularization (respectively, Landweber iterations) to set up a new search space on which the least-squares functional is minimized. We outline how these methods can be understood as rational Krylov space methods, i.e., based on the space of rational functions of the forward operator. The main result is that these methods form an optimal-order regularization schemes when combined with the discrepancy principle as stopping rule and when the underlying regularization parameters are sufficiently large.

preprint2022arXiv

A numerical comparison of some heuristic stopping rules for nonlinear Landweber iteration

The choice of a suitable regularization parameter is an important part of most regularization methods for inverse problems. In the absence of reliable estimates of the noise level, heuristic parameter choice rules can be used to accomplish this task. While they are already fairly well-understood and tested in the case of linear problems, not much is known about their behaviour for nonlinear problems and even less in the respective case of iterative regularization. Hence, in this paper, we numerically study the performance of some of these rules when used to determine a stopping index for Landweber iteration for various nonlinear inverse problems. These are chosen from different practically relevant fields such as integral equations, parameter estimation, and tomography.

preprint2016arXiv

Convergence of the gradient method for ill-posed problems

We study the convergence of the gradient descent method for solving ill-posed problems where the solution is characterized as a global minimum of a differentiable functional in a Hilbert space. The classical least-squares functional for nonlinear operator equations is a special instance of this framework and the gradient method then reduces to Landweber iteration. The main result of this article is a proof of weak and strong convergence under new nonlinearity conditions that generalize the classical tangential cone conditions.

preprint2016arXiv

Optimization of the shape (and topology) of the initial conditions for diffusion parameter identification

The design of an experiment, e.g., the setting of initial conditions, strongly influences the accuracy of the whole process of determining model parameters from data. We impose a sensitivity-based approach for choosing optimal design variables and study the optimization of the shape (and topology) of the initial conditions for an inverse problem of a diffusion parameter identification. Our approach, although case independent, is illustrated at the FRAP (Fluorescence Recovery After Photobleaching) experimental technique. The core idea resides in the maximization of a sensitivity measure, which depends on a specific experimental setting of initial conditions. By a numerical optimization, we find an interesting pattern of increasingly complicated (with respect to connectivity) optimal initial shapes. The proposed modification of the FRAP experimental protocol is rather surprising but entirely realistic and the resulting enhancement of the parameter estimate accuracy is significant.

preprint2015arXiv

Projection methods for ill-posed problems revisited

The discretization of least-squares problems for linear ill-posed operator equations in Hilbert spaces is considered. The main subject of this article concerns conditions for convergence of the associated discretized minimum-norm least-squares solution to the exact solution using exact attainable data. The two cases of global convergence (convergence for all exact solution) or local convergence (convergence for a specific exact solution) are investigated. We review the existing results and prove new equivalent condition when the discretized solution always converges to the exact solution. An important tool is to recognize the discrete solution operator as oblique projection. Hence, global convergence can be characterized by certain subspaces having uniformly bounded angles. We furthermore derive practically useful conditions when this holds and put them into the context of known results. For local convergence we generalize results on the characterization of weak or strong convergence and state some new sufficient conditions. We furthermore provide an example of a bounded sequence of discretized solutions which does not converge at all, not even weakly.

preprint2015arXiv

Towards analytical model optimization in atmospheric tomography

Modern ground-based telescopes rely on a technology called adaptive optics (AO) in order to compensate for the loss of image quality caused by atmospheric turbulence. Next-generation AO systems designed for a wide field of view require a stable and high-resolution reconstruction of the refractive index fluctuations in the atmosphere. By introducing a novel Bayesian method, we address the problem of estimating an atmospheric turbulence strength profile and reconstructing the refractive index fluctuations simultaneously, where we only use wavefront measurements of incoming light from guide stars. Most importantly, we demonstrate how this method can be used for model optimization as well. We propose two different algorithms for solving the maximum a posteriori estimate: the first approach is based on alternating minimization and has the advantage of integrability into existing atmospheric tomography methods. In the second approach, we formulate a convex non-differentiable optimization problem, which is solved by an iterative thresholding method. This approach clearly illustrates the underlying sparsity-enforcing mechanism for the strength profile. By introducing a tuning/regularization parameter, an automated model reduction of the layer structure of the atmosphere is achieved. Using numerical simulations, we demonstrate the performance of our method in practice.

preprint2011arXiv

Analysis and Approximation of the Canonical Polyadic Tensor Decomposition

We study the least-squares (LS) functional of the canonical polyadic (CP) tensor decomposition. Our approach is based on the elimination of one factor matrix which results in a reduced functional. The reduced functional is reformulated into a projection framework and into a Rayleigh quotient. An analysis of this functional leads to several conclusions: new sufficient conditions for the existence of minimizers of the LS functional, the existence of a critical point in the rank-one case, a heuristic explanation of "swamping" and computable bounds on the minimal value of the LS functional. The latter result leads to a simple algorithm -- the Centroid Projection algorithm -- to compute suboptimal solutions of tensor decompositions. These suboptimal solutions are applied to iterative CP algorithms as initial guesses, yielding a method called centroid projection for canonical polyadic (CPCP) decomposition which provides a significant speedup in our numerical experiments compared to the standard methods.