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

Fix hessian with fixed atoms #1175

Merged
merged 6 commits into from
Feb 12, 2025
Merged

Conversation

foxtran
Copy link
Contributor

@foxtran foxtran commented Feb 4, 2025

Fixes #475

This patch updates hessian.F90 to properly work with freezed atoms (and a little bit with fixed, but I did not get the difference between fix and freeze). With this patch, input (provided in #475) can be computed without segfaults (fixed in commit 53aaa53) and produced Hessian has zero values for fixed atoms and fixed-free atom crossings, so it looks fine for me :-)

Unfortunately, I did not implement all branches and therefore hessian with 1-2 fixed atoms can not be computed.

Copy link
Member

@marvinfriede marvinfriede left a comment

Choose a reason for hiding this comment

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

Confirmed. #475 works now.

@marvinfriede marvinfriede merged commit ca631a3 into grimme-lab:main Feb 12, 2025
17 checks passed
@foxtran foxtran deleted the fix/fix-hessian branch February 12, 2025 12:24
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Hessian calculation runs into segmentation fault when atoms are frozen
3 participants