20
20
# [documentation](https://MikaelSlevinsky.github.io/FastTransforms).
21
21
22
22
using FastTransforms, LinearAlgebra, Plots
23
+ const GENFIGS = joinpath (dirname (dirname (pathof (FastTransforms))), " docs/src/generated" )
24
+ ! isdir (GENFIGS) && mkdir (GENFIGS)
23
25
plotlyjs ()
24
26
25
27
# Our function $f$ on the disk:
@@ -41,8 +43,12 @@ F = [f(r*cospi(θ), r*sinpi(θ)) for r in r, θ in θ]
41
43
# We superpose a surface plot of $f$ on top of the grid:
42
44
X = [r* cospi (θ) for r in r, θ in θ]
43
45
Y = [r* sinpi (θ) for r in r, θ in θ]
44
- scatter3d (vec (X), vec (Y), vec (0 F); markersize= 0.75 , markercolor= :red , size= (1000 , 1000 ))
46
+ scatter3d (vec (X), vec (Y), vec (0 F); markersize= 0.75 , markercolor= :red , size= (800 , 600 ))
45
47
surface! (X, Y, F; legend= false , xlabel= " x" , ylabel= " y" , zlabel= " f" )
48
+ savefig (joinpath (GENFIGS, " zernike.html" ))
49
+ # ##```@raw html
50
+ # ##<object type="text/html" data="../zernike.html" style="width:100%;height:600px;"></object>
51
+ # ##```
46
52
47
53
# We precompute a (generalized) Zernike--Chebyshev×Fourier plan:
48
54
α, β = 0 , 0
@@ -95,8 +101,12 @@ F = [f(x[n], w[n]*z) for n in 1:N, z in z]
95
101
# We superpose a surface plot of $f$ on top of the grid:
96
102
X = [x for x in x, z in z]
97
103
Y = [w* z for w in w, z in z]
98
- scatter3d (vec (X), vec (Y), vec (0 F); markersize= 0.75 , markercolor= :green , size= (1000 , 1000 ))
104
+ scatter3d (vec (X), vec (Y), vec (0 F); markersize= 0.75 , markercolor= :green , size= (800 , 600 ))
99
105
surface! (X, Y, F; legend= false , xlabel= " x" , ylabel= " y" , zlabel= " f" )
106
+ savefig (joinpath (GENFIGS, " dunklxu.html" ))
107
+ # ##```@raw html
108
+ # ##<object type="text/html" data="../dunklxu.html" style="width:100%;height:600px;"></object>
109
+ # ##```
100
110
101
111
# We precompute a Dunkl-Xu--Chebyshev plan:
102
112
P = plan_rectdisk2cheb (F, β)
0 commit comments