share/odepack/fortran/dmnorm.f

   1 *DECK DMNORM
   2       DOUBLE PRECISION FUNCTION DMNORM (N, V, W)
   3 C-----------------------------------------------------------------------
   4 C This function routine computes the weighted max-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   DMNORM = MAX(i=1,...,N) ABS(V(i))*W(i)
   8 C-----------------------------------------------------------------------
   9       INTEGER N,   I
  10       DOUBLE PRECISION V, W,   VM
  11       DIMENSION V(N), W(N)
  12       VM = 0.0D0
  13       DO 10 I = 1,N
  14  10     VM = MAX(VM,ABS(V(I))*W(I))
  15       DMNORM = VM
  16       RETURN
  17 C----------------------- End of Function DMNORM ------------------------
  18       END