Require same number of functions as nodes in nodal_mass_matrix_for_face