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

Work around numerical instability of the mio N example #93

Open
wants to merge 1 commit into
base: main
Choose a base branch
from

Conversation

vanderhe
Copy link
Member

@vanderhe vanderhe commented Aug 8, 2024

This slight change is enough to work around the numerical instabilities encountered for the nitrogen STO basis of the mio example. In fact, it is even closer to the geometric series (0.5000 1.2051 2.9044 7.0000) suggested by skexp.

Difference to the original slateratom binary (+original exponents), as reported by skdiff:

*** Atomic data:
Onsite:         2.739e-08     1
Hubbards:       5.742e-07     1
Hubbard (s):    5.742e-07
Occupations:    0.000e+00     0

*** Integral tables:
Hamiltonian:    1.820e-06 ( 104,  9)
Overlap:        4.158e-06 ( 190,  9)

Difference to the original slateratom binary (+updated exponents), as reported by skdiff:

*** Atomic data:
Onsite:         8.272e-09     1
Hubbards:       3.631e-07     1
Hubbard (s):    3.631e-07
Occupations:    0.000e+00     0

*** Integral tables:
Hamiltonian:    3.706e-07 (  19,  5)
Overlap:        2.660e-11 (  45,  6)

Closes #92 (not really a fix though).

@vanderhe vanderhe requested a review from aradi August 8, 2024 14:59
Copy link
Member

@bhourahine bhourahine left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Seems OK, but would be better to switch to something algorithmically more robust.

@vanderhe
Copy link
Member Author

Seems OK, but would be better to switch to something algorithmically more robust.

Yeah, this PR does not really qualify as a fix. Strictly spoken is not even the mio nitrogen anymore, but what can you do...

@vanderhe vanderhe force-pushed the fixMioNitrogenExample branch from 0633566 to 1a62f7d Compare November 19, 2024 09:25
@vanderhe vanderhe added the documentation Improvements or additions to documentation label Dec 9, 2024
@bhourahine
Copy link
Member

@vanderhe did you get any success in trying to make the change numerically smaller?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
documentation Improvements or additions to documentation
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Problematic nitrogen basis in examples
2 participants