如何在Armadillo或RcppArmadillo中将矩阵复制到立方体中
在Armadillo或RcppArmadillo中，可以使用reshape()函数将矩阵复制到立方体中。
reshape()函数用于将一个矩阵重塑为指定维度的数组。在将矩阵复制到立方体中时，可以使用reshape()函数将矩阵的行和列重塑为立方体的高度、宽度和深度。
以下是使用reshape()函数将矩阵复制到立方体的示例代码：
#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::export]]
arma::cube copyMatrixToCube(const arma::mat& matrix, int height, int width, int depth) {
// Reshape the matrix to a cube
arma::cube cube(matrix.memptr(), height, width, depth, false);
return cube;
}在上述代码中，copyMatrixToCube()函数接受一个矩阵(matrix)和立方体的高度(height)、宽度(width)和深度(depth)作为参数。函数内部使用reshape()函数将矩阵复制到一个立方体(cube)中，并返回该立方体。
使用Armadillo和RcppArmadillo进行矩阵和立方体操作时，可以参考以下腾讯云相关产品和产品介绍链接地址：
- 腾讯云服务器（云服务器ECS）：提供高性能、可扩展的云服务器实例，适用于各种计算场景。链接地址：https://cloud.tencent.com/product/cvm
- 腾讯云容器服务（容器实例 TKE）：提供高性能、高可靠的容器化应用运行环境，支持快速部署和弹性伸缩。链接地址：https://cloud.tencent.com/product/tke
- 腾讯云对象存储（对象存储 COS）：提供安全、稳定、低成本的云端存储服务，适用于海量数据存储和访问。链接地址：https://cloud.tencent.com/product/cos
请注意，以上链接仅供参考，具体产品选择应根据实际需求进行评估和决策。
相关·内容
我有一个具有正定矩阵的三维数组，我想得到一个与所有矩阵的Cholesky因子相同大小的数组。我使用的是Armadillo库和cube类型，为此我尝试使用方便的函数each_slice。但是我没有让lambda表达式正确地工作，所以希望有人能帮助我指出我的错误。
下面是一个很小的例子：
// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
// [[Rcpp::export]]
arma::cube chol_array(arma::cube Sigma) {
arma::cube Sigma_chol =
浏览 3提问于2018-01-05得票数 2
回答已采纳
在R中，我可以根据它们的索引提取矩阵元素，如下所示
> m <- matrix(1:6, nrow = 3)
> m
[,1] [,2]
[1,] 1 4
[2,] 2 5
[3,] 3 6
> row_index <- c(1, 2)
> col_index <- c(2, 2)
> m[cbind(row_index, col_index)]
[1] 4 5
Armadillo / Rcpp::Armadillo有没有一种原生的方法来做到这一点？我能做的最多的是一个自定义函数，它使用行和列索引来计
浏览 0提问于2016-05-28得票数 2
2回答
这是我刚发现的东西。由于某种原因，在Armadillo中，将稠密矩阵乘以稀疏矩阵比乘以稀疏和稠密矩阵(即倒序)要慢得多。
// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
// [[Rcpp::export]]
arma::sp_mat mult_sp_den_to_sp(arma::sp_mat& a, arma::mat& b)
{
// sparse x dense -> sparse
arma::sp_mat result(a * b);
return
浏览 0提问于2018-04-11得票数 3
回答已采纳
我和Armadillo和RcppArmadillo有个奇怪的问题。我正在创建一个由零值填充的立方体，我希望将特定的元素转换为零值。然而，当我使用辅助功能时，其他元素长乐的值会稍微变小，并且常常与nan相等。有人知道是什么原因造成的吗？
示例：
#include <RcppArmadillo.h>
using namespace arma;
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::export]]
cube testc() {
cube tester = cube(10,10,2);
uvec indexes =
浏览 4提问于2021-01-10得票数 1
回答已采纳
在我的RcppArmadillo项目中，我有一些矩阵(例如mat A，B，C;)和一个S4对象，例如D(来自R中的外部函数)。由于我需要在这些矩阵和D之间进行一些计算，所以我希望在RcppArmadillo中将"D“转换为合适的数据类型，例如arma::mat D。有可能吗？做这件事最好的方法是什么？这是一个类似的代码：
#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::export]]
Rcpp::List func1(arma::mat A, arma::mat B){
浏览 3提问于2017-06-12得票数 1
回答已采纳
2回答
问候和问候，
我试图使用字段对象类型，而不是列表数据类型，以避免发出copy命令。我试图这样做，以减少从列表中删除一个矩阵时所需的时间，这个矩阵已经定义在armadillo的数据结构中。(例如移除arma::mat * as<arma::mat>(NumericMatrix) )。
具体来说，我试图将一个字段传递到一个函数中，并对其进行操作。编译时收到的错误消息是：
没有用于初始化“arma：：field< arma::Mat >”的匹配构造函数
函数定义上的错误：#include <RcppArmadillo.h>
注意:目标是能够在函数声明中使用
浏览 3提问于2014-11-18得票数 0
回答已采纳
我试图用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 it：“zgetri_'的未定义引用”时，这会引发一个错误。如果我将cx_mat替换为mat，那么它编译和工作都很好，但是它只适用于实际的矩阵。使用inv()会
浏览 7提问于2016-01-02得票数 3
回答已采纳
我正在尝试通过RInside在C++中使用R。我在向R传递armadillo矩阵并返回结果时遇到了问题。下面我可以从一个R库函数返回一个结果，但是我得到了错误的结果。我正在使用moments包中的skewness函数作为示例，它在R中正常工作。我检查了RInside中的示例，但我仍然不确定如何使用RcppArmadillo。我如何正确地将c++中的armadillo矩阵传递给R？
#include <RInside.h>
#include <RcppArmadillo.h>
using namespace
浏览 4提问于2012-11-14得票数 4
回答已采纳
解决方案现在在中联机。
我从RcppArmadillo中的mvtnorm包中重新实现了dmvnorm。我某种程度上喜欢Armadillo，但我想它也适用于普通的Rcpp。从and范数出发的方法是基于mahalanobis距离，所以我有一个函数，然后是多元正规密度函数。
让我给你看我的密码：
#include <RcppArmadillo.h>
#include <Rcpp.h>
// [[Rcpp::depends("RcppArmadillo")]]
// [[Rcpp::export]]
arma::vec mahalanobis_arma(
浏览 5提问于2013-07-12得票数 3
回答已采纳
我有一个很大的矩阵。我试着用动态变化的重量来做样本。由于它被迫在R中使用循环，所以我尝试在Rcpp中实现它，这样它就有机会运行得更快一些。在做了一些实验之后，我想我已经想出了如何用正确的权重随机获取一个索引。
诀窍是，我只是在任何给定时间从列的子集中取样(如果它在C中更有效--矩阵实际上是对称的)，这可以改变为行。我的索引只为这个列的子集定义。在R，我会做一些事情
large_matrix[, columns_of_interest][index]
这个很好用。我将如何使用Rcpp/Armadillo进行等效的操作？我猜
cppFunction("arma::vec element_f
浏览 3提问于2018-06-04得票数 1
回答已采纳
1回答
假设我必须用矩阵A和B编写一些RcppArmadillo代码，用OpenMP创建一个矩阵，它包含3列和3行，等于A的列数乘以B的行数。
我编写了这段代码，但当我试图运行它时，它会崩溃。我最好的猜测是，我在创建row变量时出错了，但我不知道如何修复它。
#include "RcppArmadillo.h"
#include <omp.h>
// [[Rcpp::depends(RcppArmadillo)]]
using namespace Rcpp;
// [[Rcpp::plugins(openmp)]]
// [[Rcpp::export]]
arma::ma
浏览 0提问于2019-02-16得票数 0
回答已采纳
我正在编写一个包，它的一些函数调用了RcppArmadillo中的RcppArmadillo::sample。但是，我在编译时遇到了以下错误。
包含在函数的Citrus.cpp:2: Citrus.cpp:2中的文件中：：RcppArmadillo：：sample(const&，int，bool，(Rcpp：：n数控向量)，T=arma：：subview_col: Citrus.cpp:241:从这里实例化了./R/x86_64-unknown-linux-gnu-library/3.0/RcppArmadillo/include/RcppArmadilloExtensions/sa
浏览 1提问于2017-04-03得票数 0
我所处的情况是，我对Armadillo对象有一个Rcpp::XPtr (例如，arma::Mat，它可能是支持的数据类型之一的矩阵)。现在，我想编写一个查询元素数量的函数。到目前为止，我能想到的最好方法是：(受启发)：
#define DISPATCH_DATA_TYPE(CALL) \
{ \
switch (data_type)
浏览 0提问于2019-07-31得票数 1
回答已采纳
我正在通过RcppArmadillo使用armadillo和R。
我有一个稀疏矩阵和一个行号作为输入。我想搜索矩阵的相应行，并返回非零的位置。
到目前为止，我的函数看起来像
// [[Rcpp::export]]
arma::uvec findAdjacentStates(arma::sp_mat adjacency, int row) {
arma::uvec out(adjacency.n_cols);
arma::SpSubview<double>::const_iterator start = adjacency.row(row).begin();
arma::
浏览 1提问于2017-11-07得票数 0
3回答
我是C++编程的新手(使用Rcpp无缝集成到R中)，我希望得到一些关于如何加快计算速度的建议。
请考虑以下示例：
testmat <- matrix(1:9, nrow=3)
testvec <- 1:3
testmat*testvec
# [,1] [,2] [,3]
#[1,] 1 4 7
#[2,] 4 10 16
#[3,] 9 18 27
在这里，R回收了testvec，因此，松散地说，testvec“变成了一个与testmat相同维的矩阵，用于这种乘法。然后返回Hadamard产品。
浏览 1提问于2014-07-24得票数 7
回答已采纳
我想知道如何使用RcppArmadillo从QR分解返回“紧凑”表单。使用RcppArmadillo在矩阵上运行QR分解相对简单，如下所示：
library(inline)
src <- '
using namespace Rcpp;
using namespace arma;
mat X = as<mat>(A);
mat Q, R;
qr(Q,R,X);
return Rcpp::List::create(Rcpp::Named("Q")=Q,Rcpp::Named("R")=R);
'
fun <- cxxf
浏览 4提问于2015-03-02得票数 2
新的find_finite和find_nonfinite函数在Armadillo 4.300是伟大的补充！但是，在我使用Rcpp的测试中，它们比标准循环慢2.5倍。下面是一些代码，用于计算与R的na.rm=TRUE选项相对应的按情况删除的和和平均值。R的性能基准测试表明，第一个版本(sum_arma和mean_arma)比循环快3.5倍。我做的一切都是对的？有什么办法可以提高性能吗？
C++代码
#include <numeric>
#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
// [
浏览 4提问于2014-05-07得票数 2
回答已采纳
我试图使用RCPP返回一组矩阵。我下面的代码效率极低。我想知道下面的代码是否有效。
#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::export]]
Rcpp::List hello(
const arma::rowvec& g,
const int& n,
const int& p,
const arma::mat& S,
const arma::mat& zc,
const arma::
浏览 3提问于2019-12-19得票数 1
回答已采纳
1回答
请注意，这个错误是从一个更大的上下文中提取的，我显然不能在这里完全报告。
我在文件fun.cpp中有以下函数
#include <RcppArmadilloExtensions/sample.h>
using namespace Rcpp;
// [[Rcpp::depends(RcppArmadillo)]]
arma::vec colMeans(arma::mat data){
int n_0 = data.n_rows;
arma::vec xbar(data.n_cols);
for(int i = 0; i < data.n_rows;
浏览 0提问于2017-02-15得票数 1
回答已采纳
1回答
我猜想，或者更确切地说，我希望，我有一个单一的可修复的问题，或者可能有许多较小的问题，应该放弃。无论哪种方式，我对Rcpp都比较陌生，在并行计算方面非常不知情，无法在线找到解决方案。
问题是，R或R中的“致命错误”被困在一个循环中，就像10次迭代中的5分钟，当非并行版本同时进行5K的迭代时，粗略地说。
由于这个算法适合于一个更大的项目，我调用了其他几个函数，它们都在Rcpp中，我只用'arma‘对象重写了它们，因为这似乎有助于其他人，这里。我还使用我在Rcpp中编写的“热映射”优化器运行优化部分，同样只在“arma”中运行，但没有改进--我还应该指出，返回的部分是'arma::
浏览 0提问于2018-07-05得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
