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           0 :   SUBROUTINE wann_abinv(atoms,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             : 
      22             :     COMPLEX, INTENT (INOUT) :: acof(:,0:,:)
      23             :     COMPLEX, INTENT (INOUT) :: bcof(:,0:,:)
      24             :     COMPLEX, INTENT (INOUT) :: ccof(-atoms%llod:,:,:,:)!(-llod:llod,neigd,nlod,natd)
      25             : 
      26             :     !     .. Local Scalars ..
      27             :     INTEGER :: itype,ineq,iatom,ilo,l
      28             : 
      29           0 :     iatom=0
      30           0 :     DO itype=1,atoms%ntype
      31           0 :        DO ineq=1,atoms%neq(itype)
      32           0 :           iatom=iatom+1
      33           0 :           IF(atoms%invsat(iatom).NE.2) CYCLE
      34           0 :           DO l=1,atoms%lmax(itype),2
      35             :              acof(:,l**2:l*(l+2),iatom) = (-1)**l *&
      36           0 :                   acof(:,l**2:l*(l+2),iatom)
      37             :              bcof(:,l**2:l*(l+2),iatom) = (-1)**l * &
      38           0 :                   bcof(:,l**2:l*(l+2),iatom)
      39             :           ENDDO
      40           0 :           DO ilo=1,atoms%nlo(itype)
      41           0 :              l=atoms%llo(ilo,itype)
      42           0 :              IF(l.GT.0) THEN
      43           0 :                 IF(MOD(l,2).EQ.0)CYCLE 
      44             :                 ccof(-l:l,:,ilo,iatom) = (-1)**l * &
      45           0 :                      ccof(-l:l,:,ilo,iatom)
      46             :              ENDIF
      47             :           ENDDO
      48             :        ENDDO
      49             :     ENDDO
      50             : 
      51           0 :   END SUBROUTINE wann_abinv
      52             : END MODULE m_wann_abinv
      53             :       
      54             : 
      55             :