5 research outputs found
Iteration-complexity of an inner accelerated inexact proximal augmented Lagrangian method based on the classical Lagrangian function and a full Lagrange multiplier update
This paper establishes the iteration-complexity of an inner accelerated
inexact proximal augmented Lagrangian (IAIPAL) method for solving
linearly-constrained smooth nonconvex composite optimization problems that is
based on the classical augmented Lagrangian (AL) function. More specifically,
each IAIPAL iteration consists of inexactly solving a proximal AL subproblem by
an accelerated composite gradient (ACG) method followed by a classical Lagrange
multiplier update. Under the assumption that the domain of the composite
function is bounded and the problem has a Slater point, it is shown that IAIPAL
generates an approximate stationary solution in ACG
iterations (up to a logarithmic term) where is the tolerance for both
stationarity and feasibility. Moreover, the above bound is derived without
assuming that the initial point is feasible. Finally, numerical results are
presented to demonstrate the strong practical performance of IAIPAL
On QuasiâNewton methods in fast Fourier transformâbased micromechanics
This work is devoted to investigating the computational power of QuasiâNewton methods in the context of fast Fourier transform (FFT)âbased computational micromechanics. We revisit FFTâbased NewtonâKrylov solvers as well as modern QuasiâNewton approaches such as the recently introduced Anderson accelerated basic scheme. In this context, we propose two algorithms based on the BroydenâFletcherâGoldfarbâShanno (BFGS) method, one of the most powerful QuasiâNewton schemes. To be specific, we use the BFGS update formula to approximate the global Hessian or, alternatively, the local material tangent stiffness. Both for Newton and QuasiâNewton methods, a globalization technique is necessary to ensure global convergence. Specific to the FFTâbased context, we promote a Dongâtype line search, avoiding function evaluations altogether. Furthermore, we investigate the influence of the forcing term, that is, the accuracy for solving the linear system, on the overall performance of inexact (Quasiâ)Newton methods. This work concludes with numerical experiments, comparing the convergence characteristics and runtime of the proposed techniques for complex microstructures with nonlinear material behavior and finite as well as infinite material contrast
Low-rank and sparse reconstruction in dynamic magnetic resonance imaging via proximal splitting methods
Dynamic magnetic resonance imaging (MRI) consists of collecting multiple MR images in time, resulting in a spatio-temporal signal. However, MRI intrinsically suffers from long acquisition times due to various constraints. This limits the full potential of dynamic MR imaging, such as obtaining high spatial and temporal resolutions which are crucial to observe dynamic phenomena. This dissertation addresses the problem of the reconstruction of dynamic MR images from a limited amount of samples arising from a nuclear magnetic resonance experiment. The term limited can be explained by the approach taken in this thesis to speed up scan time, which is based on violating the Nyquist criterion by skipping measurements that would be normally acquired in a standard MRI procedure. The resulting problem can be classified in the general framework of linear ill-posed inverse problems. This thesis shows how low-dimensional signal models, specifically lowrank and sparsity, can help in the reconstruction of dynamic images from partial measurements. The use of these models are justified by significant developments in signal recovery techniques from partial data that have emerged in recent years in signal processing. The major contributions of this thesis are the development and characterisation of fast and efficient computational tools using convex low-rank and sparse constraints via proximal gradient methods, the development and characterisation of a novel joint reconstructionâseparation method via the low-rank plus sparse matrix decomposition technique, and the development and characterisation of low-rank based recovery methods in the context of dynamic parallel MRI. Finally, an additional contribution of this thesis is to formulate the various MR image reconstruction problems in the context of convex optimisation to develop algorithms based on proximal splitting methods