Coupled solve

dune.mmesh.iterativeSolve(schemes, targets, callback=None, iter=100, tol=1e-08, f_tol=None, verbose=False, accelerate=False)[source]

Helper function to solve bulk and interface scheme coupled iteratively.

Parameters
  • schemes – pair of schemes

  • targets – pair of discrete functions that should be solved for AND that are used in the coupling forms

  • callback – update function that is called every time before solving a scheme

  • iter – maximum number of iterations

  • tol – objective tolerance between two iterates in two norm

  • verbose – print residuum for each iteration

  • accelerate – use a vector formulation of Aitken’s fix point acceleration proposed by Irons and Tuck.

Returns

if converged and number of iterations

Note

The targets also must be used in the coupling forms.

dune.mmesh.monolithicSolve(schemes, targets, callback=None, iter=30, tol=1e-08, f_tol=1e-05, eps=1.49012e-08, verbose=0)[source]
Helper function to solve bulk and interface scheme coupled monolithically.

A newton method assembling the underlying jacobian matrix. The coupling jacobian blocks are evaluated by finite differences. We provide a fast version with a C++ backend using UMFPACK.

Parameters
  • schemes – pair of schemes

  • targets – pair of discrete functions that should be solved for AND that are used in the coupling forms

  • callback – update function that is called every time before solving a scheme

  • iter – maximum number of iterations

  • tol – objective residual of iteration step in two norm

  • f_tol – objective residual of function value in two norm

  • eps – step size for finite difference

  • verbose – 1: print residuum for each newton iteration, 2: print details

Returns

if converged