| 查看: 888 | 回复: 10 | |||
| 当前主题已经存档。 | |||
落日照秋草木虫 (正式写手)
|
[交流]
【求助】求助LAPACK 问题！【已完结】已有1人参与
|
||
|
请问哪位知道ＺＧＥＴＲＩ这个语句如何使用啊？ [ Last edited by 余泽成 on 2010-3-28 at 19:08 ] |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
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