Test Angular momentum operator (#83)

* Test Angular momentum operator

* Update Project.toml

* Update Project.toml

* Update Project.toml

Author: dlfivefifty

Project.toml

@@ -24,23 +24,23 @@ StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
-ArrayLayouts = "0.6, 0.7"
+ArrayLayouts = "0.7"
 BandedMatrices = "0.16"
-BlockArrays = "0.14.1, 0.15"
+BlockArrays = "0.15"
 BlockBandedMatrices = "0.10.1"
-ClassicalOrthogonalPolynomials = "0.3.2"
-ContinuumArrays = "0.7.3"
+ClassicalOrthogonalPolynomials = "0.4"
+ContinuumArrays = "0.8"
 DomainSets = "0.4, 0.5"
-FastTransforms = "0.11, 0.12"
+FastTransforms = "0.12"
 FillArrays = "0.11"
-HarmonicOrthogonalPolynomials = "0.0.4, 0.0.5, 0.1"
+HarmonicOrthogonalPolynomials = "0.1"
 InfiniteArrays = "0.10"
 InfiniteLinearAlgebra = "0.5"
-LazyArrays = "0.21"
-LazyBandedMatrices = "0.5"
-QuasiArrays = "0.4, 0.5"
-SpecialFunctions = "0.10, 1"
-StaticArrays = "0.12, 1"
+LazyArrays = "0.21.5"
+LazyBandedMatrices = "0.5.8"
+QuasiArrays = "0.6"
+SpecialFunctions = "1"
+StaticArrays = "1"
 julia = "1.6"
 
 [extras]

test/test_triangle.jl

@@ -317,4 +317,20 @@ import MultivariateOrthogonalPolynomials: tri_forwardrecurrence, grid, TriangleR
             end
         end
     end
+
+    @testset "AngularMomentum" begin
+        P = JacobiTriangle()
+        P¹ = JacobiTriangle(1,1,1)
+        xy = axes(P,1)
+        x,y = first.(xy),last.(xy)
+        ∂ˣ = PartialDerivative{1}(xy)
+        ∂ʸ = PartialDerivative{2}(xy)
+        L1 = x .* ∂ʸ
+        L2 = y .* ∂ˣ
+        L = x .* ∂ʸ - y .* ∂ˣ
+        A = P¹ \ (L1 * P)
+        B = P¹ \ (L2 * P)
+        C = P¹ \ (L * P)
+        @test C[1:10,1:10] ≈ A[1:10,1:10] - B[1:10,1:10]
+    end
 end