@@ -2184,8 +2184,8 @@ rc = qmckl_check(context,
2184
2184
0,
2185
2185
&(asymp_jasb_hpc[0]) )
2186
2186
);
2187
- assert(fabs(asymp_jasb_doc[0]-asymp_jasb_hpc[0]) < 1.e-12 );
2188
- assert(fabs(asymp_jasb_doc[1]-asymp_jasb_hpc[1]) < 1.e-12 );
2187
+ assert(fabs(asymp_jasb_doc[0]-asymp_jasb_hpc[0]) < 1.e-8 );
2188
+ assert(fabs(asymp_jasb_doc[1]-asymp_jasb_hpc[1]) < 1.e-8 );
2189
2189
2190
2190
#+end_src
2191
2191
@@ -2567,7 +2567,7 @@ assert(qmckl_electron_provided(context));
2567
2567
assert(rc == QMCKL_SUCCESS);
2568
2568
2569
2569
for (int64_t i=0 ; i<walk_num*elec_num*elec_num*(cord_num+1) ; i++) {
2570
- assert(fabs(ee_distance_rescaled_hpc[i] - ee_distance_rescaled_doc[i]) < 1.e-12 );
2570
+ assert(fabs(ee_distance_rescaled_hpc[i] - ee_distance_rescaled_doc[i]) < 1.e-10 );
2571
2571
}
2572
2572
}
2573
2573
#+end_src
@@ -2989,7 +2989,7 @@ assert(qmckl_electron_provided(context));
2989
2989
assert(rc == QMCKL_SUCCESS);
2990
2990
2991
2991
for (int i = 0; i < walk_num*nucl_num*elec_num*4; i++) {
2992
- assert(fabs(ee_distance_rescaled_gl_doc[i] - ee_distance_rescaled_gl_hpc[i]) < 1.e-12 );
2992
+ assert(fabs(ee_distance_rescaled_gl_doc[i] - ee_distance_rescaled_gl_hpc[i]) < 1.e-8 );
2993
2993
}
2994
2994
}
2995
2995
@@ -3543,7 +3543,7 @@ assert(qmckl_jastrow_champ_provided(context));
3543
3543
assert (rc == QMCKL_SUCCESS);
3544
3544
3545
3545
for (int64_t i = 0; i < walk_num; i++) {
3546
- assert(fabs(factor_ee_doc[i] - factor_ee_hpc[i]) < 1.e-12 );
3546
+ assert(fabs(factor_ee_doc[i] - factor_ee_hpc[i]) < 1.e-8 );
3547
3547
}
3548
3548
}
3549
3549
#+end_src
@@ -4310,7 +4310,7 @@ assert(qmckl_jastrow_champ_provided(context));
4310
4310
4311
4311
for (int64_t i = 0 ; i < walk_num*4*elec_num ; i++) {
4312
4312
printf("%ld %f %f\n", i, factor_ee_gl_hpc[i], factor_ee_gl_doc[i]);
4313
- assert(fabs(factor_ee_gl_hpc[i] - factor_ee_gl_doc[i]) < 1.e-12 );
4313
+ assert(fabs(factor_ee_gl_hpc[i] - factor_ee_gl_doc[i]) < 1.e-8 );
4314
4314
}
4315
4315
}
4316
4316
#+end_src
@@ -5002,7 +5002,7 @@ assert(fabs(en_distance_rescaled[0][0][6] - 0.4726452953409436) < 1.e-12);
5002
5002
assert(rc == QMCKL_SUCCESS);
5003
5003
5004
5004
for (int64_t i=0 ; i<walk_num*nucl_num*elec_num ; ++i) {
5005
- assert(fabs(en_distance_rescaled_doc[i] - en_distance_rescaled_hpc[i]) < 1.e-12 );
5005
+ assert(fabs(en_distance_rescaled_doc[i] - en_distance_rescaled_hpc[i]) < 1.e-8 );
5006
5006
}
5007
5007
}
5008
5008
#+end_src
@@ -5519,7 +5519,7 @@ assert(qmckl_electron_provided(context));
5519
5519
assert(rc == QMCKL_SUCCESS);
5520
5520
5521
5521
for (int i = 0; i < walk_num*nucl_num*elec_num*4; i++) {
5522
- assert(fabs(en_distance_rescaled_gl_doc[i] - en_distance_rescaled_gl_hpc[i]) < 1.e-12 );
5522
+ assert(fabs(en_distance_rescaled_gl_doc[i] - en_distance_rescaled_gl_hpc[i]) < 1.e-8 );
5523
5523
}
5524
5524
}
5525
5525
#+end_src
@@ -5987,7 +5987,7 @@ assert(fabs(22.781375792083587 - factor_en[0]) < 1.e-12);
5987
5987
assert(rc == QMCKL_SUCCESS);
5988
5988
5989
5989
for (int64_t i = 0; i < walk_num; i++) {
5990
- assert(fabs(factor_en_doc[i] - factor_en_hpc[i]) < 1.e-12 );
5990
+ assert(fabs(factor_en_doc[i] - factor_en_hpc[i]) < 1.e-10 );
5991
5991
}
5992
5992
}
5993
5993
#+end_src
@@ -11762,7 +11762,7 @@ TODO
11762
11762
factor_een_gl_hpc);
11763
11763
11764
11764
for (int64_t i = 0; i < walk_num*4*elec_num; ++i) {
11765
- assert(fabs(factor_een_gl_doc[i] - factor_een_gl_hpc[i]) < 1e-12 );
11765
+ assert(fabs(factor_een_gl_doc[i] - factor_een_gl_hpc[i]) < 1e-8 );
11766
11766
}
11767
11767
}
11768
11768
#+end_src
0 commit comments