On 15 Sep 2020, at 09:18, Mohammad Sarraf Joshaghani <m.sarraf.j@gmail.com> wrote:
The inner(grad(sol),n)*v in last two terms give me an error: UFL:ERROR Shapes do not match: <Grad id=4986788992> and <Zero id=4998074176>.
But the projection was okay when changed them to sol*v. Do you know where is this error coming from? (FYI: my sol Function is defined on DG1 functionspace)
I always mix up the definitions of jump(foo, n). I actually prefer just to write out what you want explicitly in terms of the restricted values.
solve(F == 0, fluxes)
# Now fluxes contains your integrals.
Lawrence
Also, how this globally defined fluxes are connected to cells? In other words for a triangular mesh, is it possible to know which 3 components of fluxes.vector().array() are related to cell E?
Vt.cell_node_map().values contains for each cell the three fluxes.