Skip to content

Commit

Permalink
Fixed test, and removed ;
Browse files Browse the repository at this point in the history
  • Loading branch information
rafavzqz committed Jul 23, 2024
1 parent e593bf5 commit dc291f6
Show file tree
Hide file tree
Showing 10 changed files with 11 additions and 11 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -51,7 +51,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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 128)
%! assert (space.ndof, 430)
%! assert (sp_mul.ndof, 231)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 108)
%! assert (space.ndof, 416)
%! assert (sp_mul.ndof, 225)
Expand Down
2 changes: 1 addition & 1 deletion geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed1_cube.m
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 64)
%! assert (space.ndof, 540)
%! assert (sp_mul.ndof, 216)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -167,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 108)
%! assert (space.ndof, 923)
%! assert (sp_mul.ndof, 360)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 56)
%! assert (space.ndof, 801)
%! assert (sp_mul.ndof, 316)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 100)
%! assert (space.ndof, 312)
%! assert (sp_mul.ndof, 169)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -46,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 54)
%! assert (space.ndof, 560)
%! assert (sp_mul.ndof, 225)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -65,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 81)
%! assert (space.ndof, 820)
%! assert (sp_mul.ndof, 325)
Expand Down
4 changes: 2 additions & 2 deletions geopdes/inst/examples/maxwell/ex_maxwell_eig_mixed2_Lshaped.m
Original file line number Diff line number Diff line change
Expand Up @@ -51,8 +51,8 @@
%! 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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! assert (msh.nel, 128)
%! assert (space.ndof, 430)
%! assert (sp_mul.ndof, 200)
%! 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)
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +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)
%! eigv = eigv(eigv>0 & eigv<Inf);
%! nzeros = numel (find (eigv < 1e-4));
%! assert (msh.nel, 100)
%! assert (space.ndof, 312)
Expand Down

0 comments on commit dc291f6

Please sign in to comment.