I was wondering if the static condensation steps (involving harmonic_extension_trans, etc…) are performed internally by BVP if eliminate_internal is set for the bilinear form or if we should only use BVP if static condensation is not required.
No – as far as I know, there is no problem. I had just been looking at some code which seemed to indicate that BVP performed those steps simply based on the ‘eliminate_internal’ keyword, but I couldn’t find documentation to confirm that. That is a nice feature!