如何测试dpotrf

Question

我在一个科学应用程序中执行一些测试。此应用程序使用Lapack dpotrf Lapack函数。我对线性代数不是很了解。我必须模拟一个对dpotrf的大调用，然后生成一个随机矩阵来对其执行dpotrf。

这个矩阵必须有一个特殊的模式吗？我可以只生成一个简单的NxN双矩阵吗？

Answer 1

从dpotrf.f出发，DPOTRF计算实对称正定矩阵A的Cholesky分解。

因此，输入应该是positive-definite和对称的。

另外，dpotrf期望输入是一个数组，即您的输入是(N*N)x1。我通常将输出重塑为NxN。