Avoid fails in the tests from Inf or -Inf

rafavzqz · Jul 23, 2024 · e593bf5 · e593bf5
1 parent e349b12
commit e593bf5
Show file tree

Hide file tree

Showing 20 changed files with 21 additions and 1 deletion.
diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_Lshaped.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_Lshaped.m
@@ -52,6 +52,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 128)
 %! assert (space.ndof, 430)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_Lshaped_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_Lshaped_mp.m
@@ -88,6 +88,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = mp_solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 108)
 %! assert (space.ndof, 416)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_cube.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_cube.m
@@ -47,6 +47,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 64)
 %! assert (space.ndof, 540)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_cube_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_cube_mp.m
@@ -175,6 +175,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = mp_solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 108)
 %! assert (space.ndof, 923)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_curvedL.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_curvedL.m
@@ -51,6 +51,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 128)
 %! assert (space.ndof, 430)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_fichera_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_fichera_mp.m
@@ -46,6 +46,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = mp_solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 56)
 %! assert (space.ndof, 801)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_Lshaped.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_Lshaped.m
@@ -51,7 +51,8 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 128)
 %! assert (space.ndof, 430)
 %! assert (sp_mul.ndof, 231)
-%! assert (eigv(2:6), [1.47283540635722; 3.53399368976646; 9.86962050365007; 9.86962103171719; 11.38943951321592], 1e-13)
+%! assert (eigv(1:5), [1.47283540635722; 3.53399368976646; 9.86962050365007; 9.86962103171719; 11.38943951321592], 1e-13)
diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_Lshaped_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_Lshaped_mp.m
@@ -88,6 +88,7 @@
 %! [geometry, msh, space, sp_mul, eigv, eigf] = ...
 %!                   mp_solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 108)
 %! assert (space.ndof, 416)
 %! assert (sp_mul.ndof, 225)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_cube.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_cube.m
@@ -45,6 +45,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 64)
 %! assert (space.ndof, 540)
 %! assert (sp_mul.ndof, 216)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_cube_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_cube_mp.m
@@ -167,6 +167,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = mp_solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 108)
 %! assert (space.ndof, 923)
 %! assert (sp_mul.ndof, 360)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_fichera_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_fichera_mp.m
@@ -45,6 +45,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = mp_solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 56)
 %! assert (space.ndof, 801)
 %! assert (sp_mul.ndof, 316)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_square.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_square.m
@@ -49,6 +49,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 100)
 %! assert (space.ndof, 312)
 %! assert (sp_mul.ndof, 169)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_thick_L.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_thick_L.m
@@ -46,6 +46,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 54)
 %! assert (space.ndof, 560)
 %! assert (sp_mul.ndof, 225)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_thick_L_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_thick_L_mp.m
@@ -65,6 +65,7 @@
 %! [geometry, msh, space, sp_mul, eigv, eigf] = ...
 %!                   mp_solve_maxwell_eig_mixed1 (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 81)
 %! assert (space.ndof, 820)
 %! assert (sp_mul.ndof, 325)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed2_Lshaped.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed2_Lshaped.m
@@ -51,6 +51,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = solve_maxwell_eig_mixed2_2d (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! assert (msh.nel, 128)
 %! assert (space.ndof, 430)
 %! assert (sp_mul.ndof, 200)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed2_square.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed2_square.m
@@ -49,6 +49,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, sp_mul, eigv, eigf] = solve_maxwell_eig_mixed2_2d (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(eigv>0 & eigv<Inf)
 %! nzeros = numel (find (eigv < 1e-4));
 %! assert (msh.nel, 100)
 %! assert (space.ndof, 312)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_ring_1eighth.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_ring_1eighth.m
@@ -51,6 +51,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 64)
 %! assert (space.ndof, 220)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_square.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_square.m
@@ -51,6 +51,7 @@
 %! method_data.nquad      = [4 4];     % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 100)
 %! assert (space.ndof, 312)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_thick_L.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_thick_L.m
@@ -48,6 +48,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 54)
 %! assert (space.ndof, 560)

diff --git a/geopdes/inst/examples/maxwell/ex_maxwell_eig_thick_L_mp.m b/geopdes/inst/examples/maxwell/ex_maxwell_eig_thick_L_mp.m
@@ -65,6 +65,7 @@
 %! method_data.nquad      = [3 3 3]; % Points for the Gaussian quadrature rule
 %! [geometry, msh, space, eigv, eigf] = mp_solve_maxwell_eig (problem_data, method_data);
 %! [eigv, perm] = sort (eigv);
+%! eigv = eigv(abs(eigv)<Inf);
 %! nzeros = numel (find (eigv < 1e-10));
 %! assert (msh.nel, 81)
 %! assert (space.ndof, 820)