Add tests for Meshes.Quadrangle (#95)

* add tests for Meshes.Quadrangle

* add Quadrangle to support matrix in docs

Author: JoshuaLampert

docs/src/supportmatrix.md

@@ -39,6 +39,7 @@ Cubature integration rules are recommended (and the default).
 | `Line` | ✅ | ✅ | ✅ |
 | `ParaboloidSurface` | ✅ | ✅ | ✅ |
 | `Plane` | ✅ | ✅ | ✅ |
+| `Quadrangle` | ✅ | ✅ | ✅ |
 | `Ray` | ✅ | ✅ | ✅ |
 | `Ring` | ✅ | ✅ | ✅ |
 | `Rope` | ✅ | ✅ | ✅ |

test/Project.toml

@@ -3,6 +3,7 @@ Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
 ExplicitImports = "7d51a73a-1435-4ff3-83d9-f097790105c7"
 Meshes = "eacbb407-ea5a-433e-ab97-5258b1ca43fa"
 QuadGK = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
+SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 TestItemRunner = "f8b46487-2199-4994-9208-9a1283c18c0a"
 TestItems = "1c621080-faea-4a02-84b6-bbd5e436b8fe"
@@ -13,6 +14,7 @@ Aqua = "0.7, 0.8"
 ExplicitImports = "1.6.0"
 Meshes = "0.50, 0.51"
 QuadGK = "2.1.1"
+SpecialFunctions = "2"
 TestItemRunner = "1"
 TestItems = "1"
 Unitful = "1.19"

test/combinations.jl

@@ -203,6 +203,35 @@ end
     @test_throws "not supported" volumeintegral(f, plane)
 end
 
+@testitem "Meshes.Quadrangle" setup=[Setup] begin
+    using SpecialFunctions: erf
+    quadrangle = Quadrangle((-1.0, 0.0), (-1.0, 1.0), (1.0, 1.0), (1.0, 0.0))
+
+    function f(p::P) where {P <: Meshes.Point}
+        ur = hypot(p.coords.x, p.coords.y)
+        r = ustrip(u"m", ur)
+        exp(-r^2)
+    end
+    fv(p) = fill(f(p), 3)
+
+    # Scalar integrand
+    sol = 0.5 * π * erf(1)^2 * u"m^2"
+    @test integral(f, quadrangle, GaussLegendre(100)) ≈ sol
+    @test integral(f, quadrangle, GaussKronrod()) ≈ sol
+    @test integral(f, quadrangle, HAdaptiveCubature()) ≈ sol
+
+    # Vector integrand
+    vsol = fill(sol, 3)
+    @test integral(fv, quadrangle, GaussLegendre(100)) ≈ vsol
+    @test integral(fv, quadrangle, GaussKronrod()) ≈ vsol
+    @test integral(fv, quadrangle, HAdaptiveCubature()) ≈ vsol
+
+    # Integral aliases
+    @test_throws "not supported" lineintegral(f, quadrangle)
+    @test surfaceintegral(f, quadrangle) ≈ sol
+    @test_throws "not supported" volumeintegral(f, quadrangle)
+end
+
 @testitem "Meshes.Ray" setup=[Setup] begin
     a = Point(0.0u"m", 0.0u"m", 0.0u"m")
     v = Vec(1.0u"m", 1.0u"m", 1.0u"m")