Use meshmode.discretization.poly_element.{mass_matrix, diff_matrices}