Some more deprecation fixes

Author: dlfivefifty

src/MultivariateOrthogonalPolynomials.jl

@@ -3,7 +3,7 @@ __precompile__()
 module MultivariateOrthogonalPolynomials
 using Base, RecipesBase, ApproxFun, BandedMatrices, BlockArrays,
     FastTransforms, FastGaussQuadrature, StaticArrays, FillArrays,
-    Libdl, SpecialFunctions
+    LinearAlgebra, Libdl, SpecialFunctions
 
 # package code goes here
 import Base: values,getindex,setindex!,*,.*,+,.+,-,.-,==,<,<=,>,

src/Triangle.jl

@@ -388,21 +388,21 @@ function clenshaw2D(Jx,Jy,cfs::Vector{Vector{T}},x,y) where T
         resize!(Abk2x,Int(K))
         BLAS.blascopy!(Int(K),bk1,1,Abk2x,1)
         resize!(Abk2y,Int(K))
-        Abk2y[1:Int(K)-1]=0
-        Abk2y[end]=bk1[Int(K)+1]
+        Abk2y[1:Int(K)-1] .= 0
+        Abk2y[end] = bk1[Int(K)+1]
 
         Abk1x,Abk2x=Abk2x,Abk1x
         Abk1y,Abk2y=Abk2y,Abk1y
 
         bk2 = bk1  ::Vector{T}
 
         bk1 = (x*Abk1x) ::Vector{T}
-        Base.axpy!(y,Abk1y,bk1)
+        LinearAlgebra.axpy!(y,Abk1y,bk1)
         mul!(bk1,Bx,Abk1x,-one(T),one(T))
         mul!(bk1,By,Abk1y,-one(T),one(T))
         mul!(bk1,Cx,Abk2x,-one(T),one(T))
         mul!(bk1,Cy,Abk2y,-one(T),one(T))
-        Base.axpy!(one(T),cfs[Int(K)],bk1)
+        LinearAlgebra.axpy!(one(T),cfs[Int(K)],bk1)
     end
 
     K = Block(1)

test/TriangleTest.jl

@@ -100,6 +100,36 @@ ncoefficients(f)
 
 cjt
 
+f = Fun((x,y) -> cos(100x*y), KoornwinderTriangle(0.0,0.0,0.0))
+        plot(f)
+
+ncoefficients(f)
+
+
+@time f = Fun((x,y) -> cos(100x*y), Triangle(),100)
+    @time g = Laplacian()*f
+    34
+ncoefficients(f)
+
+@which Derivative([1,0])*f
+
+@which ApproxFun.mul_coefficients(view(L,1:ncoefficients(f), 1:ncoefficients(f)),
+    f.coefficients)
+
+using BlockBandedMatrices
+V = view(L,1:ncoefficients(f), 1:ncoefficients(f))
+    @which AbstractMatrix(V)
+
+
+L = Derivative([1,0]) : space(f); N = 100
+    @time L[Block.(1:N), Block.(1:N)]
+
+ncoefficients(f)
+plot(g)
+
+
+
+
 @time f = Fun((x,y)->cos(100x*y),KoornwinderTriangle(0.0,-0.5,-0.5)); # 1.15s
 @time f = Fun((x,y)->cos(500x*y),KoornwinderTriangle(0.0,-0.5,-0.5),40_000); # 0.2
 @test f(0.1,0.2) ≈ cos(500*0.1*0.2)