From 858ea62a0ad16ec3da33d536a0b4cf9bc49b6c55 Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Thu, 26 Jan 2023 18:19:24 +0100 Subject: [PATCH] Shadowing variables --- org/qmckl_jastrow.org | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/org/qmckl_jastrow.org b/org/qmckl_jastrow.org index 039756d..34f201d 100644 --- a/org/qmckl_jastrow.org +++ b/org/qmckl_jastrow.org @@ -5023,13 +5023,13 @@ qmckl_exit_code qmckl_compute_een_rescaled_e_hpc ( een_rescaled_e_ij[kk]= ( kk < (elec_pairs) ? 1.0 : 0.0 ); } - int64_t k = 0; + int64_t kk = 0; for (int i = 0; i < elec_num; ++i) { for (int j = 0; j < i; ++j) { - // een_rescaled_e_ij(k, 2) = dexp(-rescale_factor_ee * ee_distance(i, j, nw)); - een_rescaled_e_ij[k + elec_pairs] = exp(-rescale_factor_ee * \ + // een_rescaled_e_ij(kk, 2) = dexp(-rescale_factor_ee * ee_distance(i, j, nw)); + een_rescaled_e_ij[kk + elec_pairs] = exp(-rescale_factor_ee * \ ee_distance[j + i*elec_num + nw*(elec_num*elec_num)]); - k = k + 1; + kk += 1; } } @@ -5052,7 +5052,7 @@ qmckl_exit_code qmckl_compute_een_rescaled_e_hpc ( // Up to here it should work. for ( int l = 1; l < (cord_num+1); ++l) { - k = 0; + int k = 0; for (int i = 0; i < elec_num; ++i) { for (int j = 0; j < i; ++j) { x = een_rescaled_e_ij[k + l*elec_pairs];