BoundaryMesh returns wrong Data

[Johannes Neumann]

When running this simple code with Dolfin 1.5/Ubuntu 14.04 and Dolfin 1.3/OpenSuse

mesh = UnitSquareMesh(1,1)
print BoundaryMesh( mesh, "exterior" ).entity_map( 1 ).array()

I get

[65871248 47760800        3        4]

so the fist two values are wrong/not initialized.

[Chris Richardson]

BoundaryMesh is quite broken, generally, and will be superceded soon (?) by MeshView. See also issue #313 and issue #381

[Johannes Neumann]

Alright. Thanks for the information. For those interested, here is a workaround in 2D until MeshView becomes available.

mesh.init( 0, 1 )
edgeCount = np.zeros( mesh.num_edges(), dtype = int )
for cell in cells( mesh ):
    edgeCount[cell.entities( 1 ) ] += 1
boundaryEdges = np.where( edgeCount < 2 ) # belong to only one triangle

[Anders Logg (Chalmers)]

Removing milestone: 1.7 (automated comment)