Skip to content

Algebraic multigrid with filtering contact (AMGFContact) preconditioner for contact linear systems#1542

Open
thartland wants to merge 19 commits intodevelopfrom
hartland1/feature/amgf
Open

Algebraic multigrid with filtering contact (AMGFContact) preconditioner for contact linear systems#1542
thartland wants to merge 19 commits intodevelopfrom
hartland1/feature/amgf

Conversation

@thartland
Copy link
Contributor

@thartland thartland commented Feb 16, 2026

This PR adds the capability of using AMGF (algebraic multigrid with filtering arXiv preprint)-preconditioned Krylov-subspace linear solvers for the Newton linear solves that arise in outer loop solves of contact problems.

  • Obtain list of contact dofs via ContactData::addContactInteraction calls.
  • Create a boundary_dof --> subspace transfer operator (mfem::HypreParMatrix) call. Much of this infastructure is in the utilities of continuationsolvers but a version most suitable for usage in smith would be best for less complex code.
  • Obtain subspace transfer (displacement dofs --> contact dofs) operator via the list of contact dofs in appropriate place for passing the subspace transfer operator to the preconditioner provided that said linear solver preconditioner is AMGF.
  • Resolve usage errors conflating print_level and preconditioner_print_level from linear solver options in the buildPreconditioner method.
  • Test.

For future investigation/potential PR(s)
(1) Make AMGFContact robust to indefiniteness issues stemming from using exact Hessians in Newton solves.
(2) More careful consideration for setting outer/inner tolerances (as well as norms used) for Newton-Krylov.

Note -- AMGFContact preconditioned Krylov-subspace methods will fails unless the problem is a contact problem. The failure is indicated by mfem generated error messages indicating that the subspace transfer operator has not been set.

@thartland thartland added the WIP Work in progress label Feb 16, 2026
@thartland thartland changed the title WIP: AMGF WIP: Algebraic multigrid with filtering (AMGF) preconditioner for contact linear systems Feb 25, 2026
@thartland thartland removed WIP Work in progress labels Feb 27, 2026
@thartland thartland changed the title WIP: Algebraic multigrid with filtering (AMGF) preconditioner for contact linear systems Algebraic multigrid with filtering (AMGF) preconditioner for contact linear systems Feb 27, 2026
@thartland thartland requested review from ebchin and tupek2 February 27, 2026 22:38
@thartland thartland changed the title Algebraic multigrid with filtering (AMGF) preconditioner for contact linear systems Algebraic multigrid with filtering contact (AMGFContact) preconditioner for contact linear systems Mar 3, 2026
…n error in the gcc Release build due to a potential null pointer deference warning (error) in the body of the iterative solver GetComm method. Replacing the iterative solver GetComm call with a HypreParMatrix GetComm call, should serve the same purpose
@thartland thartland added the ready for review Ready for active inspection by reviewers label Mar 5, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

ready for review Ready for active inspection by reviewers

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant