How to define a function space using a BoundaryMesh - disagreement between geometric and toplogical dimension

Asked by Renier Marchand

Hi,

I am trying to create a FunctionSpace on the boundary of a problem domain.

Specifically using BoundaryMesh as follows:

import dolfin as D

mesh = D.UnitCube(2,2,2)
print mesh.topology().dim(), mesh.geometry().dim()
V = D.FunctionSpace(mesh, 'N1curl', 1)

boundary_mesh = D.BoundaryMesh(mesh)
print boundary_mesh.topology().dim(), boundary_mesh.geometry().dim()
W = D.FunctionSpace(boundary_mesh, 'N1curl', 1)

The last statement gives an error:
RuntimeError: *** Error: Geometric dimension of the UFC dofmap and the Mesh do not match.

This also happens when I try to use other types of elements. I assume this is because there are a difference between the topological and geometric dimension.

Is this correct? Or is this an error?

Renier

Question information

Language:
English Edit question
Status:
Answered
For:
DOLFIN Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Martin Sandve Alnæs (martinal) said :
#1

There are plans to implement function spaces on boundaries next year.

Can you help with this problem?

Provide an answer of your own, or ask Renier Marchand for more information if necessary.

To post a message you must log in.