1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
subroutine dset(n,dx,dy,incy)
c!but
c dset affecte un scalaire a tous les elements d'un vecteur
c!liste d'appel
c subroutine dset(n,dx,dy,incy)
c double precision dx,dy(n*incy)
c integer n,incy
c
c n : nombre d'elements du vecteur dy
c dx : scalaire a affecter
c dy : tableau contenant le vecteur
c incy : increment entre deux elements consecutifs du vecteur y
c dans le tableau dy
c!
c Copyright INRIA
double precision dx,dy(*)
c
if (n.le.0) return
iy = 1
if (incy.lt.0) iy = (-n+1)*incy + 1
do 10 i = 1,n
dy(iy) = dx
iy = iy + incy
10 continue
return
end
|
