| 查看: 908 | 回复: 10 | |||
| 当前主题已经存档。 | |||
落日照秋草木虫 (正式写手)
|
[交流]
【求助】求助LAPACK 问题！【已完结】 已有1人参与
|
||
|
请问哪位知道ＺＧＥＴＲＩ这个语句如何使用啊？ [ Last edited by 余泽成 on 2010-3-28 at 19:08 ] |
回复此楼» 猜你喜欢
实验室接单子
已经有3人回复
-
要不要辞职读博？
已经有5人回复
-
不自信的我
已经有10人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
dickli2008
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 393.4
- 散金: 112
- 红花: 1
- 帖子: 236
- 在线: 114.5小时
- 虫号: 556129
- 注册: 2008-05-10
- 专业: 凝聚态物性 II ：电子结构
10楼2010-03-22 11:00:36
dic213
木虫 (著名写手)
第三军团儿童团团长
- 应助: 0 (幼儿园)
- 金币: 1777.4
- 散金: 5
- 帖子: 2476
- 在线: 37.1小时
- 虫号: 42931
- 注册: 2004-04-03
- 性别: GG
- 专业: 计算数学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+2):辛苦了 1-27 19:07
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+2):辛苦了 1-27 19:07
|
SUBROUTINE ZGETRI( N, A, LDA, IPIV, WORK, LWORK, INFO ) 016: * Purpose 017: * ======= 018: * 019: * ZGETRI computes the inverse of a matrix using the LU factorization 020: * computed by ZGETRF. 021: * 022: * This method inverts U and then computes inv(A) by solving the system 023: * inv(A)*L = inv(U) for inv(A). 024: * 025: * Arguments 026: * ========= 027: * 028: * N (input) INTEGER 029: * The order of the matrix A. N >= 0. 030: * 031: * A (input/output) COMPLEX*16 array, dimension (LDA,N) 032: * On entry, the factors L and U from the factorization 033: * A = P*L*U as computed by ZGETRF. 034: * On exit, if INFO = 0, the inverse of the original matrix A. 035: * 036: * LDA (input) INTEGER 037: * The leading dimension of the array A. LDA >= max(1,N). 038: * 039: * IPIV (input) INTEGER array, dimension (N) 040: * The pivot indices from ZGETRF; for 1<=i<=N, row i of the 041: * matrix was interchanged with row IPIV(i). 042: * 043: * WORK (workspace/output) COMPLEX*16 array, dimension (MAX(1,LWORK)) 044: * On exit, if INFO=0, then WORK(1) returns the optimal LWORK. 045: * 046: * LWORK (input) INTEGER 047: * The dimension of the array WORK. LWORK >= max(1,N). 048: * For optimal performance LWORK >= N*NB, where NB is 049: * the optimal blocksize returned by ILAENV. 050: * 051: * If LWORK = -1, then a workspace query is assumed; the routine 052: * only calculates the optimal size of the WORK array, returns 053: * this value as the first entry of the WORK array, and no error 054: * message related to LWORK is issued by XERBLA. 055: * 056: * INFO (output) INTEGER 057: * = 0: successful exit 058: * < 0: if INFO = -i, the i-th argument had an illegal value 059: * > 0: if INFO = i, U(i,i) is exactly zero; the matrix is 060: * singular and its inverse could not be computed. 061: * 062: * ===================================================================== http://www.netlib.org/lapack/explore-html/zgetri.f.html |
2楼2010-01-27 17:14:25
落日照秋草
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 2388.6
- 散金: 617
- 红花: 8
- 帖子: 670
- 在线: 912小时
- 虫号: 472183
- 注册: 2007-12-04
- 专业: 药物化学
3楼2010-01-27 17:33:22
dic213
木虫 (著名写手)
第三军团儿童团团长
- 应助: 0 (幼儿园)
- 金币: 1777.4
- 散金: 5
- 帖子: 2476
- 在线: 37.1小时
- 虫号: 42931
- 注册: 2004-04-03
- 性别: GG
- 专业: 计算数学
4楼2010-01-28 08:47:03
lurencyj
木虫 (著名写手)
- 应助: 159 (高中生)
- 金币: 2869.2
- 散金: 520
- 红花: 8
- 沙发: 10
- 帖子: 1244
- 在线: 148.3小时
- 虫号: 888093
- 注册: 2009-10-29
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2010-01-28 09:04:14
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
6楼2010-01-28 10:19:16
ftp0000
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1822.1
- 帖子: 292
- 在线: 6.6小时
- 虫号: 585530
- 注册: 2008-08-08
- 专业: 凝聚态物性 II ：电子结构
7楼2010-02-11 20:23:50
yy3780
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4304.2
- 散金: 50
- 红花: 3
- 沙发: 3
- 帖子: 1079
- 在线: 401.9小时
- 虫号: 653605
- 注册: 2008-11-13
- 专业: 造血相关器官（肝脏/脾脏/
8楼2010-02-12 22:44:59
dickli2008
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 393.4
- 散金: 112
- 红花: 1
- 帖子: 236
- 在线: 114.5小时
- 虫号: 556129
- 注册: 2008-05-10
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+3):谢谢应助，辛苦了！ 2010-03-22 15:44
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+3):谢谢应助，辛苦了！ 2010-03-22 15:44
|
完整的功能提示给你吧： ?getri Computes the inverse of an LU-factored general matrix. Syntax FORTRAN 77: call sgetri( n, a, lda, ipiv, work, lwork, info ) call dgetri( n, a, lda, ipiv, work, lwork, info ) call cgetri( n, a, lda, ipiv, work, lwork, info ) call zgetri( n, a, lda, ipiv, work, lwork, info ) Fortran 95: call getri( a, ipiv [,info] ) Description This routine is declared in mkl_lapack.fi for FORTRAN 77 interface, in lapack.f90 for Fortran 95 interface, and in mkl_lapack.h for C interface. The routine computes the inverse inv(A) of a general matrix A. Before calling this routine, call ?getrf to factorize A. Input Parameters n INTEGER. The order of the matrix A; n ≥ 0. a, work REAL for sgetri DOUBLE PRECISION for dgetri COMPLEX for cgetri DOUBLE COMPLEX for zgetri. Arrays: a(lda,*), work(*). a(lda,*) contains the factorization of the matrix A, as returned by ?getrf: A = P*L*U. The second dimension of a must be at least max(1,n). work(*) is a workspace array of dimension at least max(1,lwork). lda INTEGER. The first dimension of a; lda ≥ max(1, n). ipiv INTEGER. Array, DIMENSION at least max(1, n). The ipiv array, as returned by ?getrf. lwork INTEGER. The size of the work array; lwork ≥ n. If lwork = -1, then a workspace query is assumed; the routine only calculates the optimal size of the work array, returns this value as the first entry of the work array, and no error message related to lwork is issued by xerbla. See Application Notes below for the suggested value of lwork. Output Parameters a Overwritten by the n-by-n matrix inv(A). If info = 0, on exit work(1) contains the minimum value of lwork required for optimum performance. Use this lwork for subsequent runs. work(1) info INTEGER. If info = 0, the execution is successful. If info = -i, the i-th parameter had an illegal value. If info = i, the i-th diagonal element of the factor U is zero, U is singular, and the inversion could not be completed. Fortran 95 Interface Notes Routines in Fortran 95 interface have fewer arguments in the calling sequence than their FORTRAN 77 counterparts. For general conventions applied to skip redundant or reconstructible arguments, see Fortran 95 Interface Conventions. 703 LAPACK Routines: Linear Equations 3 Specific details for the routine getri interface are as follows: a Holds the matrix A of size (n,n). ipiv Holds the vector of length n. Application Notes For better performance, try using lwork = n*blocksize, where blocksize is a machine-dependent value (typically, 16 to 64) required for optimum performance of the blocked algorithm. If you are in doubt how much workspace to supply, use a generous value of lwork for the first run or set lwork = -1. If you choose the first option and set any of admissible lwork sizes, which is no less than the minimal value described, the routine completes the task, though probably not so fast as with a recommended workspace, and provides the recommended workspace in the first element of the corresponding array work on exit. Use this value (work(1)) for subsequent runs. If you set lwork = -1, the routine returns immediately and provides the recommended workspace in the first element of the corresponding array (work). This operation is called a workspace query. Note that if you set lwork to less than the minimal required value and not -1, the routine returns immediately with an error exit and does not provide any information on the recommended workspace. The computed inverse X satisfies the following error bound: |XA - I| ≤ c(n)ε|X|P|L||U|, where c(n) is a modest linear function of n; ε is the machine precision; I denotes the identity matrix; P, L, and U are the factors of the matrix factorization A = P*L*U. The total number of floating-point operations is approximately (4/3)n3 for real flavors and (16/3)n3 for complex flavors. |
9楼2010-03-22 10:59:02