Line data    Source code

       1             : MODULE m_coreerr
       2             : 
       3             :    CONTAINS
       4             :    
       5         634 :    SUBROUTINE coreerr(err,var,s,nsol,pow,qow,piw,qiw)
       6             : !   ********************************************************************
       7             : !   *                                                                  *
       8             : !   *   CALCULATE THE MISMATCH OF THE RADIAL WAVE FUNCTIONS AT THE     *
       9             : !   *   POINT  NMATCH  FOR OUT- AND INWARD INTEGRATION                 *
      10             : !   *                                                                  *
      11             : !   ********************************************************************
      12             : 
      13             :       IMPLICIT NONE
      14             : 
      15             :       INTEGER nsol,s
      16             : 
      17             :       REAL err(4),piw(2,2),pow(2,2),qiw(2,2),qow(2,2),var(4)
      18             : 
      19             :       INTEGER t
      20             : 
      21         634 :       err(1) = pow(s,s) - piw(s,s)*var(2)
      22         634 :       err(2) = qow(s,s) - qiw(s,s)*var(2)
      23             : 
      24         634 :       IF (nsol.EQ.1) RETURN
      25             : 
      26         401 :       t = 3 - s
      27             : 
      28         401 :       err(1) = err(1) + pow(s,t)*var(3) - piw(s,t)*var(2)*var(4)
      29         401 :       err(2) = err(2) + qow(s,t)*var(3) - qiw(s,t)*var(2)*var(4)
      30         401 :       err(3) = pow(t,s) - piw(t,s)*var(2) + pow(t,t)*var(3) - piw(t,t)*var(2)*var(4)
      31         401 :       err(4) = qow(t,s) - qiw(t,s)*var(2) + qow(t,t)*var(3) - qiw(t,t)*var(2)*var(4)
      32             : 
      33             :    END SUBROUTINE coreerr
      34             : 
      35             : END MODULE m_coreerr