MatrixFree spacedim different than dim

[giuseppe orlando]

Hello everyone,

I am writing this post to ask a curiosity about the MatrixFree infrastructure. Is it possible to employ it for solving PDEs on surface (like step-38 tutorial). 

Indeed, looking at the documentation of the reinit function, it seems that only a 'DofHanlder<dim>' instance is accepted, whereas one would need a 'DofHandler<dim, spacedim>', with 'dim' different than 'spacedim' in order to handle this case. Am I correct or is there any workaround I am currently missing?

Best

Giuseppe