Question #137463 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/137463
Status: Open => Answered
Garth Wells proposed the following answer:
It's because you haven't specified a function space for Q. If a space is
not provided, FFC will interpolate the function in a continuous Lagrange
space, with a guess at the polynomial order to preserve the convergence
order. If you want to interpolate Q in a discontinuous space, try
VDG = VectorFunctionSpace(mesh, "DG", degree)
Q = Expression(("(x[1]>=0.5? 10.0 : 0.0)","(x[0]>=0.5? 100.0 : 0.0)"),
element = VDG.ufl_element())
--
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@lists.launchpad.net
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
