From: Denis Davydov Date: Wed, 15 Jun 2016 09:21:52 +0000 (+0200) Subject: parpack: change two Assert so that LinearOperator can be used for mass_matrix X-Git-Tag: v8.5.0-rc1~971^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F2683%2Fhead;p=dealii.git parpack: change two Assert so that LinearOperator can be used for mass_matrix --- diff --git a/include/deal.II/lac/parpack_solver.h b/include/deal.II/lac/parpack_solver.h index 2e5f0058b6..de67f19d17 100644 --- a/include/deal.II/lac/parpack_solver.h +++ b/include/deal.II/lac/parpack_solver.h @@ -610,12 +610,14 @@ void PArpackSolver::solve PArpackExcInvalidEigenvalueSize(n_eigenvalues, eigenvalues.size())); - Assert (n_eigenvalues < mass_matrix.m(), - PArpackExcInvalidNumberofEigenvalues(n_eigenvalues, mass_matrix.m())); + // use eigenvectors to get the problem size so that it is possible to + // employ LinearOperator for mass_matrix. + Assert (n_eigenvalues < eigenvectors[0].size(), + PArpackExcInvalidNumberofEigenvalues(n_eigenvalues, eigenvectors[0].size())); - Assert (additional_data.number_of_arnoldi_vectors < mass_matrix.m(), + Assert (additional_data.number_of_arnoldi_vectors < eigenvectors[0].size(), PArpackExcInvalidNumberofArnoldiVectors( - additional_data.number_of_arnoldi_vectors, mass_matrix.m())); + additional_data.number_of_arnoldi_vectors, eigenvectors[0].size())); Assert (additional_data.number_of_arnoldi_vectors > 2*n_eigenvalues+1, PArpackExcSmallNumberofArnoldiVectors(