From: Denis Davydov Date: Sun, 18 Feb 2018 21:37:00 +0000 (+0100) Subject: scalapack: fix a minor resize bug in eigenpairs_symmetric() X-Git-Tag: v9.0.0-rc1~394^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b85125b39ca64abb1e286a6005bc3c095a22cca7;p=dealii.git scalapack: fix a minor resize bug in eigenpairs_symmetric() --- diff --git a/source/lac/scalapack.cc b/source/lac/scalapack.cc index a3eef89667..49c5c51e6a 100644 --- a/source/lac/scalapack.cc +++ b/source/lac/scalapack.cc @@ -621,8 +621,8 @@ ScaLAPACKMatrix::eigenpairs_symmetric(const bool compute_eigenvector plamch( &(this->grid->blacs_context), &cmach, abstol); abstol *= 2; ifail.resize(n_rows); - iclustr.resize(n_local_rows * n_local_columns); - gap.resize(n_local_rows * n_local_columns); + iclustr.resize(2 * grid->n_process_rows * grid->n_process_columns); + gap.resize(grid->n_process_rows * grid->n_process_columns); psyevx(&jobz, &range, &uplo, &n_rows, A_loc, &submatrix_row, &submatrix_column, descriptor, &vl, &vu, &il, &iu, &abstol, &m, &nz, &ev[0], &orfac,