-
Notifications
You must be signed in to change notification settings - Fork 14
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 cutoff
in LODE TabulatedRadialIntegral
#250
Conversation
Here is a pre-built version of the code in this pull request: wheels.zip, you can install it locally by unzipping |
@@ -24,7 +24,7 @@ pub struct LodeRadialIntegralSplineParameters { | |||
pub max_radial: usize, | |||
/// Number of angular components | |||
pub max_angular: usize, | |||
/// cutoff radius, this is also the maximal value that can be interpolated | |||
/// k-space cutoff radius, this is also the maximal value that can be interpolated | |||
pub cutoff: f64, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Maybe we should call this also this k_cutoff
to avoid confusion with the real space version?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
yes, that sounds like a good idea!
13eeaed
to
2e4ee9a
Compare
}; | ||
|
||
let gto = LodeRadialIntegralGto::new(LodeRadialIntegralGtoParameters { | ||
max_radial: parameters.max_radial, | ||
max_angular: parameters.max_angular, | ||
cutoff: parameters.cutoff, | ||
cutoff: 5.0, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I changed this value here to a different one as k_cutoff
. Otherwise one might be confused why the real and the k-space cutoff has to be same.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
👍
Fixes #249
The
cutoff
for the LODETabulatedRadialIntegral
was set to real space environment cutoff. This is wrong since for LODE we spline in k-space from 0 until thek_cutoff
.📚 Documentation preview 📚: https://rascaline--250.org.readthedocs.build/en/250/