precompile common function calls

Author: jishnub

Project.toml

@@ -16,6 +16,7 @@ FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
 IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
+SnoopPrecompile = "66db9d55-30c0-4569-8b51-7e840670fc0c"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
@@ -34,6 +35,7 @@ FastTransforms = "0.12, 0.13, 0.14"
 FillArrays = "0.11, 0.12, 0.13"
 IntervalSets = "0.5, 0.6, 0.7"
 Reexport = "0.2, 1"
+SnoopPrecompile = "1"
 SpecialFunctions = "0.10, 1.0, 2"
 julia = "1.6"
 

src/ApproxFunOrthogonalPolynomials.jl

@@ -109,5 +109,6 @@ include("roots.jl")
 include("specialfunctions.jl")
 include("fastops.jl")
 include("show.jl")
+include("precompile.jl")
 
 end

src/precompile.jl

@@ -0,0 +1,20 @@
+using SnoopPrecompile
+
+@precompile_setup begin
+    splist = Any[Jacobi(1, 1), Chebyshev(), Ultraspherical(1)]
+    append!(splist, Any[Jacobi(1, 1, 0..1), Chebyshev(0..1), Ultraspherical(1, 0..1)])
+    push!(splist, Chebyshev(Segment(1.0+im,2.0+2im)))
+    a = Fun(ChebyshevInterval{BigFloat}(),BigFloat[1,2,3])
+    @precompile_all_calls begin
+        for S in splist
+        	v = [0.0, 1.0]
+        	f = Fun(S,v)
+        	abs(DefiniteIntegral()*f - sum(f))
+        	norm(Derivative()*f-f')
+        	norm(differentiate(integrate(f)) - f)
+        	norm(differentiate(cumsum(f))-f)
+        	norm(first(cumsum(f)))
+        end
+        Fun(sin,Interval(big"0.0", big"1.0")
+    end
+end