Stack overflow with ApproxFun type

[aml5600]

MWE:

using ApproxFun, OrdinaryDiffEq
Ω = -4..4
g = Fun(x->exp(-0.5*x^2), Ω)
fode!(du,u,p,t) = du[1] = u[1]' 
ode_prob = ODEProblem(fode!, [g], (0., 1.); saveat=0.1)
sol = solve(ode_prob, Tsit5())

causes

ERROR: LoadError: StackOverflowError:
Stacktrace:
 [1] recursive_unitless_bottom_eltype(::Type{Any}) at /Users//.julia/packages/RecursiveArrayTools/oacFD/src/utils.jl:86 (repeats 79984 times)

[aml5600]

removing the array wrappers causes:

julia> sol = solve(ode_prob, Tsit5())
ERROR: MethodError: no method matching similar(::Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}}, ::Type{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}}})
Closest candidates are:
  similar(::VSCodeServer.JuliaInterpreter.Compiled, ::Any) at /Users/aleonard31/.vscode/extensions/julialang.language-julia-1.1.25/scripts/packages/JuliaInterpreter/src/types.jl:7
  similar(::Array{T,1}, ::Type) where T at array.jl:377
  similar(::Array{T,2}, ::Type) where T at array.jl:378
  ...
Stacktrace:
 [1] alg_cache(::Tsit5, ::Function, ::Function, ::Type{T} where T, ::Type{T} where T, ::Type{T} where T, ::Function, ::Function, ::ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing}, ::Float64, ::Float64, ::Float64, ::DiffEqBase.NullParameters, ::Bool, ::Val{true}) at /Users/aleonard31/.julia/packages/OrdinaryDiffEq/Yv0iV/src/caches/low_order_rk_caches.jl:357
 [2] __init(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Tsit5, ::Tuple{}, ::Tuple{}, ::Tuple{}, ::Type{Val{true}}; saveat::Float64, tstops::Tuple{}, d_discontinuities::Tuple{}, save_idxs::Nothing, save_everystep::Bool, save_on::Bool, save_start::Bool, save_end::Bool, callback::Nothing, dense::Bool, calck::Bool, dt::Float64, dtmin::Nothing, dtmax::Float64, force_dtmin::Bool, adaptive::Bool, gamma::Rational{Int64}, abstol::Nothing, reltol::Nothing, qmin::Rational{Int64}, qmax::Int64, qsteady_min::Int64, qsteady_max::Int64, qoldinit::Rational{Int64}, fullnormalize::Bool, failfactor::Int64, beta1::Nothing, beta2::Nothing, maxiters::Int64, internalnorm::typeof(DiffEqBase.ODE_DEFAULT_NORM), internalopnorm::typeof(LinearAlgebra.opnorm), isoutofdomain::typeof(DiffEqBase.ODE_DEFAULT_ISOUTOFDOMAIN), unstable_check::typeof(DiffEqBase.ODE_DEFAULT_UNSTABLE_CHECK), verbose::Bool, timeseries_errors::Bool, dense_errors::Bool, advance_to_tstop::Bool, stop_at_next_tstop::Bool, initialize_save::Bool, progress::Bool, progress_steps::Int64, progress_name::String, progress_message::typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE), userdata::Nothing, allow_extrapolation::Bool, initialize_integrator::Bool, alias_u0::Bool, alias_du0::Bool, initializealg::OrdinaryDiffEq.DefaultInit, kwargs::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}) at /Users/aleonard31/.julia/packages/OrdinaryDiffEq/Yv0iV/src/solve.jl:296
 [3] __solve(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Tsit5; kwargs::Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}}) at /Users/aleonard31/.julia/packages/OrdinaryDiffEq/Yv0iV/src/solve.jl:4
 [4] solve_call(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Tsit5; merge_callbacks::Bool, kwargs::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}) at /Users/aleonard31/.julia/packages/DiffEqBase/ORQhu/src/solve.jl:65
 [5] solve_call(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Tsit5) at /Users/aleonard31/.julia/packages/DiffEqBase/ORQhu/src/solve.jl:52
 [6] solve_up(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Nothing, ::Function, ::DiffEqBase.NullParameters, ::Tsit5; kwargs::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}) at /Users/aleonard31/.julia/packages/DiffEqBase/ORQhu/src/solve.jl:86
 [7] solve_up(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Nothing, ::Function, ::DiffEqBase.NullParameters, ::Tsit5) at /Users/aleonard31/.julia/packages/DiffEqBase/ORQhu/src/solve.jl:79
 [8] solve(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Tsit5; sensealg::Nothing, u0::Nothing, p::Nothing, kwargs::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}) at /Users/aleonard31/.julia/packages/DiffEqBase/ORQhu/src/solve.jl:74
 [9] solve(::ODEProblem{Fun{Chebyshev{Interval{:closed,:closed,Int64},Float64},Float64,Array{Float64,1}},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(fode!),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Symbol,Float64,Tuple{Symbol},NamedTuple{(:saveat,),Tuple{Float64}}},DiffEqBase.StandardODEProblem}, ::Tsit5) at /Users/aleonard31/.julia/packages/DiffEqBase/ORQhu/src/solve.jl:72
 [10] top-level scope at REPL[9]:1