mkl里函数进行cholesky分级用法

一路向北 · 发表于 2021-4-23 20:29:53

想请教下,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:32

由于图片太大,不能上图,很抱歉,希望有人能解答下

li913 · 发表于 2021-4-25 18:28:27

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

一路向北 · 发表于 2021-4-26 14:05:53

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

好的十分感谢

[数学库] mkl里函数进行cholesky分级用法