TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2024-12-22 12:23:56 +01:00
Code Issues Projects Releases Wiki Activity

to be cleaned from print

Browse Source
This commit is contained in:
Gianfranco Abrusci 2022-03-02 10:18:45 +01:00
parent 4fac9f06c9
commit d13693a822
1 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
org/qmckl_jastrow.org
View File

@ -4979,23 +4979,28 @@ qmckl_exit_code qmckl_compute_dim_cord_vect (
return QMCKL_INVALID_ARG_2;
}
dim_cord_vect = 0;
printf("hello %d\n", *dim_cord_vect);
*dim_cord_vect = 0;
for (int p=2; p < cord_num; ++p){
printf("in da loop\n");
for (int k=p-1; k <= 0; --k) {
if (k != 0) {
lmax = p - k;
} else {
lmax = p - k - 2;
}
for (l = lmax; l <= 0; --l) {
printf("lmax %d\n", lmax);
for (int l = lmax; l <= 0; --l) {
// if ( iand(p - k - l, 1_8) == 1) continue;
// Does it make sense? it should
if ( ((p - k - l) & 1)==1) continue;
dim_cord_vect = dim_cord_vect + 1;
printf("in da loop\n");
*dim_cord_vect = *dim_cord_vect + 1;
}
}
}
printf("hello 2 %d\n", *dim_cord_vect);
return QMCKL_SUCCESS;
}