anchor.diff:           !DIR$ FORCEINLINE
            call bielec_integrals_index(i,j,k,l,idx)
            !DIR$ FORCEINLINE
  -         call map_get(mo_integrals_map,idx,mo_integrals_cache(i,j,k,l))
  +         call map_get(mo_integrals_map,idx,integral)
  +         ii = l-mo_integrals_cache_min
  +         ii = ior( ishft(ii,6), k-mo_integrals_cache_min)
  +         ii = ior( ishft(ii,6), j-mo_integrals_cache_min)
  +         ii = ior( ishft(ii,6), i-mo_integrals_cache_min)
  +         mo_integrals_cache(ii) = integral
          enddo
        enddo
      enddo
  
positive.code:   integer, intent(in)            :: i,j,k,l
    integer(key_kind)              :: idx
    type(map_type), intent(inout)  :: map
    real(integral_kind)            :: tmp
    PROVIDE ao_bielec_integrals_in_map
    !DIR$ FORCEINLINE
    if (ao_overlap_abs(i,k)*ao_overlap_abs(j,l) < ao_integrals_threshold ) then
      tmp = 0.d0
    else if (ao_bielec_integral_schwartz(i,k)*ao_bielec_integral_schwartz(j,l) < ao_integrals_threshold) then
      tmp = 0.d0
    else
      if ( (i >= ao_integrals_cache_min) .and. &
          (j >= ao_integrals_cache_min) .and. &
          (k >= ao_integrals_cache_min) .and. &
          (l >= ao_integrals_cache_min) .and. &
          (i <= ao_integrals_cache_max) .and. &
          (j <= ao_integrals_cache_max) .and. &
          (k <= ao_integrals_cache_max) .and. &
          (l <= ao_integrals_cache_max) ) then
          tmp = ao_integrals_cache(i,j,k,l)
      else
        !DIR$ FORCEINLINE
        call bielec_integrals_index(i,j,k,l,idx)
        call map_get(map,idx,tmp)
      endif
    endif
    get_ao_bielec_integral = tmp
  end
  
  
  
negative.code:          !DIR$ FORCEINLINE
           call bielec_integrals_index(i,j,k,l,idx)
           !DIR$ FORCEINLINE
           call map_get(ao_integrals_map,idx,integral)
           ii = l-ao_integrals_cache_min
           ii = ior( ishft(ii,6), k-ao_integrals_cache_min)
           ii = ior( ishft(ii,6), j-ao_integrals_cache_min)
           ii = ior( ishft(ii,6), i-ao_integrals_cache_min)
           ao_integrals_cache(ii) = integral
         enddo
       enddo
     enddo
  