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.

serrqp.f 3.3 kB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138
  1. *> \brief \b SERRQP
  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 SERRQP( PATH, NUNIT )
  12. *
  13. * .. Scalar Arguments ..
  14. * CHARACTER*3 PATH
  15. * INTEGER NUNIT
  16. * ..
  17. *
  18. *
  19. *> \par Purpose:
  20. * =============
  21. *>
  22. *> \verbatim
  23. *>
  24. *> SERRQP tests the error exits for SGEQP3.
  25. *> \endverbatim
  26. *
  27. * Arguments:
  28. * ==========
  29. *
  30. *> \param[in] PATH
  31. *> \verbatim
  32. *> PATH is CHARACTER*3
  33. *> The LAPACK path name for the routines to be tested.
  34. *> \endverbatim
  35. *>
  36. *> \param[in] NUNIT
  37. *> \verbatim
  38. *> NUNIT is INTEGER
  39. *> The unit number for output.
  40. *> \endverbatim
  41. *
  42. * Authors:
  43. * ========
  44. *
  45. *> \author Univ. of Tennessee
  46. *> \author Univ. of California Berkeley
  47. *> \author Univ. of Colorado Denver
  48. *> \author NAG Ltd.
  49. *
  50. *> \date December 2016
  51. *
  52. *> \ingroup single_lin
  53. *
  54. * =====================================================================
  55. SUBROUTINE SERRQP( PATH, NUNIT )
  56. *
  57. * -- LAPACK test routine (version 3.7.0) --
  58. * -- LAPACK is a software package provided by Univ. of Tennessee, --
  59. * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
  60. * December 2016
  61. *
  62. * .. Scalar Arguments ..
  63. CHARACTER*3 PATH
  64. INTEGER NUNIT
  65. * ..
  66. *
  67. * =====================================================================
  68. *
  69. * .. Parameters ..
  70. INTEGER NMAX
  71. PARAMETER ( NMAX = 3 )
  72. * ..
  73. * .. Local Scalars ..
  74. CHARACTER*2 C2
  75. INTEGER INFO, LW
  76. * ..
  77. * .. Local Arrays ..
  78. INTEGER IP( NMAX )
  79. REAL A( NMAX, NMAX ), TAU( NMAX ), W( 3*NMAX+1 )
  80. * ..
  81. * .. External Functions ..
  82. LOGICAL LSAMEN
  83. EXTERNAL LSAMEN
  84. * ..
  85. * .. External Subroutines ..
  86. EXTERNAL ALAESM, CHKXER, SGEQP3
  87. * ..
  88. * .. Scalars in Common ..
  89. LOGICAL LERR, OK
  90. CHARACTER*32 SRNAMT
  91. INTEGER INFOT, NOUT
  92. * ..
  93. * .. Common blocks ..
  94. COMMON / INFOC / INFOT, NOUT, OK, LERR
  95. COMMON / SRNAMC / SRNAMT
  96. * ..
  97. * .. Executable Statements ..
  98. *
  99. NOUT = NUNIT
  100. WRITE( NOUT, FMT = * )
  101. C2 = PATH( 2: 3 )
  102. LW = 3*NMAX + 1
  103. A( 1, 1 ) = 1.0E+0
  104. A( 1, 2 ) = 2.0E+0
  105. A( 2, 2 ) = 3.0E+0
  106. A( 2, 1 ) = 4.0E+0
  107. OK = .TRUE.
  108. *
  109. IF( LSAMEN( 2, C2, 'QP' ) ) THEN
  110. *
  111. * Test error exits for QR factorization with pivoting
  112. *
  113. * SGEQP3
  114. *
  115. SRNAMT = 'SGEQP3'
  116. INFOT = 1
  117. CALL SGEQP3( -1, 0, A, 1, IP, TAU, W, LW, INFO )
  118. CALL CHKXER( 'SGEQP3', INFOT, NOUT, LERR, OK )
  119. INFOT = 2
  120. CALL SGEQP3( 1, -1, A, 1, IP, TAU, W, LW, INFO )
  121. CALL CHKXER( 'SGEQP3', INFOT, NOUT, LERR, OK )
  122. INFOT = 4
  123. CALL SGEQP3( 2, 3, A, 1, IP, TAU, W, LW, INFO )
  124. CALL CHKXER( 'SGEQP3', INFOT, NOUT, LERR, OK )
  125. INFOT = 8
  126. CALL SGEQP3( 2, 2, A, 2, IP, TAU, W, LW-10, INFO )
  127. CALL CHKXER( 'SGEQP3', INFOT, NOUT, LERR, OK )
  128. END IF
  129. *
  130. * Print a summary line.
  131. *
  132. CALL ALAESM( PATH, OK, NOUT )
  133. *
  134. RETURN
  135. *
  136. * End of SERRQP
  137. *
  138. END