Jupyter Kernel crashes when doing interpolation on product space

NGSolve
1

I noticed the following bug when using a jupyter notebook. Here is the main part of the code (this also involves ngsxfem):

Whbase = VectorH1(mesh, order=FE_degree, dirichlet="bottom|top")
Wh = Restrict(Whbase, hasneg)

Qhbase = H1(mesh, order=FE_degree - 1)
Qh = Restrict(Qhbase, hasneg)
WxQ = FESpace([Wh, Qh], dgjumps=True)

gfvv_p = GridFunction(WxQ)
gfvv, gfp = gfvv_p.components

This works fine. But if I try to do this:

gfvv_p.components[0].Set(CF((0.0,0.0)))
gfvv_p.components[1].Set(CF(0))

or

gfvv.Set((0,0))
gfp.Set(0)

Then, the jupyter kernel crashes when I execute the next cell (which doesn’t do much). The same thing happens if I use the “Interpolate” method.