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