segmentation fault in create_transfer_matrix (PETScDMCollection)

[Søren Madsen]

Hi!

I need an interpolation matrix between meshes with N1curl elements, but create_transfer_matrix seg. faults when using these elements. It does not seg. fault when using Lagrange elements. Here is a MWE

from dolfin import *

N=5
meshes = [RectangleMesh.create([Point(0,0), Point(1,1)],[N,N],CellType.Type.triangle),
          RectangleMesh.create([Point(0,0), Point(1,1)],[2*N,2*N],CellType.Type.triangle)]
ele="N1curl"
#ele="Lagrange"
V = [FunctionSpace(mesh, ele, 1) for mesh in meshes]

dm=dm_collection = PETScDMCollection(V)

T=dm.create_transfer_matrix(V[0],V[1])

print(V[1].dim(),V[0].dim())
print(T.size(0),T.size(1))

When running without MPI I just get a seg. fault. When running with -n 2 I get a bit more information

Process 0: Computed global bounding box tree with 3 boxes.
Process 1: Computed global bounding box tree with 3 boxes.
corrupted double-linked list

‌

Is N1curl not supported by this function or is this a (known) bug?

‌

Best Regards,

Søren