Subroutine */
| 7123 | } /* dgetf2_ */ |
| 7124 | |
| 7125 | /* Subroutine */ int dgetrf_(integer *m, integer *n, doublereal *a, integer * |
| 7126 | lda, integer *ipiv, integer *info) |
| 7127 | { |
| 7128 | /* System generated locals */ |
| 7129 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5; |
| 7130 | |
| 7131 | /* Local variables */ |
| 7132 | static integer i__, j, jb, nb; |
| 7133 | extern /* Subroutine */ int dgemm_(char *, char *, integer *, integer *, |
| 7134 | integer *, doublereal *, doublereal *, integer *, doublereal *, |
| 7135 | integer *, doublereal *, doublereal *, integer *); |
| 7136 | static integer iinfo; |
| 7137 | extern /* Subroutine */ int dtrsm_(char *, char *, char *, char *, |
| 7138 | integer *, integer *, doublereal *, doublereal *, integer *, |
| 7139 | doublereal *, integer *), dgetf2_( |
| 7140 | integer *, integer *, doublereal *, integer *, integer *, integer |
| 7141 | *), xerbla_(char *, integer *); |
| 7142 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 7143 | integer *, integer *, ftnlen, ftnlen); |
| 7144 | extern /* Subroutine */ int dlaswp_(integer *, doublereal *, integer *, |
| 7145 | integer *, integer *, integer *, integer *); |
| 7146 | |
| 7147 | |
| 7148 | /* |
| 7149 | -- LAPACK routine (version 3.2) -- |
| 7150 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 7151 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 7152 | November 2006 |
| 7153 | |
| 7154 | |
| 7155 | Purpose |
| 7156 | ======= |
| 7157 | |
| 7158 | DGETRF computes an LU factorization of a general M-by-N matrix A |
| 7159 | using partial pivoting with row interchanges. |
| 7160 | |
| 7161 | The factorization has the form |
| 7162 | A = P * L * U |
| 7163 | where P is a permutation matrix, L is lower triangular with unit |
| 7164 | diagonal elements (lower trapezoidal if m > n), and U is upper |
| 7165 | triangular (upper trapezoidal if m < n). |
| 7166 | |
| 7167 | This is the right-looking Level 3 BLAS version of the algorithm. |
| 7168 | |
| 7169 | Arguments |
| 7170 | ========= |
| 7171 | |
| 7172 | M (input) INTEGER |
| 7173 | The number of rows of the matrix A. M >= 0. |
| 7174 | |
| 7175 | N (input) INTEGER |
| 7176 | The number of columns of the matrix A. N >= 0. |
| 7177 | |
| 7178 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 7179 | On entry, the M-by-N matrix to be factored. |
| 7180 | On exit, the factors L and U from the factorization |
| 7181 | A = P*L*U; the unit diagonal elements of L are not stored. |
| 7182 |