anchor.diff:          call check_allocate('seen_iha', size(seen_ijij), ierr)
           allocate(seen_ijij((active_basis_offset*(active_basis_offset+1))/2), stat=ierr)
           call check_allocate('seen_ijij', size(seen_ijij), ierr)
  -        allocate(seen_iaib(active_basis_offset,(nbasis*(nbasis+1))/2), stat=ierr)
  +        allocate(seen_iaib(-active_basis_offset+1:0,(nbasis*(nbasis+1))/2), stat=ierr)
           call check_allocate('seen_iaib', size(seen_iaib), ierr)
           seen_iha = .false.
           seen_ijij = 0
  
positive.code: 
          end do
  
          deallocate(seen_iha, stat=ierr)
          call check_deallocate('seen_iha', ierr)
          deallocate(seen_ijij, stat=ierr)
  
negative.code:         use point_group_symmetry, only: init_pg_symmetry
          use system, only: fcidump, uhf, ecore, nel
  
          use checking, only: check_allocate, check_deallocate
          use errors, only: stop_all
          use parallel, only: parent
          use ranking, only: insertion_rank_rp
          use utils, only: get_free_unit, tri_ind_reorder, int_fmt
  
          use, intrinsic :: iso_fortran_env
  
  