TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2024-07-18 00:43:51 +02:00
Code Issues Projects Releases Wiki Activity

still failing

Browse Source
This commit is contained in:
Gianfranco Abrusci 2022-02-02 18:21:40 +01:00
parent f6f346d511
commit 67f80532f9
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
org/qmckl_jastrow.org
View File

@ -1340,7 +1340,7 @@ qmckl_exit_code qmckl_compute_asymp_jasb (
const double rescale_factor_kappa_ee, const double rescale_factor_kappa_ee,
double* const asymp_jasb ) { double* const asymp_jasb ) {
// Put some code here // Put some code here
int64_t i, p, info; int64_t info;
double kappa_inv, x, asym_one; double kappa_inv, x, asym_one;
kappa_inv = 1.0/ rescale_factor_kappa_ee; kappa_inv = 1.0/ rescale_factor_kappa_ee;
@ -1363,7 +1363,7 @@ qmckl_exit_code qmckl_compute_asymp_jasb (
for (int i=0; i<=1; ++i) { for (int i=0; i<=1; ++i) {
x = kappa_inv; x = kappa_inv;
for (int p = 1; p < bord_num-1; ++i){ for (int p = 1; p < bord_num; ++i){
x = x * kappa_inv; x = x * kappa_inv;
asymp_jasb[i] = asymp_jasb[i] + bord_vector[p + 1] * x; asymp_jasb[i] = asymp_jasb[i] + bord_vector[p + 1] * x;
} }