HERWIG/jimmy/divon4/rlen.F

   1 *CMZ :          23/08/93  13.30.12  by  Jonathan Butterworth
   2 *-- Author :
   3       SUBROUTINE RLEN (LEN,VEC,VLEN)
   4       INTEGER   LEN
   5       DOUBLE PRECISION   VEC(LEN), VLEN
   6       INTEGER  I
   7       DOUBLE PRECISION   ABSV, RATIO, RMIN, TOL, TVMX, VMAX
   8       SAVE
   9       DATA TOL / 1.0D-20 /
  10       DATA RMIN / 1.0D-34/
  11       VMAX=0.0D+0
  12       DO 10 I=1,LEN
  13       ABSV=ABS(VEC(I))
  14       IF(ABSV.GT.VMAX) VMAX=ABSV
  15  10   CONTINUE
  16       VLEN=0.0D+0
  17       IF(VMAX.EQ.0.0D+0) RETURN
  18       TVMX=0.0D+0
  19       IF(VMAX.GT.RMIN) TVMX=TOL*VMAX
  20       DO 20 I=1,LEN
  21       IF(ABS(VEC(I)).LE.TVMX) GOTO 20
  22       RATIO=VEC(I)/VMAX
  23       VLEN=VLEN+RATIO*RATIO
  24  20   CONTINUE
  25       VLEN=VMAX*SQRT(VLEN)
  26       RETURN
  27       END