  Many problems in science and engineering can be rigorously recast into
minimizing a suitable energy functional. We have been developing efficient and
flexible solution strategies to tackle various minimization problems by
employing finite element discretization with P1 triangular elements [1,2]. An
extension to rectangular hp-finite elements in 2D is introduced in this
contribution.
