LCOV - code coverage report
Current view: top level - core - rsimp.f (source / functions) Hit Total Coverage
Test: combined.info Lines: 0 21 0.0 %
Date: 2019-09-08 04:53:50 Functions: 0 1 0.0 %

          Line data    Source code
       1           0 :       REAL FUNCTION rsimp(mrad,f,r,irn,dx)
       2             : c.....................................................rsimp
       3             : 
       4             :       IMPLICIT NONE
       5             : c===================
       6             : c     radial integration via simpson
       7             : C     ..
       8             : C     .. Scalar Arguments ..
       9             :       INTEGER,INTENT (IN) :: mrad
      10             :       REAL dx
      11             :       INTEGER irn
      12             : C     ..
      13             : C     .. Array Arguments ..
      14             :       REAL f(mrad),r(mrad)
      15             : C     ..
      16             : C     .. Local Scalars ..
      17             :       REAL s
      18             :       INTEGER i,isw,nl,np
      19             : C     ..
      20           0 :       isw = 0
      21           0 :       rsimp = 0.0
      22           0 :       IF (irn.LE.2) RETURN
      23           0 :       IF (irn/2*2.EQ.irn) isw = 1
      24           0 :       np = irn - isw
      25           0 :       s = f(1)*r(1) + f(np)*r(np)
      26           0 :       nl = np - 1
      27           0 :       DO 10 i = 2,nl,2
      28           0 :          s = s + 4.0*f(i)*r(i)
      29           0 :    10 CONTINUE
      30           0 :       nl = nl - 1
      31           0 :       IF (nl.LT.3) GO TO 30
      32           0 :       DO 20 i = 3,nl,2
      33           0 :          s = s + 2.0*f(i)*r(i)
      34           0 :    20 CONTINUE
      35           0 :    30 s = s*dx/3.0
      36           0 :       IF (isw.EQ.1) GO TO 40
      37             :       rsimp = s
      38           0 :       RETURN
      39           0 :    40 rsimp = s + (f(irn)*r(irn)+f(irn-1)*r(irn-1))*0.50*dx
      40           0 :       RETURN
      41             :       END

Generated by: LCOV version 1.13