scnc
/
phono3py
mirror of https://github.com/phonopy/phono3py.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Minor code clean-up

This commit is contained in:
Atsushi Togo 2025-06-14 18:15:03 +09:00
parent 26cac3aef5
commit 55c334e6c2
1 changed files with 10 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
c/reciprocal_to_normal.c
View File

@ -207,21 +207,19 @@ static void get_fc3_e0(lapack_complex_double *fc3_e0,
const lapack_complex_double *fc3_reciprocal, const lapack_complex_double *fc3_reciprocal,
const lapack_complex_double *e0, const lapack_complex_double *e0,
const int64_t band_index_0, const int64_t num_band) { const int64_t band_index_0, const int64_t num_band) {
int64_t j, k, ll; int64_t j, k;
lapack_complex_double fc3_elem; lapack_complex_double fc3_elem;
for (j = 0; j < num_band; j++) { for (j = 0; j < num_band; j++) {
for (k = 0; k < num_band; k++) { for (k = 0; k < num_band * num_band; k++) {
for (ll = 0; ll < num_band; ll++) { fc3_elem =
fc3_elem = phonoc_complex_prod( phonoc_complex_prod(fc3_reciprocal[j * num_band * num_band + k],
fc3_reciprocal[j * num_band * num_band + k * num_band + ll], e0[band_index_0 * num_band + j]);
e0[band_index_0 * num_band + j]); fc3_e0[k] = lapack_make_complex_double(
fc3_e0[k * num_band + ll] = lapack_make_complex_double( lapack_complex_double_real(fc3_e0[k]) +
lapack_complex_double_real(fc3_e0[k * num_band + ll]) + lapack_complex_double_real(fc3_elem),
lapack_complex_double_real(fc3_elem), lapack_complex_double_imag(fc3_e0[k]) +
lapack_complex_double_imag(fc3_e0[k * num_band + ll]) + lapack_complex_double_imag(fc3_elem));
lapack_complex_double_imag(fc3_elem));
}
} }
} }
} }