anchor.diff: -!
  -!
  -!
  +
   ! 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
  
positive.code:       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 + &
  
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)
  