Merge branch 'tests' of github.com:JuliaGeometry/MeshIntegrals.jl int…

…o tests
JuliaGeometry · Dec 5, 2024 · 4b9157e · 4b9157e
2 parents 43edac0 + 452278b
commit 4b9157e
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/test/combinations.jl b/test/combinations.jl
@@ -146,15 +146,16 @@ end
         ux = ustrip(p.coords.x)
         (1 / sqrt(1 + cos(ux)^2)) * u"Ω"
     end
-    solution = 2π * u"Ω*m"
+    fv(p) = fill(f(p), 3)
+    sol = 2π * u"Ω*m"
+    vsol = fill(sol, 3)
 
     # Scalar integrand
     @test integral(f, curve, GaussLegendre(100))≈sol rtol=0.5e-2
     @test integral(f, curve, GaussKronrod())≈sol rtol=0.5e-2
     @test integral(f, curve, HAdaptiveCubature())≈sol rtol=0.5e-2
 
     # Vector integrand
-    vsol = fill(sol, 3)
     @test integral(fv, curve, GaussLegendre(100))≈vsol rtol=0.5e-2
     @test integral(fv, curve, GaussKronrod())≈vsol rtol=0.5e-2
     @test integral(fv, curve, HAdaptiveCubature())≈vsol rtol=0.5e-2
@@ -179,8 +180,8 @@ end
         sqrt(a^2 - t^2) * u"Ω"
     end
     fv(p) = fill(f(p), 3)
-    solution = π * a^2 / 4 * u"Ω*m"
-    vsol = fill(solution, 3)
+    sol = π * a^2 / 4 * u"Ω*m"
+    vsol = fill(sol, 3)
 
     # Scalar integrand
     @test integral(f, box, GaussLegendre(100))≈sol rtol=1e-6