The FEM package Gridap.jl stalls with Julia 1.6.0-beta1, but works with Julia 1.5.3 . #39190
Description
I have seen a major regression in our FEM package Gridap.jl with Julia 1.6.0-beta1 and Julia nightly.
I am not even able to run a very simple code like this one (Julia simply stalls):
using Gridap
model = CartesianDiscreteModel((0,1,0,1),(2,2))
Ω = Triangulation(model)
dΩ = Measure(Ω,1)
reffe = ReferenceFE(lagrangian,Float64,1)
V = TestFESpace(model,reffe,conformity=:H1)I think the problem is in the last line.
Everything has worked fine in Gridap with previous Julia versions including Julia 1.5.3.
I am aware that function TestFESpace leads to deep stacks at some points, but Julia 1.5.3 seems to be able to deal with it.
We, the Gridap.jl folks, would be very grateful if this can be fixed before releasing Julia 1.6.0. We are ready to help.
As a reference, the Gridap test suite and the Tutorials work well for Julia 1.5.3, but they stall with Julia 1.6.0-beta1 and with Julia nightly, see the logs of these GH actions (this also shows that the problem can be reproduced in other machines):
https://github.com/gridap/Gridap.jl/actions/runs/476998614
https://github.com/gridap/Tutorials/actions/runs/477005492
Thanks for the help!