examples/elasticity.py

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

optimize

solve

solving 63 dof system to machine precision using arnoldi solver

residual decreased by 15.7 orders using 0 krylov vectors

solver returned with residual 5e-17

constrained 63/882 dofs

optimum value 0.00e+00

solve

solving 819 dof system to machine precision using arnoldi solver

residual decreased by 14.8 orders using 0 krylov vectors

residual decreased by 0.3 orders using 1 krylov vectors

solver returned with residual 1e-15