代码之家 › 专栏 › 技术社区 › Stéphane Laurent

使用RcppEigen将矩阵x3d转换为矩阵x3d

rcpp eigen casting r

Stéphane Laurent · 技术社区 · 7 年前

与 RcppEigen ,我需要一个 Matrix3d 作为函数的参数。但这是不可能的,它只接受 MatrixXd . 我试着投了以下几张,但没用:

Rcpp::List MtoAxisAngle(Eigen::MatrixXd & m)
{
  Eigen::Matrix3d m3 = m.cast();
  ......

有什么解决办法吗?

1 回复 | 直到 7 年前

Ralf Stubner 7 年前

可以在新矩阵的初始化中使用现有矩阵:

// [[Rcpp::depends(RcppEigen)]]
#include <RcppEigen.h>

// [[Rcpp::export]]
Rcpp::List MtoAxisAngle(Eigen::MatrixXd & m)
{
  Eigen::Matrix3d m3(m);
  return Rcpp::List::create(Rcpp::Named("size") = m3.size());
}


/*** R
MtoAxisAngle(matrix(1:9,3,3))
# MtoAxisAngle(matrix(1:16,4,4))
*/

因为带有不一致矩阵的调用会杀死R,所以应该事先检查大小是否正确。

推荐文章

P_B · 如何从矩阵中减去均值向量

1 年前

Lionnel104 · 特征矩阵与向量之间的元素运算

1 年前

Kota · 使用-O2编译时,特征::Vector2f未正确初始化

1 年前

Riobaldo Tatarana · Eigen映射的奇异行为

1 年前

atb · 生成具有特征的矩阵的列运算的无循环和

2 年前

NateM · 有效张量乘法

7 年前

avgn · 为什么我的特征代码不能扩展到更多线程

7 年前

beesleep · 保存为成员的SelfAdjointEigenSolver的结果在另一个作用域中重新初始化

7 年前

avgn · 特征稠密稀疏矩阵积是线程化的吗?

7 年前

Bubble · 特征FFT库

7 年前