Same file with your added line. I use: gmsh -2 mesh.geo ________________________________________ 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:57 To: firedrake@imperial.ac.uk Subject: Re: [firedrake] dimension with external mesh On 17/02/16 15:47, Tomasz Salwa [RPG] wrote:
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
This sounds like somehow your mesh file is corrupted (or invalid). Is it the same geo as before? Can you send the gmsh command you used to create it? Lawrence