diff --git a/org/qmckl_jastrow.org b/org/qmckl_jastrow.org index d0da070..9601b51 100644 --- a/org/qmckl_jastrow.org +++ b/org/qmckl_jastrow.org @@ -1370,7 +1370,11 @@ print("asymp_jasb[1] : ", asymp_jasb[1]) #+begin_src c :tangle (eval c_test) assert(qmckl_electron_provided(context)); +int64_t type_nucl_num = n2_type_nucl_num; int64_t* type_nucl_vector = &(n2_type_nucl_vector[0]); +int64_t aord_num = n2_aord_num; +int64_t bord_num = n2_bord_num; +int64_t cord_num = n2_cord_num; double* aord_vector = &(n2_aord_vector[0][0]); double* bord_vector = &(n2_bord_vector[0]); double* cord_vector = &(n2_cord_vector[0][0]); @@ -1380,11 +1384,11 @@ rc = qmckl_init_jastrow(context); assert(!qmckl_jastrow_provided(context)); /* Set the data */ -rc = qmckl_set_jastrow_ord_num(context, n2_aord_num, n2_bord_num, n2_cord_num); +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, n2_type_nucl_num); +rc = qmckl_set_jastrow_type_nucl_num(context, type_nucl_num); assert(rc == QMCKL_SUCCESS); -rc = qmckl_set_jastrow_type_nucl_vector(context, n2_type_nucl_vector, nucl_num); +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);