Thank you. Still something wrong, now I'm getting this: File "/home/mmtjs/work/programs/lin_coupled_fd_2d/lib/beam.py", line 49, in build_function_spaces self.V = fd.VectorFunctionSpace(self.mesh, "CG", 1) File "/opt/firedrake/1510/lib/python2.7/site-packages/PyOP2-0.11.0_388_gde0c811-py2.7-linux-x86_64.egg/pyop2/caching.py", line 161, in __new__ obj = make_obj() File "/opt/firedrake/1510/lib/python2.7/site-packages/PyOP2-0.11.0_388_gde0c811-py2.7-linux-x86_64.egg/pyop2/caching.py", line 142, in make_obj obj.__init__(*args, **kwargs) File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/functionspace.py", line 613, in __init__ mesh.init() File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/mesh.py", line 475, in init self._callback(self) File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/mesh.py", line 438, in callback self._coordinate_fs = functionspace.VectorFunctionSpace(self, "Lagrange", 1) File "/opt/firedrake/1510/lib/python2.7/site-packages/PyOP2-0.11.0_388_gde0c811-py2.7-linux-x86_64.egg/pyop2/caching.py", line 161, in __new__ obj = make_obj() File "/opt/firedrake/1510/lib/python2.7/site-packages/PyOP2-0.11.0_388_gde0c811-py2.7-linux-x86_64.egg/pyop2/caching.py", line 142, in make_obj obj.__init__(*args, **kwargs) File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/functionspace.py", line 615, in __init__ super(VectorFunctionSpace, self).__init__(mesh, element, name, dim=dim, rank=1) File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/functionspace.py", line 140, in __init__ dmplex.get_facet_nodes(mesh.exterior_facets.facet_cell, File "/opt/firedrake/1510/lib/python2.7/site-packages/PyOP2-0.11.0_388_gde0c811-py2.7-linux-x86_64.egg/pyop2/utils.py", line 64, in __get__ obj.__dict__[self.__name__] = result = self.fget(obj) File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/mesh.py", line 505, in exterior_facets boundary_ids, unique_markers=unique_ids) File "/opt/firedrake/1510/lib/python2.7/site-packages/firedrake/mesh.py", line 54, in __init__ "Every marker has to be contained in unique_markers" AssertionError: Every marker has to be contained in unique_markers ________________________________________ From: firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> on behalf of Lawrence Mitchell <lawrence.mitchell@imperial.ac.uk> Sent: 17 February 2016 15:34 To: firedrake@imperial.ac.uk Subject: Re: [firedrake] dimension with external mesh On 17/02/16 15:28, Tomasz Salwa [RPG] wrote:
Attached.
print repr(mesh.ufl_cell()) returns:
Cell('interval', 1)
OK, this is your problem. The mesh is actually a mesh of intervals. Looking at your geo file, I see you need to add a physical surface, otherwise only the physical line is meshed. Something like: Physical Surface(6) = {6}; Then when I generate and read in the mesh, I get quadrilateral cells. Cheers, Lawrence