Statistiques
| Révision :

root / src / blas / ccopy.f @ 4

Historique | Voir | Annoter | Télécharger (969 octet)

1 1 pfleura2
      SUBROUTINE CCOPY(N,CX,INCX,CY,INCY)
2 1 pfleura2
*     .. Scalar Arguments ..
3 1 pfleura2
      INTEGER INCX,INCY,N
4 1 pfleura2
*     ..
5 1 pfleura2
*     .. Array Arguments ..
6 1 pfleura2
      COMPLEX CX(*),CY(*)
7 1 pfleura2
*     ..
8 1 pfleura2
*
9 1 pfleura2
*  Purpose
10 1 pfleura2
*  =======
11 1 pfleura2
*
12 1 pfleura2
*     CCOPY copies a vector x to a vector y.
13 1 pfleura2
*
14 1 pfleura2
*  Further Details
15 1 pfleura2
*  ===============
16 1 pfleura2
*
17 1 pfleura2
*     jack dongarra, linpack, 3/11/78.
18 1 pfleura2
*     modified 12/3/93, array(1) declarations changed to array(*)
19 1 pfleura2
*
20 1 pfleura2
*     .. Local Scalars ..
21 1 pfleura2
      INTEGER I,IX,IY
22 1 pfleura2
*     ..
23 1 pfleura2
      IF (N.LE.0) RETURN
24 1 pfleura2
      IF (INCX.EQ.1 .AND. INCY.EQ.1) GO TO 20
25 1 pfleura2
*
26 1 pfleura2
*        code for unequal increments or equal increments
27 1 pfleura2
*          not equal to 1
28 1 pfleura2
*
29 1 pfleura2
      IX = 1
30 1 pfleura2
      IY = 1
31 1 pfleura2
      IF (INCX.LT.0) IX = (-N+1)*INCX + 1
32 1 pfleura2
      IF (INCY.LT.0) IY = (-N+1)*INCY + 1
33 1 pfleura2
      DO 10 I = 1,N
34 1 pfleura2
          CY(IY) = CX(IX)
35 1 pfleura2
          IX = IX + INCX
36 1 pfleura2
          IY = IY + INCY
37 1 pfleura2
   10 CONTINUE
38 1 pfleura2
      RETURN
39 1 pfleura2
*
40 1 pfleura2
*        code for both increments equal to 1
41 1 pfleura2
*
42 1 pfleura2
   20 DO 30 I = 1,N
43 1 pfleura2
          CY(I) = CX(I)
44 1 pfleura2
   30 CONTINUE
45 1 pfleura2
      RETURN
46 1 pfleura2
      END