Standardize usage to be

P=plan_paduatransform(v)
P*v == paduatransform(v)

Author: dlfivefifty

src/PaduaTransform.jl

@@ -18,7 +18,7 @@ end
 """
 Inverse Padua Transform maps the 2D Chebyshev coefficients to the values of the interpolation polynomial at the Padua points.
 """
-function ipaduatransform{T}(P::IPaduaTransformPlan,v::AbstractVector{T})
+function *{T}(P::IPaduaTransformPlan,v::AbstractVector{T})
     cfsmat=trianglecfsmat(P,v)
     n,m=size(cfsmat)
     scale!(view(cfsmat,:,2:m-1),0.5)
@@ -28,6 +28,8 @@ function ipaduatransform{T}(P::IPaduaTransformPlan,v::AbstractVector{T})
     return paduavals
 end
 
+ipaduatransform(v::AbstractVector) = plan_ipaduatransform(v)*v
+
 """
 Creates (n+2)x(n+1) Chebyshev coefficient matrix from triangle coefficients.
 """
@@ -89,7 +91,7 @@ end
 """
 Padua Transform maps from interpolant values at the Padua points to the 2D Chebyshev coefficients.
 """
-function paduatransform{T}(P::PaduaTransformPlan,v::AbstractVector{T})
+function *{T}(P::PaduaTransformPlan,v::AbstractVector{T})
     N=length(v)
     n=Int(cld(-3+sqrt(1+8N),2))
     vals=paduavalsmat(P,v)
@@ -104,6 +106,8 @@ function paduatransform{T}(P::PaduaTransformPlan,v::AbstractVector{T})
     return cfs
 end
 
+paduatransform(v::AbstractVector) = plan_paduatransform(v)*v
+
 """
 Creates (n+2)x(n+1) matrix of interpolant values on the tensor grid at the (n+1)(n+2)/2 Padua points.
 """

test/runtests.jl

@@ -213,12 +213,14 @@ N=div((n+1)*(n+2),2)
 v=rand(N)  #Length of v is the no. of Padua points
 Pl=plan_paduatransform(v)
 IPl=plan_ipaduatransform(v)
-@test_approx_eq paduatransform(Pl,ipaduatransform(IPl,v)) v
-@test_approx_eq ipaduatransform(IPl,paduatransform(Pl,v)) v
+@test_approx_eq Pl*(IPl*v) v
+@test_approx_eq IPl*(Pl*v) v
+@test_approx_eq Pl*v paduatransform(v)
+@test_approx_eq IPl*v ipaduatransform(v)
 
 println("Testing runtimes for (I)Padua Transforms")
-@time paduatransform(Pl,v)
-@time ipaduatransform(IPl,v)
+@time Pl*v
+@time IPl*v
 
 println("Accuracy of 2d function interpolation at a point")
 function trianglecfsmat{T}(cfs::AbstractVector{T})
@@ -250,8 +252,7 @@ function paduaeval(f::Function,x::AbstractFloat,y::AbstractFloat,m::Integer)
     pvals=Array(T,M)
     p=paduapoints(T,m)
     map!(f,pvals,p[:,1],p[:,2])
-    plan=plan_paduatransform(pvals)
-    coeffs=paduatransform(plan,pvals)
+    coeffs=paduatransform(pvals)
     cfs_mat=trianglecfsmat(coeffs)
     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