Dear all, while trying to run a fenics model, which used to work about two years ago, I have found the following problem: I generate diffpack meshes using netgen and convert them to .xml using dolfin-convert. The meshes produced by the recent version of dolfin-convert contain a new section , which looks like this: ... .... The problem is that I am not able to define and project an Expression on such mesh. This can be demonstrated by the following minimal example, where pipe.xml and pipe2.xml are meshes without and with the section respectively. The projected field v2 is zero everywhere, while v1 is the expected parabolic velocity profile. The example meshes can be found at http://www.roucka.eu/pub/pipe.xml http://www.roucka.eu/pub/pipe2.xml from dolfin import * #pipe radius [m] radius = 0.0094 # Load mesh and subdomains mesh1 = Mesh("pipe.xml") mesh2 = Mesh("pipe2.xml") # Define function spaces V1 = VectorFunctionSpace(mesh1, "CG", 2) V2 = VectorFunctionSpace(mesh2, "CG", 2) # Create velocity profile function sqradius = radius**2 profile = "2.0*(1.0-(x[0]*x[0]+x[1]*x[1])/(" + str(sqradius) + "))" v1 = project(Expression(("0.0", "0.0", profile)), V1) v2 = project(Expression(("0.0", "0.0", profile)), V2) # this shows that both meshes are loaded correctly? plot(mesh1) plot(mesh2) # but the projections are not identical. plot(v1) plot(v2) interactive() exit() I can solve the problem by deleting the from my meshes, but I believe that it is not the correct solution. Anyone knows what am I doing wrong? Thanks in advance for any help. Stepan

When	Who	Action	New State
2013-02-25	Stepan Roucka	Give more information	Open
2013-02-25	Johan Hake	Answer	Answered
2013-02-25	Jan Blechta	Comment	Answered
2013-02-25	Stepan Roucka	Confirm	Solved
2013-02-25	Stepan Roucka	Comment	Solved

DOLFIN

History of question #222796