From 2f246780eb0a08ee87ab5d98fe9c0e2f17685594 Mon Sep 17 00:00:00 2001
From: ydamour <yann.damour@hotmail.fr>
Date: Tue, 13 Jun 2023 14:05:13 +0200
Subject: [PATCH] fix bug in get_excitation_general

---
 src/utils_cc/org/phase.org | 2 ++
 src/utils_cc/phase.irp.f   | 2 ++
 2 files changed, 4 insertions(+)

diff --git a/src/utils_cc/org/phase.org b/src/utils_cc/org/phase.org
index 5f67859c..2156a251 100644
--- a/src/utils_cc/org/phase.org
+++ b/src/utils_cc/org/phase.org
@@ -137,6 +137,7 @@ subroutine get_excitation_general(det1,det2,degree,n,list_anni,list_crea,phase,N
   do j = 1, 2
     k = 1
     do i = 1, n1(j)
+       if (k > n_anni(j)) exit
        if (l1(i,j) /= list_anni(k,j)) cycle
        pos_anni(k,j) = i
        k = k + 1
@@ -147,6 +148,7 @@ subroutine get_excitation_general(det1,det2,degree,n,list_anni,list_crea,phase,N
   do j = 1, 2
     k = 1
     do i = 1, n2(j)
+       if (k > n_crea(j)) exit
        if (l2(i,j) /= list_crea(k,j)) cycle
        pos_crea(k,j) = i
        k = k + 1
diff --git a/src/utils_cc/phase.irp.f b/src/utils_cc/phase.irp.f
index 01b41f49..e0703fb8 100644
--- a/src/utils_cc/phase.irp.f
+++ b/src/utils_cc/phase.irp.f
@@ -96,6 +96,7 @@ subroutine get_excitation_general(det1,det2,degree,n,list_anni,list_crea,phase,N
   do j = 1, 2
     k = 1
     do i = 1, n1(j)
+       if (k > n_anni(j)) exit
        if (l1(i,j) /= list_anni(k,j)) cycle
        pos_anni(k,j) = i
        k = k + 1
@@ -106,6 +107,7 @@ subroutine get_excitation_general(det1,det2,degree,n,list_anni,list_crea,phase,N
   do j = 1, 2
     k = 1
     do i = 1, n2(j)
+       if (k > n_crea(j)) exit
        if (l2(i,j) /= list_crea(k,j)) cycle
        pos_crea(k,j) = i
        k = k + 1