pure integer function find_guga_excit_lvl_to_doubles(ilutI, ilutJ)
! make an highly optimized excitation level finder for guga up
! to doubles (for efficiency reasons)
! maybe use this and the SD version to initialize a pointer
! at startup to use the correct one and to be sure to have a
! correct GUGA excit-lvl info at all necessary stages
integer(n_int), intent(in) :: ilutI(0:GugaBits%len_tot), ilutJ(0:GugaBits%len_tot)
unused_var(ilutI)
unused_var(ilutJ)
find_guga_excit_lvl_to_doubles = nel
call stop_all("here", "todo")
end function find_guga_excit_lvl_to_doubles