Fortran—将变量传递到CGESV中

Question

我正在尝试测试LAPACK方法CGESV，但遇到了一个问题。我想在代码的其他部分重用我的"A"矩阵，但当我将其传递到方法中时，它会发生变化。"A"的定义：

(input/output) COMPLEX array, dimension (LDA,N)
      On entry, the N-by-N coefficient matrix A.
      On exit, the factors L and U from the factorization
      A = P*L*U; the unit diagonal elements of L are not stored.

在将a的值传递到CGESV之后，除了创建一个临时变量来存储值之外，有没有一种方法可以保持a的值？

Answer 1

正如您已经注意到的，A矩阵被P*L*U分解覆盖。如果矩阵的大小不是很大，可以复制A矩阵的内容，并使用该副本进行分解。

         CALL CCOPY(N*N, A, 1, A_NEW, 1)

如果矩阵大小太大，以至于无法在内存中保留两个副本，则可以使用分解后的矩阵执行数学运算。例如计算y=A*x

*        y = x
         CALL CCOPY(N, X, 1, Y, 1)
*        y = U * y
         CALL CTRMV('Upper', 'No transpose', 'Non-unit', N, A, N, Y, 1)
*        y = L * y
         CALL CTRMV('Lower', 'No transpose', 'Unit', N, A, N, Y, 1)
*        y = P * y
         CALL DLASWP( 1, Y, N, 1, N, IPIV, 1 )

所需的额外存储器是大小为IPIV的整数N。

Answer 2

例程在适当的位置执行它们的工作，因此保留原始数组的唯一方法是制作一个副本。