From 47e9b559d92a2d4a334eb672c7429ffdfb5f6687 Mon Sep 17 00:00:00 2001
From: Qiming Sun <osirpt.sun@gmail.com>
Date: Wed, 27 Sep 2023 23:12:24 -0700
Subject: [PATCH] Fix SR int1e_grids integrals

---
 src/g1e_grids.c | 10 +++++++++-
 1 file changed, 9 insertions(+), 1 deletion(-)

diff --git a/src/g1e_grids.c b/src/g1e_grids.c
index 52729a4..e40ecad 100644
--- a/src/g1e_grids.c
+++ b/src/g1e_grids.c
@@ -27,7 +27,15 @@ void CINTinit_int1e_grids_EnvVars(CINTEnvVars *envs, int *ng, int *shls,
         envs->common_factor = 2 * M_PI
                 * CINTcommon_fac_sp(envs->i_l) * CINTcommon_fac_sp(envs->j_l);
 
-        int nroots = envs->nrys_roots;
+        int rys_order = envs->nrys_roots;
+        int nroots = rys_order;
+        double omega = env[PTR_RANGE_OMEGA];
+        if (omega < 0 && rys_order <= 3) {
+                nroots *= 2;
+        }
+        envs->rys_order = rys_order;
+        envs->nrys_roots = nroots;
+
         int dli, dlj;
         int ibase = envs->li_ceil > envs->lj_ceil;
         if (ibase) {