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.7 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125
  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. * Arguments:
  30. * ==========
  31. *
  32. *> \param[in] N
  33. *> \verbatim
  34. *> N is INTEGER
  35. *> number of elements in input vector(s)
  36. *> \endverbatim
  37. *>
  38. *> \param[in] CX
  39. *> \verbatim
  40. *> CX is COMPLEX array, dimension ( 1 + ( N - 1 )*abs( INCX ) )
  41. *> \endverbatim
  42. *>
  43. *> \param[in] INCX
  44. *> \verbatim
  45. *> INCX is INTEGER
  46. *> storage spacing between elements of CX
  47. *> \endverbatim
  48. *>
  49. *> \param[out] CY
  50. *> \verbatim
  51. *> CY is COMPLEX array, dimension ( 1 + ( N - 1 )*abs( INCY ) )
  52. *> \endverbatim
  53. *>
  54. *> \param[in] INCY
  55. *> \verbatim
  56. *> INCY is INTEGER
  57. *> storage spacing between elements of CY
  58. *> \endverbatim
  59. *
  60. * Authors:
  61. * ========
  62. *
  63. *> \author Univ. of Tennessee
  64. *> \author Univ. of California Berkeley
  65. *> \author Univ. of Colorado Denver
  66. *> \author NAG Ltd.
  67. *
  68. *> \date November 2017
  69. *
  70. *> \ingroup complex_blas_level1
  71. *
  72. *> \par Further Details:
  73. * =====================
  74. *>
  75. *> \verbatim
  76. *>
  77. *> jack dongarra, linpack, 3/11/78.
  78. *> modified 12/3/93, array(1) declarations changed to array(*)
  79. *> \endverbatim
  80. *>
  81. * =====================================================================
  82. SUBROUTINE CCOPY(N,CX,INCX,CY,INCY)
  83. *
  84. * -- Reference BLAS level1 routine (version 3.8.0) --
  85. * -- Reference BLAS is a software package provided by Univ. of Tennessee, --
  86. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  87. * November 2017
  88. *
  89. * .. Scalar Arguments ..
  90. INTEGER INCX,INCY,N
  91. * ..
  92. * .. Array Arguments ..
  93. COMPLEX CX(*),CY(*)
  94. * ..
  95. *
  96. * =====================================================================
  97. *
  98. * .. Local Scalars ..
  99. INTEGER I,IX,IY
  100. * ..
  101. IF (N.LE.0) RETURN
  102. IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
  103. *
  104. * code for both increments equal to 1
  105. *
  106. DO I = 1,N
  107. CY(I) = CX(I)
  108. END DO
  109. ELSE
  110. *
  111. * code for unequal increments or equal increments
  112. * not equal to 1
  113. *
  114. IX = 1
  115. IY = 1
  116. IF (INCX.LT.0) IX = (-N+1)*INCX + 1
  117. IF (INCY.LT.0) IY = (-N+1)*INCY + 1
  118. DO I = 1,N
  119. CY(IY) = CX(IX)
  120. IX = IX + INCX
  121. IY = IY + INCY
  122. END DO
  123. END IF
  124. RETURN
  125. END