c - 在调用 lapack 子例程之前我们需要复制矩阵吗

Question

我正在解决一些问题，其中间步骤需要求解 Ax=b 线性系统。我计划使用 lapack 通过使用像 DGESV 这样的子例程来解决它。然而我发现 lapack 将矩阵“A”转换为“PLU”(一些较低和较高的矩阵格式)。我希望“A”在我 future 的运营中完好无损。我每次发送到 lapack 之前都需要复制吗？复制“b”是可以的，但是，如果“A”很大，那么复制本身就很耗时。有没有办法避免它或 lapack 的其他一些功能可以避免更改“A”。另外，如果复制不可避免，那么有没有一种有效的方法来实现它？我使用“C”作为编程语言。

Answer 1

LAPACK 通常尝试将结果存储在输入参数中，但这里您希望 A 保持不变。

LAPACK 还在内部使用工作数组。您基本上可以通过将原始元素读取到工作数组来遵循类似的路线。从固定源(可能是文件)或只是将其保存在内存中，解决问题获取系数，然后检索原始元素并将它们放置在同一个数组中。