Line data    Source code

       1             :       module m_wann_gwf_auxovlp
       2             :       implicit none
       3             :       contains
       4             : 
       5           0 :       subroutine wann_gwf_auxovlp(q1,q2,a,mmnaux)
       6             :       use m_constants, only : pimach
       7             : 
       8             :       implicit none
       9             :       real,    intent(in)  :: q1(3),q2(3),a
      10             :       complex, intent(out) :: mmnaux
      11             : 
      12             :       integer :: i
      13             :       real :: pi, a_aux, b_aux, tau(3)
      14             :       complex :: ov(3)
      15             : 
      16           0 :       pi = pimach()
      17           0 :       a_aux = a
      18           0 :       b_aux = 0.98*a
      19           0 :       tau = 2.0*pi*(q1-q2)/a_aux
      20           0 :       ov = cmplx(1.,0.)
      21             : 
      22           0 :       do i=1,3
      23           0 :          if(abs(tau(i)).lt.1e-10) cycle
      24           0 :          ov(i) = 8.0*pi*pi*sin(tau(i)*b_aux/2.0)
      25           0 :          ov(i) = ov(i) / ( 4.0*pi*pi - (tau(i)*b_aux)**2 )
      26           0 :          ov(i) = ov(i) / ( tau(i)*b_aux )
      27             :       enddo
      28           0 :       mmnaux = ov(1)*ov(2)*ov(3)
      29             : 
      30           0 :       write(*,*)'ov(1) =',ov(1)
      31           0 :       write(*,*)'ov(2) =',ov(2)
      32           0 :       write(*,*)'ov(3) =',ov(3)
      33           0 :       write(*,*)'mmnaux=',mmnaux
      34             : 
      35           0 :       end subroutine wann_gwf_auxovlp
      36             : 
      37             :       end module m_wann_gwf_auxovlp