Update test_rectdisk.jl

Author: dlfivefifty

test/test_rectdisk.jl

@@ -34,33 +34,34 @@ using ForwardDiff
         @test WP == WP
 
         x, y = coordinates(P)
+        L = WP \ WQ
+        R = Q \ P
+
+        ∂x = Derivative(P, (1,0))
+        ∂y = Derivative(P, (0,1))
+
+        Dx = Q \ (∂x * P)
+        Dy = Q \ (∂y * P)
+
+        X = P \ (x .* P)
+        Y = P \ (y .* P)
+
 
         @testset "lowering/raising" begin
-            L = WP \ WQ
             @test WP[SVector(0.1,0.2),Block.(1:6)]'L[Block.(1:6),Block.(1:4)] ≈ WQ[SVector(0.1,0.2),Block.(1:4)]'
-            R = Q \ P
             @test Q[SVector(0.1,0.2),Block.(1:4)]'R[Block.(1:4),Block.(1:4)] ≈ P[SVector(0.1,0.2),Block.(1:4)]'
 
             @test (DunklXuDisk() \ WeightedDunklXuDisk(1.0))[Block.(1:N), Block.(1:N)] ≈ (WeightedDunklXuDisk(0.0) \ WeightedDunklXuDisk(1.0))[Block.(1:N), Block.(1:N)]
         end
 
 
         @testset "jacobi" begin
-            X = P \ (x .* P)
-            Y = P \ (y .* P)
-
             @test (L * R)[Block.(1:N), Block.(1:N)] ≈ (I - X^2 - Y^2)[Block.(1:N), Block.(1:N)]
             @test P[SVector(0.1,0.2),Block.(1:5)]'X[Block.(1:5),Block.(1:4)] ≈ 0.1P[SVector(0.1,0.2),Block.(1:4)]'
             @test P[SVector(0.1,0.2),Block.(1:5)]'Y[Block.(1:5),Block.(1:4)] ≈ 0.2P[SVector(0.1,0.2),Block.(1:4)]'
         end
 
         @testset "derivatives" begin
-            ∂x = Derivative(P, (1,0))
-            ∂y = Derivative(P, (0,1))
-
-            Dx = Q \ (∂x * P)
-            Dy = Q \ (∂y * P)
-
             @test Q[SVector(0.1,0.2),Block.(1:3)]'Dx[Block.(1:3),Block.(1:4)] ≈ [ForwardDiff.gradient(𝐱 -> DunklXuDisk{eltype(𝐱)}(P.β)[𝐱,k], SVector(0.1,0.2))[1] for k=1:10]'
             Mx = Q \ (x .* Q)
             My = Q \ (y .* Q)