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

https://www.cse.buffalo.edu/~knepley/