Hi All,


I am attempting to use the complex branch of Firedrake but I'm having problems - I hope someone can help!


I installed Firedrake using the firedrake-install script from the complex branch, using the command:

python3 firedrake-install --disable-ssh --complex --minimal-petsc --adjoint --slope --install gusto --package-branch ufl complex --package-branch tsfc complex --package-branch PyOP2 complex --package-branch COFFEE complex --venv complex-firedrake


I have activated the correct virtual environment and checked out the complex branch of firedrake:


(complex-firedrake) jshipton@JSMac~/complex-firedrake/src/firedrake $ git branch

complex

  master


I'm not sure which of the tests are supposed to work on this branch so I just tried the following:


from firedrake import *

mesh = UnitSquareMesh(10, 10)

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


and got the error:

Traceback (most recent call last):

  File "test_complex.py", line 4, in <module>

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

  File "<decorator-gen-716>", line 2, in FunctionSpace

  File "/Users/jshipton/complex-firedrake/src/PyOP2/pyop2/profiling.py", line 60, in wrapper

    return f(*args, **kwargs)

  File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/functionspace.py", line 118, in FunctionSpace

    element = make_scalar_element(mesh, family, degree, vfamily, vdegree)

  File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/functionspace.py", line 40, in make_scalar_element

    mesh.init()

  File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/mesh.py", line 993, in init

    self._callback(self)

  File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/mesh.py", line 1335, in callback

    (self.num_vertices(), geometric_dim))

  File "firedrake/dmplex.pyx", line 851, in firedrake.dmplex.reordered_coords

    ELSE:

ValueError: Buffer dtype mismatch, expected 'PetscReal' but got 'complex double'


Did I miss a step or is this currently broken?


I've attached the install log just incase that has any useful information but it looks happy with everything...


Thanks!


Jemma