Fortran Coder

标题: mkl里函数进行cholesky分级用法 [打印本页]

作者: 一路向北 时间: 2021-4-23 20:29
标题: mkl里函数进行cholesky分级用法

想请教下,mkl里的函数 dpstf2(uplo, n, a, lda, piv, rank, tol, work, info)
1.这个函数work填啥
2.piv类型为什么是integer, piv不是P*A*PT里的p吗,

作者: 一路向北 时间: 2021-4-23 20:30
由于图片太大,不能上图,很抱歉,希望有人能解答下

作者: li913 时间: 2021-4-25 18:28
1、work是工作数组，实质上就是程序执行时需要的临时空间，按照手册开辟足够空间就行，可大不可小。
2、手册里可没有说 piv 是 p，而是写的 P（piv(k),k）=1

作者: 一路向北 时间: 2021-4-26 14:05

li913 发表于 2021-4-25 18:28
1、work是工作数组，实质上就是程序执行时需要的临时空间，按照手册开辟足够空间就行，可大不可小。
2、手 ...

好的十分感谢

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)