hermit.fenics.solvers

Linear / Newton solves of the shell residual (DOLFINx 0.9 / 0.11, MUMPS).

Direct solve is a port of custom_solve_direct from femo dev_coupling: assemble the tangent and residual, apply BC lifting, solve one Newton step (exact for a linear problem). direct_residual_inputs lets the caller add extra RHS terms (the direct generalized load vector) without going through UFL assembly.

Functions

linear_solve(residual, w[, bcs, direct_residual_inputs])

Solve F(w) = 0 for a linear F: one Newton step from the current w.

Module Contents

hermit.fenics.solvers.linear_solve(residual, w, bcs=(), direct_residual_inputs=None)

Solve F(w) = 0 for a linear F: one Newton step from the current w.