Commit eef27ff0 authored by Matthias Redies's avatar Matthias Redies

wfolap1 is also dead

parent 0943aebf
......@@ -307,57 +307,4 @@ CONTAINS
! wfolap = wfolap + cdum
END FUNCTION wfolap_noinv
FUNCTION wfolap1(cmt1, cpw1, cmt2, cpw2, ngpt1, ngpt2, olappw, olapmt,&
& atoms, mpbasis,hybrid)
USE m_types
IMPLICIT NONE
TYPE(t_mpbasis), intent(in) :: mpbasis
TYPE(t_hybrid), INTENT(IN) :: hybrid
TYPE(t_atoms), INTENT(IN) :: atoms
! -scalars -
COMPLEX :: wfolap1
INTEGER, INTENT(IN) :: ngpt1, ngpt2
! - arrays -
COMPLEX, INTENT(IN) :: cmt1(:,:), cmt2(:,:)
#if ( defined(CPP_INVERSION) )
REAL, INTENT(IN) :: cpw1(:), cpw2(:)
#else
COMPLEX, INTENT(IN) :: cpw1(:), cpw2(:)
#endif
#if ( defined(CPP_INVERSION) )
REAL, INTENT(IN) :: olappw(:,:)
#else
COMPLEX, INTENT(IN) :: olappw(:,:)
#endif
REAL, INTENT(IN) :: olapmt(maxval(mpbasis%num_radfun_per_l), maxval(mpbasis%num_radfun_per_l), 0:atoms%lmaxd, atoms%ntype)
! - local -
INTEGER :: itype, ieq, iatom, l, m, lm, nn
wfolap1 = 0
iatom = 0
DO itype = 1, atoms%ntype
DO ieq = 1, atoms%neq(itype)
iatom = iatom + 1
lm = 0
DO l = 0, atoms%lmax(itype)
DO M = -l, l
nn = mpbasis%num_radfun_per_l(l, itype)
wfolap1 = wfolap1 + &
& dot_product(cmt1(lm + 1:lm + nn, iatom),&
& matmul(olapmt(:nn, :nn, l, itype),&
& cmt2(lm + 1:lm + nn, iatom)))
lm = lm + nn
END DO
END DO
END DO
END DO
wfolap1 = wfolap1 + dot_product(cpw1, matmul(olappw, cpw2))
END FUNCTION wfolap1
END MODULE m_olap
......@@ -74,7 +74,7 @@ CONTAINS
USE m_types
USE m_util, ONLY: modulo1
use m_intgrf, only: intgrf, intgrf_init
USE m_olap, ONLY: wfolap_inv, wfolap_noinv, wfolap1, wfolap_init
USE m_olap, ONLY: wfolap_inv, wfolap_noinv, wfolap_init
USE m_trafo, ONLY: waveftrafo_symm
USE m_io_hybrid
IMPLICIT NONE
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment