dolfin team mailing list archive

[Martin Sandve Alnæs <question158115@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Martin Sandve Alnæs proposed the following answer:
Then I suggest you use the pure index notation.

Martin
Den 19. mai 2011 18.01 skrev "Gerard Awanou" <
question158115@xxxxxxxxxxxxxxxxxxxxx> følgende:
> 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.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help : https://help.launchpad.net/ListHelp

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