R qr分解的源代码没看懂

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，我最近在C++中利用R中的qr源代码编程一个qr分解的函数。主要是想求得矩阵qr分解后的pivot，R中qr分解的源代码如下

1. qr.default <- function(x, tol = 1e-07, LAPACK = FALSE, ...)
   
2. {
   
3.     x <- as.matrix(x)
   
4.     if(is.complex(x))
   
5.         return(structure(.Call("La_zgeqp3", x, PACKAGE = "base"), class="qr"))
   
6.     ## otherwise :
   
7.     if(!is.double(x))
   
8.         storage.mode(x) <- "double"
   
9.     if(LAPACK) {
   
10.         res <- .Call("La_dgeqp3", x, PACKAGE = "base")
    
11.         if(!is.null(cn <- colnames(x)))
    
12.             colnames(res$qr) <- cn[res$pivot]
    
13.         attr(res, "useLAPACK") <- TRUE
    
14.         class(res) <- "qr"
    
15.         return(res)
    
16.     }
    
17. 
    
18.     p <- ncol(x) # guaranteed to be integer
    
19.     n <- nrow(x)
    
20.     res <- .Fortran("dqrdc2",
    
21.              qr=x,
    
22.              n,
    
23.              n,
    
24.              p,
    
25.              as.double(tol),
    
26.              rank=integer(1L),
    
27.              qraux = double(p),
    
28.              pivot = as.integer(1L:p),
    
29.              double(2*p),
    
30.              PACKAGE="base")[c(1,6,7,8)]# c("qr", "rank", "qraux", "pivot")
    
31.     if(!is.null(cn <- colnames(x)))
    
32.         colnames(res$qr) <- cn[res$pivot]
    
33.     class(res) <- "qr"
    
34.     res
    
35. }

复制代码

C++的Eigen也有QR分解的函数，但是没有pivot 项。哪位高人知道如何将其转换为C++函数。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：源代码 Default fault Pivot 大家好 源代码 小数点