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

topology 1/2

integrating 1+1+1 distinct blocks

topology 2/2

integrating 1+1+1 distinct blocks

solve

solving 33 dof system to machine precision using direct solver

refinement iteration 1

residual increased to 3e-17 (discarding)

solver returned with residual 2e-17

constrained 33/242 dofs

optimum value -8.33e-17

topology 1/1

integrating 1+1 distinct blocks

solve

solving 209 dof system to machine precision using direct solver

refinement iteration 1

residual decreased to 1e-15

refinement iteration 2

residual decreased to 8e-16

refinement iteration 3

residual increased to 9e-16 (discarding)

solver returned with residual 8e-16