Skip to content
Snippets Groups Projects
Commit 1a3fa0f0 authored by Mikhail Svechnikov's avatar Mikhail Svechnikov
Browse files

cleanup

parent 94ec0ac8
No related branches found
No related tags found
1 merge request!2358Optimize exponents summation in "exact" mesocrystal FF algorithm (#892)
Pipeline #129753 passed
...@@ -109,12 +109,11 @@ complex_t ReMesocrystal::structureFactor(const WavevectorInfo& wavevectors) cons ...@@ -109,12 +109,11 @@ complex_t ReMesocrystal::structureFactor(const WavevectorInfo& wavevectors) cons
} }
} else { } else {
std::vector<complex_t> exponents_c(range.z() + 1); std::vector<complex_t> exponents_c(range.z() + 1);
for (int k = 0; k <= range.z(); k++)
exponents_c[k] = pow(exp_c, k + min.z());
std::vector<complex_t> exp_c_sum_factor(range.z() + 1); std::vector<complex_t> exp_c_sum_factor(range.z() + 1);
for (int k = 0; k <= range.z(); k++) for (int k = 0; k <= range.z(); k++) {
exponents_c[k] = pow(exp_c, k + min.z());
exp_c_sum_factor[k] = (pow(exp_c, k + 1) - 1.) / (exp_c - 1.); exp_c_sum_factor[k] = (pow(exp_c, k + 1) - 1.) / (exp_c - 1.);
}
for (int i = 0; i <= range.x(); i++) { for (int i = 0; i <= range.x(); i++) {
const complex_t e_a = pow(exp_a, i + min.x()); const complex_t e_a = pow(exp_a, i + min.x());
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment