anchor.diff:        Real (8) :: gaunt
         Complex (8) zfmtinp, zdotc
         External gaunt, zfmtinp, zdotc
  -!
  -!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
  -      input%groundstate%lmaxvr = 3
  -      lmmaxvr = 16
  -!
   ! check if q-vector is zero
         t1 = input%properties%elnes%vecql(1) ** 2 + &
        & input%properties%elnes%vecql(2) ** 2 + &
  
positive.code: !
  !
  !
  ! Copyright (C) 2008 J. K. Dewhurst, S. Sharma and C. Ambrosch-Draxl.
  ! This file is distributed under the terms of the GNU General Public License.
  ! See the file COPYING for license details.
  !
  !
  Subroutine genexpiqr (ik, emat)
        Use modmain
        Use modinput
  
negative.code:       Else
           emat (:, :) = em (:, :)
        End If
        Deallocate (igkqig, gnt, jlqr, vgkql, vgkqc, gkqc, tpgkqc)
        Deallocate (sfacgkq, apwalm1, apwalm2, evecfv1, evecfv2)
        If (input%groundstate%tevecsv) deallocate (evecsv1, evecsv2)
  