anchor.diff:      call map_unique(map_c)
   
     enddo
  +  deallocate(key,value)
   
   end
  
positive.code:       if (dabs(matrix_B(l,d)) < 1.d-10) then
          cycle
        endif
        do k=k_start,k_end
          do j=j_start,j_end
            do i=i_start,i_end
              call bielec_integrals_index(i,j,k,l,idx)
              call map_get(map_a,idx,tmp)
              T(i,j,k) = tmp
            enddo
          enddo
        enddo
  
        V = 0.d0
        do a=a_start,a_end
          do k=k_start,k_end
            do j=j_start,j_end
              do i=i_start,i_end
                V(j,k,a) = V(j,k,a) + T(i,j,k)*matrix_B(i,a)
              enddo
            enddo
          enddo
        enddo
  !        call DGEMM('T','N', (j_end-j_start+1),(k_end-k_start+1),     &
  !            (i_end-i_start+1), 1.d0,                                 &
  !          T, size(T,1)*
  
        T = 0.d0
        do b=b_start,b_end
          do a=a_start,a_end
            do k=k_start,k_end
              do j=j_start,j_end
                T(k,a,b) = T(k,a,b) + V(j,k,a)*matrix_B(j,b)
              enddo
            enddo
          enddo
        enddo
  
        V = 0.d0
        do c=c_start,c_end
          do b=b_start,b_end
            do a=a_start,a_end
              do k=k_start,k_end
                V(a,b,c) = V(a,b,c) + T(k,a,b)*matrix_B(k,c)
              enddo
            enddo
          enddo
        enddo
  
        do c=c_start,c_end
          do b=b_start,b_end
            do a=a_start,a_end
  !    do c=c_start,c_end
  !      do b=b_start,d
  !        do a=a_start,min(b,c)
              U(a,b,c) = U(a,b,c) + V(a,b,c) * matrix_B(l,d)
            enddo
          enddo
        enddo
  
      enddo
  
  
negative.code: !    do c=c_start,c_end
  !      do b=b_start,d
  !        do a=a_start,min(b,c)
            if (dabs(U(a,c,b)) < 1.d-15) then
              cycle
            endif
            idx = idx+1_8
            call bielec_integrals_index(a,b,c,d,key(idx))
            value(idx) = U(a,c,b)
          enddo
        enddo
      enddo
  