Dear all, I made a periodic mesh in Gmsh with the opposite sides being identified with one another except for the top and bottom. The script I used for the .geo file is this: Box(1) = {0, 0, 0, 1, 1, 0.67}; Periodic Surface {2} = {1} Translate {1, 0, 0}; Periodic Surface {4} = {3} Translate {0, 1, 0}; But when I tried to define a vector function space over this mesh, I got the following error: File "main.py", line 59, in <module> H1_vec = VectorFunctionSpace(mesh, "CG", 1, 5) File "/home/andrewh/firedrake/src/firedrake/firedrake/functionspace.py", line 177, in VectorFunctionSpace sub_element = make_scalar_element(mesh, family, degree, vfamily, vdegree) File "/home/andrewh/firedrake/src/firedrake/firedrake/functionspace.py", line 40, in make_scalar_element mesh.init() File "/home/andrewh/firedrake/src/firedrake/firedrake/mesh.py", line 1295, in init self._callback(self) File "/home/andrewh/firedrake/src/firedrake/firedrake/mesh.py", line 1678, in callback (self.num_vertices(), geometric_dim)) File "firedrake/cython/dmcommon.pyx", line 1077, in firedrake.cython.dmcommon.reordered_coords ValueError: cannot reshape array of size 1172043 into shape (130337,3) Does anyone know what's going on here? The mesh I made seems very simple so I'm not sure what could have gone wrong. Andrew Hicks