On Mon, May 18, 2020 at 4:56 PM Rodrigo Piccinini <rbpiccinini@gmail.com> wrote:
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.

The error message is saying that dS is not alright because integrating over a vertical face is different
than over a horizontal face. There must be different measures for them, like dH and dV, but I don't
know what they are.

  Thanks,

     Matt
 
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
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake


--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener