Test for bivariate points (#312)

Author: jishnub

Project.toml

@@ -30,7 +30,7 @@ ApproxFunOrthogonalPolynomialsPolynomialsExt = "Polynomials"
 ApproxFunOrthogonalPolynomialsStaticExt = "Static"
 
 [compat]
-ApproxFunBase = "0.8.58, 0.9.15"
+ApproxFunBase = "0.8.60, 0.9.16"
 ApproxFunBaseTest = "0.1"
 Aqua = "0.7"
 BandedMatrices = "0.16, 0.17"

test/MultivariateTest.jl

@@ -72,6 +72,12 @@ include("testutils.jl")
         end
     end
 
+    @testset "points" begin
+        S = Chebyshev(1..2) ⊗ Chebyshev(3..4)
+        @test all(p ∈ domain(factor(S, 1)) for p in points(S, 10, 10)[1])
+        @test all(p ∈ domain(factor(S, 2)) for p in points(S, 10, 10)[2])
+    end
+
     @testset "Evaluation" begin
 
         @testset "2D" begin