From 89ec220f4893d4edf4b04d88051095d330761ebc Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Fri, 8 Jul 2022 17:01:09 +0200
Subject: [PATCH] davidson assumend converged if zero krylov vector encountered

---
 davidson.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/davidson.h b/davidson.h
index e340545..9d3ad95 100644
--- a/davidson.h
+++ b/davidson.h
@@ -215,7 +215,7 @@ for(i=0; i<n; ++i)
 
 //orthogonalise to previous vectors
 typename LA_traits<T>::normtype vnorm= vec1.norm();
-if(vnorm==0.) laerror("Zero Krylov vector in Davidson - perhaps try different initial guess");
+if(vnorm==0.) goto converged; //Zero Krylov vector - for tiny matrices probably means converged
 else
 {
 vec1 *= (1./vnorm);