ConceptPythonBeginner Finite elements: hat functions, the weak form, and a heated bar u_fem = solve(*held) # solves, then puts the held zeros back in
ToolPythonBeginner scikit-fem from the ground up: the sag of a membrane under pressure u = solve(*condense(A, f, D=boundary))u_ramp = solve(*condense(A, asm(ramp, basis), D=boundary))return basis, solve(*condense(A, f, D=mesh.boundary_nodes()))u_free = solve(A, f)u_bad = solve(*condense(A, f_bad, D=boundary))