From: Denis Davydov <davydden@gmail.com>
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<NumberType>::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,