[firedrake] DG0 elements in 3D mesh ; TensorProductCell error

Hi all, Could somebody help me with this? I'm having the following error msg when trying to assemble the bilinear form bellow. a = jump(u)*jump(v)*dS petsc_a = assemble(a).M.handle ValueError: interior_facet integral cannot be used with a

TensorProductCell; need to distinguish between vertical and horizontal contributions.

I'm using a 3D mesh built as an extruded quadrilateral mesh. u and v are in DG0 function space - FunctionSpace(mesh, "DG", 0). Ultimately, I'd like to solve a heterogeneous diffusion equation in 3D. Thanks in advance -- Rodrigo