the spatial configuration $\h{J}_{\square}$.

Parallelization of key tasks
----------------------------

The library TBB [@TBB] is used to perform as many computationally intensive distributed tasks as possible. These include:

1.  the assembly of the tangent and residual contributions from the cells in the volume and on the surface;

2.  updating the state of the `ContinuumPoint` after solving the linearized problem in the Newton scheme.

The main tool for using TBB in [deal.II]{} is the [WorkStream]{} class. The bottleneck, however, for large computations is the linear solver. The code could, in theory, be parallelized in a distributed sense relatively easily. The primary challenge would be to ensure that the independent distributed volume and surface meshes can communicate, via a map of the common degrees of freedom.

The choice of linear solver and preconditioner
----------------------------------------------

The choice of linear solver and preconditioner are specified in the parameter file `parameter.prm`. The default choice is the conjugate gradient solver (the matrix problem is symmetric) with Jacobi preconditioning as provided by the [deal.II]{} library. Both this choice of preconditioner and solver are multithreaded. The choice works well for the majority of problems investigated.

Numerical results {#sec_numerical_results}
=================

The objective of the current section is to elucidate key 