Branch data     Line data    Source code

       1                 :             : !!****m* ABINIT/m_qplusg
       2                 :             : !! NAME
       3                 :             : !!  m_qplusg
       4                 :             : !!
       5                 :             : !! FUNCTION
       6                 :             : !!
       7                 :             : !! COPYRIGHT
       8                 :             : !! Copyright (C) 1999-2026 ABINIT group (MG, FB, BG)
       9                 :             : !! This file is distributed under the terms of the
      10                 :             : !! GNU General Public License, see ~abinit/COPYING
      11                 :             : !! or http://www.gnu.org/copyleft/gpl.txt .
      12                 :             : !! For the initials of contributors, see ~abinit/doc/developers/contributors.txt .
      13                 :             : !!
      14                 :             : !! SOURCE
      15                 :             : 
      16                 :             : #if defined HAVE_CONFIG_H
      17                 :             : #include "config.h"
      18                 :             : #endif
      19                 :             : 
      20                 :             : #include "abi_common.h"
      21                 :             : 
      22                 :             : module m_qplusg
      23                 :             : 
      24                 :             :  use defs_basis
      25                 :             : 
      26                 :             :  implicit none
      27                 :             : 
      28                 :             :  private
      29                 :             : !!***
      30                 :             : 
      31                 :             :  public :: cmod_qpg
      32                 :             : !!***
      33                 :             : 
      34                 :             : CONTAINS
      35                 :             : !!***
      36                 :             : 
      37                 :             : !!****f* m_qplusg/cmod_qpg
      38                 :             : !! NAME
      39                 :             : !! cmod_qpg
      40                 :             : !!
      41                 :             : !! FUNCTION
      42                 :             : !! Set up table of lengths |q+G| for the Coulomb potential.
      43                 :             : !!
      44                 :             : !! INPUTS
      45                 :             : !! gvec(3,npwvec)=Reduced coordinates of the G vectors.
      46                 :             : !! gprimd(3,3)=Dimensional primitive translations for reciprocal space ($\textrm{bohr}^{-1}$)
      47                 :             : !! iq=Index specifying the q point.
      48                 :             : !! npwvec=Number of planewaves
      49                 :             : !! nq=Number of q points.
      50                 :             : !! q=Coordinates of q points.
      51                 :             : !!
      52                 :             : !! OUTPUT
      53                 :             : !! qplusg(npwvec)=Norm of q+G vector
      54                 :             : !!
      55                 :             : !! SOURCE
      56                 :             : 
      57                 :        2903 : subroutine cmod_qpg(nq, iq, q, npwvec, gvec, gprimd, qplusg)
      58                 :             : 
      59                 :             : !Arguments ------------------------------------
      60                 :             : !scalars
      61                 :             :  integer,intent(in) :: iq,npwvec,nq
      62                 :             : !arrays
      63                 :             :  integer,intent(in) :: gvec(3,npwvec)
      64                 :             :  real(dp),intent(in) :: gprimd(3,3),q(3,nq)
      65                 :             :  real(dp),intent(out) :: qplusg(npwvec)
      66                 :             : 
      67                 :             : !Local variables ------------------------------
      68                 :             : !scalars
      69                 :             :  integer :: ig
      70                 :             : !arrays
      71                 :             :  real(dp) :: gmet(3,3),gpq(3)
      72                 :             : !************************************************************************
      73                 :             : 
      74                 :             :  ! Compute reciprocal space metrics
      75   [ +  +  +  +  :      116120 :  gmet = MATMUL(TRANSPOSE(gprimd),gprimd)
                   +  + ]
      76                 :             : 
      77   [ +  +  +  +  :        5454 :  if (ALL(ABS(q(:,iq)) < tol3)) then
                   +  + ]
      78                 :             :    ! Treat q as if it were zero except when G=0
      79   [ +  +  +  +  :       11664 :    qplusg(1)=two_pi*SQRT(DOT_PRODUCT(q(:,iq),MATMUL(gmet,q(:,iq))))
                   +  + ]
      80         [ +  + ]:       95983 :    do ig=2,npwvec
      81         [ +  + ]:      381016 :      gpq(:)=gvec(:,ig)
      82   [ +  +  +  +  :     1524793 :      qplusg(ig)=two_pi*SQRT(DOT_PRODUCT(gpq,MATMUL(gmet,gpq)))
                   +  + ]
      83                 :             :    end do
      84                 :             :  else
      85         [ +  + ]:      233340 :    do ig=1,npwvec
      86         [ +  + ]:      924664 :      gpq(:)=gvec(:,ig)+q(:,iq)
      87   [ +  +  +  +  :     3700830 :      qplusg(ig)=two_pi*SQRT(DOT_PRODUCT(gpq,MATMUL(gmet,gpq)))
                   +  + ]
      88                 :             :    end do
      89                 :             :  end if
      90                 :             : 
      91                 :        2903 : end subroutine cmod_qpg
      92                 :             : !!***
      93                 :             : 
      94                 :             : end module m_qplusg
      95                 :             : !!***