language-agnostic - 如何测试 dpotrf

Question

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

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

Answer 1

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

因此，输入预计是正定的和对称的。

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