anchor.diff:              else
                 call mqc_error_l('errorMsg',6,'found',found)
               endIf
  -          endIf
  -          if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
  -            call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
  -            tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  -          endIf
  -          nBasis = fileInfo%getVal('nBasis')
  -          call mqc_matrix_spinBlockGHF(tmpMatrixAlpha)
  -          tmpMatrixBeta = tmpMatrixAlpha%mat([nBasis+1,-1],[nBasis+1,-1])
  -          tmpMatrixBetaAlpha = tmpMatrixAlpha%mat([1,nBasis],[nBasis+1,-1])
  -          tmpMatrixAlphaBeta = tmpMatrixAlpha%mat([nBasis+1,-1],[1,nBasis])
  -          tmpMatrixAlpha = tmpMatrixAlpha%mat([1,nBasis],[1,nBasis])
  -          call mqc_integral_allocate(est_integral,'fock','general',tmpMatrixAlpha, &
  -            tmpMatrixBeta,tmpMatrixAlphaBeta,tmpMatrixBetaAlpha)
  +          else
  +            if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
  +              call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
  +              tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  +            endIf
  +            nBasis = fileInfo%getVal('nBasis')
  +            call mqc_matrix_spinBlockGHF(tmpMatrixAlpha)
  +            tmpMatrixBeta = tmpMatrixAlpha%mat([nBasis+1,-1],[nBasis+1,-1])
  +            tmpMatrixBetaAlpha = tmpMatrixAlpha%mat([1,nBasis],[nBasis+1,-1])
  +            tmpMatrixAlphaBeta = tmpMatrixAlpha%mat([nBasis+1,-1],[1,nBasis])
  +            tmpMatrixAlpha = tmpMatrixAlpha%mat([1,nBasis],[1,nBasis])
  +            call mqc_integral_allocate(est_integral,'fock','general',tmpMatrixAlpha, &
  +              tmpMatrixBeta,tmpMatrixAlphaBeta,tmpMatrixBetaAlpha)
  +          endIf
           else
             call mqc_error_L('Unknown wavefunction type in getESTObj', 6, &
                  'fileinfo%isRestricted()', fileinfo%isRestricted(), &
  
positive.code:             else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            endIf
            call mqc_integral_allocate(est_integral,'mo coefficients','space',tmpMatrixAlpha)
          elseIf(fileinfo%isUnrestricted()) then
            call fileInfo%getArray('ALPHA MO COEFFICIENTS',tmpMatrixAlpha,foundOut=found)
            if(present(foundObj)) foundObj = found
  
negative.code:             else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            else
              call mqc_eigenvalues_allocate(est_eigenvalues,'mo energies','space',tmpVectorAlpha)
            endIf
          elseIf(fileinfo%isUnrestricted()) then
            call fileInfo%getArray('ALPHA ORBITAL ENERGIES',vectorOut=tmpVectorAlpha,foundOut=found)
            if(present(foundObj)) foundObj = found
  