anchor.diff:          cycle
         endif
         
  -      ll = l-l_start+1
  -!      T2d = 0.d0
  -!      do a=2,a_array(1,ll),3
  -!        ik = a_array(a,ll)
  -!        j = a_array(a+1,ll)
  -!        T2d(ik,j) = transfer(a_array(a+2,ll), 1.d0)
  -!      enddo
  -
  -      a=2
  +      a=l_pointer(l)
         do j=j_start,j_end
           ik=0
           do k=k_start,k_end
             do i=i_start,k
               ik = ik+1
  -            if ( (ik /= a_array(a,ll)).or.(j /= a_array(a+1,ll)) ) then
  +            if ( (ik /= a_array(a)).or.(j /= a_array(a+1))  &
  +                 .or.(a >= l_pointer(l+1)) ) then
                 T2d(ik,j) = 0.d0
               else
  -              T2d(ik,j) = transfer(a_array(a+2,ll), 1.d0)
  +              T2d(ik,j) = transfer(a_array(a+2), 1.d0)
                 a=a+3
               endif
             enddo
  
positive.code:   ! Create a temporary memory-mapped file
    integer                        :: fd
    type(c_ptr)                    :: c_pointer
    integer*8, pointer             :: a_array(:,:)
    call mmap(trim(ezfio_filename)//'/work/four_idx',                  &
        (/ (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,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) )
  
  
    !$OMP DO SCHEDULE(dynamic,4)
    do l=l_start,l_end
      a = 1
      ll = l-l_start+1
      do j=j_start,j_end
        ik=0
        do k=k_start,k_end
          do i=i_start,k
            ik = ik+1
            call bielec_integrals_index(i,j,k,l,idx)
            call map_get(map_a,idx,tmp)
            if (tmp /= 0.d0) then
              a = a+1
              a_array(a,ll) = ik
              a = a+1
              a_array(a,ll) = j
              a = a+1
              a_array(a,ll) = transfer(dble(tmp), 1_8)
            endif
          enddo
        enddo
      enddo
      a_array(a+1,ll) = 0
      a_array(1,ll) = a
    enddo
    !$OMP END DO
  
    allocate( T2d((i_end-i_start+1)*(k_end-k_start+2)/2, j_start:j_end), &
              V2d((i_end-i_start+1)*(k_end-k_start+2)/2, b_start:b_end), &
  
negative.code:   call map_sort(map_c)
  
    call munmap( &
        (/ 12_8 * map_a % n_elements /), 8, fd, c_pointer)
    deallocate(l_pointer)
  
  end
  