For post-processing purpose (e.g., monitoring residuals) I find the function get_dirichlet_dof_values very useful.

However, I noticed a possible inconsistency in the API: some related functions such as get_free_dof_values or get_cell_is_dirichlet provide methods that accept a ::FEFunction, while get_dirichlet_dof_values must be called on the finite element space associated with the function. The code below illustrates the difference:

  fv = get_free_dof_values(xh)
  dv = get_dirichlet_dof_values(get_fe_space(xh))

I'm wondering if it makes sense to add an overload accepting a FEFunction, following the pattern of get_cell_is_dirichlet defined in:

See #1186 for a possible implementation. Thanks in advance