Subroutine */
| 6947 | } /* sgesv_ */ |
| 6948 | |
| 6949 | /* Subroutine */ int sgetf2_(integer *m, integer *n, real *a, integer *lda, |
| 6950 | integer *ipiv, integer *info) |
| 6951 | { |
| 6952 | /* System generated locals */ |
| 6953 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 6954 | real r__1; |
| 6955 | |
| 6956 | /* Local variables */ |
| 6957 | static integer i__, j, jp; |
| 6958 | extern /* Subroutine */ int sger_(integer *, integer *, real *, real *, |
| 6959 | integer *, real *, integer *, real *, integer *), sscal_(integer * |
| 6960 | , real *, real *, integer *); |
| 6961 | static real sfmin; |
| 6962 | extern /* Subroutine */ int sswap_(integer *, real *, integer *, real *, |
| 6963 | integer *); |
| 6964 | extern doublereal slamch_(char *); |
| 6965 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 6966 | extern integer isamax_(integer *, real *, integer *); |
| 6967 | |
| 6968 | |
| 6969 | /* |
| 6970 | -- LAPACK routine (version 3.2) -- |
| 6971 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 6972 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 6973 | November 2006 |
| 6974 | |
| 6975 | |
| 6976 | Purpose |
| 6977 | ======= |
| 6978 | |
| 6979 | SGETF2 computes an LU factorization of a general m-by-n matrix A |
| 6980 | using partial pivoting with row interchanges. |
| 6981 | |
| 6982 | The factorization has the form |
| 6983 | A = P * L * U |
| 6984 | where P is a permutation matrix, L is lower triangular with unit |
| 6985 | diagonal elements (lower trapezoidal if m > n), and U is upper |
| 6986 | triangular (upper trapezoidal if m < n). |
| 6987 | |
| 6988 | This is the right-looking Level 2 BLAS version of the algorithm. |
| 6989 | |
| 6990 | Arguments |
| 6991 | ========= |
| 6992 | |
| 6993 | M (input) INTEGER |
| 6994 | The number of rows of the matrix A. M >= 0. |
| 6995 | |
| 6996 | N (input) INTEGER |
| 6997 | The number of columns of the matrix A. N >= 0. |
| 6998 | |
| 6999 | A (input/output) REAL array, dimension (LDA,N) |
| 7000 | On entry, the m by n matrix to be factored. |
| 7001 | On exit, the factors L and U from the factorization |
| 7002 | A = P*L*U; the unit diagonal elements of L are not stored. |
| 7003 | |
| 7004 | LDA (input) INTEGER |
| 7005 | The leading dimension of the array A. LDA >= max(1,M). |
| 7006 |