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.

csscal.f 2.1 kB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394
  1. *> \brief \b CSSCAL
  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 CSSCAL(N,SA,CX,INCX)
  12. *
  13. * .. Scalar Arguments ..
  14. * REAL SA
  15. * INTEGER INCX,N
  16. * ..
  17. * .. Array Arguments ..
  18. * COMPLEX CX(*)
  19. * ..
  20. *
  21. *
  22. *> \par Purpose:
  23. * =============
  24. *>
  25. *> \verbatim
  26. *>
  27. *> CSSCAL scales a complex vector by a real 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 CSSCAL(N,SA,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. REAL SA
  62. INTEGER INCX,N
  63. * ..
  64. * .. Array Arguments ..
  65. COMPLEX CX(*)
  66. * ..
  67. *
  68. * =====================================================================
  69. *
  70. * .. Local Scalars ..
  71. INTEGER I,NINCX
  72. * ..
  73. * .. Intrinsic Functions ..
  74. INTRINSIC AIMAG,CMPLX,REAL
  75. * ..
  76. IF (N.LE.0 .OR. INCX.LE.0) RETURN
  77. IF (INCX.EQ.1) THEN
  78. *
  79. * code for increment equal to 1
  80. *
  81. DO I = 1,N
  82. CX(I) = CMPLX(SA*REAL(CX(I)),SA*AIMAG(CX(I)))
  83. END DO
  84. ELSE
  85. *
  86. * code for increment not equal to 1
  87. *
  88. NINCX = N*INCX
  89. DO I = 1,NINCX,INCX
  90. CX(I) = CMPLX(SA*REAL(CX(I)),SA*AIMAG(CX(I)))
  91. END DO
  92. END IF
  93. RETURN
  94. END