Line data    Source code

       1              : !--------------------------------------------------------------------------------
       2              : ! Copyright (c) 2016 Peter Grünberg Institut, Forschungszentrum Jülich, Germany
       3              : ! This file is part of FLEUR and available as free software under the conditions
       4              : ! of the MIT license as expressed in the LICENSE file in more detail.
       5              : !--------------------------------------------------------------------------------
       6              : 
       7              : MODULE m_wann_abinv
       8              : CONTAINS
       9           40 :   SUBROUTINE wann_abinv(atoms,sym,acof,bcof,ccof)
      10              :     !     ***************************************************************
      11              :     !     Transform acof,bcof,ccof in case of atoms related by inversion
      12              :     !     symmetry to obtain the coefficients in the global frame.
      13              :     !     Based on abcrot.
      14              :     !     Frank Freimuth
      15              :     !     ***************************************************************
      16              :     USE m_types
      17              :     IMPLICIT NONE
      18              :     !     ..
      19              :     !     .. Scalar Arguments ..
      20              :     TYPE(t_atoms),INTENT(IN) :: atoms
      21              :     TYPE(t_sym),INTENT(IN)   :: sym
      22              : 
      23              :     COMPLEX, INTENT (INOUT) :: acof(:,0:,:)
      24              :     COMPLEX, INTENT (INOUT) :: bcof(:,0:,:)
      25              :     COMPLEX, INTENT (INOUT) :: ccof(-atoms%llod:,:,:,:)!(-llod:llod,neigd,nlod,natd)
      26              : 
      27              :     !     .. Local Scalars ..
      28              :     INTEGER :: itype,ineq,iatom,ilo,l
      29              : 
      30           40 :     iatom=0
      31           80 :     DO itype=1,atoms%ntype
      32          160 :        DO ineq=1,atoms%neq(itype)
      33           80 :           iatom=iatom+1
      34           80 :           IF(sym%invsat(iatom).NE.2) CYCLE
      35          160 :           DO l=1,atoms%lmax(itype),2
      36              :              acof(:,l**2:l*(l+2),iatom) = (-1)**l *&
      37         7680 :                   acof(:,l**2:l*(l+2),iatom)
      38              :              bcof(:,l**2:l*(l+2),iatom) = (-1)**l * &
      39         7720 :                   bcof(:,l**2:l*(l+2),iatom)
      40              :           ENDDO
      41           80 :           DO ilo=1,atoms%nlo(itype)
      42            0 :              l=atoms%llo(ilo,itype)
      43           80 :              IF(l.GT.0) THEN
      44            0 :                 IF(MOD(l,2).EQ.0)CYCLE 
      45              :                 ccof(-l:l,:,ilo,iatom) = (-1)**l * &
      46            0 :                      ccof(-l:l,:,ilo,iatom)
      47              :              ENDIF
      48              :           ENDDO
      49              :        ENDDO
      50              :     ENDDO
      51              : 
      52           40 :   END SUBROUTINE wann_abinv
      53              : END MODULE m_wann_abinv
      54              :       
      55              : 
      56              :