1
0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2024-07-18 00:43:51 +02:00

Testing - II. #22

This commit is contained in:
vijay gopal chilkuri 2021-07-07 20:45:55 +05:30
parent 22c8f01bab
commit 478ce7cc38

View File

@ -1402,28 +1402,28 @@ double* aord_vector = &(n2_aord_vector[0][0]);
double* bord_vector = &(n2_bord_vector[0]);
double* cord_vector = &(n2_cord_vector[0][0]);
///* Initialize the Jastrow data */
//rc = qmckl_init_jastrow(context);
//assert(!qmckl_jastrow_provided(context));
//
///* Set the data */
//rc = qmckl_set_jastrow_ord_num(context, aord_num, bord_num, cord_num);
//assert(rc == QMCKL_SUCCESS);
//rc = qmckl_set_jastrow_type_nucl_num(context, type_nucl_num);
//assert(rc == QMCKL_SUCCESS);
//rc = qmckl_set_jastrow_type_nucl_vector(context, type_nucl_vector, nucl_num);
//assert(rc == QMCKL_SUCCESS);
//rc = qmckl_set_jastrow_aord_vector(context, aord_vector);
//assert(rc == QMCKL_SUCCESS);
//rc = qmckl_set_jastrow_bord_vector(context, bord_vector);
//assert(rc == QMCKL_SUCCESS);
//rc = qmckl_set_jastrow_cord_vector(context, cord_vector);
//assert(rc == QMCKL_SUCCESS);
//rc = qmckl_set_jastrow_dependencies(context);
//assert(rc == QMCKL_SUCCESS);
//
///* Check if Jastrow is properly initialized */
//assert(qmckl_jastrow_provided(context));
/* Initialize the Jastrow data */
rc = qmckl_init_jastrow(context);
assert(!qmckl_jastrow_provided(context));
/* Set the data */
rc = qmckl_set_jastrow_ord_num(context, aord_num, bord_num, cord_num);
assert(rc == QMCKL_SUCCESS);
rc = qmckl_set_jastrow_type_nucl_num(context, type_nucl_num);
assert(rc == QMCKL_SUCCESS);
rc = qmckl_set_jastrow_type_nucl_vector(context, type_nucl_vector, nucl_num);
assert(rc == QMCKL_SUCCESS);
rc = qmckl_set_jastrow_aord_vector(context, aord_vector);
assert(rc == QMCKL_SUCCESS);
rc = qmckl_set_jastrow_bord_vector(context, bord_vector);
assert(rc == QMCKL_SUCCESS);
rc = qmckl_set_jastrow_cord_vector(context, cord_vector);
assert(rc == QMCKL_SUCCESS);
rc = qmckl_set_jastrow_dependencies(context);
assert(rc == QMCKL_SUCCESS);
/* Check if Jastrow is properly initialized */
assert(qmckl_jastrow_provided(context));
//
//double asymp_jasb[2];
//rc = qmckl_get_jastrow_asymp_jasb(context, asymp_jasb);