examples/laplace.py

nelems=10 Number of elements along edge.
etype=square Type of elements (square/triangle/mixed).
btype=std Type of basis function (std/spline), availability depending on the selected element type.
degree=1 Polynomial degree.

optimize

solve

solving 21 dof system to machine precision using arnoldi solver

residual decreased by 16.6 orders using 0 krylov vectors

solver returned with residual 1e-17

constrained 21/121 dofs

optimum value 4.32e-10

solve

solving 100 dof system to machine precision using arnoldi solver

residual decreased by 15.0 orders using 0 krylov vectors

residual decreased by 0.2 orders using 1 krylov vectors

solver returned with residual 1e-15

L2 error: 2.56e-04