@hild here is the code for mixed laplacian with computation of the divergence