matrix - Lapack + c + 矩阵

Question

你能帮我理解如何调用 Lapack 例程吗？特别是，我想调用 DTRTRI (Uplo, DIAG, N, DA, LDA, INFO)。

Answer 1

查看Armadillo C++ Library的源代码，它调用了许多 Lapack 函数。

Armadillo 直接从 C++ 调用 Fortran 函数，无需使用任何 C 语言包装器或转换 Fortran 代码。

Answer 2

只需为 LAPACK 使用现有的C 语言包装器之一（它本身是用 FORTRAN 编写的）。

具体来说，尝试lapwrapc。

它有一个在 src/lapack_dtrtri.c 中定义的 DTRTRI 的 C 接口

void lapack_dtrtri( const enum lapack_uplo_type uplo, 
                    const enum lapack_diag_type diag, 
                    const int n, 
                    double * a, 
                    const int lda, 
                    int * info )