我正在尝试使用RcppArmadillo反转一个复数方阵:
# include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
using namespace Rcpp;
using namespace arma;
// [[Rcpp::export]]
cx_mat fn(cx_mat x) {
return x.i();
}
然而,当我使用
sourceCpp
时会出现错误:"undefined reference to zgetri_'
"。如果我只是将cx_mat
替换为mat
,那么它就可以编译并正常工作,但这只适用于实矩阵。使用inv()
会导致相同的编译错误。有趣的是,使用伪逆pinv()
可以通过编译,但与R中的solve()
结果略有不同。# include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
using namespace Rcpp;
using namespace arma;
// [[Rcpp::export]]
cx_mat fn(cx_mat x) {
return pinv(x);
}
> a<-matrix(c(3+0.1i,7+2i,4+0i,2+0.5i),ncol=2)
> a
[,1] [,2]
[1,] 3+0.1i 4+0.0i
[2,] 7+2.0i 2+0.5i
> identical(solve(a), fn(a))
[1] FALSE
> solve(a) - fn(a)
[,1] [,2]
[1,] -6.938894e-17-7.80626e-18i 8.326673e-17-6.93889e-18i
[2,] 1.665335e-16+1.95156e-17i -1.665335e-16+4.16334e-17i
我知道在这种情况下,差异是机器精度级别的调整,但我仍然想知道是否有办法使inv()
或.i()
适用于复杂矩阵。谢谢。