dolfin team mailing list archive

[dbeacham <question91944@xxxxxxxxxxxxxxxxxxxxx>]

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

dbeacham posted a new comment:
I assume the above discussion regarding dolfin assembly over meshes embedded in higher degree space (ie 1d mesh in R^2),
is why the below code only prints the x[0] correctly, with x[1] incorrectly taking the value 0,0.5 or 1?

from dolfin import *

class MyExp(Expression):
    def eval(self, values, x):
        print x # <-----------------prints x[0], x[1]:  but x[1] is not true value
        values[0] = 0

mesh = UnitSquare(10,10)
bmesh = Mesh(BoundaryMesh(mesh))
for i,j in [(i,j) for i in range(2) for j in range(2)]:
    bmesh.init(i,j)

V = VectorFunctionSpace(bmesh, "CG", 2)

exp = MyExp()
f = Function(V)
f.interpolate(exp)

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