dolfin team mailing list archive

[Gerard Awanou <question158115@xxxxxxxxxxxxxxxxxxxxx>]

Question #158115 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/158115

    Status: Answered => Open

Gerard Awanou is still having a problem:
Thanks Mr Wells and Alnaes

The entire bilinear form a(u,v) reads as follows

\sum_T \int_T D^2u:D^2v dx + \sum_e \int_e average(\partial^2 u/\partial n^2)*jump(\partial v/\partial n)
average(\partial^2 v/\partial n^2)*jump(\partial u/\partial n) ds + penalty term = F(v)

The first part can be handled as: inner(grad(grad(u)),
grad(grad(v)))*dx. The jump term is encoded as in jump(grad(v), n)).

The suggestion dot(n, grad(dot(grad(u), n))) does not seem right. Once a
dot product is taken, the second grad operator will act only on the
edge. Perhaps dot(n,dot(grad(grad(u)),n ) but it does not 'compile'.

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.