anchor.diff:    ! Create a temporary memory-mapped file
     integer                        :: fd
     type(c_ptr)                    :: c_pointer
  -  integer*8, pointer             :: a_array(:,:,:)
  +  integer*8, pointer             :: a_array(:,:)
     call mmap(trim(ezfio_filename)//'/work/four_idx',                  &
  -      (/ 2_8,int(i_end-i_start+1,8),int(j_end-j_start+1,8),int(k_end-k_start+1,8), int(l_end-l_start+1,8) /), 8, fd, .False., c_pointer)
  -  call c_f_pointer(c_pointer, a_array, (/ 4, (i_end-i_start+1)*(j_end-j_start+1)*(k_end-k_start+1)/2, l_end-l_start+1 /))
  -
  +      (/ (int(i_end-i_start+1,8)*int(j_end-j_start+2,8)*int(k_end-k_start+1,8)),int(l_end-l_start+1,8) /), 16, fd, .False., c_pointer)
  +  call c_f_pointer(c_pointer, a_array, (/ ((i_end-i_start+1)*(j_end-j_start+1)*(k_end-k_start+1)*3_8)/2_8, l_end-l_start+1_8 /))
   
     !$OMP PARALLEL DEFAULT(NONE) SHARED(a_array,c_pointer,fd,          &
         !$OMP  a_start,a_end,b_start,b_end,c_start,c_end,d_start,d_end,&
         !$OMP  i_start,i_end,j_start,j_end,k_start,k_end,l_start,l_end,&
         !$OMP  i_min,i_max,j_min,j_max,k_min,k_max,l_min,l_max,        &
         !$OMP  map_a,map_c,matrix_B)                                   &
  -      !$OMP  PRIVATE(key,value,T,U,V,i,j,k,l,idx,ik,   &
  +      !$OMP  PRIVATE(key,value,T,U,V,i,j,k,l,idx,ik,ll,   &
         !$OMP  a,b,c,d,tmp,T2d,V2d)
     allocate( key(i_max*j_max*k_max), value(i_max*j_max*k_max) )
     allocate( U(a_start:a_end, c_start:c_end, b_start:b_end) )
  
positive.code: 
      !$OMP CRITICAL
      call map_append(map_c, key, value, idx) 
      call map_sort(map_c)
      !$OMP END CRITICAL
  
  
  
negative.code:   double precision, allocatable  :: T2d(:,:), V2d(:,:)
    integer                        :: i_max, j_max, k_max, l_max
    integer                        :: i_min, j_min, k_min, l_min
    integer                        :: i, j, k, l, ik, ll
    integer                        :: a, b, c, d
    double precision, external     :: get_ao_bielec_integral
    integer(key_kind)              :: idx
  