ArpackSolver::AdditionalData additional_data(
num_arnoldi_vectors, ArpackSolver::largest_magnitude, true);
ArpackSolver eigensolver(solver_control, additional_data);
- eigensolver.solve(stiffness_matrix,
- mass_matrix,
- inverse,
- eigenvalues,
- eigenfunctions,
- eigenvalues.size());
+ check_solver_within_range(eigensolver.solve(stiffness_matrix,
+ mass_matrix,
+ inverse,
+ eigenvalues,
+ eigenfunctions,
+ eigenvalues.size()),
+ solver_control.last_step(),
+ 2,
+ 10);
// make sure that we have eigenvectors and they are mass-orthonormal:
// a) (A*x_i-\lambda*B*x_i).L2() == 0
-DEAL::Convergence step 8 value 0.00000
+DEAL::Solver stopped within 2 - 10 iterations
DEAL:: Eigenvalue 0 : (4.93877,0.00000)
DEAL:: Eigenvalue 1 : (12.3707,0.00000)
DEAL:: Eigenvalue 2 : (12.3707,0.00000)
ArpackSolver::AdditionalData additional_data(
num_arnoldi_vectors, ArpackSolver::largest_magnitude, true);
ArpackSolver eigensolver(solver_control, additional_data);
- eigensolver.solve(stiffness_matrix,
- mass_matrix,
- inverse,
- eigenvalues,
- eigenfunctions,
- eigenvalues.size());
+ check_solver_within_range(eigensolver.solve(stiffness_matrix,
+ mass_matrix,
+ inverse,
+ eigenvalues,
+ eigenfunctions,
+ eigenvalues.size()),
+ solver_control.last_step(),
+ 2,
+ 10);
// make sure that we have eigenvectors and they are mass-orthonormal:
// a) (A*x_i-\lambda*B*x_i).L2() == 0
-DEAL::Convergence step 8 value 0.00000
+DEAL::Solver stopped within 2 - 10 iterations
DEAL::8 iterations used
DEAL:: Eigenvalue 0 : (4.93877,0.00000)
DEAL:: Eigenvalue 1 : (12.3707,0.00000)