davidson assumend converged if zero krylov vector encountered

2022-07-08 17:01:09 +02:00 · 2022-07-08 17:01:09 +02:00 · 89ec220f48
commit 89ec220f48
parent 3f0fd3f8b0
1 changed files with 1 additions and 1 deletions
--- 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);