Compute the Laplacian of trial functions

See: Computing element-wise spatial derivatives of the solution - #4 by dokken