examples/elasticity.py

nelems=10 Number of elements along edge.
etype=square Type of elements (square/triangle/mixed).
btype=std Type of basis function (std/spline), with availability depending on the configured element type.
degree=1 Polynomial degree.
poisson=.25 Poisson's ratio, nonnegative and strictly smaller than 1/2.

optimize

solve

solving 33 dof system to machine precision using arnoldi solver

residual decreased by 16.2 orders using 0 krylov vectors

solver returned with residual 2e-17

constrained 33/242 dofs

optimum value 0.00e+00

solve

solving 209 dof system to machine precision using arnoldi solver

residual decreased by 15.1 orders using 0 krylov vectors

residual decreased by 0.2 orders using 1 krylov vectors

solver returned with residual 1e-15