Firedrake CUDA bug

cuda_broken_109e128133d6622d99d80fe99144cf269e1039bf.py

#!/usr/bin/env python
# Run with PYOP2_BACKEND=cuda
from firedrake import *

degree = 1
mesh = UnitSquareMesh(10, 10)

def p():
    V = FunctionSpace(mesh, "CG", degree)
    u = TrialFunction(V)
    v = TestFunction(V)
    a = u * v * dx
    A = assemble(a)
    A.M

# This should work.but fails with:
# AttributeError: 'Dat' object has no attribute 'vec_ro'
p()

Traceback

Traceback (most recent call last):
  File "../project/bugs/cuda_broken_109e128133d6622d99d80fe99144cf269e1039bf.py", line 17, in <module>
    p()
  File "../project/bugs/cuda_broken_109e128133d6622d99d80fe99144cf269e1039bf.py", line 8, in p
    V = FunctionSpace(mesh, "CG", degree)
  File "/homes/hjd11/proj/Firedrake/PyOP2/pyop2/caching.py", line 161, in __new__
    obj = make_obj()
  File "/homes/hjd11/proj/Firedrake/PyOP2/pyop2/caching.py", line 142, in make_obj
    obj.__init__(*args, **kwargs)
  File "/homes/hjd11/proj/Firedrake/firedrake/firedrake/functionspace.py", line 518, in __init__
    mesh.init()
  File "/homes/hjd11/proj/Firedrake/firedrake/firedrake/mesh.py", line 421, in init
    self._callback(self)
  File "/homes/hjd11/proj/Firedrake/firedrake/firedrake/mesh.py", line 388, in callback
    self._coordinate_fs = functionspace.VectorFunctionSpace(self, "Lagrange", 1)
  File "/homes/hjd11/proj/Firedrake/PyOP2/pyop2/caching.py", line 161, in __new__
    obj = make_obj()
  File "/homes/hjd11/proj/Firedrake/PyOP2/pyop2/caching.py", line 142, in make_obj
    obj.__init__(*args, **kwargs)
  File "/homes/hjd11/proj/Firedrake/firedrake/firedrake/functionspace.py", line 585, in __init__
    super(VectorFunctionSpace, self).__init__(mesh, element, name, dim=dim, rank=1)
  File "/homes/hjd11/proj/Firedrake/firedrake/firedrake/functionspace.py", line 121, in __init__
    with function.Function(self).dat.vec_ro as v:
AttributeError: 'Dat' object has no attribute 'vec_ro'