journal/nme Int J Numer Methods Eng. 2019;120:880-900. † Inhomogeneous Dirichlet BVPs can be traced back to homogeneous ones: If u = u 0 at Γ D , let u 0 =û 0 | Γ D hold with a functionû 0 defined on Ω, then v ∶= u −û 0 equally solves −∇ · (v) = f + ∇ · (û 0 ) in Ω with v = 0 at Γ D and the solution of the primary problem is u = v +û 0 .