home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / octave-1.1.1p1-src.tgz / tar.out / fsf / octave / libcruft / odepack / vnorm.f < prev    next >
Text File  |  1996-09-28  |  726b  |  19 lines
  1.       DOUBLE PRECISION FUNCTION VNORM (N, V, W)
  2. CLLL. OPTIMIZE
  3. C-----------------------------------------------------------------------
  4. C THIS FUNCTION ROUTINE COMPUTES THE WEIGHTED ROOT-MEAN-SQUARE NORM
  5. C OF THE VECTOR OF LENGTH N CONTAINED IN THE ARRAY V, WITH WEIGHTS
  6. C CONTAINED IN THE ARRAY W OF LENGTH N..
  7. C   VNORM = SQRT( (1/N) * SUM( V(I)*W(I) )**2 )
  8. C-----------------------------------------------------------------------
  9.       INTEGER N,   I
  10.       DOUBLE PRECISION V, W,   SUM
  11.       DIMENSION V(N), W(N)
  12.       SUM = 0.0D0
  13.       DO 10 I = 1,N 
  14.  10     SUM = SUM + (V(I)*W(I))**2
  15.       VNORM = DSQRT(SUM/DBLE(N))      
  16.       RETURN
  17. C----------------------- END OF FUNCTION VNORM -------------------------
  18.       END 
  19.