reduce the grid size

Author: MikaelSlevinsky

examples/disk.jl

@@ -26,7 +26,7 @@ plotlyjs()
 f = (x,y) -> (x^2-y^2+1)/((x^2-y^2+1)^2+(2x*y+1)^2)
 
 # The Zernike polynomial degree:
-N = 20
+N = 15
 M = 4N-3
 
 # The radial grid:
@@ -41,7 +41,7 @@ F = [f(r*cospi(θ), r*sinpi(θ)) for r in r, θ in θ]
 # We superpose a surface plot of $f$ on top of the grid:
 X = [r*cospi(θ) for r in r, θ in θ]
 Y = [r*sinpi(θ) for r in r, θ in θ]
-scatter3d(vec(X), vec(Y), vec(0X); markersize=0.75, markercolor=:red, size=(1000, 1000))
+scatter3d(vec(X), vec(Y), vec(0F); markersize=0.75, markercolor=:red, size=(1000, 1000))
 surface!(X, Y, F; legend=false, xlabel="x", ylabel="y", zlabel="f")
 
 # We precompute a (generalized) Zernike--Chebyshev×Fourier plan:
@@ -95,7 +95,7 @@ F = [f(x[n], w[n]*z) for n in 1:N, z in z]
 # We superpose a surface plot of $f$ on top of the grid:
 X = [x for x in x, z in z]
 Y = [w*z for w in w, z in z]
-scatter3d(vec(X), vec(Y), vec(0X); markersize=0.75, markercolor=:green, size=(1000, 1000))
+scatter3d(vec(X), vec(Y), vec(0F); markersize=0.75, markercolor=:green, size=(1000, 1000))
 surface!(X, Y, F; legend=false, xlabel="x", ylabel="y", zlabel="f")
 
 # We precompute a Dunkl-Xu--Chebyshev plan: