dvnrms.f 1.1 KB

123456789101112131415161718192021222324252627282930313233343536
  1. *DECK DVNRMS
  2. DOUBLE PRECISION FUNCTION DVNRMS (N, V, W)
  3. C***BEGIN PROLOGUE DVNRMS
  4. C***SUBSIDIARY
  5. C***PURPOSE Subsidiary to DDEBDF
  6. C***LIBRARY SLATEC
  7. C***TYPE DOUBLE PRECISION (VNWRMS-S, DVNRMS-D)
  8. C***AUTHOR (UNKNOWN)
  9. C***DESCRIPTION
  10. C
  11. C DVNRMS computes a weighted root-mean-square vector norm for the
  12. C integrator package DDEBDF.
  13. C
  14. C***SEE ALSO DDEBDF
  15. C***ROUTINES CALLED (NONE)
  16. C***REVISION HISTORY (YYMMDD)
  17. C 820301 DATE WRITTEN
  18. C 890531 Changed all specific intrinsics to generic. (WRB)
  19. C 890831 Modified array declarations. (WRB)
  20. C 890911 Removed unnecessary intrinsics. (WRB)
  21. C 891214 Prologue converted to Version 4.0 format. (BAB)
  22. C 900328 Added TYPE section. (WRB)
  23. C***END PROLOGUE DVNRMS
  24. INTEGER I, N
  25. DOUBLE PRECISION SUM, V, W
  26. DIMENSION V(*),W(*)
  27. C***FIRST EXECUTABLE STATEMENT DVNRMS
  28. SUM = 0.0D0
  29. DO 10 I = 1, N
  30. SUM = SUM + (V(I)/W(I))**2
  31. 10 CONTINUE
  32. DVNRMS = SQRT(SUM/N)
  33. RETURN
  34. C ----------------------- END OF FUNCTION DVNRMS
  35. C ------------------------
  36. END