有没有办法找到哪个直接求解器正在使用 solve()？

Question

正如标题所说，对于R中给定的稀疏矩阵，有没有办法找出哪个矩阵分解正在应用函数solve()？例如，在 Matlab 中有 spparms('spumoni', 2); ，其中 return 一些关于用于求解稀疏矩阵分解的算法的信息。

Answer 1

好吧，运行 R 可以到达

> methods(solve)

这将产生

[1] solve.default solve.qr

如果你输入

> solve.default

你会回来的

function (a, b, tol = .Machine$double.eps, LINPACK = FALSE, ...)
{
    if (!missing(LINPACK))
        warning("the LINPACK argument has been defunct since R 3.1.0")
    if (is.complex(a) || (!missing(b) && is.complex(b))) {
        a <- as.matrix(a)
        if (missing(b)) {
            b <- diag(1 + (0+0i), nrow(a))
            colnames(b) <- rownames(a)
        }
        return(.Internal(La_solve_cmplx(a, b)))
    }
    if (inherits(a, "qr")) {
        warning("solve.default called with a \"qr\" object: use 'qr.solve'")
        return(solve.qr(a, b, tol))
    }
    a <- as.matrix(a)
    if (missing(b)) {
        b <- diag(1, nrow(a))
        colnames(b) <- rownames(a)
    }
    .Internal(La_solve(a, b, tol))
}

这意味着它是 La_solve 或 La_solve_cmplx。查看他们的实现，例如here 可以发现 La_solve 将调用 LAPACK 例程 DGESV，而 La_solve_cmplx 将调用 LAPACK 例程 ZGESV。

简单吧？

有没有办法找到哪个直接求解器正在使用 solve()？

Is there a way to find which direct solver is using solve()?

math

r

algebra

sparse-matrix

matrix-decomposition