LCOV - code coverage report
Current view: top level - xc-pot - pbecor2.f90 (source / functions) Hit Total Coverage
Test: FLEUR test coverage Lines: 8 8 100.0 %
Date: 2024-03-28 04:22:06 Functions: 1 1 100.0 %

          Line data    Source code
       1             : MODULE m_pbecor2
       2             : !---------------------
       3             : ! slimmed down version of gcor used in pw91 routines, to interpolate
       4             : ! lsd correlation energy, as given by (10) of
       5             : ! j. p. perdew and y. wang, phys. rev. b {\bf 45}, 13244 (1992).
       6             : ! k. burke, may 11, 1996.
       7             : !---------------------
       8             : CONTAINS
       9   460796259 :    SUBROUTINE pbecor2( &
      10             :       a,a1,b1,b2,b3,b4,rtrs, &
      11             :       gg,ggrs)
      12             :       IMPLICIT NONE
      13             : 
      14             :       REAL, INTENT (IN)  :: a,a1,b1,b2,b3,b4,rtrs
      15             :       REAL, INTENT (OUT) :: gg,ggrs
      16             : 
      17             :       REAL :: q0,q1,q2,q3
      18             : !     ..
      19   460796259 :       q0 = -2.e0*a* (1.e0+a1*rtrs*rtrs)
      20   460796259 :       q1 = 2.e0*a*rtrs* (b1+rtrs* (b2+rtrs* (b3+b4*rtrs)))
      21   460796259 :       q2 = log(1.e0+1.e0/q1)
      22   460796259 :       gg = q0*q2
      23   460796259 :       q3 = a* (b1/rtrs+2.e0*b2+rtrs* (3.e0*b3+4.e0*b4*rtrs))
      24   460796259 :       ggrs = -2.e0*a*a1*q2 - q0*q3/ (q1* (1.e0+q1))
      25             : 
      26   460796259 :    END SUBROUTINE pbecor2
      27             : END MODULE m_pbecor2

Generated by: LCOV version 1.14