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.

ccopy.f 2.1 kB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394
  1. *> \brief \b CCOPY
  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 CCOPY(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. *> CCOPY copies a vector x to a vector y.
  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 CCOPY(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. INTEGER I,IX,IY
  69. * ..
  70. IF (N.LE.0) RETURN
  71. IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
  72. *
  73. * code for both increments equal to 1
  74. *
  75. DO I = 1,N
  76. CY(I) = CX(I)
  77. END DO
  78. ELSE
  79. *
  80. * code for unequal increments or equal increments
  81. * not equal to 1
  82. *
  83. IX = 1
  84. IY = 1
  85. IF (INCX.LT.0) IX = (-N+1)*INCX + 1
  86. IF (INCY.LT.0) IY = (-N+1)*INCY + 1
  87. DO I = 1,N
  88. CY(IY) = CX(IX)
  89. IX = IX + INCX
  90. IY = IY + INCY
  91. END DO
  92. END IF
  93. RETURN
  94. END