You can not select more than 25 topics Topics must start with a chinese character,a letter or number, can include dashes ('-') and can be up to 35 characters long.

cswap.f 2.2 kB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798
  1. *> \brief \b CSWAP
  2. *
  3. * =========== DOCUMENTATION ===========
  4. *
  5. * Online html documentation available at
  6. * http://www.netlib.org/lapack/explore-html/
  7. *
  8. * Definition:
  9. * ===========
  10. *
  11. * SUBROUTINE CSWAP(N,CX,INCX,CY,INCY)
  12. *
  13. * .. Scalar Arguments ..
  14. * INTEGER INCX,INCY,N
  15. * ..
  16. * .. Array Arguments ..
  17. * COMPLEX CX(*),CY(*)
  18. * ..
  19. *
  20. *
  21. *> \par Purpose:
  22. * =============
  23. *>
  24. *> \verbatim
  25. *>
  26. *> CSWAP interchanges two vectors.
  27. *> \endverbatim
  28. *
  29. * Authors:
  30. * ========
  31. *
  32. *> \author Univ. of Tennessee
  33. *> \author Univ. of California Berkeley
  34. *> \author Univ. of Colorado Denver
  35. *> \author NAG Ltd.
  36. *
  37. *> \date November 2011
  38. *
  39. *> \ingroup complex_blas_level1
  40. *
  41. *> \par Further Details:
  42. * =====================
  43. *>
  44. *> \verbatim
  45. *>
  46. *> jack dongarra, linpack, 3/11/78.
  47. *> modified 12/3/93, array(1) declarations changed to array(*)
  48. *> \endverbatim
  49. *>
  50. * =====================================================================
  51. SUBROUTINE CSWAP(N,CX,INCX,CY,INCY)
  52. *
  53. * -- Reference BLAS level1 routine (version 3.4.0) --
  54. * -- Reference BLAS is a software package provided by Univ. of Tennessee, --
  55. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  56. * November 2011
  57. *
  58. * .. Scalar Arguments ..
  59. INTEGER INCX,INCY,N
  60. * ..
  61. * .. Array Arguments ..
  62. COMPLEX CX(*),CY(*)
  63. * ..
  64. *
  65. * =====================================================================
  66. *
  67. * .. Local Scalars ..
  68. COMPLEX CTEMP
  69. INTEGER I,IX,IY
  70. * ..
  71. IF (N.LE.0) RETURN
  72. IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
  73. *
  74. * code for both increments equal to 1
  75. DO I = 1,N
  76. CTEMP = CX(I)
  77. CX(I) = CY(I)
  78. CY(I) = CTEMP
  79. END DO
  80. ELSE
  81. *
  82. * code for unequal increments or equal increments not equal
  83. * to 1
  84. *
  85. IX = 1
  86. IY = 1
  87. IF (INCX.LT.0) IX = (-N+1)*INCX + 1
  88. IF (INCY.LT.0) IY = (-N+1)*INCY + 1
  89. DO I = 1,N
  90. CTEMP = CX(IX)
  91. CX(IX) = CY(IY)
  92. CY(IY) = CTEMP
  93. IX = IX + INCX
  94. IY = IY + INCY
  95. END DO
  96. END IF
  97. RETURN
  98. END