dolfin team mailing list archive

["B. Emek Abali" <question155482@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

B. Emek Abali proposed the following answer:
K. B. Ølgaard, A. Logg, and G. N. Wells. Automated code generation for discontinuous Galerkin methods.
SIAM J. Sci. Comput., 31(2):849–864, 2008. URL http://dx.doi.org//10.1137/070710032.

has a good explanation for you,

the restriction is needed while the compiler don't know where to evaluate on an interior facet, while a facet has two cells, so you need to change 
in
bound_mem 
conc*vel[i]*normal[i] -> conc('+')*vel[i]*normal('+')[i] - conc('-')*vel[i]*normal('-')[i]

but I have not done it yet so this is just a guess !

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