Esope source files

   1 : C VECUNI    SOURCE    CHAT      05/01/13    04:05:26     5004           
   2 : c---------------------------------------------------------------------
   3 : c
   4 :       SUBROUTINE VECUNI (N, V, MODUL)
   5 : c
   6 : c=====================================================================
   7 : c                                                                    =
   8 : c   This routine computes the length of the vector 'v' and converts  =
   9 : c   it to a unit one.                                                =
  10 : c                                                                    =
  11 : c=====================================================================
  12 :       IMPLICIT INTEGER(I-N)
  13 :       integer*4 n
  14 :       real*8    v (n), modul
  15 : c
  16 :       integer*4 i
  17 :       real*8    moduli
  18 : c
  19 :       modul = 0.0 d0
  20 :       do i = 1, n
  21 :          modul = modul + v (i) **2
  22 :       end do
  23 : c
  24 :       if (modul .eq. 1.0 d0) return
  25 : c
  26 :       modul = SQRT (modul)
  27 :       moduli = 1.0 d0 / modul
  28 :       do i = 1, n
  29 :          v (i) = v (i) * moduli
  30 :       end do
  31 : c
  32 :       return
  33 :       end
  34 :  
  35 :