adjust LAPACKFullMatrix::print_formatted() to output lower triangular part for Cholesky

diff --git a/source/lac/lapack_full_matrix.cc b/source/lac/lapack_full_matrix.cc
index bed672572a4cf73e5207c901afcc07a2a146c093..3e8246ce8a9a7f9b48491a224860c7b5928a50dc 100644 (file)
--- a/source/lac/lapack_full_matrix.cc
+++ b/source/lac/lapack_full_matrix.cc
@@ -1440,7 +1440,9 @@ LAPACKFullMatrix<number>::print_formatted (
 
   for (size_type i=0; i<this->n_rows(); ++i)
     {
-      for (size_type j=0; j<this->n_cols(); ++j)
+      // Cholesky is stored in lower triangular, so just output this part:
+      const size_type nc = state == LAPACKSupport::cholesky ? i+1 : this->n_cols();
+      for (size_type j=0; j<nc; ++j)
         // we might have complex numbers, so use abs also to check for nan
         // since there is no isnan on complex numbers
         if (std::isnan(std::abs((*this)(i,j))))