anchor.diff:    implicit none
     
     integer :: i
  -  print *, "pro"
  -  print *, "pross"
     do i=1, N_det
       call get_mask_phase(psi_det_sorted(1,1,i), psi_phasemask(1,1,i))
     end do
  
positive.code:   integer(bit_kind)              :: hole_mask(N_int,2), particle_mask(N_int,2)
    double precision               :: fock_diag_tmp(2,mo_tot_num+1)
    
    print *, "BEG"
    
    
    call build_fock_tmp(fock_diag_tmp,psi_det_generators(1,1,i_generator),N_int)
  
  
negative.code:       particle_mask(k,:) = hole_mask(k,:)
      enddo
      
      call select_doubles(i_generator,hole_mask,particle_mask,fock_diag_tmp,E0,pt2,b)
      call select_singles(i_generator,hole_mask,particle_mask,fock_diag_tmp,E0,pt2,b)
    enddo
  end subroutine
  
  
  