From 18b9f4a0d0724a23daf113f62aedebc7c8fe5d95 Mon Sep 17 00:00:00 2001 From: Jukka Aho Date: Sat, 9 May 2026 18:29:39 +0300 Subject: [PATCH] chore(test): delete wedge quadrature regression Remove wedge-specific quadrature tests pending new coverage. - Drop `test/quadrature/test_wedges.jl`. --- test/quadrature/test_wedges.jl | 54 ---------------------------------- 1 file changed, 54 deletions(-) delete mode 100644 test/quadrature/test_wedges.jl diff --git a/test/quadrature/test_wedges.jl b/test/quadrature/test_wedges.jl deleted file mode 100644 index ac056d3..0000000 --- a/test/quadrature/test_wedges.jl +++ /dev/null @@ -1,54 +0,0 @@ -# This file is a part of JuliaFEM. -# License is MIT: see https://github.com/JuliaFEM/JuliaFEM.jl/blob/master/LICENSE.md - -@testset "Wedge Quadrature Rules" begin - - @testset "GaussLegendre{2} - 6 points (default)" begin - points = get_quadrature_points(Wedge, GaussLegendre{2}()) - @test length(points) == 6 - - # Check weights sum to 1 (volume of reference wedge) - weight_sum = sum(p.weight for p in points) - @test weight_sum ≈ 1.0 - - # All weights should be equal - @test all(p.weight ≈ 1/6 for p in points) - - # Check z coordinates are at ±1/sqrt(3) - a = 1/sqrt(3) - @test sum(abs(p.coords[3]) ≈ a for p in points) == 6 - end - - @testset "GaussLegendre{2,:B} - 6 points (variant B)" begin - points = get_quadrature_points(Wedge, GaussLegendre{2,:B}()) - @test length(points) == 6 - - # Check weights sum to 1 - weight_sum = sum(p.weight for p in points) - @test weight_sum ≈ 1.0 - - # All weights should be equal - @test all(p.weight ≈ 1/6 for p in points) - end - - @testset "GaussLegendre{5} - 21 points" begin - points = get_quadrature_points(Wedge, GaussLegendre{5}()) - @test length(points) == 21 - - # Check weights sum to 1 - weight_sum = sum(p.weight for p in points) - @test weight_sum ≈ 1.0 rtol=1e-10 - - # Should have 3 layers (z = -alpha, 0, alpha) - a = sqrt(3/5) - z_values = [p.coords[3] for p in points] - @test sum(abs(z) ≈ a for z in z_values) == 14 # 7 points at each ±alpha - @test sum(abs(z) ≈ 0.0 for z in z_values) == 7 # 7 points at z=0 - end - - @testset "Point coordinates are Vec type" begin - points = get_quadrature_points(Wedge, GaussLegendre{2}()) - @test points[1].coords isa Vec{3,Float64} - @test !(points[1].coords isa SVector) - end -end