Subroutine */
| 6960 | } /* dgesv_ */ |
| 6961 | |
| 6962 | /* Subroutine */ int dgetf2_(integer *m, integer *n, doublereal *a, integer * |
| 6963 | lda, integer *ipiv, integer *info) |
| 6964 | { |
| 6965 | /* System generated locals */ |
| 6966 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 6967 | doublereal d__1; |
| 6968 | |
| 6969 | /* Local variables */ |
| 6970 | static integer i__, j, jp; |
| 6971 | extern /* Subroutine */ int dger_(integer *, integer *, doublereal *, |
| 6972 | doublereal *, integer *, doublereal *, integer *, doublereal *, |
| 6973 | integer *), dscal_(integer *, doublereal *, doublereal *, integer |
| 6974 | *); |
| 6975 | static doublereal sfmin; |
| 6976 | extern /* Subroutine */ int dswap_(integer *, doublereal *, integer *, |
| 6977 | doublereal *, integer *); |
| 6978 | |
| 6979 | extern integer idamax_(integer *, doublereal *, integer *); |
| 6980 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 6981 | |
| 6982 | |
| 6983 | /* |
| 6984 | -- LAPACK routine (version 3.2) -- |
| 6985 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 6986 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 6987 | November 2006 |
| 6988 | |
| 6989 | |
| 6990 | Purpose |
| 6991 | ======= |
| 6992 | |
| 6993 | DGETF2 computes an LU factorization of a general m-by-n matrix A |
| 6994 | using partial pivoting with row interchanges. |
| 6995 | |
| 6996 | The factorization has the form |
| 6997 | A = P * L * U |
| 6998 | where P is a permutation matrix, L is lower triangular with unit |
| 6999 | diagonal elements (lower trapezoidal if m > n), and U is upper |
| 7000 | triangular (upper trapezoidal if m < n). |
| 7001 | |
| 7002 | This is the right-looking Level 2 BLAS version of the algorithm. |
| 7003 | |
| 7004 | Arguments |
| 7005 | ========= |
| 7006 | |
| 7007 | M (input) INTEGER |
| 7008 | The number of rows of the matrix A. M >= 0. |
| 7009 | |
| 7010 | N (input) INTEGER |
| 7011 | The number of columns of the matrix A. N >= 0. |
| 7012 | |
| 7013 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 7014 | On entry, the m by n matrix to be factored. |
| 7015 | On exit, the factors L and U from the factorization |
| 7016 | A = P*L*U; the unit diagonal elements of L are not stored. |
| 7017 | |
| 7018 | LDA (input) INTEGER |
| 7019 | The leading dimension of the array A. LDA >= max(1,M). |