From 50a5aa11a6670f80c8fd3354eb91b991a3844c97 Mon Sep 17 00:00:00 2001
From: Kristian Arjas <kristian.arjas@protonmail.com>
Date: Wed, 15 Sep 2021 13:40:20 +0300
Subject: [PATCH] Add residual to solver (fix regression from 087cd1c)

---
 qpms/beyn.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/qpms/beyn.c b/qpms/beyn.c
index 8d98ff2..b5b8177 100644
--- a/qpms/beyn.c
+++ b/qpms/beyn.c
@@ -485,10 +485,10 @@ static int beyn_process_matrices(BeynSolver *solver, beyn_function_M_t M_functio
     if (res_tol > 0 && residual > res_tol) continue;
 
     eigenvalues[KRetained] = z;
-    if(eigenvectors) {
+    if(eigenvectors) { // eigenvectors is NULL when calculating the "coarse" contour for error estimates
       for(int j = 0; j < m; ++j) 
         MAT(eigenvectors, l, m, KRetained, j) = MAT(V0S, m, K, j, k);
-
+      solver->residuals[KRetained] = residual;
     }
     ++KRetained;
   }