我正在科学应用中执行一些测试。此应用程序使用Lapack dpotrf Lapack函数。我不太了解线性代数。我必须模拟对dpotrf的大调用,然后生成一个随机矩阵以对其执行dpotrf。
此矩阵必须具有特殊模式吗?我可以只生成一个简单的NxN双矩阵吗?
从dpotrf.f,DPOTRF计算实对称正定矩阵A的Cholesky分解。
因此,输入应为positive-definite并且是对称的。
[此外,dpotrf希望输入为数组,即您的输入为(N * N)x1。我通常将输出调整为NxN。
这里是一个C代码示例:https://gist.github.com/Imxset21/6786ed18ff92908e84c0