dolfin team mailing list archive

[Garth Wells <question137463@xxxxxxxxxxxxxxxxxxxxx>]

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.