1
C/MEMBR ADD NAME=WASUM,SSI=0
3
double precision function wasum(n,xr,xi,incx)
6
c cette fonction determine la addition des normes l1 des
7
c composantes d'un vecteur complexe dont les parties reelles
8
c sont rangees dans le vecteur double precision xr et les
9
c parties imaginaires dans le vecteur double precision xi.
13
c double precision function wasum(n,xr,xi,incx)
15
c n: entier, taille du vecteur traite
17
c xr, xi: vecteurs double precision contenant,
18
c respectivement, les parties reelles et imaginaires du
21
c incx: increment entre deux composantes consecutives des
26
c cleve moler.- mathlab.
29
double precision xr(*),xi(*),s
32
if (n .le. 0) go to 20
35
s = s + abs(xr(ix)) + abs(xi(ix))