TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2025-01-03 10:06:09 +01:00
Code Issues Projects Releases Wiki Activity

compute_dim_cord_vect done

Browse Source
This commit is contained in:
Gianfranco Abrusci 2022-03-09 11:15:15 +01:00
parent d13693a822
commit a4ba5deac4
1 changed files with 6 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
org/qmckl_jastrow.org
View File

@ -4979,28 +4979,21 @@ qmckl_exit_code qmckl_compute_dim_cord_vect (
return QMCKL_INVALID_ARG_2;
}
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) {
for (int p=2; p <= cord_num; ++p){
for (int k=p-1; k >= 0; --k) {
if (k != 0) {
lmax = p - k;
} else {
lmax = p - k - 2;
}
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
for (int l = lmax; l >= 0; --l) {
if ( ((p - k - l) & 1)==1) continue;
printf("in da loop\n");
*dim_cord_vect = *dim_cord_vect + 1;
*dim_cord_vect=*dim_cord_vect+1;
}
}
}
printf("hello 2 %d\n", *dim_cord_vect);
return QMCKL_SUCCESS;
}