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

          Line data    Source code
       1           0 :       SUBROUTINE cfnorm(mrad,is,it,nsol,nmatch,jtop,var,gc,fc,rc,rc2,
       2           0 :      +                  dx,gck,fck)
       3             : c...........................................................cfnorm
       4             : c wavefunctions normalization
       5             : c
       6             :       IMPLICIT NONE
       7             : c
       8             : C     .. Scalar Arguments ..
       9             :       INTEGER, INTENT (IN) :: mrad
      10             :       REAL dx
      11             :       INTEGER is,it,jtop,nmatch,nsol
      12             : C     ..
      13             : C     .. Array Arguments ..
      14             :       REAL fc(2,2,mrad),fck(2,2,mrad),gc(2,2,mrad),gck(2,2,mrad),
      15             :      +     rc(mrad),rc2(mrad),var(4)
      16             : C     ..
      17             : C     .. Local Scalars ..
      18             :       REAL xnorm
      19             :       INTEGER i,j,k,n
      20             : C     ..
      21             : C     .. Local Arrays ..
      22           0 :       REAL rint(mrad)
      23             : C     ..
      24             : C     .. External Functions ..
      25             :       REAL rsimp
      26             :       EXTERNAL rsimp
      27             : C     ..
      28             : C     .. Intrinsic Functions ..
      29             :       INTRINSIC sqrt
      30             : C     ..
      31           0 :       DO n = 1,mrad
      32           0 :          DO k = 1,nsol
      33           0 :             gck(k,is,n) = 0.0
      34           0 :             fck(k,is,n) = 0.0
      35             :          END DO
      36             :       END DO
      37             : C                         ---------------------------------
      38             : C--->                     NORMALIZE WAVEFUNCTIONS ACCORDING
      39             : C                               TO MATCHING CONDITIONS
      40             : C                         ---------------------------------
      41             : C                                    INWARD - SOLUTION
      42           0 :       DO 30 n = nmatch,jtop
      43           0 :          DO 20 j = 1,nsol
      44           0 :             DO 10 i = 1,nsol
      45           0 :                gc(i,j,n) = gc(i,j,n)*var(2)
      46           0 :                fc(i,j,n) = fc(i,j,n)*var(2)
      47           0 :    10       CONTINUE
      48           0 :    20    CONTINUE
      49           0 :    30 CONTINUE
      50             : c
      51           0 :       IF (nsol.EQ.2) THEN
      52             : C                                   OUTWARD - SOLUTION
      53           0 :          DO 50 n = 1, (nmatch-1)
      54           0 :             DO 40 i = 1,nsol
      55           0 :                gc(i,it,n) = gc(i,it,n)*var(3)
      56           0 :                fc(i,it,n) = fc(i,it,n)*var(3)
      57           0 :    40       CONTINUE
      58           0 :    50    CONTINUE
      59             : C                                    INWARD - SOLUTION
      60           0 :          DO 70 n = nmatch,jtop
      61           0 :             DO 60 i = 1,nsol
      62           0 :                gc(i,it,n) = gc(i,it,n)*var(4)
      63           0 :                fc(i,it,n) = fc(i,it,n)*var(4)
      64           0 :    60       CONTINUE
      65           0 :    70    CONTINUE
      66             :       END IF
      67             : C                                    SUM FOR EACH KAPPA
      68           0 :       DO n = 1,jtop
      69           0 :          DO k = 1,nsol
      70           0 :             gck(k,is,n) = 0.00
      71           0 :             fck(k,is,n) = 0.00
      72             :          END DO
      73             :       END DO
      74           0 :       DO n = 1,jtop
      75           0 :          DO k = 1,nsol
      76           0 :             DO j = 1,nsol
      77           0 :                gck(k,is,n) = gck(k,is,n) + gc(k,j,n)
      78           0 :                fck(k,is,n) = fck(k,is,n) + fc(k,j,n)
      79             :             END DO
      80             :          END DO
      81             :       END DO
      82             : C                       -----------------------------------
      83             : C                       CALCULATE  NORM  AND NORMALIZE TO 1
      84             : C                       -----------------------------------
      85           0 :       DO n = 1,jtop
      86           0 :          rint(n) = 0.00
      87           0 :          DO k = 1,nsol
      88           0 :             rint(n) = rint(n) + rc2(n)* (gck(k,is,n)**2+fck(k,is,n)**2)
      89             :          END DO
      90             :       END DO
      91           0 :       xnorm = rsimp(mrad,rint,rc,jtop,dx)
      92           0 :       xnorm = 1.00/sqrt(xnorm)
      93           0 :       DO n = 1,jtop
      94           0 :          DO k = 1,nsol
      95           0 :             gck(k,is,n) = gck(k,is,n)*xnorm
      96           0 :             fck(k,is,n) = fck(k,is,n)*xnorm
      97             :          END DO
      98             :       END DO
      99           0 :       RETURN
     100             :       END

Generated by: LCOV version 1.13