Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

--cosmo infinity calculation fails for HS- #1086

Open
TyBalduf opened this issue Aug 12, 2024 · 0 comments
Open

--cosmo infinity calculation fails for HS- #1086

TyBalduf opened this issue Aug 12, 2024 · 0 comments
Labels
unconfirmed This report has not yet been confirmed by the developers

Comments

@TyBalduf
Copy link
Contributor

Describe the bug
Running HS- with --cosmo infinity (or --cpcmx <any_solvent> since this will have do a cosmo infinity solvent run first), I get this output in the SCC.

 iter      E             dE          RMSdq      gap      omega  full diag
   1     -4.0645745 -0.406457E+01  0.682E+00    5.96       0.0  T
   2     -4.0657628 -0.118828E-02  0.374E+00    5.60       1.0  T
Upps!

I found the error message here in the code and looking at the subroutine gjinv, it seems like can only fail if the matrix its working with is singular. This is based on the (commented out) warning messages in the subroutine.

There isn't likely much that can be done about a singular matrix (maybe some kind of pseudo-inverse?), but it would be useful to enable these warning messages to clarify what has gone wrong with the calculation.

Additional context
Input: hs_debug_sdf.txt
Output: new_hs.txt

@TyBalduf TyBalduf added the unconfirmed This report has not yet been confirmed by the developers label Aug 12, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
unconfirmed This report has not yet been confirmed by the developers
Projects
None yet
Development

No branches or pull requests

1 participant