Fixes eval errors

Author: MikeAClarke

src/PaduaTransform.jl

@@ -96,10 +96,10 @@ function paduatransform{T}(P::PaduaTransformPlan,v::AbstractVector{T})
     tensorcfs=P.dctplan*vals
     m,l=size(tensorcfs)
     scale!(tensorcfs,T(2)/(n*(n+1)))
-    scale!(sub(tensorcfs,1,:),0.5)
-    scale!(sub(tensorcfs,:,1),0.5)
-    scale!(sub(tensorcfs,m,:),0.5)
-    scale!(sub(tensorcfs,:,l),0.5)
+    scale!(view(tensorcfs,1,:),0.5)
+    scale!(view(tensorcfs,:,1),0.5)
+    scale!(view(tensorcfs,m,:),0.5)
+    scale!(view(tensorcfs,:,l),0.5)
     cfs=trianglecfsvec(P,tensorcfs)
     return cfs
 end
@@ -115,11 +115,11 @@ function paduavalsmat(P::PaduaTransformPlan,v::AbstractVector)
         d=div(n+2,2)
         m=0
         @inbounds for i=1:n+1
-            vals[1+mod(i,2):2:n+1+mod(i,2),i]=sub(v,m+1:m+d)
+            vals[1+mod(i,2):2:n+1+mod(i,2),i]=view(v,m+1:m+d)
             m+=d
         end
     else
-        @inbounds vals[1:2:end]=sub(v,:)
+        @inbounds vals[1:2:end]=view(v,:)
     end
     return vals
 end

test/runtests.jl

@@ -245,17 +245,16 @@ Interpolates a 2d function at a given point using 2d Chebyshev series.
 function paduaeval(f::Function,x::AbstractFloat,y::AbstractFloat,m::Integer)
     T=promote_type(typeof(x),typeof(y))
     M=div((m+1)*(m+2),2)
-    pvals=Array(T,M)
-    p=paduapoints(m)
-    pvals=map!(f,p[:,1],p[:,2])
-    plan=plan_paduatransform(pvals)
-    coeffs=paduatransform(plan,pvals)
+    p=paduapoints(T,m)
+    map!(f,p,p[:,1],p[:,2])
+    plan=plan_paduatransform(p)
+    coeffs=paduatransform(plan,p)
     cfs_mat=trianglecfsmat(coeffs)
     cfs_mat=view(cfs_mat,1:m+1,:)
     f_x=sum([cfs_mat[k,j]*cos((j-1)*acos(x))*cos((k-1)*acos(y)) for k=1:m+1, j=1:m+1])
     return f_x
 end
-f_xy = (x,y) -> x^2*y+x^3
+f_xy = (x,y) ->x^2*y+x^3
 g_xy = (x,y) ->cos(exp(2*x+y))*sin(y)
 x=0.1;y=0.2
 m=130