Fixed broken filter_integrals in h_apply

scripts/generate_h_apply.py

@@ -11,6 +11,7 @@ declarations
 decls_main
 deinit_thread
 init_main
+filter_integrals
 filter2p
 filter2h2p_double
 filter2h2p_single
@@ -104,6 +105,11 @@ class H_apply(object):
     s["do_mono_excitations"] = d[do_mono_exc]
     s["do_double_excitations"] = d[do_double_exc]
     s["keys_work"]  += "call fill_H_apply_buffer_no_selection(key_idx,keys_out,N_int,iproc)"
+    s["filter_integrals"] = "array_pairs = .True."
+    if SingleRef:
+      s["filter_integrals"] = """
+      call get_mo_bielec_integrals_existing_ik(i_a,j_a,mo_num,array_pairs,mo_integrals_map)
+      """
 
     s["generate_psi_guess"]  = """
   ! Sort H_jj to find the N_states lowest states

src/determinants/h_apply.template.f

@@ -240,7 +240,10 @@ subroutine $subroutine_diexcOrg(key_in,key_mask,hole_1,particl_1,hole_2, particl
       !   hole = a^(+)_j_a(ispin) a_i_a(ispin)|key_in> : single exc :: orb(i_a,ispin) --> orb(j_a,ispin)
       hole_save = hole
 
-      if (ispin == 1) then
+      ! Build array of the non-zero integrals of second excitation
+      $filter_integrals
+
+        if (ispin == 1) then
         integer                        :: jjj
 
         i=0