ewset.f

Go to the documentation of this file.

      SUBROUTINE ewset (N, ITOL, RTOL, ATOL, YCUR, EWT)
CLLL. OPTIMIZE
C-----------------------------------------------------------------------
C THIS SUBROUTINE SETS THE ERROR WEIGHT VECTOR EWT ACCORDING TO
C     EWT(I) = RTOL(I)*ABS(YCUR(I)) + ATOL(I),  I = 1,...,N,
C WITH THE SUBSCRIPT ON RTOL AND/OR ATOL POSSIBLY REPLACED BY 1 ABOVE,
C DEPENDING ON THE VALUE OF ITOL.
C-----------------------------------------------------------------------
      INTEGER n, itol
      INTEGER i
      DOUBLE PRECISION rtol, atol, ycur, ewt
      dimension rtol(*), atol(*), ycur(n), ewt(n) 
C
      go to(10, 20, 30, 40), itol
 10   CONTINUE
      DO 15 i = 1,n 
 15     ewt(i) = rtol(1)*dabs(ycur(i)) + atol(1)
      RETURN
 20   CONTINUE
      DO 25 i = 1,n 
 25     ewt(i) = rtol(1)*dabs(ycur(i)) + atol(i)
      RETURN
 30   CONTINUE
      DO 35 i = 1,n 
 35     ewt(i) = rtol(i)*dabs(ycur(i)) + atol(1)
      RETURN
 40   CONTINUE
      DO 45 i = 1,n 
 45     ewt(i) = rtol(i)*dabs(ycur(i)) + atol(i)
      RETURN
C----------------------- END OF SUBROUTINE EWSET -----------------------
      END