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.

cscal.f 1.9 kB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091
  1. *> \brief \b CSCAL
  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 CSCAL(N,CA,CX,INCX)
  12. *
  13. * .. Scalar Arguments ..
  14. * COMPLEX CA
  15. * INTEGER INCX,N
  16. * ..
  17. * .. Array Arguments ..
  18. * COMPLEX CX(*)
  19. * ..
  20. *
  21. *
  22. *> \par Purpose:
  23. * =============
  24. *>
  25. *> \verbatim
  26. *>
  27. *> CSCAL scales a vector by a constant.
  28. *> \endverbatim
  29. *
  30. * Authors:
  31. * ========
  32. *
  33. *> \author Univ. of Tennessee
  34. *> \author Univ. of California Berkeley
  35. *> \author Univ. of Colorado Denver
  36. *> \author NAG Ltd.
  37. *
  38. *> \date November 2011
  39. *
  40. *> \ingroup complex_blas_level1
  41. *
  42. *> \par Further Details:
  43. * =====================
  44. *>
  45. *> \verbatim
  46. *>
  47. *> jack dongarra, linpack, 3/11/78.
  48. *> modified 3/93 to return if incx .le. 0.
  49. *> modified 12/3/93, array(1) declarations changed to array(*)
  50. *> \endverbatim
  51. *>
  52. * =====================================================================
  53. SUBROUTINE CSCAL(N,CA,CX,INCX)
  54. *
  55. * -- Reference BLAS level1 routine (version 3.4.0) --
  56. * -- Reference BLAS is a software package provided by Univ. of Tennessee, --
  57. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  58. * November 2011
  59. *
  60. * .. Scalar Arguments ..
  61. COMPLEX CA
  62. INTEGER INCX,N
  63. * ..
  64. * .. Array Arguments ..
  65. COMPLEX CX(*)
  66. * ..
  67. *
  68. * =====================================================================
  69. *
  70. * .. Local Scalars ..
  71. INTEGER I,NINCX
  72. * ..
  73. IF (N.LE.0 .OR. INCX.LE.0) RETURN
  74. IF (INCX.EQ.1) THEN
  75. *
  76. * code for increment equal to 1
  77. *
  78. DO I = 1,N
  79. CX(I) = CA*CX(I)
  80. END DO
  81. ELSE
  82. *
  83. * code for increment not equal to 1
  84. *
  85. NINCX = N*INCX
  86. DO I = 1,NINCX,INCX
  87. CX(I) = CA*CX(I)
  88. END DO
  89. END IF
  90. RETURN
  91. END