It would be good to have a more efficient interface than ```julia weights = diag(mass_matrix(D)) ```