From 6d4d051833013b54de0ca97d3839289d3608e455 Mon Sep 17 00:00:00 2001
From: Denis Davydov <davydden@gmail.com>
Date: Mon, 18 Dec 2017 15:19:34 +0100
Subject: [PATCH] adjust LAPACKFullMatrix::print_formatted() to output lower
 triangular part for Cholesky

---
 source/lac/lapack_full_matrix.cc | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/source/lac/lapack_full_matrix.cc b/source/lac/lapack_full_matrix.cc
index bed672572a..3e8246ce8a 100644
--- 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))))
-- 
2.39.5