Move over PointSpace

Author: dlfivefifty

src/ApproxFunBase.jl

@@ -111,6 +111,8 @@ include("Operators/Operator.jl")
 include("Caching/caching.jl")
 include("PDE/PDE.jl")
 include("Spaces/Spaces.jl")
+include("hacks.jl")
 include("testing.jl")
+include("specialfunctions.jl")
 
 end #module

src/Fun.jl

@@ -457,6 +457,8 @@ function differentiate(f::Fun,k::Integer)
     (k==0) ? f : differentiate(differentiate(f),k-1)
 end
 
+# use conj(transpose(f)) for ArraySpace
+adjoint(f::Fun) = differentiate(f)
 
 
 

src/LinearAlgebra/LinearAlgebra.jl

@@ -9,7 +9,6 @@ include("lyap.jl")
 include("bary.jl")
 include("clenshaw.jl")
 include("ultraspherical.jl")
-include("fourier.jl")
 
 include("hesseneigs.jl")
 

src/LinearAlgebra/fourier.jl

@@ -129,6 +129,10 @@ bvp(d,k) = vcat(Operator{prectype(d)}[ldiffbc(d,i) for i=0:div(k,2)-1],
 
 periodic(d,k) = Operator{prectype(d)}[Evaluation(d,leftendpoint,i)-Evaluation(d,rightendpoint,i) for i=0:k]
 
+# shorthand for second order
+ivp(d) = ivp(d,2)
+bvp(d) = bvp(d,2)
+
 
 
 for op in (:rdirichlet,:ldirichlet,:lneumann,:rneumann,:ivp,:bvp)

src/Operators/functionals/Evaluation.jl

@@ -245,64 +245,4 @@ Base.isfinite(f::Fun{CS}) where {CS<:ConstantSpace} = isfinite(Number(f))
 
 
 
-## ConstantSpace and PointSpace default overrides
-
-for SP in (:ConstantSpace,)
-    for OP in (:abs,:sign,:exp,:sqrt,:angle)
-        @eval begin
-            $OP(z::Fun{<:$SP,<:Complex}) = Fun(space(z),$OP.(coefficients(z)))
-            $OP(z::Fun{<:$SP,<:Real}) = Fun(space(z),$OP.(coefficients(z)))
-            $OP(z::Fun{<:$SP}) = Fun(space(z),$OP.(coefficients(z)))
-        end
-    end
-
-    # we need to pad coefficients since 0^0 == 1
-    for OP in (:^,)
-        @eval begin
-            function $OP(z::Fun{<:$SP},k::Integer)
-                k ≠ 0 && return Fun(space(z),$OP.(coefficients(z),k))
-                Fun(space(z),$OP.(pad(coefficients(z),dimension(space(z))),k))
-            end
-            function $OP(z::Fun{<:$SP},k::Number)
-                k ≠ 0 && return Fun(space(z),$OP.(coefficients(z),k))
-                Fun(space(z),$OP.(pad(coefficients(z),dimension(space(z))),k))
-            end
-        end
-    end
-end
 
-
-for OP in (:(Base.max),:(Base.min))
-    @eval begin
-        $OP(a::Fun{CS1,T},b::Fun{CS2,V}) where {CS1<:ConstantSpace,CS2<:ConstantSpace,T<:Real,V<:Real} =
-            Fun($OP(Number(a),Number(b)),space(a) ∪ space(b))
-        $OP(a::Fun{CS,T},b::Real) where {CS<:ConstantSpace,T<:Real} =
-            Fun($OP(Number(a),b),space(a))
-        $OP(a::Real,b::Fun{CS,T}) where {CS<:ConstantSpace,T<:Real} =
-            Fun($OP(a,Number(b)),space(b))
-    end
-end
-
-for OP in (:<,:(Base.isless),:(<=),:>,:(>=))
-    @eval begin
-        $OP(a::Fun{CS},b::Fun{CS}) where {CS<:ConstantSpace} = $OP(convert(Number,a),Number(b))
-        $OP(a::Fun{CS},b::Number) where {CS<:ConstantSpace} = $OP(convert(Number,a),b)
-        $OP(a::Number,b::Fun{CS}) where {CS<:ConstantSpace} = $OP(a,convert(Number,b))
-    end
-end
-
-# from DualNumbers
-for (funsym, exp) in Calculus.symbolic_derivatives_1arg()
-    funsym == :abs && continue
-    funsym == :sign && continue
-    funsym == :exp && continue
-    funsym == :sqrt && continue
-    @eval begin
-        $(funsym)(z::Fun{CS,T}) where {CS<:ConstantSpace,T<:Real} =
-            Fun($(funsym)(Number(z)),space(z))
-        $(funsym)(z::Fun{CS,T}) where {CS<:ConstantSpace,T<:Complex} =
-            Fun($(funsym)(Number(z)),space(z))
-        $(funsym)(z::Fun{CS}) where {CS<:ConstantSpace} =
-            Fun($(funsym)(Number(z)),space(z))
-    end
-end