Mercurial > octave-nkf
view libcruft/odepack/vnorm.f @ 14626:f947d2922feb stable rc-3-6-2-0
3.6.2-rc0 release candidate
* configure.ac (AC_INIT): Version is now 3.6.2-rc0.
(OCTAVE_RELEASE_DATE): Now 2012-05-11.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 11 May 2012 13:46:18 -0400 |
parents | 30c606bec7a8 |
children |
line wrap: on
line source
DOUBLE PRECISION FUNCTION VNORM (N, V, W) CLLL. OPTIMIZE C----------------------------------------------------------------------- C THIS FUNCTION ROUTINE COMPUTES THE WEIGHTED ROOT-MEAN-SQUARE NORM C OF THE VECTOR OF LENGTH N CONTAINED IN THE ARRAY V, WITH WEIGHTS C CONTAINED IN THE ARRAY W OF LENGTH N.. C VNORM = SQRT( (1/N) * SUM( V(I)*W(I) )**2 ) C----------------------------------------------------------------------- INTEGER N, I DOUBLE PRECISION V, W, SUM DIMENSION V(N), W(N) SUM = 0.0D0 DO 10 I = 1,N 10 SUM = SUM + (V(I)*W(I))**2 VNORM = DSQRT(SUM/DBLE(N)) RETURN C----------------------- END OF FUNCTION VNORM ------------------------- END