Mercurial > octave-nkf
view libcruft/odepack/svnorm.f @ 16088:b29b10fbb744 stable release-3-6-4
Version 3.6.4 released.
* configure.ac (AC_INIT): Version is now 3.6.4.
(OCTAVE_RELEASE_DATE): Now 2013-02-21.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Thu, 21 Feb 2013 15:17:54 -0500 |
parents | 96ba591be50f |
children |
line wrap: on
line source
REAL FUNCTION SVNORM (N, V, W) C***BEGIN PROLOGUE SVNORM C***SUBSIDIARY C***PURPOSE Weighted root-mean-square vector norm. C***TYPE SINGLE PRECISION (SVNORM-S, DVNORM-D) C***AUTHOR Hindmarsh, Alan C., (LLNL) C***DESCRIPTION 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 SVNORM = SQRT( (1/N) * SUM( V(i)*W(i) )**2 ) C C***SEE ALSO SLSODE C***ROUTINES CALLED (NONE) C***REVISION HISTORY (YYMMDD) C 791129 DATE WRITTEN C 890501 Modified prologue to SLATEC/LDOC format. (FNF) C 890503 Minor cosmetic changes. (FNF) C 930809 Renamed to allow single/double precision versions. (ACH) C***END PROLOGUE SVNORM C**End INTEGER N, I REAL V, W, SUM DIMENSION V(N), W(N) C C***FIRST EXECUTABLE STATEMENT SVNORM SUM = 0.0E0 DO 10 I = 1,N 10 SUM = SUM + (V(I)*W(I))**2 SVNORM = SQRT(SUM/N) RETURN C----------------------- END OF FUNCTION SVNORM ------------------------ END