small 2-body fixes

Author: EmielSlootman

org/qmckl_forces.org

@@ -3085,6 +3085,8 @@ integer(qmckl_exit_code) function qmckl_compute_forces_jastrow_een( &
 
   forces_jastrow_een = 0.d0
 
+  if (cord_num == 0) return
+
   do nw =1, walk_num
     do n = 1, dim_c_vector
       l = lkpm_combined_index(n, 1)
@@ -4961,34 +4963,36 @@ qmckl_exit_code qmckl_provide_forces_jastrow_single_een(qmckl_context context)
                            NULL);
   }
 
-  /* Check if en rescaled distance is provided */
-  rc = qmckl_provide_een_rescaled_n(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+  if (ctx->jastrow_champ.cord_num > 0) {
 
-  /* Check if en rescaled distance derivatives is provided */
-  rc = qmckl_provide_een_rescaled_n_gl(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    /* Check if en rescaled distance is provided */
+    rc = qmckl_provide_een_rescaled_n(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
-  /* Check if en rescaled distance is provided */
-  rc = qmckl_provide_een_rescaled_single_n(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    /* Check if en rescaled distance derivatives is provided */
+    rc = qmckl_provide_een_rescaled_n_gl(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
-  /* Check if en rescaled distance derivatives is provided */
-  rc = qmckl_provide_een_rescaled_single_n_gl(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    /* Check if en rescaled distance is provided */
+    rc = qmckl_provide_een_rescaled_single_n(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
-  rc = qmckl_provide_tmp_c(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    /* Check if en rescaled distance derivatives is provided */
+    rc = qmckl_provide_een_rescaled_single_n_gl(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
-  rc = qmckl_provide_forces_tmp_c(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    rc = qmckl_provide_tmp_c(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
-  rc = qmckl_provide_jastrow_champ_delta_p(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    rc = qmckl_provide_forces_tmp_c(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
-  rc = qmckl_provide_forces_jastrow_delta_p(context);
-  if(rc != QMCKL_SUCCESS) return rc;
+    rc = qmckl_provide_jastrow_champ_delta_p(context);
+    if(rc != QMCKL_SUCCESS) return rc;
 
+    rc = qmckl_provide_forces_jastrow_delta_p(context);
+    if(rc != QMCKL_SUCCESS) return rc;
+  }
 
 
   /* Compute if necessary */