@inproceedings{KrauseWeiser2021, author = {Krause, Rolf and Weiser, Martin}, title = {Multilevel augmented Lagrangian solvers for overconstrained contact formulations}, volume = {71}, booktitle = {ESAIM: ProcS}, doi = {10.1051/proc/202171175}, pages = {175 -- 184}, year = {2021}, abstract = {Multigrid methods for two-body contact problems are mostly based on special mortar discretizations, nonlinear Gauss-Seidel solvers, and solution-adapted coarse grid spaces. Their high computational efficiency comes at the cost of a complex implementation and a nonsymmetric master-slave discretization of the nonpenetration condition. Here we investigate an alternative symmetric and overconstrained segment-to-segment contact formulation that allows for a simple implementation based on standard multigrid and a symmetric treatment of contact boundaries, but leads to nonunique multipliers. For the solution of the arising quadratic programs, we propose augmented Lagrangian multigrid with overlapping block Gauss-Seidel smoothers. Approximation and convergence properties are studied numerically at standard test problems.}, language = {en} } @book{DeuflhardWeiser2020, author = {Deuflhard, Peter and Weiser, Martin}, title = {Numerische Mathematik 3. Adaptive L{\"o}sung partieller Differentialgleichungen}, edition = {2}, publisher = {de Gruyter}, isbn = {978-3-11-069168-9}, doi = {10.1515/9783110689655}, pages = {456}, year = {2020}, language = {de} } @article{SchielaStoeckleinWeiser2021, author = {Schiela, Anton and St{\"o}cklein, Matthias and Weiser, Martin}, title = {A primal dual projection algorithm for efficient constraint preconditioning}, volume = {43}, journal = {SIAM Journal on Scientific Computing}, number = {6}, doi = {10.1137/20M1380739}, pages = {A4095 -- A4120}, year = {2021}, abstract = {We consider a linear iterative solver for large scale linearly constrained quadratic minimization problems that arise, for example, in optimization with PDEs. By a primal-dual projection (PDP) iteration, which can be interpreted and analysed as a gradient method on a quotient space, the given problem can be solved by computing sulutions for a sequence of constrained surrogate problems, projections onto the feasible subspaces, and Lagrange multiplier updates. As a major application we consider a class of optimization problems with PDEs, where PDP can be applied together with a projected cg method using a block triangular constraint preconditioner. Numerical experiments show reliable and competitive performance for an optimal control problem in elasticity.}, language = {en} } @misc{BaumannDudaSchielaetal.2026, author = {Baumann, Felix and Duda, Georg and Schiela, Anton and Weiser, Martin}, title = {Identification of Stress in Heterogeneous Contact Models}, volume = {173}, journal = {Non-Smooth and Complementarity-Based Distributed Parameter Systems}, editor = {Hinterm{\"u}ller, Michael and Herzog, Roland and Kanzow, Christian and Ulbrich, Michael and Ulbrich, Stefan}, publisher = {Springer Nature}, year = {2026}, abstract = {We develop a heterogeneous model of the lower limb system to simulate muscle forces and stresses acting on the knee joint. The modelling of the bone dynamics leads to an index-3 DAE, which we discretize by higher order collocation methods. Furthermore, we present an elastomechanical contact knee joint model of the articular cartilage. For the solution of the contact problem we develop an efficient multigrid solver, based on an Augmented-Lagrangian relaxation of the contact constraints. We formulate the identification of joint forces and resulting stresses with respect to different knee joint models as an inverse problem based on medical gait data.}, language = {en} }