Merge pull request #2 from JuliaApproximation/master

merge updates

Author: TSGut

.github/workflows/CompatHelper.yml

@@ -1,25 +1,26 @@
 name: CompatHelper
-
 on:
   schedule:
-    - cron: '00 * * * *'
-
+    - cron: 0 0 * * *
+  workflow_dispatch:
 jobs:
-  build:
-    runs-on: ${{ matrix.os }}
-    strategy:
-      matrix:
-        julia-version: [1.2.0]
-        julia-arch: [x86]
-        os: [ubuntu-latest]
+  CompatHelper:
+    runs-on: ubuntu-latest
     steps:
-      - uses: julia-actions/setup-julia@latest
-        with:
-          version: ${{ matrix.julia-version }}
-      - name: Install dependencies
-        run: julia -e 'using Pkg; Pkg.add(Pkg.PackageSpec(name = "CompatHelper", url = "https://github.com/bcbi/CompatHelper.jl.git"))'
-      - name: CompatHelper.main
+      - name: "Install CompatHelper"
+        run: |
+          import Pkg
+          name = "CompatHelper"
+          uuid = "aa819f21-2bde-4658-8897-bab36330d9b7"
+          version = "2"
+          Pkg.add(; name, uuid, version)
+        shell: julia --color=yes {0}
+      - name: "Run CompatHelper"
+        run: |
+          import CompatHelper
+          CompatHelper.main()
+        shell: julia --color=yes {0}
         env:
           GITHUB_TOKEN: ${{ secrets.GITHUB_TOKEN }}
-          JULIA_DEBUG: CompatHelper
-        run: julia -e 'using CompatHelper; CompatHelper.main()'
+          COMPATHELPER_PRIV: ${{ secrets.DOCUMENTER_KEY }}
+          # COMPATHELPER_PRIV: ${{ secrets.COMPATHELPER_PRIV }}

Project.toml

@@ -1,6 +1,6 @@
 name = "MultivariateOrthogonalPolynomials"
 uuid = "4f6956fd-4f93-5457-9149-7bfc4b2ce06d"
-version = "0.1.0"
+version = "0.1.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -33,7 +33,7 @@ ContinuumArrays = "0.7.3"
 DomainSets = "0.4"
 FastTransforms = "0.11, 0.12"
 FillArrays = "0.11"
-HarmonicOrthogonalPolynomials = "0.0.4, 0.0.5"
+HarmonicOrthogonalPolynomials = "0.0.4, 0.0.5, 0.1"
 InfiniteArrays = "0.10"
 InfiniteLinearAlgebra = "0.5"
 LazyArrays = "0.21"

examples/diskhelmholtz.jl

@@ -1,80 +1,18 @@
 using MultivariateOrthogonalPolynomials, FastTransforms, BlockBandedMatrices, Plots
 plotly()
-import MultivariateOrthogonalPolynomials: ZernikeITransform, grid, plotgrid
 
+Z = Zernike(1)
+W = Weighted(Z)
+xy = axes(Z,1); x,y = first.(xy),last.(xy)
+Δ = Z \ (Laplacian(xy) * W)
+S = Z \ W
+k = 2
+f = @.(cos(x*exp(y)))
+F = factorize(Δ + k^2 * S)
+c = (Z \ f)
+F \ c
 
-Z = Zernike()[:,Block.(Base.OneTo(10))]
-g = grid(Z)
-θ = getproperty.(g[1,:],:θ)
-[permutedims(RadialCoordinate.(1,θ)); g; permutedims(RadialCoordinate.(0,θ))]
-import Mul
+u = W * ((Δ + k^2 * S) \ (Z \ f))
 
-plotgrid(Z)
-
-
-Z = Zernike()
-xy = axes(Z,1)
-x,y = first.(xy),last.(xy)
-u = Z * (Z \ @.(cos(10x*y)))
-surface(u)
-
-@.(cos(10x*y))
-
-plot(u)
-
-g = plotgrid(Z[:,Block.(Base.OneTo(10))])
-surface(first.(g), last.(g), u[g])
-plot(u)
-
-
-G = grid(Z)
-
-contourf(first.(G), last.(G), ones(size(G)...))
-
-
-scatter(vec(first.(G)), vec(last.(G)))
-G[1,1].θ
-
-WZ = Weighted(Zernike(1))
-xy = axes(WZ,1)
-x,y = first.(xy),last.(xy)
-f = Zernike(1) \ @. exp(x*cos(y))
-
-N = 50
-KR = Block.(Base.OneTo(N))
-Δ = (Zernike(1) \ (Laplacian(xy) * WZ))[KR,KR]
-C = (Zernike(1) \ WZ)[KR,KR]
-k = 5
-L = Δ - k^2 * C
-
-v = f[KR]
-@time u = (L \ v);
-
-g = MultivariateOrthogonalPolynomials.grid(Zernike(1)[:,KR])
-U = ZernikeITransform{Float64}(N, 0, 1) * u
-
-plot(first.(g), last.(g), U)
-
-
-
-F = factorize(Zernike(1)[:,KR]).plan
-F \ u
-u
-
-
-F*u
-
-
-F.plan \ v
-F |> typeof |> fieldnames
-
-grid(WZ[:,KR])
-
-F |>typeof |> fieldnames
-F.F * v
-
-m = DiskTrav(v).matrix
-
-plan_disk2cxf(m, 0, 0) * m