Support expansion of non-type inferred functions (#105)

* Support expansion of non-type inferred functions

* Update adaptivetransform.jl

* Update adaptivetransform.jl

* Update adaptivetransform.jl

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.6.4"
+version = "0.6.5"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/adaptivetransform.jl

@@ -1,5 +1,6 @@
+
 transform_ldiv(A::AbstractQuasiArray{T}, f::AbstractQuasiArray{V}, ::Tuple{<:Any,InfiniteCardinal{0}}) where {T,V}  =
-    adaptivetransform_ldiv(convert(AbstractQuasiArray{promote_type(T,V)}, A), f)
+    adaptivetransform_ldiv(A, f)
 
 pad(c::AbstractVector{T}, ax::Union{OneTo,OneToInf}) where T = [c; Zeros(length(ax)-length(c))]
 pad(c, ax...) = PaddedArray(c, ax)
@@ -17,13 +18,15 @@ padresize!(cfs::PseudoBlockVector, m, ax...) = padresize!(cfs.blocks, m, ax...)
 
 increasingtruncations(::OneToInf) = oneto.(2 .^ (4:∞))
 increasingtruncations(::BlockedUnitRange) = broadcast(n -> Block.(oneto(n)), (2 .^ (4:∞)))
-function adaptivetransform_ldiv(A::AbstractQuasiArray{U}, f::AbstractQuasiVector{V}) where {U,V}
-    T = promote_type(eltype(U),eltype(V))
 
+
+function adaptivetransform_ldiv(A::AbstractQuasiArray{U}, f::AbstractQuasiVector{V}) where {U,V}
     r = checkpoints(A)
     fr = f[r]
     maxabsfr = norm(fr,Inf)
 
+    # don't use V as eltype might be Any
+    T = promote_type(eltype(U), eltype(typeof(first(fr))))
     tol = 20eps(real(T))
     ax = axes(A,2)
 

test/test_chebyshev.jl

@@ -473,6 +473,12 @@ import ContinuumArrays: MappedWeightedBasisLayout, Map, WeightedBasisLayout
         @test h'h ≈ cosh(2) - sinh(2)/2
         @test f'h ≈ h'f ≈ 2.498566528528904 # 1/2 π BesselI[1, 2]
     end
+
+    @testset "non-type inferred" begin
+        T = Chebyshev()
+        f = x -> abs(x) ≤ 1 ? 1 : "hi"
+        @test T \ f.(axes(T,1)) ≈ [1; zeros(∞)]
+    end
 end
 
 struct QuadraticMap{T} <: Map{T} end