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

Bugfix: Isotope and boundary scattering terms #208

Merged
merged 11 commits into from
Jan 3, 2024
Merged

Conversation

jcoulter12
Copy link
Collaborator

@jcoulter12 jcoulter12 commented Jan 2, 2024

Summary

Bugfix of two scattering rates:

  • There was an error in the calculation of the off-diagonal isotope scattering rates when symmetries without symmetries.
  • Boundary scattering was previously overestimated due to a missing square root factor.

This is an important bug fix for those running phonon relaxon or variational BTE calculations (where the full matrix is needed and symmetry cannot be used) and also for those using the boundary scattering feature.

Other minor changes:

  • The test suite has been updated to reflect these fixes.
  • Unrelated, the symmetrizeMatrix variable again has been set to true, as this should have been amended in PR Switch manual tranpose of distributed matrix for ScaLAPACK pdtran #201.
  • Minor fix to handle elph matrix element files of up to 200 GB. May also handle larger ones, however, this has not yet been tested.

Major thanks goes to Michele Simoncelli (@MSimoncelli) for his help in the debugging process.

@jcoulter12 jcoulter12 added the bug Something isn't working label Jan 2, 2024
@jcoulter12 jcoulter12 merged commit 08f5b78 into develop Jan 3, 2024
32 checks passed
@jcoulter12 jcoulter12 deleted the isotopeFix branch January 3, 2024 03:09
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant