]> https://gitweb.dealii.org/ - dealii.git/commitdiff
scalapack: fix a minor resize bug in eigenpairs_symmetric() 5921/head
authorDenis Davydov <davydden@gmail.com>
Sun, 18 Feb 2018 21:37:00 +0000 (22:37 +0100)
committerDenis Davydov <davydden@gmail.com>
Mon, 19 Feb 2018 07:29:04 +0000 (08:29 +0100)
source/lac/scalapack.cc

index a3eef89667e4e338562b013eb9909a7756637491..49c5c51e6ae72b49a5a826d8854d17d8197a414e 100644 (file)
@@ -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,

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.