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