minor performance improvements

Author: jishnub

README.md

@@ -145,7 +145,7 @@ julia> dPl!(P,0.5)
  0.0  0.375  0.0
 ```
 
-Note that the column number that will be populated depends on the order of the derivative, assuming 0-based indexing. Therefore `Pl!` will fill column 0, `dPl!` will fill columnn 1 and `d2Pl!` will fill column 2. Combinations of these will fill multiple columns as expected.
+Note that the column number that will be populated depends on the order of the derivative, assuming 0-based indexing. Therefore `Pl!` will fill column 0, `dPl!` will fill column 1 and `d2Pl!` will fill column 2. Combinations of these will fill multiple columns as expected.
 
 # License
 

src/LegendrePolynomials.jl

@@ -309,15 +309,15 @@ It defaults to lmax = maximum(axes(arr,1))
 function dPl!(arr::AbstractVector{<:Real},x::Real;lmax::Integer=maximum(axes(arr,1)))
 	P = zeros(axes(arr,1),0:1)
 	Pl_dPl_d2Pl!(P,x;lmax=lmax)
-	@inbounds @. arr[:] = P[:,1]
+	@views @. arr = P[:,1]
 	arr
 end
 
 function dPl!(arr::AbstractMatrix{<:Real},x::Real;
 	lmax::Integer=maximum(axes(arr,1)))
 	P = zeros(axes(arr,1),0:1)
 	Pl_dPl_d2Pl!(P,x;lmax=lmax)
-	@inbounds @. arr[:,1] = P[:,1]
+	@views @. arr[:,1] = P[:,1]
 	arr
 end
 
@@ -358,15 +358,15 @@ It defaults to lmax = maximum(axes(arr,1))
 function d2Pl!(arr::AbstractVector{<:Real},x::Real;lmax::Integer=maximum(axes(arr,1)))
 	P = zeros(axes(arr,1),0:2)
 	Pl_dPl_d2Pl!(P,x;lmax=lmax)
-	@. arr[:] = P[:,2]
+	@views @. arr = P[:,2]
 	arr
 end
 
 function d2Pl!(arr::AbstractMatrix{<:Real},
 	x::Real;lmax::Integer=maximum(axes(arr,1)))
 	P = zeros(axes(arr,1),0:2)
 	Pl_dPl_d2Pl!(P,x;lmax=lmax)
-	@. arr[:,2] = P[:,2]
+	@views @. arr[:,2] = P[:,2]
 	arr
 end
 
@@ -391,8 +391,7 @@ function dPl_d2Pl!(arr::AbstractMatrix{<:Real},x::Real;
 	P = zeros(axes(arr,1),0:2)
 	Pl_dPl_d2Pl!(P,x;lmax=lmax)
 
-	@. arr[:,1] = P[:,1]
-	@. arr[:,2] = P[:,2]
+	@views @. arr[:,1:2] = P[:,1:2]
 
 	arr
 end
@@ -418,8 +417,8 @@ function Pl_d2Pl!(arr::AbstractMatrix{<:Real},x::Real;
 	P = zeros(axes(arr,1),0:2)
 	Pl_dPl_d2Pl!(P,x;lmax=lmax)
 
-	@. arr[:,0] = P[:,0]
-	@. arr[:,2] = P[:,2]
+	@views @. arr[:,0] = P[:,0]
+	@views @. arr[:,2] = P[:,2]
 
 	arr
 end