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
  -          call mqc_integral_allocate(est_integral,'overlap','spin',tmpMatrixAlpha, &
  -            tmpMatrixAlpha)
  +          else
  +            if(MQC_Matrix_HaveComplex(tmpMatrixAlpha)) then
  +              call mqc_matrix_symm2full(tmpMatrixAlpha,'hermitian')
  +              tmpMatrixAlpha = transpose(tmpMatrixAlpha)
  +            endIf
  +            call mqc_integral_allocate(est_integral,'overlap','spin',tmpMatrixAlpha, &
  +              tmpMatrixAlpha)
  +          endIf
           elseIf(fileinfo%isGeneral()) then
             call fileInfo%getArray('OVERLAP',tmpMatrixAlpha,foundOut=found)
             if(present(foundObj)) foundObj = found
  
positive.code:             else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            endIf
            call fileInfo%getArray('BETA FOCK MATRIX',tmpMatrixBeta,foundOut=found)
            if(present(foundObj)) foundObj = found
            if(.not.found) then
              errorMsg = 'BETA FOCK MATRIX not present on file'
              if(present(foundObj)) then
                write(6,'(A)') errorMsg
              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
            if(MQC_Matrix_HaveComplex(tmpMatrixBeta)) then
              call mqc_matrix_symm2full(tmpMatrixBeta,'hermitian')
              tmpMatrixBeta = transpose(tmpMatrixBeta)
            endIf
            call mqc_integral_allocate(est_integral,'fock','spin',tmpMatrixAlpha, &
              tmpMatrixBeta)
          elseIf(fileinfo%isGeneral()) then
            call fileInfo%getArray('ALPHA FOCK MATRIX',tmpMatrixAlpha,foundOut=found)
            if(present(foundObj)) foundObj = found
  
negative.code:             else
                call mqc_error_l('errorMsg',6,'found',found)
              endIf
            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(), &
  