Dear All,

How do I instruct Firedrake about the dimension of the problem when using external mesh from gmsh?

I had a working code with utility mesh. I changed:

mesh = RectangleMesh(nx, nz, Lx, Lz)

V = VectorFunctionSpace(mesh, "CG", 1)

to:

mesh = Mesh("mesh.msh")

V = VectorFunctionSpace(mesh, "CG", 1, dim=2)

yet I get the error at further expressions that worked previously:

...

  File "/home/mmtjs/work/programs/lin_coupled_fd_2d/lib/beam.py", line 73, in initialize_solvers
    + self.mu * ( fd.inner( delX, delv_B + fd.transpose(delv_B) ) )
  File "/opt/firedrake/1510/lib/python2.7/site-packages/ufl/exproperators.py", line 177, in _add
    return Sum(self, o)
  File "/opt/firedrake/1510/lib/python2.7/site-packages/ufl/algebra.py", line 55, in __new__
    error("Can't add expressions with different shapes.")
  File "/opt/firedrake/1510/lib/python2.7/site-packages/ufl/log.py", line 151, in error
    raise self._exception_type(self._format_raw(*message))
ufl.log.UFLException: Can't add expressions with different shapes.

Thank you,

Tomasz