Allows Lower|Upper as a template argument of CG and MINRES: in this case the full matrix will be considered.

1 files changed, 4 insertions, 1 deletions

diff --git a/test/sparse_solver.h b/test/sparse_solver.h
index 887ff02a9..6cf99e0b0 100644
--- a/test/sparse_solver.h
+++ b/test/sparse_solver.h
@@ -195,7 +195,10 @@ int generate_sparse_spd_problem(Solver& , typename Solver::MatrixType& A, typena
   dA = dM * dM.adjoint();
   
   halfA.resize(size,size);
-  halfA.template selfadjointView<Solver::UpLo>().rankUpdate(M);
+  if(Solver::UpLo==(Lower|Upper))
+    halfA = A;
+  else
+    halfA.template selfadjointView<Solver::UpLo>().rankUpdate(M);
   
   return size;
 }