Page 1 of 1

envmc matrix inversion not converged

Posted: Tue Jul 11, 2017 2:59 pm
by Katharina Doblhoff
Dear Casino users!
I have a system for which I constantly get the message
Warning: [EMIN_MATRIX_ALGEBRA] Could not invert S (stage 1).
when performing an emin optimization of the Jastrow.

Looking at the code, this error is thrown if the divide and conquer algorithm in the singular value decomposition of the S matrix (I think it is eqn 226 in the manual) does not converge or if the input values are forbidden.

Does anybody have a good idea of what to DO against this (except for exchanging gesdd by gesdv in the code, which would - I guess - be slower but maybe more stable)?

Thanks,
Katharina

Re: envmc matrix inversion not converged

Posted: Fri Aug 11, 2017 8:31 am
by Katharina Doblhoff
I solved the issue: It seems to be a problem with early lapack versions...
So in case anybody else stumbles into this issue: Try using newer libraries!