In our optimization codes, we have been differentiating with respect to the geometry degrees of freedom. This was accomplished by:
c_el = domain.ufl_domain().ufl_coordinate_element()
meshFunctionSpace = fem.FunctionSpace(domain, c_el)
However this code does not run with newer versions, giving the error:
FunctionSpace.__init__() missing 1 required positional argument: 'cppV'
What is the current accepted method/best practice to get such function space?